通过PK内核函数获取面的面积

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
{
5 o7 i, f6 s! B1 q' g/ f6 Q8 ?        tag_t       ps_tag = NULL_TAG;
& C4 h+ J- L; X: K6 t7 V' a( H' ?        double     amount[3] = { 0.0, 0.0, 0.0 };
        double     mass[3] = { 0.0, 0.0, 0.0 };
% H! \% [7 y, d6 a2 ?0 s  S& o        double     c_of_g[3] = { 0.0, 0.0, 0.0 };
3 |- f/ J; R+ q% V        double     m_of_i[6] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
        double     periphery[3] = { 0.0, 0.0, 0.0 };
        double     areas = 0.0;
2 |& @0 Y/ Z" O; a+ o# g" n5 d        PK_FACE_t  pkFace = 0;
1 C, h  K# L: l; @6 G! ~0 x+ h3 W        PK_TOPOL_eval_mass_props_o_t props_o_t;       
        if (NULL_TAG == face_tag || UF_PS_ask_ps_tag_of_object(face_tag, &ps_tag))
                return 0;
        pkFace = ps_tag;
        if (tol < 0.99)
                tol = 0.99;
1 v* m- w$ I5 x3 r' I2 J5 {& F7 x        if (tol > 1.0)
                tol = 1.0;
        PK_TOPOL_eval_mass_props_o_m(props_o_t);
! n% e5 W, t9 q; P        props_o_t.mass = PK_mass_c_of_g_c;
        PK_TOPOL_eval_mass_props(1, &pkFace, tol, &props_o_t, amount, mass, c_of_g, m_of_i, periphery);
  C( F' h) N4 R9 U/ `) x        areas = amount[0] * 1000.0 * 1000.0;
$ E6 ]7 W% N: c7 s3 [        return areas;
}
0 f* Q: ~) z3 x' q( \