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

f1405602

double AskFaceAreasByPK(tag_t face_tag, double tol)
2 i; ~/ |& w8 J6 N{
2 a8 `$ ?4 N" Z' B        tag_t       ps_tag = NULL_TAG;
+ q4 m: r0 F. [  ?3 S$ k$ o        double     amount[3] = { 0.0, 0.0, 0.0 };
        double     mass[3] = { 0.0, 0.0, 0.0 };
: x* |  W( [' E0 o  M$ K        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 };
  X, f# L& `- N1 R        double     periphery[3] = { 0.0, 0.0, 0.0 };
0 H. R+ v( V  w8 e: z/ {& y        double     areas = 0.0;
& b" @) |/ h% S, z( B+ ^        PK_FACE_t  pkFace = 0;
. s, _: x* C. E( c* J        PK_TOPOL_eval_mass_props_o_t props_o_t;       
6 B7 r4 b+ C; R9 _( d        if (NULL_TAG == face_tag || UF_PS_ask_ps_tag_of_object(face_tag, &ps_tag))
                return 0;
        pkFace = ps_tag;
& L  {# F2 P9 F' n0 i        if (tol < 0.99)
! o/ D# n6 F4 Z+ k- W+ Y$ R* r5 V                tol = 0.99;
$ R8 A  Y) J1 F( r) k        if (tol > 1.0)
                tol = 1.0;
        PK_TOPOL_eval_mass_props_o_m(props_o_t);
: h" e" W7 Q5 ~! D% d- [" f! L        props_o_t.mass = PK_mass_c_of_g_c;
0 @  y9 `5 v0 A4 \+ H/ \        PK_TOPOL_eval_mass_props(1, &pkFace, tol, &props_o_t, amount, mass, c_of_g, m_of_i, periphery);
: o6 R+ P: }" r        areas = amount[0] * 1000.0 * 1000.0;
        return areas;
}
2 b$ ^; v6 j4 T) ]: e7 _