New approach for external direct solvers in PETSc

Lisandro Dalcin dalcinl at
Fri Jul 25 17:47:56 CDT 2008

Perhaps yor are missing in petscmat.h a line like the one below?? In
the docsyou mention it, but it is not listed alongside the other

#define MAT_SOLVER_PETSC      "petsc"

BTW, Have you ever considered the posibility of not using the 'Mat'
class for factored matrices?? I believe that normal matrices and
factored matrices share very few in common as to being instances of
the same class. Looking at the source code in src/mat, there are many
checks like 'if (mat-->factor)' or  'if (!mat->factor)' that seems to
support my concern. What do you think?

On Fri, Jul 25, 2008 at 7:02 PM, Barry Smith <bsmith at> wrote:
>   petsc-dev users,
>     I have finally completed my changes to PETSc for a new approach to
> accessing external direct solvers
> in PETSc like Spooles, MUMPS etc. I will be pushing it to petsc-dev sometime
> the middle of next week.
> If you are using the direct solvers you might consider cloning from
>  or
> ssh://
> and trying it out before then. Please report problem to
> petsc-maint at as you hit them.
>   Barry
> From the changes file
> The "parallel direct solver" matrix types like
> MATAIJSPOOLES are ALL gone. Now you use -pc_factor_mat_solver_package
> spooles etc or PCFactorSetMatSolverPackage() or if working directly with
> matrices, MatGetFactor(A,MATSPOOLES,...)

Lisandro Dalcín
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594

More information about the petsc-dev mailing list