作文网站哪个平台好,广东省建设集团有限公司,石家庄网络seo推广,sae wordpress 主题最近在用fenics框架跑有限元代码#xff0c;其中有一个部分是把在矩阵里定义的初始值#xff0c;赋值到有限元空间里#xff0c;这就涉及到了初始矩阵和有限元空间坐标的转化#xff0c;部分代码如下
for i in range(len(dof_coordinates)):# x, y dof_coordinates[i…最近在用fenics框架跑有限元代码其中有一个部分是把在矩阵里定义的初始值赋值到有限元空间里这就涉及到了初始矩阵和有限元空间坐标的转化部分代码如下
for i in range(len(dof_coordinates)):# x, y dof_coordinates[i]#原本的区间是0到physical_length的,根据num_grids_per_axis长度进行缩放.# ix, iy min(int(x / physical_length * (num_grids_per_axis - 1)), num_grids_per_axis - 1), min(int(y / physical_length * (num_grids_per_axis - 1)), num_grids_per_axis - 1)ix, iy min(int(x * (num_grids_per_axis - 1) / physical_length), num_grids_per_axis - 1), min(int(y * (num_grids_per_axis - 1) / physical_length), num_grids_per_axis - 1)
这里有一点需要注意
x要先和num_grids_per_axis相乘然后再除以physical_length。先除再乘的话数学逻辑上一样但是如果一个很小的数除以一个很大的数会有舍入误差。所以先乘再除。
