在FIPY中无法再现球面网格的结果

Question

我正在尝试从this post中复制一个结果(考虑到jeguyer在答复中提供的修复)。但是，当我试图运行完全相同的代码时，我会收到以下警告：

/home/dmitry/~./conda/envs/FIPY/lib/python3.9/site-packages/fipy/viewers/matplotlibViewer/matplotlibViewer.py:195: UserWarning: Matplotlib is currently using module://matplotlib_inline.backend_inline, which is a non-GUI backend, so cannot show the figure.
  self.fig.show()
/home/dmitry/~./conda/envs/FIPY/lib/python3.9/site-packages/fipy/variables/arithmeticCellToFaceVariable.py:32: RuntimeWarning: invalid value encountered in subtract
  return (cell2 - cell1) * alpha + cell1
/home/dmitry/~./conda/envs/FIPY/lib/python3.9/site-packages/fipy/variables/faceGradVariable.py:124: RuntimeWarning: invalid value encountered in subtract
  N = (N2 - numerix.take(self.var, id1, axis=-1)) / dAP

由此产生的X_ca数组都由-inf组成。有什么方法可以修复它吗？提前谢谢你。

Answer 1

正如在my answer中所指出的，您所链接的问题：

有些求解者不喜欢球形网格，这可能是因为细胞体积的变化范围很大。SciPy LinearLUSolver似乎起作用了。明智的预处理可能会帮助其他的解决者。

SciPy LU解决。其他求解者，例如，PETSc LU没有。在Windows上解决问题时，默认情况下会得到SciPy解决程序。Ubuntu和macOS默认为PETSc。

我不知道要让其他的解决者去工作需要些什么。