[petsc-users] PaStiX does not work

[Gong Ding]

Dear Barry
The mainly petsc related file is attached.
The whole project is very large, it is a 3D/parallel code for semiconductor simulation.
You can download the source code from http://www.cogenda.com/downloads/category/7-genius-open-source-edition.html
It can be compiled with petsc 3.0. I am now adding 3.1 support.

I had use MUMPS solver for quite a long time by the same code: 
PCFactorSetMatSolverPackage (pc, MAT_SOLVER_MUMPS);
It works well for most of the time. However, it may crash (when other solver like KSP type works).  

I'd like to try PASTIX if it is faster and more stable. 

Yours
Gong Ding


> 
> On Apr 6, 2010, at 7:04 PM, Gong Ding wrote:
> 
>> Dear Petsc developer,
>>
>> I found a problem that PASTIX solver can not be loaded from PETSC,  
>> even for 3.1 version.
>>
>> When I use following code, the ksp_view told me that PC is ilu. And  
>> the snes solver can not convergence at all.
>>
>> ierr = KSPSetType (ksp, (char*) KSPPREONLY); assert(!ierr);//it  
>> seems PaStiX don't work with KSPPREONLY
>> ierr = PCSetType  (pc, (char*) PCLU);
>> ierr = PCFactorSetMatSolverPackage (pc, MAT_SOLVER_PASTIX); assert(! 
>> ierr);
>>
>> Is it a bug or I made something wrong?
> 
>   You are doing something wrong but one cannot tell from just this  
> code fragment.  Send the entire code to petsc-maint at mcs.anl.gov and  
> the options you used to run it and the output it produced.
> 
>    Barry
> 
>>
>> Yours
>> Gong Ding
>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: fvm_nonlinear_solver.cc
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100407/c47cb2ce/attachment-0001.diff>