New approach for external direct solvers in PETSc

[Lisandro Dalcin]

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
packages.

#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 mcs.anl.gov> 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
>
> http://petsc.cs.iit.edu/hg/petsc/petsc-dev-new-solvers  or
>
> ssh://petsc@petsc.cs.iit.edu//hg/petsc/petsc-dev-new-solvers
>
> and trying it out before then. Please report problem to
> petsc-maint at mcs.anl.gov 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