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 <[email protected]> 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://[email protected]//hg/petsc/petsc-dev-new-solvers

and trying it out before then. Please report problem to [email protected] 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