[petsc-users] ILUTP in PETSc

[Barry Smith]

   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 <lu_qin_2000 at yahoo.com> 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 <xsli at lbl.gov>
> To: Barry Smith <bsmith at mcs.anl.gov> 
> Cc: Qin Lu <lu_qin_2000 at yahoo.com>; "petsc-users at mcs.anl.gov" <petsc-users at mcs.anl.gov> 
> 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 <bsmith at mcs.anl.gov> 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 <lu_qin_2000 at yahoo.com> 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.pdf that 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
>> 
>>