Hi,Are there any appropriate solutions to the following problems? why?
load “MUMPS_seq”
…
problem kg2(psai2,u,solver=sparsesolver,master=-1)=int2d(Th1)(psai2*u/(dt^2))-
…
0 : Allocation error for id%IS( 101830817 )
** ERROR RETURN ** FROM ZMUMPS INFO(1)= -13
** INFO(2)= 101830817
On return from ZMUMPS, INFOG(1)= -13
On return from ZMUMPS, INFOG(2)= 101830817
Erreur Mumps seq: number -13
Fatal Erreur MUMPS_seq Factorize
current line = 57
Exec error : MUMPS_seq Factorize
– number :0
** Instance Error 1 in ZMUMPS_F77 1
** MPI_ABORT called
or
load “UMFPACK64”
…
problem kg2(psai2,u,solver=“UMFPACK64”)=int2d(Th1)(psai2*u/(dt^2))-
…
Error umpfack umfpack_zl_numeric status -1
Error umfpack_di_solve status -3
– Solve :
min (-0.155751,-0.155755) max (0.155751,0.155755)
Hi, Does this compile right?
load “PETSc-complex”
…
problem kg2(psai2,u,sparams=“-pc_type gamg”)=int2d(Th1)(psai2u/(dt^2))-int2d(Th1)(2psai1u/(dt^2))+int2d(Th1)(psai0u/(dt^2))
…
Hi, professor, the programming language is as follows
…
varf rhs(psai1,u)=int2d(Th1)(2psai1u/(dt^2))-int2d(Th1)(lambda*(psai1conj(psai1))psai1u);
…
but,When it’s running
varf rhs(psai1,u)=int2d(Th1)(2psai1u/(dt^2))-int2d(Th1)(lambda(psai1*conj(psai1)) error operator * <10LinearCombI7MGauche4C_F0E>, <10LinearCombI7MGauche4C_F0E>
why?