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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
{
        tag_t       ps_tag = NULL_TAG;
        double     amount[3] = { 0.0, 0.0, 0.0 };
  H- a3 {& S; J! \$ q( F        double     mass[3] = { 0.0, 0.0, 0.0 };
        double     c_of_g[3] = { 0.0, 0.0, 0.0 };
        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 };
6 D4 G5 i/ ~6 N3 w1 Q! @: l        double     areas = 0.0;
) C! f; S% j8 W0 J! `" @        PK_FACE_t  pkFace = 0;
6 Q% }3 N. w6 d1 z7 D& J        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))
# l# D! M( h# m6 O) ?% @! u                return 0;
$ t6 y+ q) t3 j' B) m; z        pkFace = ps_tag;
        if (tol < 0.99)
                tol = 0.99;
        if (tol > 1.0)
                tol = 1.0;
        PK_TOPOL_eval_mass_props_o_m(props_o_t);
+ n# }& z3 V& T8 c8 U$ g        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);
( t" A* i! P0 e  o5 N% F4 I1 w        areas = amount[0] * 1000.0 * 1000.0;
9 S& N9 z5 N+ R7 c, _& M        return areas;
}
7 t3 q" H/ [% d( y! q- c$ |
& S+ `& C: t2 G- a, w: }