I tried to use command line options as the example suggested ('-ksp_type
preonly -pc_type ilu -pc_factor_mat_solver_package superlu
-mat_superlu_ilu_droptol 1.e-8') without changing my source code, but then the
call to KSPSetUp returned error number 56.
Does this mean I still need to change the source code (such as adding calls to
PCFactorSetMatSolverPackage, PCFactorGetMatrix, etc.)in addition to the command
line options?
I ask this since the use of SuperLU seems to be different from using Hypre,
which can be invoked with command line options without changing source code.
Thanks a lot,
Qin
----- Original Message -----
From: Barry Smith <[email protected]>
To: Qin Lu <[email protected]>
Cc: Xiaoye S. Li <[email protected]>; "[email protected]"
<[email protected]>
Sent: Monday, May 12, 2014 5:11 PM
Subject: Re: [petsc-users] ILUTP in PETSc
See for example:
http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MATSOLVERSUPERLU.html
On May 12, 2014, at 4:54 PM, Qin Lu <[email protected]> wrote:
> Hello,
>
> I have built PETSc with SuperLU, but what are PETSc's command line options to
> invoke SuperLU's ILUTP preconditioner and to set the dropping tolerance?
> (-mat_superlu_ilu_droptol for the latter?)
>
> Do I need to do some programming in order to call SuperLU's preconditioner,
> or the command line options would work?
>
> Many thanks,
> Qin
>
>
> From: Xiaoye S. Li <[email protected]>
> To: Barry Smith <[email protected]>
> Cc: Qin Lu <[email protected]>; "[email protected]"
> <[email protected]>
> Sent: Friday, May 2, 2014 3:40 PM
> Subject: Re: [petsc-users] ILUTP in PETSc
>
>
>
> The sequential SuperLU has ILUTP implementation, not in parallel versions.
> PETSc already supports the option of using SuperLU, so you should be able to
> try easily.
>
> In SuperLU distribution:
>
> EXAMPLE/zitersol.c : an example to use GMRES with ILUTP preconditioner
>(returned from driver SRC/zgsisx.c)
>
> SRC/zgsitrf.c : the actual ILUTP factorization routine
>
>
> Sherry Li
>
>
>
> On Fri, May 2, 2014 at 12:25 PM, Barry Smith <[email protected]> wrote:
>
>
>> At http://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html there
>> are two listed. ./configure —download-hypre
>>
>> mpiexec -n 23 ./yourprogram -pc_type hypre -pc_hypre_type ilupt or euclid
>>
>> you can also add -help to see what options are available.
>>
>> Both pretty much suck and I can’t image much reason for using them.
>>
>> Barry
>>
>>
>>
>> On May 2, 2014, at 10:27 AM, Qin Lu <[email protected]> wrote:
>>
>>> Hello,
>>>
>>> I am interested in using ILUTP preconditioner with PETSc linear solver.
>>> There is an online doc
>>> https://fs.hlrs.de/projects/par/par_prog_ws/pdf/petsc_nersc01_short.pdfthat
>>> mentioned it is available in PETSc with other packages (page 62-63). Is
>>> there any instructions or examples on how to use it?
>>>
>>> Many thanks,
>>> Qin
>>
>>
