在使用直接分解方法求解稀疏线性系统时，使用的排序策略会显着影响因子中非零元素的填充因子。一种这样的排序策略是嵌套剖析。我想知道是否可以在仅给定网格参数的情况下提前提出嵌套解剖排序（假设具有一阶差分的 M x N 平方有限差分网格）。
编辑 我刚刚发现有这样的代码： http ://www.cise.ufl.edu/research/sparse/meshnd/
常规网格上的嵌套解剖
计算科学
线性代数
2021-12-01 08:09:02
1个回答
是的。我最近编写了代码来做到这一点。
假设你有一个网格，并且可以接受具有 100 个顶点的叶节点。然后可以定义一个递归函数，其中参数是：
- 矩形子域的尺寸和偏移量
- 指向将存储重新排序的数组的指针
然后，该例程只需计算局部维度的乘积，以确定域是否小到可以接受的小叶，然后，如果是，则编写叶节点自然索引（例如为格），否则，切割最大的子域维度，在左右两块上递归，然后写分隔符自然索引。
其它你可能感兴趣的问题