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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
{
        tag_t       ps_tag = NULL_TAG;
$ W4 H5 I$ v, _2 f3 A7 y* k* G        double     amount[3] = { 0.0, 0.0, 0.0 };
7 A/ @1 i3 S" X0 X        double     mass[3] = { 0.0, 0.0, 0.0 };
& {6 J  B* Y* o6 _# S        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 };
- i* s( y- q# y. F        double     periphery[3] = { 0.0, 0.0, 0.0 };
        double     areas = 0.0;
% b% Z0 G/ e1 k' c        PK_FACE_t  pkFace = 0;
1 n! L3 j8 `1 K2 W5 z        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))
! a( G0 H. _% W0 v' a4 N                return 0;
% H. b' p9 Q, W, _8 {9 o        pkFace = ps_tag;
6 \2 I& Z* z, T8 `8 O% t- P4 b        if (tol < 0.99)
                tol = 0.99;
        if (tol > 1.0)
                tol = 1.0;
1 Z6 c$ [2 h# {3 G/ c% i* L9 p0 r        PK_TOPOL_eval_mass_props_o_m(props_o_t);
6 b& T# h% k! C, W/ q. l        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);
        areas = amount[0] * 1000.0 * 1000.0;
        return areas;
}
% J$ }1 B; c9 u
( Z- m+ o6 N. z2 V9 a+ m9 |# \