[petsc-users] Recommendations of possible preconditioner/ksp solver/snes solver for three phase flow reservoir simulator
Chukwudi Chukwudozie
chdozie at gmail.com
Sun Dec 16 21:00:59 CST 2018
Barry,
Thanks for the reply. I have attached two files with the information.
On Sat, Dec 15, 2018 at 8:58 PM Smith, Barry F. <bsmith at mcs.anl.gov> wrote:
>
> Chuks,
>
> 1) Please run a successful case with -snes_view and send the output so
> we can see the exact solver options being used.
>
> 2) Run a less successful case with -snes_monitor
> -ksp_monitor_true_residual -ksp_converged_reason
> -fieldsplit_ksp_monitor_true_residual -fieldsplit_ksp_converged reason so
> we can see how the linear solver (and the sub linear solves) are (not)
> converging.
>
>
> Barry
>
>
> > On Dec 15, 2018, at 7:54 PM, Chukwudi Chukwudozie via petsc-users <
> petsc-users at mcs.anl.gov> wrote:
> >
> > Hi,
> >
> > I am developing a fully implicit control volume finite element black oil
> reservoir simulator using petscdmplex and snes. So far, for 2 phase flow
> problems (3 by 3 block for pressure, saturation and bottomhole pressure),
> my simulator runs fine for small time step sizes in the order of less than
> 5 days. However, increasing the time step size, it doesn't converge with -3
> (-5 sometimes) diverge reason. Considering that this is a fully implicit
> but nonlinear problem, I expect unconditional stability for reasonable time
> step sizes.
> >
> > The background DM for the 3 block problem is a composite one with the
> first two being DMPLES's (bag->plexScalNode for pressure and saturation)
> and the last one, bag->WellRedun, a redundant DM for the wellbore
> pressures. See below.
> >
> > ierr =
> DMRedundantCreate(PETSC_COMM_WORLD,0,bag->numWells,&bag->WellRedun);CHKERRQ(ierr);
> > ierr =
> DMCompositeCreate(PETSC_COMM_WORLD,&bag->MultiPhasePacker);CHKERRQ(ierr);
> > ierr =
> DMCompositeAddDM(bag->MultiPhasePacker,bag->plexScalNode);CHKERRQ(ierr);
> > ierr =
> DMCompositeAddDM(bag->MultiPhasePacker,bag->plexScalNode);CHKERRQ(ierr);
> > ierr =
> DMCompositeAddDM(bag->MultiPhasePacker,bag->WellRedun);CHKERRQ(ierr);
> >
> > I have used nested matrices as shown below. In addition, I used the
> fieldsplit (composite multiplicative) preconditioner with FGRES and least
> square for newton solver. See below.
> >
> > ierr =
> MatCreateNest(PETSC_COMM_WORLD,Msize,is,Msize,is,&bK[0],K);CHKERRQ(ierr);
> >
> >
> > ierr = SNESGetKSP(bag->snesP,&kspP);CHKERRQ(ierr);
> > ierr =
> KSPSetTolerances(kspP,1.e-8,1.e-8,PETSC_DEFAULT,PETSC_DEFAULT);CHKERRQ(ierr);
> > ierr = KSPSetType(kspP,KSPFGMRES);CHKERRQ(ierr);
> > ierr = KSPSetFromOptions(kspP);CHKERRQ(ierr);
> > ierr = KSPGetPC(kspP,&pcP);CHKERRQ(ierr);
> > ierr = PCSetType(pcP, PCFIELDSPLIT);CHKERRQ(ierr);
> > ierr =
> DMCompositeGetGlobalISs(bag->MultiPhasePacker,&is);CHKERRQ(ierr);
> > for(i = 0; i < Msize; i++) ierr =
> PCFieldSplitSetIS(pcP,NULL,is[i]);CHKERRQ(ierr);
> > ierr =
> PCFieldSplitSetType(pcP,PC_COMPOSITE_MULTIPLICATIVE);CHKERRQ(ierr);
> >
> > Any ideas on possible precoditioner options/solver options/combinations
> to improve stability will be appreciated.
> >
> > Chuks
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20181216/27c4696e/attachment-0001.html>
-------------- next part --------------
FracD:
Numerical implementation of multiphysics models.
(c) 2016-2017 Chukwudi Chukwudozie chdozie at gmail.com
Well object "well":
top: 3.000000e+03 3.000000e+03
bottom: 3.000000e+03 3.000000e+03
wellnodes: 1.500000e+03 1.500000e+03
well water rate: 1.000000e+02
well oil rate: 0.000000e+00
well gas rate: 0.000000e+00
total well rate: 0.000000e+00
well bottomhole pressure: 5.000000e+03
well height/production interval: 3.000000e+01
well radius: 2.920000e-01
well block equivalent radius: 0.000000e+00
well skin: 0.000000e+00
Well Type: "INJECTORWATER" well under "WATERRATE" condition
Well object "well":
top: 3.000000e+03 3.000000e+03
bottom: 3.000000e+03 3.000000e+03
wellnodes: 4.500000e+03 4.500000e+03
well water rate: -0.000000e+00
well oil rate: -5.000000e+02
well gas rate: -0.000000e+00
total well rate: -0.000000e+00
well bottomhole pressure: 2.000000e+03
well height/production interval: 3.000000e+01
well radius: 2.920000e-01
well block equivalent radius: 0.000000e+00
well skin: 0.000000e+00
Well Type: "PRODUCER" well under "PRESSURE" condition
Option table:
#PETSc Option Table entries:
-Cr 4e-6
-PcowInrelPermData
-RsBubblepointfilename solutiongas.txt
-bubblepointfilename bubblepoint2phase.txt
-dim 2
-elementtype QUADRILATERAL
-fluid OILWATER
-initialSw 0.2
-initialpressure 4000
-maxtimestep 10
-meshfilename 2phasemesh_quad.msh
-nw 2
-oil_FVFmodel INTERPOLATION
-oil_PVTfilename oil_pvt_1.txt
-oil_compress 1e-5
-oil_denscoeff 0.
-oil_densmodel INTERPOLATION
-oil_prefdens 0.
-oil_preffvf -3000
-oil_prefvisc -3000
-oil_refFVF 1.11
-oil_refdens 1.0
-oil_refvisc 1.2
-oil_visccoeff 4e-6
-oil_viscmodel INTERPOLATION
-oilwaterrelpermfilename owrelperm1.txt
-p_snes_view
-perm 300,300
-phi 0.3
-pref_phi 4000
-refpressure 4000
-timevalue 5
-units FIELDUNITS
-verbose
-w0_Pbh 5000
-w0_Qws 100
-w0_bottom 3000,3000
-w0_constraint waterrate
-w0_coords 1500,1500
-w0_h 30
-w0_rw 0.292
-w0_top 3000,3000
-w0_type injectorwater
-w1_Pbh 2000
-w1_Qos 500
-w1_bottom 3000,3000
-w1_constraint pressure
-w1_coords 4500,4500
-w1_h 30
-w1_rw 0.292
-w1_top 3000,3000
-w1_type producer
-w2_Pbh 2000
-w2_Qos 500
-w2_bottom 3000,3000
-w2_constraint pressure
-w2_coords 2250,3750
-w2_h 30
-w2_rw 0.292
-w2_top 3000,3000
-w2_type producer
-w3_Pbh 2000
-w3_Qos 500
-w3_bottom 3000,3000
-w3_constraint pressure
-w3_coords 3750,2250
-w3_h 30
-w3_rw 0.292
-w3_top 3000,3000
-w3_type producer
-w4_Pbh 2000
-w4_Qos 500
-w4_bottom 3000,3000
-w4_constraint pressure
-w4_coords 3750,3750
-w4_h 30
-w4_rw 0.292
-w4_top 3000,3000
-w4_type producer
-water_FVFmodel ANALYTICAL
-water_compress 3e-6
-water_denscoeff 0.
-water_densmodel ANALYTICAL
-water_prefdens 0.
-water_preffvf 14.7
-water_prefvisc 14.7
-water_refFVF 1.0142
-water_refdens 1.0
-water_refvisc 0.96
-water_visccoeff 0.
-water_viscmodel ANALYTICAL
#End of PETSc Option Table entries
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=4583
total number of function evaluations=8
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 7 iterations 3.
Checking 0
Time Step: 0
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=16048
total number of function evaluations=24
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 17 iterations 4.
Checking 1
Time Step: 1
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=15778
total number of function evaluations=27
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 18 iterations 4.
Checking 2
Time Step: 2
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=18413
total number of function evaluations=32
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 20 iterations 4.
Checking 3
Time Step: 3
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=21146
total number of function evaluations=29
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 22 iterations 4.
Checking 4
Time Step: 4
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=44142
total number of function evaluations=54
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
[ERROR] snesP diverged with reason -5
Checking 5
Time Step: 5
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=16223
total number of function evaluations=23
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 17 iterations 4.
Checking 6
Time Step: 6
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=46100
total number of function evaluations=52
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
[ERROR] snesP diverged with reason -5
Checking 7
Time Step: 7
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=27735
total number of function evaluations=34
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 33 iterations 4.
Checking 8
Time Step: 8
SNES Object: (P_) 1 MPI processes
type: newtonls
maximum iterations=50, maximum function evaluations=10000
tolerances: relative=1e-08, absolute=1e-08, solution=1e-08
total number of linear solver iterations=33027
total number of function evaluations=37
norm schedule ALWAYS
SNESLineSearch Object: (P_) 1 MPI processes
type: bt
interpolation: cubic
alpha=1.000000e-04
maxstep=1.000000e+08, minlambda=1.000000e-12
tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
maximum iterations=40
KSP Object: (P_) 1 MPI processes
type: fgmres
restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-08, absolute=1e-08, divergence=10000.
right preconditioning
using UNPRECONDITIONED norm type for convergence test
PC Object: (P_) 1 MPI processes
type: fieldsplit
FieldSplit with MULTIPLICATIVE composition: total splits = 3
Solver info for each split is in the following KSP objects:
Split number 0 Defined by IS
KSP Object: (P_fieldsplit_0_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_0_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_0_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 1 Defined by IS
KSP Object: (P_fieldsplit_1_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_1_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2500, cols=2500
package used to perform factorization: petsc
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_1_) 1 MPI processes
type: seqaij
rows=2500, cols=2500
total: nonzeros=21904, allocated nonzeros=21904
total number of mallocs used during MatSetValues calls =0
not using I-node routines
Split number 2 Defined by IS
KSP Object: (P_fieldsplit_2_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
left preconditioning
using NONE norm type for convergence test
PC Object: (P_fieldsplit_2_) 1 MPI processes
type: ilu
out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 1., needed 1.
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=2, cols=2
package used to perform factorization: petsc
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: (P_fieldsplit_2_) 1 MPI processes
type: seqaij
rows=2, cols=2
total: nonzeros=4, allocated nonzeros=4
total number of mallocs used during MatSetValues calls =0
using I-node routines: found 1 nodes, limit used is 5
linear system matrix followed by preconditioner matrix:
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
Mat Object: 1 MPI processes
type: nest
rows=5002, cols=5002
Matrix object:
type=nest, rows=3, cols=3
MatNest structure:
(0,0) : prefix="P_fieldsplit_0_", type=seqaij, rows=2500, cols=2500
(0,1) : type=seqaij, rows=2500, cols=2500
(0,2) : type=seqaij, rows=2500, cols=2
(1,0) : type=seqaij, rows=2500, cols=2500
(1,1) : prefix="P_fieldsplit_1_", type=seqaij, rows=2500, cols=2500
(1,2) : type=seqaij, rows=2500, cols=2
(2,0) : type=seqaij, rows=2, cols=2500
(2,1) : type=seqaij, rows=2, cols=2500
(2,2) : prefix="P_fieldsplit_2_", type=seqaij, rows=2, cols=2
snesP converged in 36 iterations 4.
Checking 9
Time Step: 9
rank: 0/1 NEW TIME STEP
-------------- next part --------------
FracD:
Numerical implementation of multiphysics models.
(c) 2016-2017 Chukwudi Chukwudozie chdozie at gmail.com
Well object "well":
top: 3.000000e+03 3.000000e+03
bottom: 3.000000e+03 3.000000e+03
wellnodes: 1.500000e+03 1.500000e+03
well water rate: 1.000000e+02
well oil rate: 0.000000e+00
well gas rate: 0.000000e+00
total well rate: 0.000000e+00
well bottomhole pressure: 5.000000e+03
well height/production interval: 3.000000e+01
well radius: 2.920000e-01
well block equivalent radius: 0.000000e+00
well skin: 0.000000e+00
Well Type: "INJECTORWATER" well under "WATERRATE" condition
Well object "well":
top: 3.000000e+03 3.000000e+03
bottom: 3.000000e+03 3.000000e+03
wellnodes: 4.500000e+03 4.500000e+03
well water rate: -0.000000e+00
well oil rate: -5.000000e+02
well gas rate: -0.000000e+00
total well rate: -0.000000e+00
well bottomhole pressure: 2.000000e+03
well height/production interval: 3.000000e+01
well radius: 2.920000e-01
well block equivalent radius: 0.000000e+00
well skin: 0.000000e+00
Well Type: "PRODUCER" well under "PRESSURE"