[Getfem-users] mumps and gmm

[mikael gueguen]

hi

what is the difference between MUMPS_solve and  
MUMPS_distributed_matrix_solve?
i want to use the mumps solver in parallel but i think i don't use  
correctly the routine.
in a first step, i want to assembly the matrix and vector on one  
process, and solve with mumps in parallel. Is it possible?

	i call the mumps solver like that (the assembly is made with the  
dense matrix Mwrite on one process) :
        gmm::csc_matrix< double >   Mvide;
        gmm::copy(Mwrite,Mcalcul);
        gmm::MUMPS_solve(Mcalcul,Ddl,SM);
        Mcalcul.swap(Mvide);

best regards