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 at 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
