我正在使用有限元、Newton Raphson 技术和 MUMPS 直接求解器来求解$P$在这个等式中：

$$\frac{\delta K}{\delta x} \frac{\delta P}{\delta x}= \frac{\delta K B}{\delta x} + A$$

$A$和$B$是常数参数。$K$是随空间变化的参数。$P$是解变量。

$K$分布如下：

因此，$K$变为零！所以我的矩阵的对角线上会有零，求解器最终会打印出矩阵是奇异的。

是否有任何技术仍然可以解决这个问题？