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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
1 V2 a- r) r- S0 n7 M8 Y{
        tag_t       ps_tag = NULL_TAG;
        double     amount[3] = { 0.0, 0.0, 0.0 };
        double     mass[3] = { 0.0, 0.0, 0.0 };
$ _5 e: h) M- G/ ?        double     c_of_g[3] = { 0.0, 0.0, 0.0 };
' S/ s) ]* ?6 T4 G  j# g- ^        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 };
/ i8 P' P. _& _9 N% w; J        double     areas = 0.0;
5 Y- W1 {, o. c1 E        PK_FACE_t  pkFace = 0;
        PK_TOPOL_eval_mass_props_o_t props_o_t;       
) ]* a! m# ~5 w: e. G        if (NULL_TAG == face_tag || UF_PS_ask_ps_tag_of_object(face_tag, &ps_tag))
                return 0;
        pkFace = ps_tag;
/ u; W4 N9 |( F5 I- R- @* E, d        if (tol < 0.99)
                tol = 0.99;
9 }+ V+ H$ v" c& w% C1 L) b& u        if (tol > 1.0)
                tol = 1.0;
        PK_TOPOL_eval_mass_props_o_m(props_o_t);
        props_o_t.mass = PK_mass_c_of_g_c;
& ]$ P$ b2 `) O, ^        PK_TOPOL_eval_mass_props(1, &pkFace, tol, &props_o_t, amount, mass, c_of_g, m_of_i, periphery);
2 u4 q  G3 e1 ]5 o2 f6 V        areas = amount[0] * 1000.0 * 1000.0;
& ]' a9 r/ x) L: }        return areas;
5 `1 j+ r6 V# R+ m! d& O}
6 @+ X: ?( `/ I' ]% x2 F5 j- g7 e% R
: w6 \/ ]0 a  n