# MUMPS or UMFPACK64
**URL:** https://community.freefem.org/t/mumps-or-umfpack64/1437
**Category:** General Discussion
**Created:** [January 16, 2022, 7:19am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437 "2022-01-16T07:19:18Z")
**Posts on this page:** 17
**Page:** 1
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [January 16, 2022, 7:19am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/1 "2022-01-16T07:19:18Z")
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)
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [January 16, 2022, 7:43am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/2 "2022-01-16T07:43:06Z")
The problem is too big to be factorized by a direct solver, switch to an iterative method.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [January 16, 2022, 7:50am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/3 "2022-01-16T07:50:37Z")
Thank you, I haven’t used iterative methods. Is there an example that works?
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [January 16, 2022, 7:52am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/4 "2022-01-16T07:52:50Z")
You can use PETSc with the parameter `-pc_type gamg`, for example.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [January 16, 2022, 8:23am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/5 "2022-01-16T08:23:15Z")
> [@prj](#):
>
> -pc\_type gamg
Hi, Does this compile right?
load “PETSc-complex”
…
problem kg2(psai2,u,sparams=“-pc\_type gamg”)=int2d(Th1)(psai2_u/(dt^2))-int2d(Th1)(2_psai1_u/(dt^2))+int2d(Th1)(psai0_u/(dt^2))
…
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [January 16, 2022, 8:43am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/6 "2022-01-16T08:43:01Z")
No, you need to use a `varf` instead of a `problem`.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [January 17, 2022, 11:09am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/7 "2022-01-17T11:09:51Z")
Hi, professor, the programming language is as follows
…
varf rhs(psai1,u)=int2d(Th1)(2_psai1_u/(dt^2))-int2d(Th1)(lambda\*(psai1_conj(psai1))psai1u);
…
but,When it’s running
varf rhs(psai1,u)=int2d(Th1)(2_psai1_u/(dt^2))-int2d(Th1)(lambda_(psai1\*conj(psai1)) error operator \* \<10LinearCombI7MGauche4C\_F0E\>, \<10LinearCombI7MGauche4C\_F0E\>
why?
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [January 17, 2022, 1:01pm UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/8 "2022-01-17T13:01:39Z")
The term `psai1*conj(psai1)` is not linear, so you cannot put that in a `varf` (or in a `problem` for that matter).
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [January 17, 2022, 1:23pm UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/9 "2022-01-17T13:23:46Z")
I see. Thank you very much for your reply to my question in your busy schedule.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [February 19, 2022, 8:27am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/10 "2022-02-19T08:27:44Z")
> [@Qianrui](#):
>
> load “PETSc-complex”
Hi,professor
load “PETSc-complex”
…
matrix A = va(Xhf, Xhf,sparams = “-pc\_type asm -sub\_pc\_type lu -ksp\_converged\_reason”);
…
The operation procedure is very slow. What solution can we do?
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [February 19, 2022, 8:42am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/11 "2022-02-19T08:42:52Z")
You need to use a `Mat`, not a `matrix`.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [February 20, 2022, 8:01am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/12 "2022-02-20T08:01:45Z")
Pestc can only be applied to FreeFem + + - mpi.exe version? Because FreeFem + + - mpi.exe has been debugging failed, FreeFem+ +. exe apply?
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [February 20, 2022, 8:19am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/13 "2022-02-20T08:19:07Z")
No, you need to use `FreeFem++-mpi.exe` to use PETSc/SLEPc.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [February 20, 2022, 8:25am UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/14 "2022-02-20T08:25:50Z")
Microsoft MPI was successfully installed on windows10, but examples, using freefem + + -mpi.exe failed to run. Why?
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [February 21, 2022, 1:40pm UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/15 "2022-02-21T13:40:34Z")
Hello, PETSc needs to use freefem + + - mpi .exe to run, I need to download PETSc separately, or freefem + + - mpi. exe can be called directly?
---
### Author:  [@prj](https://community.freefem.org/u/prj)
#### Post date: [February 21, 2022, 1:49pm UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/16 "2022-02-21T13:49:52Z")
It can be called directly. Please look at this screencast [Pierre Jolivet introduction to parallel FreeFEM part 1 - YouTube](https://www.youtube.com/watch?v=-Aw2O46V2bo), it answers most of your recent questions.
---
### Author:  [@Qianrui](https://community.freefem.org/u/Qianrui)
#### Post date: [February 21, 2022, 1:53pm UTC](https://community.freefem.org/t/mumps-or-umfpack64/1437/17 "2022-02-21T13:53:39Z")
Hello, can you help me solve this problem？
[https://community.freefem.org/t/freefem-mpi-exe/1539/8](https://community.freefem.org/t/freefem-mpi-exe/1539/8)