On Thu, Jan 21, 2010 at 2:07 PM, Jack Lee <span dir="ltr"><<a href="mailto:jack.lee@sjc.ox.ac.uk">jack.lee@sjc.ox.ac.uk</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Thanks for the reply. Perhaps I should have explained a little further:<br>
<br>
> 1) MUMPS is a linear solver<br>
I am aware of that, what I meant was that instead of the default iterative solver which snes employs, I'd like the option of being able to use an external solver for each linear solution of the non-linear loop.<br></blockquote>
<div><br>Okay<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
> 2)<br>
> <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html</a><br>
Currently I am applying these steps:<br>
1. set ksp context extracted from snes to pre only<br>
2. set pc type to lu<br>
3. set pc factor solver package to MUMPS, as you've explained above.<br>
However, having no examples to compare to, I'm not 100% sure if this is the way it's supposed to be done.<br></blockquote><div><br>Yep, that is it<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
> 3) There is a manual section on this<br>
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.<br>
<br>
> 4) You can use command line options in Fortran<br>
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....<br></blockquote><div><br>Well, when you call PetscInitialize(), in Fortran, it is reading the command line arguments. You, of course, do not have<br>
to use them, but they are there.<br><br> Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Cheers,<br>
Jack<br>
<br>
<br>
<br>
<a href="mailto:petsc-users-request@mcs.anl.gov" target="_blank">petsc-users-request@mcs.anl.gov</a> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Send petsc-users mailing list submissions to<br>
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/petsc-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/petsc-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:petsc-users-request@mcs.anl.gov" target="_blank">petsc-users-request@mcs.anl.gov</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:petsc-users-owner@mcs.anl.gov" target="_blank">petsc-users-owner@mcs.anl.gov</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of petsc-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. using external direct linear solvers (Jack Lee)<br>
2. Re: using external direct linear solvers (Matthew Knepley)<br>
3. using PETSC with visual studio 2005? (Zhifeng Sheng - EWI)<br>
4. Re: using PETSC with visual studio 2005? (Satish Balay)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 20 Jan 2010 22:50:47 +0000<br>
From: Jack Lee <<a href="mailto:jack.lee@sjc.ox.ac.uk" target="_blank">jack.lee@sjc.ox.ac.uk</a>><br>
Subject: [petsc-users] using external direct linear solvers<br>
To: <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
Message-ID: <<a href="mailto:4B5788C7.20902@sjc.ox.ac.uk" target="_blank">4B5788C7.20902@sjc.ox.ac.uk</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Hi there,<br>
<br>
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.<br>
<br>
Jack<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 20 Jan 2010 17:14:16 -0600<br>
From: Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
Subject: Re: [petsc-users] using external direct linear solvers<br>
To: PETSc users list <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
Message-ID:<br>
<<a href="mailto:a9f269831001201514x473b849ci76cc9ab80d754586@mail.gmail.com" target="_blank">a9f269831001201514x473b849ci76cc9ab80d754586@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
On Wed, Jan 20, 2010 at 4:50 PM, Jack Lee <<a href="mailto:jack.lee@sjc.ox.ac.uk" target="_blank">jack.lee@sjc.ox.ac.uk</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi there,<br>
<br>
Apologies if this is an obvious question, but is there an example showing<br>
how to call an external solver? I am trying to get the snes solver to use<br>
MUMPS. As the calls are made from a fortran application, I cannot use the<br>
commandline options.<br>
<br>
</blockquote>
<br>
1) MUMPS is a linear solver<br>
<br>
2)<br>
<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCFactorSetMatSolverPackage.html</a><br>
<br>
3) There is a manual section on this<br>
<br>
4) You can use command line options in Fortran<br>
<br>
Matt<br>
<br>
Jack<br>
</blockquote>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>