0

我正在使用 petsc4py，但遇到了我不明白的异常。我定义了以下函数：

def tsIJacobian(self, ts, t, u, udot, shift, A, B):
    self.setup_problem()
    psol = fe.as_backend_type(self.sol.vector()).vec()
    pA = fe.as_backend_type(self.A).mat()
    u.copy(psol)
    JU = fe.assemble(fe.derivative(ksdg.U_terms, ksdg.U))
    Jrho = fe.assemble(fe.derivative(ksdg.rho_terms, ksdg.rho))
    pJU = fe.as_backend_type(JU).mat()
    pJrho = fe.as_backend_type(Jrho).mat()
    pA.copy(A)
    A.scale(shift)
    A.axpy(1.0, pJU)
    A.axpy(1.0, pJrho)
    A.assemble()
    if not (A is B):
        A.copy(B)
        B.assemble()

然后尝试以下。（pA 已经在别处定义为 48x48 PETSc.Mat，并已组装。ksdg 是我正在研究的一个类的实例，如果我能让它工作，它最终将成为它的成员函数）：

J = pA.duplicate()
B = pA.duplicate()
tsIJacobian(ksdg, ts, 0, psol, pdsol, 0.1, J, B)

这会引发以下异常：

---------------------------------------------------------------------------
Error                                     Traceback (most recent call last)
<ipython-input-24-2a71f7ccf0af> in <module>()
----> 1 tsIJacobian(ksdg, ts, 0, psol, pdsol, 0.1, J, B)

<ipython-input-22-14579c08d6ae> in tsIJacobian(self, ts, t, u, udot, shift, A, B)
     15     A.assemble()
     16     if not (A is B):
---> 17         A.copy(B)
     18         B.assemble()

PETSc/Mat.pyx in petsc4py.PETSc.Mat.copy (src/petsc4py.PETSc.c:118071)()

Error: error code 63

看着petscerror.h.html，63是PETSC_ERR_ARG_OUTOFRANGE 63 /* input argument, out of range */。

如果有人理解为什么 PETSc 不允许我将矩阵 A 复制到 B，我将不胜感激。谢谢。

4

1 回答 1

0

好的，我有点想通了。当尝试在稀疏矩阵中创建新的非零元素时引发错误 63 异常，如果选项

PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR

是True。因此，如果我执行

B.setOption(PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR, False)

在创建 之后B，在 tsIJacobian 中不会立即引发异常。我不明白为什么我没有得到相同的例外J。

于 2017-07-01T12:40:57.373 回答