[petsc-users] petsc-users Digest, Vol 13, Issue 28

[Hong Zhang]

Jack,
>
> > 2)
> > http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html
> Currently I am applying these steps:
> 1. set ksp context extracted from snes to pre only
> 2. set pc type to lu
> 3. set pc factor solver package to MUMPS, as you've explained above.
> However, having no examples to compare to, I'm not 100% sure if this  
> is the way it's supposed to be done.

I just wrote a Fortran example
petsc-dev/src/snes/examples/tests/ex14f.F few days ago which
uses MUMPS as a linear solver within snes non-linear loop.
The example is attached, or can be obtained through petsc-dev.

petsc-dev/src/snes/examples/tests/ex12f.F is a simpler example
which you can run it with options
./ex12f -pc_type lu -pc_factor_mat_solver_package mumps
to invoke mumps' LU. Use '-snes_view' to check which solvers are used
and '-ksp_monitor -snes_monitor' for monitoring iterations.

Hong

>
> > 3) There is a manual section on this
> if you're talking about page 83, it doesn't spell out the case which  
> I'm trying to apply to, so thought I'd better double-check.
>
> > 4) You can use command line options in Fortran
> Not in the custom fortran application we're building, which is  
> compiled as a library. We're linking in petsc to exploit its  
> fantastic solver capabilities....
>
> Cheers,
> Jack
>
>
>
> petsc-users-request at mcs.anl.gov wrote:
>> Send petsc-users mailing list submissions to
>> 	petsc-users at mcs.anl.gov
>> To subscribe or unsubscribe via the World Wide Web, visit
>> 	https://lists.mcs.anl.gov/mailman/listinfo/petsc-users
>> or, via email, send a message with subject or body 'help' to
>> 	petsc-users-request at mcs.anl.gov
>> You can reach the person managing the list at
>> 	petsc-users-owner at mcs.anl.gov
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of petsc-users digest..."
>> Today's Topics:
>>   1.  using external direct linear solvers (Jack Lee)
>>   2. Re:  using external direct linear solvers (Matthew Knepley)
>>   3.  using PETSC with visual studio 2005? (Zhifeng Sheng - EWI)
>>   4. Re:  using PETSC with visual studio 2005? (Satish Balay)
>> ----------------------------------------------------------------------
>> Message: 1
>> Date: Wed, 20 Jan 2010 22:50:47 +0000
>> From: Jack Lee <jack.lee at sjc.ox.ac.uk>
>> Subject: [petsc-users] using external direct linear solvers
>> To: petsc-users at mcs.anl.gov
>> Message-ID: <4B5788C7.20902 at sjc.ox.ac.uk>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>> Hi there,
>> Apologies if this is an obvious question, but is there an example  
>> showing how to call an external solver? I am trying to get the snes  
>> solver to use MUMPS. As the calls are made from a fortran  
>> application, I cannot use the commandline options.
>> Jack
>> ------------------------------
>> Message: 2
>> Date: Wed, 20 Jan 2010 17:14:16 -0600
>> From: Matthew Knepley <knepley at gmail.com>
>> Subject: Re: [petsc-users] using external direct linear solvers
>> To: PETSc users list <petsc-users at mcs.anl.gov>
>> Message-ID:
>> 	<a9f269831001201514x473b849ci76cc9ab80d754586 at mail.gmail.com>
>> Content-Type: text/plain; charset="iso-8859-1"
>> On Wed, Jan 20, 2010 at 4:50 PM, Jack Lee <jack.lee at sjc.ox.ac.uk>  
>> wrote:
>>> Hi there,
>>>
>>> Apologies if this is an obvious question, but is there an example  
>>> showing
>>> how to call an external solver? I am trying to get the snes solver  
>>> to use
>>> MUMPS. As the calls are made from a fortran application, I cannot  
>>> use the
>>> commandline options.
>>>
>> 1) MUMPS is a linear solver
>> 2)
>> http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html
>> 3) There is a manual section on this
>> 4) You can use command line options in Fortran
>>  Matt
>> Jack

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100121/1906c039/attachment-0002.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex14f.F
Type: application/octet-stream
Size: 12434 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100121/1906c039/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100121/1906c039/attachment-0003.htm>