I did a quick test:
1) ilu with any shift type does not converge
2) gmres/4 bjacobi/ilu + shift_type nonzero does not converge
3) mpiexec -n 4 ./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor -sub_pc_type lu
-sub_pc_factor_shift_type nonzero
...
0 KSP Residual norm 3.284826093129e+03
1 KSP Residual norm
> On Oct 27, 2015, at 12:40 PM, Hong wrote:
>
> Here is the reason why it does not converge:
> ./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
> -ksp_converged_reason
> Linear solve did not converge due to DIVERGED_NANORINF iterations 0
This means it found a zero pivot eit
Here is the reason why it does not converge:
./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
-ksp_converged_reason
Linear solve did not converge due to DIVERGED_NANORINF iterations 0
KSP Object: 1 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gra
On Tue, Oct 27, 2015 at 12:36 PM, Hong wrote:
> Matt:
>
>> On Tue, Oct 27, 2015 at 11:13 AM, Hong wrote:
>>
>>> Gary :
>>> I tested your mat.bin using
>>> petsc/src/ksp/ksp/examples/tutorials/ex10.c
>>> ./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
>>> ...
>>> Mat Object:
Matt:
> On Tue, Oct 27, 2015 at 11:13 AM, Hong wrote:
>
>> Gary :
>> I tested your mat.bin using
>> petsc/src/ksp/ksp/examples/tutorials/ex10.c
>> ./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
>> ...
>> Mat Object: 1 MPI processes
>> type: seqaij
>> rows=588, cols=
On Tue, Oct 27, 2015 at 11:13 AM, Hong wrote:
> Gary :
> I tested your mat.bin using
> petsc/src/ksp/ksp/examples/tutorials/ex10.c
> ./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
> ...
> Mat Object: 1 MPI processes
> type: seqaij
> rows=588, cols=588
> total: n
Gary :
I tested your mat.bin using
petsc/src/ksp/ksp/examples/tutorials/ex10.c
./ex10 -f0 $D/mat.bin -rhs 0 -ksp_monitor_true_residual -ksp_view
...
Mat Object: 1 MPI processes
type: seqaij
rows=588, cols=588
total: nonzeros=11274, allocated nonzeros=11274
total number of malloc
On Tue, Oct 27, 2015 at 9:10 AM, Matthew Knepley wrote:
> On Tue, Oct 27, 2015 at 9:06 AM, Gary Rebt wrote:
>
>> Dear petsc-users,
>>
>> While using the FEniCS package to Solve a simple Stokes' flow problem, I
>> have run into problems with PETSc preconditioners. In particular, I would
>> like t
Dear petsc-users,
While using the FEniCS package to Solve a simple Stokes' flow problem, I have
run into problems with PETSc preconditioners. In particular, I would like to
use ILU (no parallel version) along with GMRES to solve my linear system but
the solver just hangs indefinitely at MatLUF
On Tue, Oct 27, 2015 at 9:06 AM, Gary Rebt wrote:
> Dear petsc-users,
>
> While using the FEniCS package to Solve a simple Stokes' flow problem, I
> have run into problems with PETSc preconditioners. In particular, I would
> like to use ILU (no parallel version) along with GMRES to solve my linea
Dear petsc-users,
While using the FEniCS package to Solve a simple Stokes' flow problem, I have run into problems with PETSc preconditioners. In particular, I would like to use ILU (no parallel version) along with GMRES to solve my linear system but the solver just hangs indefinitely at MatLUFa
11 matches
Mail list logo
