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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
{
        tag_t       ps_tag = NULL_TAG;
7 t3 w  a* a, |0 Z5 z; n+ k        double     amount[3] = { 0.0, 0.0, 0.0 };
        double     mass[3] = { 0.0, 0.0, 0.0 };
        double     c_of_g[3] = { 0.0, 0.0, 0.0 };
; B- u% S- n. |+ G  h, w        double     m_of_i[6] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
" c: [  B% n- p4 L/ @0 [        double     periphery[3] = { 0.0, 0.0, 0.0 };
        double     areas = 0.0;
- [2 G* Z3 Q' m        PK_FACE_t  pkFace = 0;
        PK_TOPOL_eval_mass_props_o_t props_o_t;       
+ M0 q9 N: G* ]' M  B+ r        if (NULL_TAG == face_tag || UF_PS_ask_ps_tag_of_object(face_tag, &ps_tag))
  W" b+ _9 d' }1 `                return 0;
        pkFace = ps_tag;
        if (tol < 0.99)
1 E' G( P0 \) @8 w- b                tol = 0.99;
        if (tol > 1.0)
+ }1 f5 n9 ~) m7 X                tol = 1.0;
        PK_TOPOL_eval_mass_props_o_m(props_o_t);
' c8 o7 a6 ^" s  T2 g' y        props_o_t.mass = PK_mass_c_of_g_c;
6 u  v% t. S, h: n- R" t        PK_TOPOL_eval_mass_props(1, &pkFace, tol, &props_o_t, amount, mass, c_of_g, m_of_i, periphery);
7 H; J# m& D: s* |) Y; @        areas = amount[0] * 1000.0 * 1000.0;
        return areas;
}

% g2 {, q1 N: o4 L9 C