2

问题

python - 如何计算python中网格外表面的指定函数的表面积分？

代码

我试图创建一个最小的可行示例。在此示例中，我尝试计算使用trimesh定义的框的表面积。

在此处输入图像描述

mesh.area我只计算表面积，因为可以使用trimesh包（

这就是为什么我转向quadpy进行表面集成的原因。

这是代码：

import numpy as np
import quadpy
import trimesh


# CREATE THE BOX
# object (box)
box = trimesh.creation.box(extents=[0.3, 0.3, 0.3])

# show
box.show(viewer='gl')


# SURFACE INTEGRATION 

# 1 GET TRIANGLES of trimesh mesh object
triangles = box.triangles
# 2 SHAPE the triangles array into the shape accepted by quadpy
# trimesh -> (num_traingles, corners, xyz coordinates)
# quadpy ->  (corners, num_triangles, xyz_coords)
triangles_re = np.reshape(triangles,(
                       np.shape(triangles)[1],
                       np.shape(triangles)[0],
                       np.shape(triangles)[2])
                 )
# 3 DEFINE INTEGRATION FUNCTION 
def f(x):
    # -> get area of surface Integral[1, Area] = Area
    return np.ones(np.shape(x))

sol, error_estimate = quadpy.t2.integrate_adaptive(f, triangles_re,  1.0e-10)

问题

我收到以下错误：

assert all(areas > minimum_triangle_area)

断言错误

4

0 回答 0