Thank you for your answers. Barry's answer allowed me to get up to 50x
faster code, so it was a huge help. I decided not to use ksp context
altogether for this small operation.
Thank you again,
Ali
On 11-03-2018 18:53, Matthew Knepley wrote:
On Sun, Mar 11, 2018 at 1:14 AM, Smith, Barry F.
On Sun, Mar 11, 2018 at 1:14 AM, Smith, Barry F. wrote:
>
> 1) Run the problem with -ksp_view_mat and -ksp_view_rhs and mail
> petsc-ma...@mcs.anl.gov the resulting file produced called binaryoutput
>
>2) By default PCLU does a reordering to reduce fill that could
> introduce a zero pivoit
1) Run the problem with -ksp_view_mat and -ksp_view_rhs and mail
petsc-ma...@mcs.anl.gov the resulting file produced called binaryoutput
2) By default PCLU does a reordering to reduce fill that could introduce a
zero pivoit, PCILU does not do a reordering by default. You can use
-pc_fact