On Thu, Jan 21, 2010 at 2:07 PM, Jack Lee <jack.lee at sjc.ox.ac.uk> wrote:
> Thanks for the reply. Perhaps I should have explained a little further: > > > 1) MUMPS is a linear solver > 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. > Okay > > 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. > Yep, that is it > > 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.... > Well, when you call PetscInitialize(), in Fortran, it is reading the command line arguments. You, of course, do not have to use them, but they are there. Matt > 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 >> > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100121/01089ea5/attachment.htm>
