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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
' P8 Z) r" \, o  q5 E{
' t) [$ _/ p8 `, t0 r" V( m- _        tag_t       ps_tag = NULL_TAG;
        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 };
        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 };
! h  s" f: B* j1 r! S  t2 p        double     areas = 0.0;
' N9 H! P( {+ S* H5 Y& _( }4 X        PK_FACE_t  pkFace = 0;
        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;
' V' k& E/ @8 J! X% v$ n5 l( Q2 _/ `        if (tol > 1.0)
                tol = 1.0;
8 I$ h8 Z3 z7 Y8 J" o; t        PK_TOPOL_eval_mass_props_o_m(props_o_t);
  H: [( c- x* m4 S9 |" Z        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);
# z+ o; f7 w# a. R2 |- E! x        areas = amount[0] * 1000.0 * 1000.0;
( J! b- i. A% \1 W8 r+ F        return areas;
& Q+ t. N9 O2 c7 A0 z5 S" `}