Any examples of how to set Spooles LU and CHOLESKY direct solvers using the Fortran API?

[Barry Smith]

    Stephen,

1) You no longer need to set particular Spooles matrix types. Just use  
AIJ or SBAIJ (for symmetric case)

2) call KSPGetPC(ksp,pc,ierr)
     Call PCSetType(pc,PCLU,ierr)
      call  PCFactorSetMatSolverPackage(pc,MAT_SOLVER_SPOOLES,ierr)
      call KSPSolve()

    Barry



On Jul 24, 2009, at 5:22 AM, Stephen Ball wrote:

> Hi
>
> I have recently moved from using PETSc v2.3.3 to v3.0.0 and am  
> trying to
> update my Fortran code accordingly.
>
> Do you have any examples of how to set Spooles LU and CHOLESKY direct
> solvers using the Fortran API?
>
> I am struggling somewhat to understand the correct sequence of calls  
> for
> your new API, including the matrix and PC creation and set up stages
> when using Spooles LU and CHOLESKY direct solvers.
>
> What calls are required or optional and in what circumstances?
>
> Kindest regards
>
> Stephen
>
> This e-mail and any attachments may contain confidential and
> privileged information. If you are not the intended recipient,
> please notify the sender immediately by return e-mail, delete this
> e-mail and destroy any copies. Any dissemination or use of this
> information by a person other than the intended recipient is
> unauthorized and may be illegal.