PETSc 矩阵副本引发异常:为什么?
PETSc Matrix copy raises exception: why?
我正在使用 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已经在别处定义为a 48x48 PETSc.Mat,并组装。ksdg是我正在处理的a class的一个实例,如果我能得到,它最终将成为一个成员函数它可以工作):
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,我将不胜感激。谢谢。
好的,我有点明白了。如果选项
试图在稀疏矩阵中创建一个新的非零元素,则会引发错误 63 异常
PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR
是True
。因此,如果我执行
B.setOption(PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR, False)
创建 B
后,tsIJacobian 中没有出现异常。我不明白的是为什么我没有得到 J
.
的相同异常
我正在使用 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已经在别处定义为a 48x48 PETSc.Mat,并组装。ksdg是我正在处理的a class的一个实例,如果我能得到,它最终将成为一个成员函数它可以工作):
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,我将不胜感激。谢谢。
好的,我有点明白了。如果选项
试图在稀疏矩阵中创建一个新的非零元素,则会引发错误 63 异常PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR
是True
。因此,如果我执行
B.setOption(PETSc.Mat.Option.NEW_NONZERO_ALLOCATION_ERR, False)
创建 B
后,tsIJacobian 中没有出现异常。我不明白的是为什么我没有得到 J
.