** ====================================================================================== 
**
**             ___________                          _______
**     _______/          /_____ ________ __ _   ___/       \ ____
**    /      /___    ___/      /__   __/  /  \ /  /\__ _   /     \
**   /  //  /   /   /  /  //  /  /  / /  /    /  /___/_   /  //  /
**  /  ___ /   /   /  /  _   /  /  / /  /  /    //       /  //  /
** /__/       /___/  /__//__/  /__/ /__/__/ \__//_______/______/
**
** Authors:  Dave A. May          (dave.may@erdw.ethz.ch)           
**           Laetitia Le Pourhiet (laetitia.le_pourhiet@upmc.fr)    
**           Jed Brown            (jedbrown@mcs.anl.gov)            
**
** git url: https://bitbucket.org/jedbrown/ptatin3d.git 
** commit hash: [out-of-date] Execute "make releaseinfo" to update to the most recent revision 
** log: [out-of-date] Execute "make releaseinfo" to update to the most recent revision 
**                                                                       
** TATIN_CFLAGS = -std=gnu99 -O0 -Wall -g
**                                                                       
** WARNING pTatin3d appears to have been compiled with debug options 
** For significant performance improvements, please consult the file makefile.arch  
** Adjust TATIN_CFLAGS to include aggressive compiler optimizations 
**                                                                       
** AVX not detected - optimized kernels will not be used.
** If your system supports AVX, consider adding options,
** e.g. -march=native, to TATIN_CFLAGS in makefile.arch
** ====================================================================================== 
DataBucketView <Material Constants>:
  L                  = 800 
  buffer (max)       = 0 
  allocated          = 800 
  nfields registered = 17 
    [  0]: field name  ==>>     MaterialConst_MaterialType : Mem. usage = 3.20e-03 (MB) : rank0
    [  1]: field name  ==>>   MaterialConst_ViscosityConst : Mem. usage = 1.60e-03 (MB) : rank0
    [  2]: field name  ==>>       MaterialConst_ViscosityZ : Mem. usage = 4.80e-03 (MB) : rank0
    [  3]: field name  ==>>    MaterialConst_ViscosityArrh : Mem. usage = 1.28e-02 (MB) : rank0
    [  4]: field name  ==>>      MaterialConst_ViscosityFK : Mem. usage = 3.20e-03 (MB) : rank0
    [  5]: field name  ==>>     MaterialConst_PlasticMises : Mem. usage = 3.20e-03 (MB) : rank0
    [  6]: field name  ==>>        MaterialConst_PlasticDP : Mem. usage = 9.60e-03 (MB) : rank0
    [  7]: field name  ==>>     MaterialConst_DensityConst : Mem. usage = 1.60e-03 (MB) : rank0
    [  8]: field name  ==>> MaterialConst_DensityBoussinesq : Mem. usage = 4.80e-03 (MB) : rank0
    [  9]: field name  ==>>          MaterialConst_SoftLin : Mem. usage = 3.20e-03 (MB) : rank0
    [ 10]: field name  ==>>         MaterialConst_SoftExpo : Mem. usage = 3.20e-03 (MB) : rank0
    [ 11]: field name  ==>>        EnergyMaterialConstants : Mem. usage = 1.44e-02 (MB) : rank0
    [ 12]: field name  ==>>        EnergyConductivityConst : Mem. usage = 1.60e-03 (MB) : rank0
    [ 13]: field name  ==>>    EnergyConductivityThreshold : Mem. usage = 6.40e-03 (MB) : rank0
    [ 14]: field name  ==>>              EnergySourceConst : Mem. usage = 1.60e-03 (MB) : rank0
    [ 15]: field name  ==>>              EnergySourceDecay : Mem. usage = 3.20e-03 (MB) : rank0
    [ 16]: field name  ==>> EnergySourceAdiabaticAdvection : Mem. usage = 1.60e-03 (MB) : rank0
  Total mem. usage                                                      = 3.20e-01 (MB) : <collective over 4 ranks>
[pTatin] Created output directory: sd3d_test 
[pTatin] Created log file: sd3d_test/ptatin.log-2018.03.15_08:17:07 
[pTatin] Created options file: sd3d_test/ptatin.options-2018.03.15_08:17:07 
[pTatin] Created options file: sd3d_test/ptatin.options 
  [pTatinModel]: Registering model [0] with name "template"
  [pTatinModel]: Registering model [1] with name "viscous_sinker"
  [pTatinModel]: Registering model [2] with name "Gene3D"
  [pTatinModel]: Registering model [3] with name "indentor"
  [pTatinModel]: Registering model [4] with name "rift3D_T"
  [pTatinModel]: Registering model [5] with name "advdiff_example"
  [pTatinModel]: Registering model [6] with name "delamination"
  [pTatinModel]: Registering model [7] with name "Riftrh"
  [pTatinModel]: Registering model [8] with name "Rift_oblique3d"
  [pTatinModel]: Registering model [9] with name "geomod2008"
  [pTatinModel]: Registering model [10] with name "multilayer_folding"
  [pTatinModel]: Registering model [11] with name "submarinelavaflow"
  [pTatinModel]: Registering model [12] with name "ex_subduction"
  [pTatinModel]: Registering model [13] with name "iplus"
  [pTatinModel]: Registering model [14] with name "subduction_initiation2d"
  [pTatinModel]: Registering model [15] with name "convection2d"
  [pTatinModel]: Registering model [16] with name "thermal_sb"
  [pTatinModel]: Registering model [17] with name "sd3d"
  [pTatinModel]: Registering model [18] with name "pas"
  [pTatinModel]: Registering model [19] with name "pd"
  [pTatinModelDynamic]: Dynamically registering model with name "static_box"
  [pTatinModelDynamic]: Dynamically registering model with name "static_box_thermomech"
  [pTatinModelDynamic]: Dynamically registering model with name "analytics_vv"
  [pTatinModel]: -ptatin_model "sd3d" was detected
  [model sd3d]: CASE 1A
Current Value -preexpA_0   :  1.0000e+00  
Current Value -Ascale_0   :  1.0000e+21  
Current Value -entalpy_0   :  0.0000e+00  
Current Value -Vmol_0   :  0.0000e+00  
Current Value -nexp_0   :  1.0000e+00  
Current Value -Tref_0   :  1.0000e+00  
Current Value -Eta_scale_0   :  1.0000e+24  
Current Value -P_scale_0   :  1.5152e+08  
Current Value -density_0   :  1.3721e+01  
Current Value -preexpA_1   :  1.0000e+00  
Current Value -Ascale_1   :  4.7500e+11  
Current Value -entalpy_1   :  0.0000e+00  
Current Value -Vmol_1   :  0.0000e+00  
Current Value -nexp_1   :  4.0000e+00  
Current Value -Tref_1   :  1.0000e+00  
Current Value -Eta_scale_1   :  1.0000e+24  
Current Value -P_scale_1   :  1.5152e+08  
Current Value -density_1   :  1.4375e+01  
Current Value -preexpA_2   :  1.0000e+00  
Current Value -Ascale_2   :  4.7500e+11  
Current Value -entalpy_2   :  0.0000e+00  
Current Value -Vmol_2   :  0.0000e+00  
Current Value -nexp_2   :  4.0000e+00  
Current Value -Tref_2   :  1.0000e+00  
Current Value -Eta_scale_2   :  1.0000e+24  
Current Value -P_scale_2   :  1.5152e+08  
Current Value -density_2   :  1.4375e+01  
Current Value -preexpA_3   :  1.0000e+00  
Current Value -Ascale_3   :  4.7500e+11  
Current Value -entalpy_3   :  0.0000e+00  
Current Value -Vmol_3   :  0.0000e+00  
Current Value -nexp_3   :  4.0000e+00  
Current Value -Tref_3   :  1.0000e+00  
Current Value -Eta_scale_3   :  1.0000e+24  
Current Value -P_scale_3   :  1.5152e+08  
Current Value -density_3   :  1.4375e+01  
Current Value -preexpA_4   :  1.0000e+00  
Current Value -Ascale_4   :  4.7500e+11  
Current Value -entalpy_4   :  0.0000e+00  
Current Value -Vmol_4   :  0.0000e+00  
Current Value -nexp_4   :  4.0000e+00  
Current Value -Tref_4   :  1.0000e+00  
Current Value -Eta_scale_4   :  1.0000e+24  
Current Value -P_scale_4   :  1.5152e+08  
Current Value -density_4   :  1.4375e+01  
Generated vel/pressure mesh --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
  MaterialPointsStokes: Using Q1 projection
[[Swarm initialization: 0.0001 (sec)]]
[[Swarm->coordinate assignment: 3456 points : 0.0016 (sec)]]
[[SwarmDMDA3dDataExchangerCreate: time = 2.8019e-03 (sec)]]
Generated material points --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
[[ModelApplyInitialMeshGeometry_SD3D]]
[[ModelApplyInitialMaterialGeometry_SD3D]]
DataBucketView <MaterialPoints StokesCoefficients>:
  L                  = 19860 
  buffer (max)       = 4000 
  allocated          = 23860 
  nfields registered = 3 
    [  0]: field name  ==>>                        MPntStd : Mem. usage = 3.00e-01 (MB) : rank0
    [  1]: field name  ==>>                    MPntPStokes : Mem. usage = 7.51e-02 (MB) : rank0
    [  2]: field name  ==>>                  MPntPStokesPl : Mem. usage = 3.76e-02 (MB) : rank0
  Total mem. usage                                                      = 2.10e+00 (MB) : <collective over 4 ranks>
Mesh size (8 x 8 x 8) : MG levels 2  
         level [ 0]: global Q2 elements (4 x 4 x 4) 
         level [ 1]: global Q2 elements (8 x 8 x 8) 
level [ 0]: [total cores    4]: np-I [   1]: element range I [    4 ]
                                np-J [   2]: element range J [    2,    2 ]
                                np-K [   2]: element range K [    2,    2 ]
level [ 1]: [total cores    4]: np-I [   1]: element range I [    8 ]
                                np-J [   2]: element range J [    4,    4 ]
                                np-K [   2]: element range K [    4,    4 ]
Generated velocity mesh hierarchy --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
Generated quadrature point hierarchy --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
Level [1]: Coarse grid type :: Re-discretisation :: matrix free operator 
Level [0]: Coarse grid type :: Re-discretisation :: assembled operator 
Generated stokes operators --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
   [[ COMPUTING FLOW FIELD FOR STEP : 0 ]]
PreSolve(0) --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
   --------- LINEAR STAGE ---------
*** Rheology update for RHEOLOGY_VISCOUS selected ***
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e-02; log10(max/min) 1.00e+00; cpu time 2.42e-04 (sec)
  0 SNES Function norm 1.386601539724e+00 
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e-02; log10(max/min) 1.00e+00; cpu time 3.72e-04 (sec)
  0 KSP Component U,V,W,P residual norm [ 0.000000000000e+00, 1.386601539724e+00, 0.000000000000e+00, 0.000000000000e+00 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000000000000e+00 
      1 KSP Residual norm 3.545327890629e-01 
      2 KSP Residual norm 7.021883413814e-02 
      3 KSP Residual norm 1.371789719984e-02 
      4 KSP Residual norm 2.748262359186e-03 
  1 KSP Component U,V,W,P residual norm [ 3.959819892258e-08, 1.386547893157e+00, 1.077529522696e-07, 8.624589010438e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 6.068006612263e-03 
      1 KSP Residual norm 1.461217594662e-03 
      2 KSP Residual norm 2.372703772723e-04 
      3 KSP Residual norm 4.903251703257e-05 
  2 KSP Component U,V,W,P residual norm [ 4.801883326243e-06, 1.386007552382e+00, 7.813325636904e-06, 2.869269665969e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.878573416490e-02 
      1 KSP Residual norm 4.486511306807e-03 
      2 KSP Residual norm 7.798903572489e-04 
      3 KSP Residual norm 1.612846609179e-04 
  3 KSP Component U,V,W,P residual norm [ 3.856613810554e-05, 1.383171499908e+00, 4.293390447425e-05, 6.887909344760e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.953491780703e-02 
      1 KSP Residual norm 4.371989513665e-03 
      2 KSP Residual norm 7.996408619699e-04 
      3 KSP Residual norm 1.463360312058e-04 
  4 KSP Component U,V,W,P residual norm [ 4.341691959835e-04, 1.346135305076e+00, 4.754767123636e-04, 2.333935682561e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.861363781858e-02 
      1 KSP Residual norm 3.589139473608e-03 
      2 KSP Residual norm 6.296198609134e-04 
      3 KSP Residual norm 1.261751559598e-04 
  5 KSP Component U,V,W,P residual norm [ 2.227785769248e-03, 1.134509530245e+00, 2.374791881002e-03, 5.347790105253e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.649457504265e-02 
      1 KSP Residual norm 4.537813642321e-03 
      2 KSP Residual norm 9.573583592331e-04 
      3 KSP Residual norm 1.672440469229e-04 
  6 KSP Component U,V,W,P residual norm [ 4.529727175717e-03, 5.348741187505e-01, 4.734597488711e-03, 6.749153245123e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.491435679367e-02 
      1 KSP Residual norm 6.598290398577e-03 
      2 KSP Residual norm 1.404200832620e-03 
      3 KSP Residual norm 2.108221462854e-04 
  7 KSP Component U,V,W,P residual norm [ 2.190447852254e-03, 1.077580981085e-01, 2.199221021805e-03, 3.711887506873e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.899072917342e-02 
      1 KSP Residual norm 6.366618222711e-03 
      2 KSP Residual norm 1.294785111067e-03 
      3 KSP Residual norm 1.768794073109e-04 
  8 KSP Component U,V,W,P residual norm [ 5.518129107560e-04, 2.558713596904e-02, 4.669519798595e-04, 1.866014074291e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.466818424240e-02 
      1 KSP Residual norm 5.746366398100e-03 
      2 KSP Residual norm 1.139513148914e-03 
      3 KSP Residual norm 1.928658808343e-04 
  9 KSP Component U,V,W,P residual norm [ 1.491182703390e-04, 3.552560792946e-03, 8.575357786354e-05, 7.007783070455e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.525182701861e-02 
      1 KSP Residual norm 3.459174167212e-03 
      2 KSP Residual norm 5.756313774359e-04 
      3 KSP Residual norm 1.356213627450e-04 
 10 KSP Component U,V,W,P residual norm [ 1.788476670538e-04, 3.369266068804e-04, 1.782043583048e-04, 2.046103541568e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 5.540610433500e-02 
      1 KSP Residual norm 5.921774209134e-03 
      2 KSP Residual norm 1.338907440813e-03 
      3 KSP Residual norm 1.913790494300e-04 
 11 KSP Component U,V,W,P residual norm [ 1.046696932722e-04, 9.448470945497e-05, 9.590497119224e-05, 7.432168535445e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 5.225554166262e-02 
      1 KSP Residual norm 9.897820934874e-03 
      2 KSP Residual norm 2.226988805145e-03 
      3 KSP Residual norm 2.788162902839e-04 
 12 KSP Component U,V,W,P residual norm [ 3.431150342539e-05, 3.014674110976e-05, 3.991350763519e-05, 3.368655358849e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 5.502127187292e-02 
      1 KSP Residual norm 1.070394283281e-02 
      2 KSP Residual norm 2.278240656281e-03 
      3 KSP Residual norm 2.570554996570e-04 
 13 KSP Component U,V,W,P residual norm [ 2.733836004496e-06, 3.008520514646e-06, 3.633930076539e-06, 1.179859424685e-03 ]
  Linear solve converged due to CONVERGED_RTOL iterations 13
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e-02; log10(max/min) 1.00e+00; cpu time 2.82e-04 (sec)
  1 SNES Function norm 1.179872023790e-03 
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 1
SNES Object: 4 MPI processes
  type: newtonls
  maximum iterations=1, maximum function evaluations=10000
  tolerances: relative=1e-05, absolute=1e-08, solution=1e-08
  total number of linear solver iterations=13
  total number of function evaluations=2
  norm schedule ALWAYS
  SNESLineSearch Object: 4 MPI processes
    type: basic
    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: 4 MPI processes
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=2000, nonzero initial guess
    tolerances:  relative=0.001, absolute=1e-08, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 4 MPI processes
    type: fieldsplit
      FieldSplit with Schur preconditioner, factorization UPPER
      Preconditioner for the Schur complement formed from A11
      Split info:
      Split number 0 Defined by IS
      Split number 1 Defined by IS
      KSP solver for A00 block
        KSP Object: (fieldsplit_u_) 4 MPI processes
          type: fgmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
          right preconditioning
          using UNPRECONDITIONED norm type for convergence test
        PC Object: (fieldsplit_u_) 4 MPI processes
          type: mg
            type is MULTIPLICATIVE, levels=2 cycles=v
              Cycles per PCApply=1
              Not using Galerkin computed coarse grid matrices
          Coarse grid solver -- level -------------------------------
            KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=4, initial guess is zero
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using UNPRECONDITIONED norm type for convergence test
            PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: bjacobi
                number of blocks = 4
                Local solve is same for all blocks, in the following KSP and PC objects:
              KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                        package used to perform factorization: petsc
                        total: nonzeros=126711, allocated nonzeros=126711
                        total number of mallocs used during MatSetValues calls =0
                          using I-node routines: found 225 nodes, limit used is 5
                linear system matrix = precond matrix:
                Mat Object: (Buu_) 1 MPI processes
                  type: seqaij
                  rows=675, cols=675, bs=3
                  total: nonzeros=126711, allocated nonzeros=126711
                  total number of mallocs used during MatSetValues calls =0
                    using I-node routines: found 225 nodes, limit used is 5
              linear system matrix = precond matrix:
              Mat Object: (Buu_) 4 MPI processes
                type: mpiaij
                rows=2187, cols=2187, bs=3
                total: nonzeros=533871, allocated nonzeros=533871
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
          Down solver (pre-smoother) on level 1 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=2, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: (fieldsplit_u_) 4 MPI processes
                type: shell
                rows=14739, cols=14739, bs=3
          Up solver (post-smoother) same as down solver (pre-smoother)
          linear system matrix = precond matrix:
          Mat Object: (fieldsplit_u_) 4 MPI processes
            type: shell
            rows=14739, cols=14739, bs=3
      KSP solver for S = A11 - A10 inv(A00) A01 
        KSP Object: (fieldsplit_p_) 4 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: (fieldsplit_p_) 4 MPI processes
          type: jacobi
          linear system matrix followed by preconditioner matrix:
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: schurcomplement
            rows=2048, cols=2048
              Schur complement A11 - A10 inv(A00) A01
              A11
                Mat Object: (fieldsplit_p_) 4 MPI processes
                  type: mpisbaij
                  rows=2048, cols=2048, bs=4
                  total: nonzeros=8192, allocated nonzeros=8192
                  total number of mallocs used during MatSetValues calls =0
              A10
                Mat Object: (Bpu_) 4 MPI processes
                  type: shell
                  rows=2048, cols=14739
              KSP of A00
                KSP Object: (fieldsplit_u_) 4 MPI processes
                  type: fgmres
                    restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                    happy breakdown tolerance 1e-30
                  maximum iterations=10, initial guess is zero
                  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
                  right preconditioning
                  using UNPRECONDITIONED norm type for convergence test
                PC Object: (fieldsplit_u_) 4 MPI processes
                  type: mg
                    type is MULTIPLICATIVE, levels=2 cycles=v
                      Cycles per PCApply=1
                      Not using Galerkin computed coarse grid matrices
                  Coarse grid solver -- level -------------------------------
                    KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=4, initial guess is zero
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using UNPRECONDITIONED norm type for convergence test
                    PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: bjacobi
                        number of blocks = 4
                        Local solve is same for all blocks, in the following KSP and PC objects:
                      KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                                package used to perform factorization: petsc
                                total: nonzeros=126711, allocated nonzeros=126711
                                total number of mallocs used during MatSetValues calls =0
                                  using I-node routines: found 225 nodes, limit used is 5
                        linear system matrix = precond matrix:
                        Mat Object: (Buu_) 1 MPI processes
                          type: seqaij
                          rows=675, cols=675, bs=3
                          total: nonzeros=126711, allocated nonzeros=126711
                          total number of mallocs used during MatSetValues calls =0
                            using I-node routines: found 225 nodes, limit used is 5
                      linear system matrix = precond matrix:
                      Mat Object: (Buu_) 4 MPI processes
                        type: mpiaij
                        rows=2187, cols=2187, bs=3
                        total: nonzeros=533871, allocated nonzeros=533871
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                  Down solver (pre-smoother) on level 1 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=2, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: (fieldsplit_u_) 4 MPI processes
                        type: shell
                        rows=14739, cols=14739, bs=3
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  linear system matrix = precond matrix:
                  Mat Object: (fieldsplit_u_) 4 MPI processes
                    type: shell
                    rows=14739, cols=14739, bs=3
              A01
                Mat Object: (Bup_) 4 MPI processes
                  type: shell
                  rows=14739, cols=2048
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: mpisbaij
            rows=2048, cols=2048, bs=4
            total: nonzeros=8192, allocated nonzeros=8192
            total number of mallocs used during MatSetValues calls =0
    linear system matrix followed by preconditioner matrix:
    Mat Object: (stokes_Amf_) 4 MPI processes
      type: shell
      rows=16787, cols=16787
    Mat Object: 4 MPI processes
      type: nest
      rows=16787, cols=16787
        Matrix object: 
          type=nest, rows=2, cols=2 
          MatNest structure: 
          (0,0) : prefix="fieldsplit_u_", type=shell, rows=14739, cols=14739 
          (0,1) : prefix="Bup_", type=shell, rows=14739, cols=2048 
          (1,0) : prefix="Bpu_", type=shell, rows=2048, cols=14739 
          (1,1) : prefix="fieldsplit_p_", type=mpisbaij, rows=2048, cols=2048 
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e-02; log10(max/min) 1.00e+00; cpu time 2.28e-04 (sec)
[[ModelOutput_SD3D]]
  writing pvdfilename sd3d_test/timeseries_v.pvd 
  writing pvdfilename sd3d_test/timeseries_vp.pvd 
  writing pvdfilename sd3d_test/timeseries_mpoints_cell.pvd 
slab base range: gmin 6.9697e-01 5.0000e-01 0.0000e+00
slab base range: gmax 7.5758e-01 5.0000e-01 3.7879e-01
slab edge range: gmax(x) 6.9697e-01
slab edge identified: pidx_w,rank_w 5551 1
[step 0] slab edge range: W,Dn 6.9697e-01 8.1566e-01
slab edge (back face) range: gmax(x) 6.9697e-01
slab edge (back face) identified: pidx_w,rank_w 5117 1
[step 0] slab edge (back face) range: W,Dn 6.9697e-01 5.7891e-01
slab edge (front face) range: gmin(z) 3.7879e-01
slab edge (front face) identified: pidx_w,rank_w 4843 1
[step 0] slab edge (front face) range: W,Dn 3.7879e-01 5.1578e-01
topo range: 1.0000e+00 1.0000e+00
vol,int v.v,vrms: 5.7392e-01 7.5860e+00 3.6356e+00
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e-02; log10(max/min) 1.00e+00; cpu time 4.52e-04 (sec)
Phi: 0.0000e+00
   --------- PICARD STAGE ---------
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; npoints_yielded 0; cpu time 2.12e-02 (sec)
  0 SNES Function norm 4.775270692598e-01 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; npoints_yielded 0; cpu time 2.21e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 4.323463059691e-01, 1.631895869650e+00, 2.796259607959e-01, 2.359718849371e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.999337514723e-01 
      1 KSP Residual norm 7.377663916584e-01 
      2 KSP Residual norm 3.247753649000e-01 
      3 KSP Residual norm 1.180000143558e-01 
      4 KSP Residual norm 4.476193682316e-02 
      5 KSP Residual norm 2.266441997007e-02 
      6 KSP Residual norm 7.126643874169e-03 
  1 KSP Component U,V,W,P residual norm [ 4.323113329898e-01, 1.631763922576e+00, 2.796033191056e-01, 1.561211320214e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.608872633067e-02 
      1 KSP Residual norm 7.610552612554e-03 
      2 KSP Residual norm 2.591614508146e-03 
      3 KSP Residual norm 8.134818099332e-04 
      4 KSP Residual norm 3.286406878188e-04 
      5 KSP Residual norm 1.345728859649e-04 
  2 KSP Component U,V,W,P residual norm [ 4.318722425688e-01, 1.630088176617e+00, 2.793100125193e-01, 5.697187913210e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.294086120970e-01 
      1 KSP Residual norm 4.630962292006e-02 
      2 KSP Residual norm 1.298539442426e-02 
      3 KSP Residual norm 4.741242682937e-03 
      4 KSP Residual norm 1.566799361099e-03 
      5 KSP Residual norm 6.182545100751e-04 
  3 KSP Component U,V,W,P residual norm [ 4.310636145363e-01, 1.627105399860e+00, 2.788427362297e-01, 9.249484097897e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.743232294349e-01 
      1 KSP Residual norm 9.034002744114e-02 
      2 KSP Residual norm 2.277628530902e-02 
      3 KSP Residual norm 8.689672373806e-03 
      4 KSP Residual norm 2.898199410230e-03 
      5 KSP Residual norm 1.014690899332e-03 
  4 KSP Component U,V,W,P residual norm [ 4.181022672362e-01, 1.577871786375e+00, 2.700981356665e-01, 3.062270765656e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.546074195955e-01 
      1 KSP Residual norm 7.714030011324e-02 
      2 KSP Residual norm 1.775904334489e-02 
      3 KSP Residual norm 6.204353801260e-03 
      4 KSP Residual norm 2.340059915739e-03 
  5 KSP Component U,V,W,P residual norm [ 3.887250048232e-01, 1.466456623288e+00, 2.516492917045e-01, 5.164728390805e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.476424992973e-01 
      1 KSP Residual norm 1.038753365913e-01 
      2 KSP Residual norm 2.380139338460e-02 
      3 KSP Residual norm 7.530804879961e-03 
      4 KSP Residual norm 3.307860114344e-03 
  6 KSP Component U,V,W,P residual norm [ 2.154124821157e-01, 8.169447268091e-01, 1.416703730772e-01, 8.555365402711e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.819403044579e-01 
      1 KSP Residual norm 1.085730980137e-01 
      2 KSP Residual norm 2.310334687462e-02 
      3 KSP Residual norm 7.321495060811e-03 
      4 KSP Residual norm 2.770113831632e-03 
  7 KSP Component U,V,W,P residual norm [ 9.831040678220e-02, 3.671706782828e-01, 6.308899148922e-02, 7.142478918077e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 7.166818587634e-01 
      1 KSP Residual norm 1.842646368334e-01 
      2 KSP Residual norm 4.274495865966e-02 
      3 KSP Residual norm 1.452104764682e-02 
      4 KSP Residual norm 6.756055844838e-03 
  8 KSP Component U,V,W,P residual norm [ 1.829237640260e-02, 6.452476009699e-02, 1.253464896631e-02, 3.325627046165e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.777216640745e-01 
      1 KSP Residual norm 1.950452077816e-01 
      2 KSP Residual norm 4.495350478930e-02 
      3 KSP Residual norm 1.124623888507e-02 
      4 KSP Residual norm 3.667226874458e-03 
  9 KSP Component U,V,W,P residual norm [ 3.424445899770e-03, 1.290170698799e-02, 2.685385105095e-03, 1.509279964586e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.715601642689e-01 
      1 KSP Residual norm 1.636403410059e-01 
      2 KSP Residual norm 4.131527676694e-02 
      3 KSP Residual norm 1.066949864305e-02 
      4 KSP Residual norm 3.834077584361e-03 
 10 KSP Component U,V,W,P residual norm [ 1.441490215607e-03, 2.835456666080e-03, 1.201124859330e-03, 6.996547816850e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.272613540573e+00 
      1 KSP Residual norm 2.153028286848e-01 
      2 KSP Residual norm 4.987195190209e-02 
      3 KSP Residual norm 1.175210046255e-02 
 11 KSP Component U,V,W,P residual norm [ 6.503941849258e-04, 9.165844858627e-04, 5.703092731550e-04, 3.248193633218e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.064086402635e-01 
      1 KSP Residual norm 1.367346726616e-01 
      2 KSP Residual norm 3.792511352879e-02 
      3 KSP Residual norm 9.336292182979e-03 
      4 KSP Residual norm 3.282048453045e-03 
 12 KSP Component U,V,W,P residual norm [ 3.461452438667e-04, 4.266637125579e-04, 3.492848154632e-04, 1.890428643838e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 8.634442722723e-01 
      1 KSP Residual norm 1.998760119833e-01 
      2 KSP Residual norm 4.878436346496e-02 
      3 KSP Residual norm 1.143474229067e-02 
      4 KSP Residual norm 4.236429204690e-03 
 13 KSP Component U,V,W,P residual norm [ 8.323351068219e-05, 9.991496467406e-05, 1.004660704907e-04, 7.136005277841e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 7.817142750668e-01 
      1 KSP Residual norm 1.696181825091e-01 
      2 KSP Residual norm 3.645052559000e-02 
      3 KSP Residual norm 9.476119768737e-03 
      4 KSP Residual norm 3.598975795253e-03 
 14 KSP Component U,V,W,P residual norm [ 4.724619117585e-05, 4.605591434139e-05, 4.964569171786e-05, 3.569129770909e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.180405756570e+00 
      1 KSP Residual norm 2.303125845581e-01 
      2 KSP Residual norm 5.616214474426e-02 
      3 KSP Residual norm 1.372883004927e-02 
      4 KSP Residual norm 4.185940545579e-03 
 15 KSP Component U,V,W,P residual norm [ 1.824177508403e-05, 1.489367870953e-05, 1.143303356730e-05, 1.616757747387e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.193482508461e+00 
      1 KSP Residual norm 1.775685117735e-01 
      2 KSP Residual norm 4.323695171164e-02 
      3 KSP Residual norm 9.530144746013e-03 
 16 KSP Component U,V,W,P residual norm [ 2.855169223328e-05, 2.251626862826e-05, 1.751882310573e-05, 8.772662819644e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.260267444807e+00 
      1 KSP Residual norm 2.349957847476e-01 
      2 KSP Residual norm 6.198660150927e-02 
      3 KSP Residual norm 1.275108677933e-02 
      4 KSP Residual norm 3.908411132993e-03 
 17 KSP Component U,V,W,P residual norm [ 4.319761131848e-06, 6.693352903722e-06, 7.600429730381e-06, 4.081751945950e-04 ]
  Linear solve converged due to CONVERGED_RTOL iterations 17
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.89e-02 (sec)
  1 SNES Function norm 2.585390859938e-02 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.91e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.611774152150e-01, 3.959729186057e-01, 1.610310482639e-01, 1.315550690850e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000402566235e+00 
      1 KSP Residual norm 3.255598463171e-01 
      2 KSP Residual norm 9.025152801920e-02 
      3 KSP Residual norm 2.520344417865e-02 
      4 KSP Residual norm 7.559005570874e-03 
  1 KSP Component U,V,W,P residual norm [ 1.632934288045e-02, 2.531630340631e-02, 1.056716680898e-02, 1.219141494502e-01 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.051964279873e+00 
      1 KSP Residual norm 3.815412132807e-01 
      2 KSP Residual norm 6.209131892160e-02 
      3 KSP Residual norm 1.310286814228e-02 
      4 KSP Residual norm 3.982309957017e-03 
  2 KSP Component U,V,W,P residual norm [ 6.909967394073e-03, 9.699591298120e-03, 4.124616327610e-03, 7.718551136992e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.307040508487e+00 
      1 KSP Residual norm 3.637791046713e-01 
      2 KSP Residual norm 7.933241979169e-02 
      3 KSP Residual norm 2.468891687751e-02 
      4 KSP Residual norm 9.672013322752e-03 
  3 KSP Component U,V,W,P residual norm [ 9.597075650735e-04, 1.625655515558e-03, 6.609139333696e-04, 2.934273508946e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.479634252733e+00 
      1 KSP Residual norm 2.400478604870e-01 
      2 KSP Residual norm 5.867558918112e-02 
      3 KSP Residual norm 1.984904407879e-02 
      4 KSP Residual norm 6.865889815805e-03 
  4 KSP Component U,V,W,P residual norm [ 4.021543191942e-04, 5.159370099396e-04, 2.798937528305e-04, 1.836766843996e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.032548582400e+00 
      1 KSP Residual norm 2.691954719620e-01 
      2 KSP Residual norm 5.035222907839e-02 
      3 KSP Residual norm 1.699085885972e-02 
      4 KSP Residual norm 5.536553096601e-03 
  5 KSP Component U,V,W,P residual norm [ 1.692395915481e-04, 1.609927883130e-04, 1.772253551945e-04, 5.781137027744e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.063517752669e+00 
      1 KSP Residual norm 2.604915783026e-01 
      2 KSP Residual norm 5.834639668165e-02 
      3 KSP Residual norm 1.823158181320e-02 
  6 KSP Component U,V,W,P residual norm [ 8.412834313132e-05, 9.659261075920e-05, 1.083867814409e-04, 3.342289252924e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.655534649312e+00 
      1 KSP Residual norm 2.571188710351e-01 
      2 KSP Residual norm 4.916329437885e-02 
      3 KSP Residual norm 1.423703889188e-02 
  7 KSP Component U,V,W,P residual norm [ 7.231846223840e-05, 5.590733736026e-05, 7.459007722738e-05, 1.847767878904e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.200311960727e+00 
      1 KSP Residual norm 2.431866593725e-01 
      2 KSP Residual norm 6.810097629754e-02 
      3 KSP Residual norm 2.838790837561e-02 
      4 KSP Residual norm 1.316009733144e-02 
  8 KSP Component U,V,W,P residual norm [ 1.457436284724e-05, 1.327395588587e-05, 1.346028797059e-05, 9.520941788204e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.183272909429e+00 
      1 KSP Residual norm 2.156399606469e-01 
      2 KSP Residual norm 6.593847066819e-02 
      3 KSP Residual norm 2.635127542301e-02 
      4 KSP Residual norm 1.150628432661e-02 
  9 KSP Component U,V,W,P residual norm [ 2.823651309201e-05, 3.123789657966e-05, 3.439092078397e-05, 4.585618086269e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.798359131144e+00 
      1 KSP Residual norm 2.274868821696e-01 
      2 KSP Residual norm 5.651471902318e-02 
      3 KSP Residual norm 2.750314178338e-02 
      4 KSP Residual norm 1.296492480066e-02 
 10 KSP Component U,V,W,P residual norm [ 1.038892913110e-05, 1.133799133400e-05, 1.172699062875e-05, 2.697379895880e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.760208179289e+00 
      1 KSP Residual norm 2.558899315718e-01 
      2 KSP Residual norm 7.239776613464e-02 
      3 KSP Residual norm 3.407423663598e-02 
      4 KSP Residual norm 1.607895858886e-02 
 11 KSP Component U,V,W,P residual norm [ 7.524650795477e-06, 8.220620269767e-06, 7.797348235897e-06, 1.412688768452e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.073491043496e+00 
      1 KSP Residual norm 2.002217477296e-01 
      2 KSP Residual norm 5.318876955884e-02 
      3 KSP Residual norm 1.774500484351e-02 
 12 KSP Component U,V,W,P residual norm [ 2.668927684192e-06, 2.933795369981e-06, 2.161693510953e-06, 6.435461536192e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.424003247138e+00 
      1 KSP Residual norm 2.322892342351e-01 
      2 KSP Residual norm 5.442839429043e-02 
      3 KSP Residual norm 1.412814165848e-02 
 13 KSP Component U,V,W,P residual norm [ 2.868380897807e-06, 2.960083574403e-06, 2.465929330034e-06, 4.040524935716e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.550577994904e+00 
      1 KSP Residual norm 2.600229091053e-01 
      2 KSP Residual norm 7.229804673247e-02 
      3 KSP Residual norm 3.039008073085e-02 
      4 KSP Residual norm 1.445269197182e-02 
 14 KSP Component U,V,W,P residual norm [ 8.312162698920e-07, 9.538834089967e-07, 1.125378347460e-06, 2.159289777112e-05 ]
  Linear solve converged due to CONVERGED_RTOL iterations 14
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  2 SNES Function norm 1.033310105618e-02 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.30e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.590919608326e-02, 2.055553424214e-02, 1.162685527460e-02, 4.204990350437e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.997058670739e-01 
      1 KSP Residual norm 2.538807722153e-01 
      2 KSP Residual norm 7.179550351349e-02 
      3 KSP Residual norm 2.184410412618e-02 
      4 KSP Residual norm 9.388674461237e-03 
  1 KSP Component U,V,W,P residual norm [ 3.494555577660e-03, 2.850953908795e-03, 1.566080646805e-03, 1.204886442182e-02 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.715149783021e-01 
      1 KSP Residual norm 2.643272170199e-02 
      2 KSP Residual norm 1.107558629379e-02 
      3 KSP Residual norm 5.114679790829e-03 
      4 KSP Residual norm 1.891388833410e-03 
  2 KSP Component U,V,W,P residual norm [ 1.407137068646e-03, 9.866306168380e-04, 7.195311223761e-04, 7.696860425404e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 5.475833734452e-01 
      1 KSP Residual norm 1.300898361266e-01 
      2 KSP Residual norm 6.067738620878e-02 
      3 KSP Residual norm 2.381918133631e-02 
      4 KSP Residual norm 9.281601997251e-03 
      5 KSP Residual norm 3.379920045830e-03 
  3 KSP Component U,V,W,P residual norm [ 1.393365196404e-04, 9.948279001172e-05, 9.957265033724e-05, 2.285318952031e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.112519676611e+00 
      1 KSP Residual norm 1.499992678738e-01 
      2 KSP Residual norm 4.486281842494e-02 
      3 KSP Residual norm 1.702099635824e-02 
      4 KSP Residual norm 7.528831124530e-03 
  4 KSP Component U,V,W,P residual norm [ 4.125556102741e-05, 2.806550626560e-05, 2.448748539857e-05, 1.369991706019e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.542778959873e-01 
      1 KSP Residual norm 2.161133284930e-01 
      2 KSP Residual norm 7.350058562922e-02 
      3 KSP Residual norm 2.474469764754e-02 
      4 KSP Residual norm 1.134187483115e-02 
      5 KSP Residual norm 4.966683963810e-03 
  5 KSP Component U,V,W,P residual norm [ 8.906467141300e-06, 9.069183000740e-06, 6.853420890041e-06, 4.701923627995e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.385348945026e+00 
      1 KSP Residual norm 1.737819592555e-01 
      2 KSP Residual norm 4.162730957021e-02 
      3 KSP Residual norm 1.266123007130e-02 
  6 KSP Component U,V,W,P residual norm [ 4.966590055748e-06, 4.338616707130e-06, 4.678324949808e-06, 3.066812733747e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 8.809430112226e-01 
      1 KSP Residual norm 1.947523182176e-01 
      2 KSP Residual norm 5.312993808474e-02 
      3 KSP Residual norm 1.716441367275e-02 
      4 KSP Residual norm 7.560266316506e-03 
  7 KSP Component U,V,W,P residual norm [ 3.296102902782e-06, 3.027926585170e-06, 4.441842414092e-06, 1.304220171221e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.702984604915e+00 
      1 KSP Residual norm 1.838105198304e-01 
      2 KSP Residual norm 4.097092440274e-02 
      3 KSP Residual norm 1.345608372049e-02 
  8 KSP Component U,V,W,P residual norm [ 2.067598956942e-06, 1.813234207609e-06, 2.167027577511e-06, 6.872669821619e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.440457001408e+00 
      1 KSP Residual norm 2.684170850878e-01 
      2 KSP Residual norm 7.168564386077e-02 
      3 KSP Residual norm 2.351107415828e-02 
      4 KSP Residual norm 1.092322151634e-02 
  9 KSP Component U,V,W,P residual norm [ 9.027027122065e-07, 1.704834549146e-06, 1.971821245137e-06, 3.451428259876e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.339883902459e+00 
      1 KSP Residual norm 1.903959480104e-01 
      2 KSP Residual norm 6.526258665194e-02 
      3 KSP Residual norm 2.971567409000e-02 
      4 KSP Residual norm 1.326114555326e-02 
 10 KSP Component U,V,W,P residual norm [ 6.116446858130e-07, 5.026312522287e-07, 6.588500707840e-07, 1.661589459972e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.591355222174e+00 
      1 KSP Residual norm 2.565764315846e-01 
      2 KSP Residual norm 8.165743272885e-02 
      3 KSP Residual norm 3.029696861203e-02 
      4 KSP Residual norm 1.212596249754e-02 
 11 KSP Component U,V,W,P residual norm [ 3.514219613496e-07, 4.400508705686e-07, 4.270027514169e-07, 8.139058202205e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 11
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.92e-02 (sec)
  3 SNES Function norm 6.376576964794e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 6.783122084541e-03, 6.438517077404e-03, 3.055297076739e-03, 2.736194196830e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.998744666348e-01 
      1 KSP Residual norm 2.459543302050e-01 
      2 KSP Residual norm 7.779913286055e-02 
      3 KSP Residual norm 2.860677918133e-02 
      4 KSP Residual norm 1.342147490862e-02 
      5 KSP Residual norm 6.380053381438e-03 
  1 KSP Component U,V,W,P residual norm [ 8.119968639900e-04, 7.717553377682e-04, 3.601972645802e-04, 3.189677634093e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.678128592623e-01 
      1 KSP Residual norm 2.089740732946e-02 
      2 KSP Residual norm 9.133497654095e-03 
      3 KSP Residual norm 4.156715843697e-03 
      4 KSP Residual norm 1.499623947917e-03 
  2 KSP Component U,V,W,P residual norm [ 9.623871324335e-05, 1.193072738089e-04, 1.360814529269e-04, 1.043159987508e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.126904214967e+00 
      1 KSP Residual norm 2.356283170753e-01 
      2 KSP Residual norm 1.097168088250e-01 
      3 KSP Residual norm 3.997593113355e-02 
      4 KSP Residual norm 1.651787955105e-02 
      5 KSP Residual norm 6.846821085269e-03 
  3 KSP Component U,V,W,P residual norm [ 2.192800180723e-05, 2.045451625283e-05, 1.508848852196e-05, 5.494502709858e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.188607921852e+00 
      1 KSP Residual norm 2.655217552980e-01 
      2 KSP Residual norm 1.186846261979e-01 
      3 KSP Residual norm 4.261705928326e-02 
      4 KSP Residual norm 1.875467907259e-02 
      5 KSP Residual norm 8.876111665468e-03 
  4 KSP Component U,V,W,P residual norm [ 1.208979797593e-05, 1.855673909994e-05, 2.223564804505e-05, 2.822713987106e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.421745676288e+00 
      1 KSP Residual norm 3.055310785106e-01 
      2 KSP Residual norm 1.202432638791e-01 
      3 KSP Residual norm 3.831240825845e-02 
      4 KSP Residual norm 1.496387048111e-02 
      5 KSP Residual norm 7.286339909808e-03 
  5 KSP Component U,V,W,P residual norm [ 1.378024960459e-06, 1.295021303219e-06, 1.138747936989e-06, 1.020314304677e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.432461235709e+00 
      1 KSP Residual norm 1.680518712330e-01 
      2 KSP Residual norm 5.505546737786e-02 
      3 KSP Residual norm 2.123468631714e-02 
      4 KSP Residual norm 1.099249631630e-02 
  6 KSP Component U,V,W,P residual norm [ 7.306353564904e-07, 7.114111184496e-07, 7.132638583899e-07, 7.131353869475e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.188128340658e+00 
      1 KSP Residual norm 2.788406410284e-01 
      2 KSP Residual norm 9.563102255661e-02 
      3 KSP Residual norm 3.209664145412e-02 
      4 KSP Residual norm 1.364941548210e-02 
      5 KSP Residual norm 6.885853349583e-03 
  7 KSP Component U,V,W,P residual norm [ 7.295637206752e-07, 4.578572717655e-07, 6.065473555965e-07, 2.903158903486e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.876065541325e+00 
      1 KSP Residual norm 2.394762860618e-01 
      2 KSP Residual norm 7.567355946517e-02 
      3 KSP Residual norm 2.367490149786e-02 
      4 KSP Residual norm 7.342734806868e-03 
  8 KSP Component U,V,W,P residual norm [ 4.937724344144e-07, 7.215647223426e-07, 9.016858301491e-07, 1.687772136365e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.283650793924e+00 
      1 KSP Residual norm 2.745472854606e-01 
      2 KSP Residual norm 1.036785450796e-01 
      3 KSP Residual norm 4.569879162358e-02 
      4 KSP Residual norm 2.313953742979e-02 
      5 KSP Residual norm 8.849270514125e-03 
  9 KSP Component U,V,W,P residual norm [ 2.306047448620e-07, 2.881435857252e-07, 4.309921579200e-07, 8.374862233795e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.122155064942e+00 
      1 KSP Residual norm 2.861206876030e-01 
      2 KSP Residual norm 8.168719097706e-02 
      3 KSP Residual norm 3.466550974624e-02 
      4 KSP Residual norm 1.620224427892e-02 
 10 KSP Component U,V,W,P residual norm [ 2.109417903380e-07, 2.716888176218e-07, 3.551593139552e-07, 4.371493506245e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 10
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  4 SNES Function norm 4.783461266623e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.89e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.925748731819e-03, 3.065266290166e-03, 1.311889446602e-03, 8.550911105029e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000466329111e+00 
      1 KSP Residual norm 2.315972178395e-01 
      2 KSP Residual norm 7.913254971719e-02 
      3 KSP Residual norm 3.101071309288e-02 
      4 KSP Residual norm 1.498631708590e-02 
      5 KSP Residual norm 7.491369839858e-03 
  1 KSP Component U,V,W,P residual norm [ 2.407358774441e-04, 2.539803656646e-04, 1.060497992797e-04, 1.217612980413e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.755258045437e-01 
      1 KSP Residual norm 2.962402976239e-02 
      2 KSP Residual norm 1.276975313662e-02 
      3 KSP Residual norm 4.933847459159e-03 
      4 KSP Residual norm 1.925681278796e-03 
  2 KSP Component U,V,W,P residual norm [ 4.215318313100e-05, 5.344668771923e-05, 6.775413765950e-05, 4.049008223594e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.352287231385e+00 
      1 KSP Residual norm 3.482442385354e-01 
      2 KSP Residual norm 1.563417494106e-01 
      3 KSP Residual norm 5.060447136528e-02 
      4 KSP Residual norm 2.177076728436e-02 
      5 KSP Residual norm 1.040270848168e-02 
  3 KSP Component U,V,W,P residual norm [ 6.923371702074e-06, 7.722030634457e-06, 7.738927177874e-06, 2.042722989147e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.685859817628e+00 
      1 KSP Residual norm 3.941509793464e-01 
      2 KSP Residual norm 1.710268821482e-01 
      3 KSP Residual norm 5.621033448454e-02 
      4 KSP Residual norm 2.404084402560e-02 
      5 KSP Residual norm 1.218010258731e-02 
  4 KSP Component U,V,W,P residual norm [ 5.349038672557e-06, 9.233924073745e-06, 1.162949397735e-05, 1.073242028870e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.973162325897e+00 
      1 KSP Residual norm 3.965427996039e-01 
      2 KSP Residual norm 1.532395660158e-01 
      3 KSP Residual norm 4.671596383114e-02 
      4 KSP Residual norm 1.789739801318e-02 
  5 KSP Component U,V,W,P residual norm [ 1.178681100573e-06, 9.362868882225e-07, 1.380868436608e-06, 3.738879085154e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.710376241749e+00 
      1 KSP Residual norm 2.736338105265e-01 
      2 KSP Residual norm 1.089175085651e-01 
      3 KSP Residual norm 4.864875970084e-02 
      4 KSP Residual norm 2.526160076617e-02 
      5 KSP Residual norm 9.839396389185e-03 
  6 KSP Component U,V,W,P residual norm [ 3.304236391430e-07, 3.979056597241e-07, 4.914212817281e-07, 2.641247677683e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.422127911726e+00 
      1 KSP Residual norm 3.511180624809e-01 
      2 KSP Residual norm 1.434498313889e-01 
      3 KSP Residual norm 5.724265709216e-02 
      4 KSP Residual norm 2.800574261589e-02 
      5 KSP Residual norm 1.267031034600e-02 
  7 KSP Component U,V,W,P residual norm [ 4.900182044861e-07, 7.493207175320e-07, 1.084546886951e-06, 1.098428941900e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.265682812171e+00 
      1 KSP Residual norm 2.759292581545e-01 
      2 KSP Residual norm 9.703744585989e-02 
      3 KSP Residual norm 4.006226844117e-02 
      4 KSP Residual norm 1.992597816553e-02 
  8 KSP Component U,V,W,P residual norm [ 3.120892591225e-07, 4.507166507682e-07, 6.190610943473e-07, 5.965001647212e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.427812925874e+00 
      1 KSP Residual norm 3.773977333025e-01 
      2 KSP Residual norm 1.680965967221e-01 
      3 KSP Residual norm 7.869595281623e-02 
      4 KSP Residual norm 3.434554861181e-02 
      5 KSP Residual norm 1.114757714995e-02 
  9 KSP Component U,V,W,P residual norm [ 2.118611092488e-07, 3.106589880524e-07, 4.586327165921e-07, 3.239168832014e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  5 SNES Function norm 3.799090156380e-03 
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 5
SNES Object: 4 MPI processes
  type: newtonls
  maximum iterations=5, maximum function evaluations=10000
  tolerances: relative=1e-05, absolute=1e-08, solution=1e-08
  total number of linear solver iterations=61
  total number of function evaluations=6
  norm schedule ALWAYS
  SNESLineSearch Object: 4 MPI processes
    type: basic
    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: 4 MPI processes
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=2000, nonzero initial guess
    tolerances:  relative=0.001, absolute=1e-08, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 4 MPI processes
    type: fieldsplit
      FieldSplit with Schur preconditioner, factorization UPPER
      Preconditioner for the Schur complement formed from A11
      Split info:
      Split number 0 Defined by IS
      Split number 1 Defined by IS
      KSP solver for A00 block
        KSP Object: (fieldsplit_u_) 4 MPI processes
          type: fgmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
          right preconditioning
          using UNPRECONDITIONED norm type for convergence test
        PC Object: (fieldsplit_u_) 4 MPI processes
          type: mg
            type is MULTIPLICATIVE, levels=2 cycles=v
              Cycles per PCApply=1
              Not using Galerkin computed coarse grid matrices
          Coarse grid solver -- level -------------------------------
            KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=4, initial guess is zero
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using UNPRECONDITIONED norm type for convergence test
            PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: bjacobi
                number of blocks = 4
                Local solve is same for all blocks, in the following KSP and PC objects:
              KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                        package used to perform factorization: petsc
                        total: nonzeros=126711, allocated nonzeros=126711
                        total number of mallocs used during MatSetValues calls =0
                          using I-node routines: found 225 nodes, limit used is 5
                linear system matrix = precond matrix:
                Mat Object: (Buu_) 1 MPI processes
                  type: seqaij
                  rows=675, cols=675, bs=3
                  total: nonzeros=126711, allocated nonzeros=126711
                  total number of mallocs used during MatSetValues calls =0
                    using I-node routines: found 225 nodes, limit used is 5
              linear system matrix = precond matrix:
              Mat Object: (Buu_) 4 MPI processes
                type: mpiaij
                rows=2187, cols=2187, bs=3
                total: nonzeros=533871, allocated nonzeros=533871
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
          Down solver (pre-smoother) on level 1 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=2, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: (fieldsplit_u_) 4 MPI processes
                type: shell
                rows=14739, cols=14739, bs=3
          Up solver (post-smoother) same as down solver (pre-smoother)
          linear system matrix = precond matrix:
          Mat Object: (fieldsplit_u_) 4 MPI processes
            type: shell
            rows=14739, cols=14739, bs=3
      KSP solver for S = A11 - A10 inv(A00) A01 
        KSP Object: (fieldsplit_p_) 4 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: (fieldsplit_p_) 4 MPI processes
          type: jacobi
          linear system matrix followed by preconditioner matrix:
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: schurcomplement
            rows=2048, cols=2048
              Schur complement A11 - A10 inv(A00) A01
              A11
                Mat Object: (fieldsplit_p_) 4 MPI processes
                  type: mpisbaij
                  rows=2048, cols=2048, bs=4
                  total: nonzeros=8192, allocated nonzeros=8192
                  total number of mallocs used during MatSetValues calls =0
              A10
                Mat Object: (Bpu_) 4 MPI processes
                  type: shell
                  rows=2048, cols=14739
              KSP of A00
                KSP Object: (fieldsplit_u_) 4 MPI processes
                  type: fgmres
                    restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                    happy breakdown tolerance 1e-30
                  maximum iterations=10, initial guess is zero
                  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
                  right preconditioning
                  using UNPRECONDITIONED norm type for convergence test
                PC Object: (fieldsplit_u_) 4 MPI processes
                  type: mg
                    type is MULTIPLICATIVE, levels=2 cycles=v
                      Cycles per PCApply=1
                      Not using Galerkin computed coarse grid matrices
                  Coarse grid solver -- level -------------------------------
                    KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=4, initial guess is zero
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using UNPRECONDITIONED norm type for convergence test
                    PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: bjacobi
                        number of blocks = 4
                        Local solve is same for all blocks, in the following KSP and PC objects:
                      KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                                package used to perform factorization: petsc
                                total: nonzeros=126711, allocated nonzeros=126711
                                total number of mallocs used during MatSetValues calls =0
                                  using I-node routines: found 225 nodes, limit used is 5
                        linear system matrix = precond matrix:
                        Mat Object: (Buu_) 1 MPI processes
                          type: seqaij
                          rows=675, cols=675, bs=3
                          total: nonzeros=126711, allocated nonzeros=126711
                          total number of mallocs used during MatSetValues calls =0
                            using I-node routines: found 225 nodes, limit used is 5
                      linear system matrix = precond matrix:
                      Mat Object: (Buu_) 4 MPI processes
                        type: mpiaij
                        rows=2187, cols=2187, bs=3
                        total: nonzeros=533871, allocated nonzeros=533871
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                  Down solver (pre-smoother) on level 1 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=2, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: (fieldsplit_u_) 4 MPI processes
                        type: shell
                        rows=14739, cols=14739, bs=3
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  linear system matrix = precond matrix:
                  Mat Object: (fieldsplit_u_) 4 MPI processes
                    type: shell
                    rows=14739, cols=14739, bs=3
              A01
                Mat Object: (Bup_) 4 MPI processes
                  type: shell
                  rows=14739, cols=2048
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: mpisbaij
            rows=2048, cols=2048, bs=4
            total: nonzeros=8192, allocated nonzeros=8192
            total number of mallocs used during MatSetValues calls =0
    linear system matrix followed by preconditioner matrix:
    Mat Object: (stokes_Amf_) 4 MPI processes
      type: shell
      rows=16787, cols=16787
    Mat Object: 4 MPI processes
      type: nest
      rows=16787, cols=16787
        Matrix object: 
          type=nest, rows=2, cols=2 
          MatNest structure: 
          (0,0) : prefix="fieldsplit_u_", type=shell, rows=14739, cols=14739 
          (0,1) : prefix="Bup_", type=shell, rows=14739, cols=2048 
          (1,0) : prefix="Bpu_", type=shell, rows=2048, cols=14739 
          (1,1) : prefix="fieldsplit_p_", type=mpisbaij, rows=2048, cols=2048 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.91e-02 (sec)
[[ModelOutput_SD3D]]
slab base range: gmin 6.9697e-01 5.0000e-01 0.0000e+00
slab base range: gmax 7.5758e-01 5.0000e-01 3.7879e-01
slab edge range: gmax(x) 6.9697e-01
slab edge identified: pidx_w,rank_w 5551 1
[step 0] slab edge range: W,Dn 6.9697e-01 8.1566e-01
slab edge (back face) range: gmax(x) 6.9697e-01
slab edge (back face) identified: pidx_w,rank_w 5117 1
[step 0] slab edge (back face) range: W,Dn 6.9697e-01 5.7891e-01
slab edge (front face) range: gmin(z) 3.7879e-01
slab edge (front face) identified: pidx_w,rank_w 4843 1
[step 0] slab edge (front face) range: W,Dn 3.7879e-01 5.1578e-01
topo range: 1.0000e+00 1.0000e+00
vol,int v.v,vrms: 5.7392e-01 2.0679e-02 1.8982e-01
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.19e-02 (sec)
Phi: 0.0000e+00
PostSolve(0) --> pTatin3dCurrentMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
[[ModelOutput_SD3D]]
slab base range: gmin 6.9697e-01 5.0000e-01 0.0000e+00
slab base range: gmax 7.5758e-01 5.0000e-01 3.7879e-01
slab edge range: gmax(x) 6.9697e-01
slab edge identified: pidx_w,rank_w 5551 1
[step 0] slab edge range: W,Dn 6.9697e-01 8.1566e-01
slab edge (back face) range: gmax(x) 6.9697e-01
slab edge (back face) identified: pidx_w,rank_w 5117 1
[step 0] slab edge (back face) range: W,Dn 6.9697e-01 5.7891e-01
slab edge (front face) range: gmin(z) 3.7879e-01
slab edge (front face) identified: pidx_w,rank_w 4843 1
[step 0] slab edge (front face) range: W,Dn 3.7879e-01 5.1578e-01
topo range: 1.0000e+00 1.0000e+00
vol,int v.v,vrms: 5.7392e-01 2.0679e-02 1.8982e-01
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.70e-02 (sec)
Phi: 0.0000e+00
  TimeStep control(StkCourant): | current = 1.0000e+32 : trial = 2.4662e-01 [accepted] | dt > dt_max : restricted to = 2.3891e-03 | ==>> dt used = 2.3891e-03 |
  TimeStep control(StkSurfaceCourant): | current = 2.3891e-03 : trial = 1.0000e+32 | ==>> dt used = 2.3891e-03 |
  timestep[stokes] dt_courant = 2.3891e-03 
<<----------------------------------------------------------------------------------------------->>
   [[ EXECUTING TIME STEP : 1 ]]
     dt    : 2.3891e-13 
     time  : 2.3891e-13 
=========== MaterialPointStd_UpdateCoordinates ============= 
      Number of material points                     19860 (init) / 19860 (final)
    	MaterialPointStd_UpdateLocalCoordinates     5.13e-03 (sec) efficiency 39.3%
    	SwarmUpdatePosition_Communication_Generic   3.66e-03 (sec) efficiency 15.6%
    	MaterialPointStd_Removal                    2.60e-05 (sec) efficiency 46.8%
=========== ================================== ============= 
   [[ COMPUTING FLOW FIELD FOR STEP : 1 ]]
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 SNES Function norm 3.799090156380e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.329199543321e-03, 2.852981335902e-03, 9.318811664980e-04, 3.239168832014e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.999115427678e-01 
      1 KSP Residual norm 2.181557241208e-01 
      2 KSP Residual norm 7.646882494740e-02 
      3 KSP Residual norm 3.229313601440e-02 
      4 KSP Residual norm 1.558383709326e-02 
      5 KSP Residual norm 7.669508934503e-03 
  1 KSP Component U,V,W,P residual norm [ 1.877000847945e-04, 2.315861653976e-04, 7.571890631944e-05, 1.034074921980e-03 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.196342905920e-01 
      1 KSP Residual norm 1.259857289538e-01 
      2 KSP Residual norm 3.364625564119e-02 
      3 KSP Residual norm 1.124351232544e-02 
      4 KSP Residual norm 3.921804586906e-03 
  2 KSP Component U,V,W,P residual norm [ 3.203551847554e-05, 4.325536443543e-05, 4.210240074988e-05, 3.910767825614e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.131943864949e+00 
      1 KSP Residual norm 5.244775386040e-01 
      2 KSP Residual norm 1.630184244701e-01 
      3 KSP Residual norm 5.175204978117e-02 
      4 KSP Residual norm 1.732479388698e-02 
  3 KSP Component U,V,W,P residual norm [ 8.092449373569e-06, 9.562175715944e-06, 5.671859346418e-06, 1.995870382365e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.632624381865e+00 
      1 KSP Residual norm 4.075206156238e-01 
      2 KSP Residual norm 1.343689501345e-01 
      3 KSP Residual norm 5.179288819545e-02 
      4 KSP Residual norm 1.546535641667e-02 
  4 KSP Component U,V,W,P residual norm [ 3.339136568412e-06, 4.200900163815e-06, 4.057024081394e-06, 9.640315048049e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.834780040750e+00 
      1 KSP Residual norm 5.566274863295e-01 
      2 KSP Residual norm 2.163833499234e-01 
      3 KSP Residual norm 8.045574564093e-02 
      4 KSP Residual norm 2.942349933447e-02 
      5 KSP Residual norm 1.439155581342e-02 
  5 KSP Component U,V,W,P residual norm [ 1.480463228200e-06, 2.099378459224e-06, 2.540681690210e-06, 4.100663486862e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.285120673373e+00 
      1 KSP Residual norm 3.206632279938e-01 
      2 KSP Residual norm 1.137498383491e-01 
      3 KSP Residual norm 3.532925121140e-02 
      4 KSP Residual norm 1.196052174789e-02 
  6 KSP Component U,V,W,P residual norm [ 6.956791588759e-07, 5.560390934085e-07, 6.245496273027e-07, 2.517230337349e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.073081314794e+00 
      1 KSP Residual norm 4.118215129878e-01 
      2 KSP Residual norm 1.269770665004e-01 
      3 KSP Residual norm 4.456271498145e-02 
      4 KSP Residual norm 1.961275125903e-02 
  7 KSP Component U,V,W,P residual norm [ 3.142856057890e-07, 7.473928541279e-07, 8.064778628017e-07, 1.242231005627e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.430399655012e+00 
      1 KSP Residual norm 4.993930130662e-01 
      2 KSP Residual norm 1.750757495741e-01 
      3 KSP Residual norm 5.962178070734e-02 
      4 KSP Residual norm 1.847989109994e-02 
  8 KSP Component U,V,W,P residual norm [ 4.793045109534e-07, 6.068281764811e-07, 7.623600626414e-07, 6.446911752977e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.643033736122e+00 
      1 KSP Residual norm 5.095681902371e-01 
      2 KSP Residual norm 1.677953158476e-01 
      3 KSP Residual norm 6.309892426297e-02 
      4 KSP Residual norm 3.209035500212e-02 
      5 KSP Residual norm 1.345091024010e-02 
  9 KSP Component U,V,W,P residual norm [ 2.212870254416e-07, 3.815608109321e-07, 4.542859084322e-07, 4.282859052668e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.391113921009e+00 
      1 KSP Residual norm 4.073388782640e-01 
      2 KSP Residual norm 1.301665739752e-01 
      3 KSP Residual norm 4.986505716706e-02 
      4 KSP Residual norm 2.368716771725e-02 
 10 KSP Component U,V,W,P residual norm [ 1.228793239231e-07, 2.117085842306e-07, 2.479764781355e-07, 2.424630254131e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 10
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.08e-02 (sec)
  1 SNES Function norm 2.835206361053e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.08e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 1.817217926037e-03, 2.073972283139e-03, 6.593536296897e-04, 2.424630254131e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000231492783e+00 
      1 KSP Residual norm 2.038074123900e-01 
      2 KSP Residual norm 7.101139283566e-02 
      3 KSP Residual norm 3.124206195041e-02 
      4 KSP Residual norm 1.460912161259e-02 
      5 KSP Residual norm 7.044139871745e-03 
  1 KSP Component U,V,W,P residual norm [ 1.139094493573e-04, 1.312650321727e-04, 4.222299667600e-05, 6.872392192366e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.107447191476e+00 
      1 KSP Residual norm 1.496761278527e-01 
      2 KSP Residual norm 3.700878583246e-02 
      3 KSP Residual norm 1.353964541743e-02 
      4 KSP Residual norm 5.001053195955e-03 
  2 KSP Component U,V,W,P residual norm [ 2.080337199790e-05, 2.819506164860e-05, 2.877138714911e-05, 2.636539241172e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.574658449688e+00 
      1 KSP Residual norm 5.746898092976e-01 
      2 KSP Residual norm 1.711227352448e-01 
      3 KSP Residual norm 6.095806728922e-02 
      4 KSP Residual norm 2.086575954781e-02 
  3 KSP Component U,V,W,P residual norm [ 5.736848101292e-06, 6.540697855058e-06, 4.456815139885e-06, 1.365041218188e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.202754165524e+00 
      1 KSP Residual norm 5.029049798964e-01 
      2 KSP Residual norm 1.720829716465e-01 
      3 KSP Residual norm 6.579965844874e-02 
      4 KSP Residual norm 1.876315949938e-02 
  4 KSP Component U,V,W,P residual norm [ 2.583447493162e-06, 3.489009962178e-06, 3.671190114265e-06, 6.339815996642e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.531993592553e+00 
      1 KSP Residual norm 6.452197760098e-01 
      2 KSP Residual norm 2.614249623509e-01 
      3 KSP Residual norm 9.252353327816e-02 
      4 KSP Residual norm 3.458358276692e-02 
  5 KSP Component U,V,W,P residual norm [ 1.765013537463e-06, 1.945995584607e-06, 2.219114584298e-06, 2.614502192360e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.005835857328e+00 
      1 KSP Residual norm 3.339937938798e-01 
      2 KSP Residual norm 1.131677084576e-01 
      3 KSP Residual norm 3.402233517508e-02 
      4 KSP Residual norm 1.128378287590e-02 
  6 KSP Component U,V,W,P residual norm [ 2.991172521811e-07, 3.093491173822e-07, 2.706077746992e-07, 1.745840372364e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.036756431089e+00 
      1 KSP Residual norm 3.900490685129e-01 
      2 KSP Residual norm 1.275193352259e-01 
      3 KSP Residual norm 4.991288232689e-02 
      4 KSP Residual norm 2.231937955160e-02 
      5 KSP Residual norm 9.963806560955e-03 
  7 KSP Component U,V,W,P residual norm [ 4.555746154673e-07, 6.388594377711e-07, 6.886509096881e-07, 8.572035735325e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.961008581682e+00 
      1 KSP Residual norm 5.353010774848e-01 
      2 KSP Residual norm 1.997687551058e-01 
      3 KSP Residual norm 6.376134217095e-02 
      4 KSP Residual norm 1.924612973474e-02 
  8 KSP Component U,V,W,P residual norm [ 4.588968129766e-07, 4.665213509386e-07, 5.693640404573e-07, 4.375932747513e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.944456411282e+00 
      1 KSP Residual norm 5.639201157427e-01 
      2 KSP Residual norm 1.985824216624e-01 
      3 KSP Residual norm 8.022312130677e-02 
      4 KSP Residual norm 3.932799166716e-02 
      5 KSP Residual norm 1.486331689835e-02 
  9 KSP Component U,V,W,P residual norm [ 2.177524378331e-07, 2.818562156994e-07, 3.197802733229e-07, 2.816659505730e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.753740210642e+00 
      1 KSP Residual norm 4.237074039859e-01 
      2 KSP Residual norm 1.358861227531e-01 
      3 KSP Residual norm 4.542465964276e-02 
      4 KSP Residual norm 2.090367195239e-02 
 10 KSP Component U,V,W,P residual norm [ 1.064884321266e-07, 1.335920805763e-07, 1.442076108484e-07, 1.606054007840e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 10
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  2 SNES Function norm 1.838203468011e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 1.188702969967e-03, 1.336153662158e-03, 4.250506448546e-04, 1.606054007840e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000195592856e+00 
      1 KSP Residual norm 2.010856867221e-01 
      2 KSP Residual norm 7.129302686596e-02 
      3 KSP Residual norm 3.178520238108e-02 
      4 KSP Residual norm 1.453507232110e-02 
      5 KSP Residual norm 6.825283181439e-03 
  1 KSP Component U,V,W,P residual norm [ 6.742012075830e-05, 7.641667928959e-05, 2.483224256170e-05, 4.249257953135e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.184961155820e+00 
      1 KSP Residual norm 1.620926664707e-01 
      2 KSP Residual norm 3.970873246646e-02 
      3 KSP Residual norm 1.567294847245e-02 
      4 KSP Residual norm 6.215402708846e-03 
  2 KSP Component U,V,W,P residual norm [ 1.288985326026e-05, 1.782202283035e-05, 1.834465378368e-05, 1.640453270684e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.823488205372e+00 
      1 KSP Residual norm 6.250584849515e-01 
      2 KSP Residual norm 1.844540475789e-01 
      3 KSP Residual norm 6.991202303570e-02 
      4 KSP Residual norm 2.451210178398e-02 
  3 KSP Component U,V,W,P residual norm [ 3.808640908616e-06, 4.387290359747e-06, 3.149958675617e-06, 8.656269877721e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.515484857324e+00 
      1 KSP Residual norm 5.524833375420e-01 
      2 KSP Residual norm 1.940248658553e-01 
      3 KSP Residual norm 7.007785404223e-02 
      4 KSP Residual norm 1.990522953385e-02 
  4 KSP Component U,V,W,P residual norm [ 1.709033453651e-06, 2.494995545052e-06, 2.764876360724e-06, 3.875286741607e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.808037470653e+00 
      1 KSP Residual norm 6.585543822471e-01 
      2 KSP Residual norm 2.690692060847e-01 
      3 KSP Residual norm 9.335246965272e-02 
      4 KSP Residual norm 3.613818172850e-02 
  5 KSP Component U,V,W,P residual norm [ 1.052524909348e-06, 1.183345573501e-06, 1.241449031914e-06, 1.605108952138e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.483598948328e+00 
      1 KSP Residual norm 3.363601483470e-01 
      2 KSP Residual norm 1.083026951648e-01 
      3 KSP Residual norm 3.465657597102e-02 
  6 KSP Component U,V,W,P residual norm [ 4.152293896484e-07, 4.135106124715e-07, 3.341907442869e-07, 1.079083896491e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.919756173598e+00 
      1 KSP Residual norm 3.779605281056e-01 
      2 KSP Residual norm 1.143205413888e-01 
      3 KSP Residual norm 3.979566269425e-02 
      4 KSP Residual norm 1.626348247052e-02 
  7 KSP Component U,V,W,P residual norm [ 4.292854233804e-07, 4.997471540542e-07, 5.951646743644e-07, 5.390077189155e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 4.135005787327e+00 
      1 KSP Residual norm 5.507340599513e-01 
      2 KSP Residual norm 2.032727261394e-01 
      3 KSP Residual norm 6.280254976387e-02 
      4 KSP Residual norm 1.761020126346e-02 
  8 KSP Component U,V,W,P residual norm [ 2.916333027197e-07, 2.620756392563e-07, 3.194085169141e-07, 2.733032784675e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.134780750959e+00 
      1 KSP Residual norm 5.553734036250e-01 
      2 KSP Residual norm 1.882268569496e-01 
      3 KSP Residual norm 7.192897882356e-02 
      4 KSP Residual norm 3.356497051234e-02 
      5 KSP Residual norm 1.170749717473e-02 
  9 KSP Component U,V,W,P residual norm [ 1.502307555136e-07, 1.741747536345e-07, 2.008479290702e-07, 1.687162667704e-06 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  3 SNES Function norm 1.066796162461e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 6.866638359531e-04, 7.761853375981e-04, 2.531408774088e-04, 1.687162667703e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000231217125e+00 
      1 KSP Residual norm 2.103148089651e-01 
      2 KSP Residual norm 7.316731371452e-02 
      3 KSP Residual norm 3.264880949167e-02 
      4 KSP Residual norm 1.475717416652e-02 
      5 KSP Residual norm 6.697276640044e-03 
  1 KSP Component U,V,W,P residual norm [ 3.835119777751e-05, 4.359842918260e-05, 1.454155778061e-05, 2.446814712554e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.164384541988e+00 
      1 KSP Residual norm 1.642661624736e-01 
      2 KSP Residual norm 4.127384532400e-02 
      3 KSP Residual norm 1.690006711775e-02 
      4 KSP Residual norm 7.019707074609e-03 
  2 KSP Component U,V,W,P residual norm [ 7.445644889709e-06, 1.038582297379e-05, 1.062808317466e-05, 9.412188645682e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.884313927299e+00 
      1 KSP Residual norm 6.692613499986e-01 
      2 KSP Residual norm 1.966839749216e-01 
      3 KSP Residual norm 7.610483933718e-02 
      4 KSP Residual norm 2.756744624225e-02 
  3 KSP Component U,V,W,P residual norm [ 2.246776765443e-06, 2.750204538850e-06, 1.918026527611e-06, 5.023160966346e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.639109993118e+00 
      1 KSP Residual norm 5.714528417425e-01 
      2 KSP Residual norm 1.991525853747e-01 
      3 KSP Residual norm 7.048071867708e-02 
      4 KSP Residual norm 2.000310475263e-02 
  4 KSP Component U,V,W,P residual norm [ 1.030701073191e-06, 1.486533889804e-06, 1.714177726690e-06, 2.224164630970e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.898101173622e+00 
      1 KSP Residual norm 6.656518419348e-01 
      2 KSP Residual norm 2.724815633902e-01 
      3 KSP Residual norm 9.494819133935e-02 
      4 KSP Residual norm 3.788033244028e-02 
  5 KSP Component U,V,W,P residual norm [ 5.939281343606e-07, 6.939729512135e-07, 6.527515682006e-07, 9.310050421685e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.660894065494e+00 
      1 KSP Residual norm 3.387490200269e-01 
      2 KSP Residual norm 1.069850520345e-01 
      3 KSP Residual norm 3.463972165507e-02 
  6 KSP Component U,V,W,P residual norm [ 2.633456166150e-07, 2.265140043686e-07, 1.732538766534e-07, 6.146056833758e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.887052869219e+00 
      1 KSP Residual norm 3.574404890191e-01 
      2 KSP Residual norm 1.110100181392e-01 
      3 KSP Residual norm 3.777177773974e-02 
      4 KSP Residual norm 1.536959832214e-02 
  7 KSP Component U,V,W,P residual norm [ 2.296648572656e-07, 2.859083731592e-07, 3.445118350945e-07, 3.031779272453e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 4.144000576836e+00 
      1 KSP Residual norm 5.697348022811e-01 
      2 KSP Residual norm 2.146984341904e-01 
      3 KSP Residual norm 6.335862258715e-02 
      4 KSP Residual norm 1.884338707685e-02 
  8 KSP Component U,V,W,P residual norm [ 1.749271185019e-07, 1.730905969531e-07, 1.974352587917e-07, 1.515961608259e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.369667656646e+00 
      1 KSP Residual norm 5.622601955432e-01 
      2 KSP Residual norm 1.933299274395e-01 
      3 KSP Residual norm 7.084173033950e-02 
      4 KSP Residual norm 3.198460072832e-02 
  9 KSP Component U,V,W,P residual norm [ 6.680085725879e-08, 8.669669712208e-08, 1.019941400321e-07, 9.438174888667e-07 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  4 SNES Function norm 5.864896219903e-04 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.88e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 3.748204292312e-04, 4.279599442316e-04, 1.425802155809e-04, 9.438174888667e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000337608927e+00 
      1 KSP Residual norm 2.174289700126e-01 
      2 KSP Residual norm 7.952045740061e-02 
      3 KSP Residual norm 3.545475415675e-02 
      4 KSP Residual norm 1.596829738166e-02 
      5 KSP Residual norm 7.254448170425e-03 
  1 KSP Component U,V,W,P residual norm [ 2.108509427729e-05, 2.423737004429e-05, 8.271212961376e-06, 1.349544276303e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.123481821829e+00 
      1 KSP Residual norm 1.605037000251e-01 
      2 KSP Residual norm 4.118780795124e-02 
      3 KSP Residual norm 1.697339380205e-02 
      4 KSP Residual norm 7.069694146714e-03 
  2 KSP Component U,V,W,P residual norm [ 4.357341111898e-06, 6.007151387230e-06, 6.264953228777e-06, 5.233884270774e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.815292337824e+00 
      1 KSP Residual norm 6.828176520039e-01 
      2 KSP Residual norm 2.022902062331e-01 
      3 KSP Residual norm 7.797883205953e-02 
      4 KSP Residual norm 2.811298725643e-02 
  3 KSP Component U,V,W,P residual norm [ 1.236117981060e-06, 1.564818668577e-06, 1.086258320957e-06, 2.764424941375e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.665230653242e+00 
      1 KSP Residual norm 5.773201378393e-01 
      2 KSP Residual norm 2.037314224824e-01 
      3 KSP Residual norm 7.062648767767e-02 
      4 KSP Residual norm 2.058115886743e-02 
  4 KSP Component U,V,W,P residual norm [ 6.068897977151e-07, 8.726269341677e-07, 1.009231913699e-06, 1.238298540184e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.906712521262e+00 
      1 KSP Residual norm 6.765599447074e-01 
      2 KSP Residual norm 2.779396096005e-01 
      3 KSP Residual norm 9.569275145768e-02 
      4 KSP Residual norm 3.816330417196e-02 
  5 KSP Component U,V,W,P residual norm [ 3.279458580563e-07, 3.874263167965e-07, 3.476633862774e-07, 5.101600005656e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.794511381433e+00 
      1 KSP Residual norm 3.473541963280e-01 
      2 KSP Residual norm 1.064914792753e-01 
      3 KSP Residual norm 3.634533578365e-02 
  6 KSP Component U,V,W,P residual norm [ 1.491958103412e-07, 1.301318550913e-07, 1.002600927180e-07, 3.358321106954e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.871998397110e+00 
      1 KSP Residual norm 3.526412398060e-01 
      2 KSP Residual norm 1.064561297819e-01 
      3 KSP Residual norm 3.498240106199e-02 
      4 KSP Residual norm 1.380643662074e-02 
  7 KSP Component U,V,W,P residual norm [ 1.268242920386e-07, 1.547564589265e-07, 1.892241964337e-07, 1.667286436806e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 4.101881955865e+00 
      1 KSP Residual norm 5.706355811853e-01 
      2 KSP Residual norm 2.142837568107e-01 
      3 KSP Residual norm 6.325010966828e-02 
      4 KSP Residual norm 1.892986121686e-02 
  8 KSP Component U,V,W,P residual norm [ 9.358005590436e-08, 9.187168866324e-08, 1.026979446212e-07, 8.419562635800e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.390936062425e+00 
      1 KSP Residual norm 5.570309868722e-01 
      2 KSP Residual norm 1.869283097704e-01 
      3 KSP Residual norm 6.696902681308e-02 
      4 KSP Residual norm 2.978689084518e-02 
  9 KSP Component U,V,W,P residual norm [ 3.482604369864e-08, 4.525198799353e-08, 5.425609794481e-08, 5.183929373087e-07 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  5 SNES Function norm 3.129309046637e-04 
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 5
SNES Object: 4 MPI processes
  type: newtonls
  maximum iterations=5, maximum function evaluations=10000
  tolerances: relative=1e-05, absolute=1e-08, solution=1e-08
  total number of linear solver iterations=47
  total number of function evaluations=6
  norm schedule ALWAYS
  SNESLineSearch Object: 4 MPI processes
    type: basic
    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: 4 MPI processes
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=2000, initial guess is zero
    tolerances:  relative=0.001, absolute=1e-08, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 4 MPI processes
    type: fieldsplit
      FieldSplit with Schur preconditioner, factorization UPPER
      Preconditioner for the Schur complement formed from A11
      Split info:
      Split number 0 Defined by IS
      Split number 1 Defined by IS
      KSP solver for A00 block
        KSP Object: (fieldsplit_u_) 4 MPI processes
          type: fgmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
          right preconditioning
          using UNPRECONDITIONED norm type for convergence test
        PC Object: (fieldsplit_u_) 4 MPI processes
          type: mg
            type is MULTIPLICATIVE, levels=2 cycles=v
              Cycles per PCApply=1
              Not using Galerkin computed coarse grid matrices
          Coarse grid solver -- level -------------------------------
            KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=4, initial guess is zero
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using UNPRECONDITIONED norm type for convergence test
            PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: bjacobi
                number of blocks = 4
                Local solve is same for all blocks, in the following KSP and PC objects:
              KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                        package used to perform factorization: petsc
                        total: nonzeros=126711, allocated nonzeros=126711
                        total number of mallocs used during MatSetValues calls =0
                          using I-node routines: found 225 nodes, limit used is 5
                linear system matrix = precond matrix:
                Mat Object: (Buu_) 1 MPI processes
                  type: seqaij
                  rows=675, cols=675, bs=3
                  total: nonzeros=126711, allocated nonzeros=126711
                  total number of mallocs used during MatSetValues calls =0
                    using I-node routines: found 225 nodes, limit used is 5
              linear system matrix = precond matrix:
              Mat Object: (Buu_) 4 MPI processes
                type: mpiaij
                rows=2187, cols=2187, bs=3
                total: nonzeros=533871, allocated nonzeros=533871
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
          Down solver (pre-smoother) on level 1 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=2, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: (fieldsplit_u_) 4 MPI processes
                type: shell
                rows=14739, cols=14739, bs=3
          Up solver (post-smoother) same as down solver (pre-smoother)
          linear system matrix = precond matrix:
          Mat Object: (fieldsplit_u_) 4 MPI processes
            type: shell
            rows=14739, cols=14739, bs=3
      KSP solver for S = A11 - A10 inv(A00) A01 
        KSP Object: (fieldsplit_p_) 4 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: (fieldsplit_p_) 4 MPI processes
          type: jacobi
          linear system matrix followed by preconditioner matrix:
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: schurcomplement
            rows=2048, cols=2048
              Schur complement A11 - A10 inv(A00) A01
              A11
                Mat Object: (fieldsplit_p_) 4 MPI processes
                  type: mpisbaij
                  rows=2048, cols=2048, bs=4
                  total: nonzeros=8192, allocated nonzeros=8192
                  total number of mallocs used during MatSetValues calls =0
              A10
                Mat Object: (Bpu_) 4 MPI processes
                  type: shell
                  rows=2048, cols=14739
              KSP of A00
                KSP Object: (fieldsplit_u_) 4 MPI processes
                  type: fgmres
                    restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                    happy breakdown tolerance 1e-30
                  maximum iterations=10, initial guess is zero
                  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
                  right preconditioning
                  using UNPRECONDITIONED norm type for convergence test
                PC Object: (fieldsplit_u_) 4 MPI processes
                  type: mg
                    type is MULTIPLICATIVE, levels=2 cycles=v
                      Cycles per PCApply=1
                      Not using Galerkin computed coarse grid matrices
                  Coarse grid solver -- level -------------------------------
                    KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=4, initial guess is zero
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using UNPRECONDITIONED norm type for convergence test
                    PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: bjacobi
                        number of blocks = 4
                        Local solve is same for all blocks, in the following KSP and PC objects:
                      KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                                package used to perform factorization: petsc
                                total: nonzeros=126711, allocated nonzeros=126711
                                total number of mallocs used during MatSetValues calls =0
                                  using I-node routines: found 225 nodes, limit used is 5
                        linear system matrix = precond matrix:
                        Mat Object: (Buu_) 1 MPI processes
                          type: seqaij
                          rows=675, cols=675, bs=3
                          total: nonzeros=126711, allocated nonzeros=126711
                          total number of mallocs used during MatSetValues calls =0
                            using I-node routines: found 225 nodes, limit used is 5
                      linear system matrix = precond matrix:
                      Mat Object: (Buu_) 4 MPI processes
                        type: mpiaij
                        rows=2187, cols=2187, bs=3
                        total: nonzeros=533871, allocated nonzeros=533871
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                  Down solver (pre-smoother) on level 1 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=2, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: (fieldsplit_u_) 4 MPI processes
                        type: shell
                        rows=14739, cols=14739, bs=3
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  linear system matrix = precond matrix:
                  Mat Object: (fieldsplit_u_) 4 MPI processes
                    type: shell
                    rows=14739, cols=14739, bs=3
              A01
                Mat Object: (Bup_) 4 MPI processes
                  type: shell
                  rows=14739, cols=2048
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: mpisbaij
            rows=2048, cols=2048, bs=4
            total: nonzeros=8192, allocated nonzeros=8192
            total number of mallocs used during MatSetValues calls =0
    linear system matrix followed by preconditioner matrix:
    Mat Object: (stokes_Amf_) 4 MPI processes
      type: shell
      rows=16787, cols=16787
    Mat Object: 4 MPI processes
      type: nest
      rows=16787, cols=16787
        Matrix object: 
          type=nest, rows=2, cols=2 
          MatNest structure: 
          (0,0) : prefix="fieldsplit_u_", type=shell, rows=14739, cols=14739 
          (0,1) : prefix="Bup_", type=shell, rows=14739, cols=2048 
          (1,0) : prefix="Bpu_", type=shell, rows=2048, cols=14739 
          (1,1) : prefix="fieldsplit_p_", type=mpisbaij, rows=2048, cols=2048 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.10e-02 (sec)
[[ModelOutput_SD3D]]
slab base range: gmin 6.9697e-01 5.0000e-01 6.0261e-22
slab base range: gmax 7.5758e-01 5.0000e-01 3.7879e-01
slab edge range: gmax(x) 6.9697e-01
slab edge identified: pidx_w,rank_w 5180 1
[step 1] slab edge range: W,Dn 6.9697e-01 8.7879e-01
slab edge (back face) range: gmax(x) 6.9697e-01
slab edge (back face) identified: pidx_w,rank_w 5155 1
[step 1] slab edge (back face) range: W,Dn 6.9697e-01 8.7879e-01
slab edge (front face) range: gmin(z) 3.7879e-01
slab edge (front face) identified: pidx_w,rank_w 5104 1
[step 1] slab edge (front face) range: W,Dn 3.7879e-01 8.7879e-01
topo range: 1.0000e+00 1.0000e+00
vol,int v.v,vrms: 5.7392e-01 1.4128e-02 1.5690e-01
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.20e-02 (sec)
Phi: 0.0000e+00
  TimeStep control(StkCourant): | current = 1.0000e+32 : trial = 3.1012e-02 [accepted] | dt > dt_max : restricted to = 2.3891e-03 | ==>> dt used = 2.3891e-03 |
  TimeStep control(StkSurfaceCourant): | current = 2.3891e-03 : trial = 1.0000e+32 | ==>> dt used = 2.3891e-03 |
  timestep_stokes[1] dt_courant = 2.3891e-03 
<<----------------------------------------------------------------------------------------------->>
   [[ EXECUTING TIME STEP : 2 ]]
     dt    : 2.3891e-03 
     time  : 2.3891e-03 
=========== MaterialPointStd_UpdateCoordinates ============= 
      Number of material points                     19860 (init) / 19860 (final)
    	MaterialPointStd_UpdateLocalCoordinates     1.59e-02 (sec) efficiency 40.2%
    	SwarmUpdatePosition_Communication_Generic   1.80e-02 (sec) efficiency 47.2%
    	MaterialPointStd_Removal                    2.80e-04 (sec) efficiency 3.9%
=========== ================================== ============= 
   [[ COMPUTING FLOW FIELD FOR STEP : 2 ]]
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.94e-02 (sec)
  0 SNES Function norm 1.198663390409e-03 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 4.355067200063e-04, 1.094883686344e-03, 2.199028535120e-04, 5.183929373087e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.999654994591e-01 
      1 KSP Residual norm 1.878147928072e-01 
      2 KSP Residual norm 4.316082648998e-02 
      3 KSP Residual norm 1.479838633500e-02 
      4 KSP Residual norm 4.749094954698e-03 
  1 KSP Component U,V,W,P residual norm [ 3.472673428988e-04, 8.730237147232e-04, 1.755301783718e-04, 4.817629181221e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.179476435100e-01 
      1 KSP Residual norm 2.767264750444e-02 
      2 KSP Residual norm 9.255034019856e-03 
      3 KSP Residual norm 3.236212235355e-03 
      4 KSP Residual norm 1.111653423595e-03 
  2 KSP Component U,V,W,P residual norm [ 1.631526418603e-04, 4.106053343333e-04, 8.200965108609e-05, 5.802095132239e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 7.702374494266e-01 
      1 KSP Residual norm 1.553582678458e-01 
      2 KSP Residual norm 5.275658842854e-02 
      3 KSP Residual norm 1.628464955204e-02 
      4 KSP Residual norm 5.125441092758e-03 
  3 KSP Component U,V,W,P residual norm [ 2.073612025467e-05, 5.144037183832e-05, 1.041508970681e-05, 2.532174761946e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.255701129799e+00 
      1 KSP Residual norm 2.372502385274e-01 
      2 KSP Residual norm 7.968044964182e-02 
      3 KSP Residual norm 2.319883290607e-02 
      4 KSP Residual norm 7.054771416327e-03 
  4 KSP Component U,V,W,P residual norm [ 7.167296288752e-06, 1.640989891679e-05, 5.240572480024e-06, 1.445819091515e-04 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.897345412302e+00 
      1 KSP Residual norm 5.761651951236e-01 
      2 KSP Residual norm 1.989156800442e-01 
      3 KSP Residual norm 5.605418348885e-02 
      4 KSP Residual norm 1.638063361324e-02 
  5 KSP Component U,V,W,P residual norm [ 2.489064956578e-06, 3.383194374192e-06, 2.768477154871e-06, 5.615469846819e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.411480436871e+00 
      1 KSP Residual norm 4.966884141885e-01 
      2 KSP Residual norm 1.599353842686e-01 
      3 KSP Residual norm 4.738703502064e-02 
      4 KSP Residual norm 1.377808120368e-02 
  6 KSP Component U,V,W,P residual norm [ 5.610456079894e-07, 6.712798333237e-07, 5.548047218783e-07, 2.645580893037e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.611338791395e-01 
      1 KSP Residual norm 1.807112894193e-01 
      2 KSP Residual norm 4.294071235686e-02 
      3 KSP Residual norm 1.392692203623e-02 
      4 KSP Residual norm 4.947966715191e-03 
  7 KSP Component U,V,W,P residual norm [ 4.764040350666e-07, 4.724436206597e-07, 4.904084939244e-07, 1.572229712753e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.657717590935e+00 
      1 KSP Residual norm 5.146346881304e-01 
      2 KSP Residual norm 1.740838069111e-01 
      3 KSP Residual norm 5.110242846753e-02 
      4 KSP Residual norm 1.347358353709e-02 
  8 KSP Component U,V,W,P residual norm [ 1.959069883025e-07, 2.662355167439e-07, 3.739235379685e-07, 6.249455337803e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.745147181619e+00 
      1 KSP Residual norm 3.117201043066e-01 
      2 KSP Residual norm 7.681665852565e-02 
      3 KSP Residual norm 2.261854511191e-02 
      4 KSP Residual norm 8.043806371006e-03 
  9 KSP Component U,V,W,P residual norm [ 1.030567747960e-07, 6.851665911360e-08, 9.744610182289e-08, 3.178192121753e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.220982893810e+00 
      1 KSP Residual norm 4.454017504530e-01 
      2 KSP Residual norm 1.125125898500e-01 
      3 KSP Residual norm 3.490615658791e-02 
      4 KSP Residual norm 1.071953975665e-02 
 10 KSP Component U,V,W,P residual norm [ 6.108266481562e-08, 8.328234247875e-08, 9.921625452384e-08, 1.427759593989e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.379010903471e+00 
      1 KSP Residual norm 3.306240924927e-01 
      2 KSP Residual norm 8.538786545990e-02 
      3 KSP Residual norm 2.273412386366e-02 
 11 KSP Component U,V,W,P residual norm [ 3.060480520198e-08, 4.645836558406e-08, 3.994381320546e-08, 8.029813126090e-07 ]
  Linear solve converged due to CONVERGED_RTOL iterations 11
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.02e-02 (sec)
  1 SNES Function norm 1.243645099618e-04 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.98e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 6.938362872082e-05, 9.619375099570e-05, 3.739733805983e-05, 8.029813126091e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 9.983757453406e-01 
      1 KSP Residual norm 1.857676345768e-01 
      2 KSP Residual norm 6.084042999353e-02 
      3 KSP Residual norm 1.782071960603e-02 
      4 KSP Residual norm 5.983508206576e-03 
  1 KSP Component U,V,W,P residual norm [ 1.605835137477e-06, 2.169215347428e-06, 8.827446962734e-07, 1.821122425727e-05 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.292535342463e+00 
      1 KSP Residual norm 3.404748374172e-01 
      2 KSP Residual norm 8.235828885308e-02 
      3 KSP Residual norm 3.710711350280e-02 
      4 KSP Residual norm 1.397643744266e-02 
  2 KSP Component U,V,W,P residual norm [ 4.035193707788e-07, 6.377316182337e-07, 3.818893796924e-07, 8.108403795159e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.696344258002e+00 
      1 KSP Residual norm 4.392272828753e-01 
      2 KSP Residual norm 1.115779582671e-01 
      3 KSP Residual norm 4.034889427071e-02 
      4 KSP Residual norm 1.484984809474e-02 
  3 KSP Component U,V,W,P residual norm [ 1.971514648368e-07, 2.079979377470e-07, 1.684666305697e-07, 4.959633422724e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.654306691572e+00 
      1 KSP Residual norm 5.212913793718e-01 
      2 KSP Residual norm 1.300807875822e-01 
      3 KSP Residual norm 4.652935242590e-02 
      4 KSP Residual norm 1.694724084368e-02 
  4 KSP Component U,V,W,P residual norm [ 7.339901076844e-08, 1.054147527558e-07, 7.184778843251e-08, 2.735493059989e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.091245521047e+00 
      1 KSP Residual norm 5.476833146150e-01 
      2 KSP Residual norm 1.793076938487e-01 
      3 KSP Residual norm 5.757384644020e-02 
      4 KSP Residual norm 1.975360965399e-02 
  5 KSP Component U,V,W,P residual norm [ 8.265336437963e-08, 1.130115185358e-07, 1.205668626272e-07, 1.337872335261e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.348676131414e+00 
      1 KSP Residual norm 4.255725927120e-01 
      2 KSP Residual norm 1.143826639878e-01 
      3 KSP Residual norm 3.294906864013e-02 
  6 KSP Component U,V,W,P residual norm [ 4.180048754574e-08, 3.542552537261e-08, 5.004271964844e-08, 7.593441536023e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.021805424233e+00 
      1 KSP Residual norm 4.756715276069e-01 
      2 KSP Residual norm 1.395755569200e-01 
      3 KSP Residual norm 4.452875402624e-02 
      4 KSP Residual norm 2.153127915700e-02 
  7 KSP Component U,V,W,P residual norm [ 1.535648157650e-08, 2.339367890527e-08, 2.836695005356e-08, 3.793930269395e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.052112730056e+00 
      1 KSP Residual norm 5.102067917699e-01 
      2 KSP Residual norm 1.493114912554e-01 
      3 KSP Residual norm 4.458424980829e-02 
      4 KSP Residual norm 1.759502533923e-02 
  8 KSP Component U,V,W,P residual norm [ 1.167497409910e-08, 1.801796272662e-08, 2.133624608978e-08, 1.980787096620e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.142703927287e+00 
      1 KSP Residual norm 4.960862549337e-01 
      2 KSP Residual norm 1.414135078774e-01 
      3 KSP Residual norm 4.710942411384e-02 
      4 KSP Residual norm 2.241258379340e-02 
  9 KSP Component U,V,W,P residual norm [ 5.030554445018e-09, 8.632856130726e-09, 9.438067318206e-09, 1.129166593470e-07 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  2 SNES Function norm 3.791041046837e-05 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.90e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.164851545821e-05, 2.828206528877e-05, 1.298664826528e-05, 1.129166593468e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000521580424e+00 
      1 KSP Residual norm 2.053237364498e-01 
      2 KSP Residual norm 6.773579945591e-02 
      3 KSP Residual norm 2.027996430859e-02 
      4 KSP Residual norm 6.095232980948e-03 
  1 KSP Component U,V,W,P residual norm [ 5.475646281052e-07, 7.011658477802e-07, 3.271884812136e-07, 5.831252186494e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.253164227466e+00 
      1 KSP Residual norm 3.394590121541e-01 
      2 KSP Residual norm 8.080079191324e-02 
      3 KSP Residual norm 2.682853934546e-02 
      4 KSP Residual norm 9.976137419235e-03 
  2 KSP Component U,V,W,P residual norm [ 1.500845255985e-07, 2.043840052682e-07, 1.057464624469e-07, 2.876394706155e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.711625352775e+00 
      1 KSP Residual norm 4.185189801215e-01 
      2 KSP Residual norm 1.078767701193e-01 
      3 KSP Residual norm 2.992589063896e-02 
      4 KSP Residual norm 1.065788691231e-02 
  3 KSP Component U,V,W,P residual norm [ 5.649884553387e-08, 7.087286509775e-08, 5.821457439699e-08, 1.695060870502e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.312066099926e+00 
      1 KSP Residual norm 4.575809469381e-01 
      2 KSP Residual norm 1.133515395012e-01 
      3 KSP Residual norm 3.936447829872e-02 
      4 KSP Residual norm 1.518510411477e-02 
  4 KSP Component U,V,W,P residual norm [ 2.138594445488e-08, 2.443909092113e-08, 1.922601636816e-08, 9.590477741772e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.543668796621e+00 
      1 KSP Residual norm 4.584581327300e-01 
      2 KSP Residual norm 1.433330750754e-01 
      3 KSP Residual norm 4.783698187245e-02 
      4 KSP Residual norm 1.819796084203e-02 
  5 KSP Component U,V,W,P residual norm [ 2.705247909840e-08, 3.791225476103e-08, 3.263821210000e-08, 4.576262564196e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.022259864038e+00 
      1 KSP Residual norm 4.467239261800e-01 
      2 KSP Residual norm 1.071133805975e-01 
      3 KSP Residual norm 2.711013331792e-02 
  6 KSP Component U,V,W,P residual norm [ 1.114007983412e-08, 8.536298024994e-09, 1.612567856027e-08, 2.438677252150e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.700968709575e+00 
      1 KSP Residual norm 4.119249491366e-01 
      2 KSP Residual norm 1.211140531549e-01 
      3 KSP Residual norm 4.567276523637e-02 
      4 KSP Residual norm 2.016971091599e-02 
  7 KSP Component U,V,W,P residual norm [ 3.181208562669e-09, 5.329656480307e-09, 5.184448271597e-09, 1.222622842540e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.701521166120e+00 
      1 KSP Residual norm 4.279447848709e-01 
      2 KSP Residual norm 1.304992486693e-01 
      3 KSP Residual norm 4.514989408932e-02 
      4 KSP Residual norm 1.925025526735e-02 
  8 KSP Component U,V,W,P residual norm [ 4.366669423258e-09, 5.943722938070e-09, 6.967575305443e-09, 6.280637309211e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.619005456593e+00 
      1 KSP Residual norm 4.298720243718e-01 
      2 KSP Residual norm 1.304237110974e-01 
      3 KSP Residual norm 5.016818606861e-02 
      4 KSP Residual norm 2.276374981551e-02 
  9 KSP Component U,V,W,P residual norm [ 1.294541275112e-09, 2.276824081969e-09, 1.982732848558e-09, 3.713808702581e-08 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.97e-02 (sec)
  3 SNES Function norm 1.327603524141e-05 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.28e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 7.671662163562e-06, 9.966094295976e-06, 4.251387611766e-06, 3.713808702569e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000455177552e+00 
      1 KSP Residual norm 2.109364862058e-01 
      2 KSP Residual norm 6.993548075928e-02 
      3 KSP Residual norm 2.130139118459e-02 
      4 KSP Residual norm 6.351622901283e-03 
  1 KSP Component U,V,W,P residual norm [ 1.926199399489e-07, 2.501723965911e-07, 1.075878530961e-07, 2.044246466422e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.216685097678e+00 
      1 KSP Residual norm 3.440716266818e-01 
      2 KSP Residual norm 8.178708519805e-02 
      3 KSP Residual norm 2.488615041450e-02 
      4 KSP Residual norm 8.613250836938e-03 
  2 KSP Component U,V,W,P residual norm [ 6.022852749131e-08, 7.520242545331e-08, 3.964359981548e-08, 1.064122729534e-06 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.825979174680e+00 
      1 KSP Residual norm 4.417183895927e-01 
      2 KSP Residual norm 1.150422739172e-01 
      3 KSP Residual norm 3.022483243184e-02 
      4 KSP Residual norm 9.617835406645e-03 
  3 KSP Component U,V,W,P residual norm [ 2.062961279795e-08, 2.399213581703e-08, 2.035579595767e-08, 6.033525625907e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.169035971726e+00 
      1 KSP Residual norm 4.528127433104e-01 
      2 KSP Residual norm 1.106000753859e-01 
      3 KSP Residual norm 3.816155556565e-02 
      4 KSP Residual norm 1.372912793278e-02 
  4 KSP Component U,V,W,P residual norm [ 6.970096356843e-09, 8.356902001121e-09, 6.630173053478e-09, 3.294424170443e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.558135312357e+00 
      1 KSP Residual norm 4.484146579065e-01 
      2 KSP Residual norm 1.348466599563e-01 
      3 KSP Residual norm 4.667722202217e-02 
      4 KSP Residual norm 1.714527891893e-02 
  5 KSP Component U,V,W,P residual norm [ 9.399935155403e-09, 1.065790208981e-08, 1.046851494240e-08, 1.598431190622e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.837804768980e+00 
      1 KSP Residual norm 4.549628826284e-01 
      2 KSP Residual norm 1.054942841604e-01 
      3 KSP Residual norm 2.745891249446e-02 
  6 KSP Component U,V,W,P residual norm [ 3.885372669031e-09, 3.469616195597e-09, 5.197014820990e-09, 8.541602804828e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.717537837949e+00 
      1 KSP Residual norm 3.918318103463e-01 
      2 KSP Residual norm 1.107418588130e-01 
      3 KSP Residual norm 3.666088453049e-02 
      4 KSP Residual norm 1.619977834942e-02 
  7 KSP Component U,V,W,P residual norm [ 9.992962308735e-10, 2.056427736760e-09, 1.719903215536e-09, 4.456383237968e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.431502660391e+00 
      1 KSP Residual norm 4.181725339425e-01 
      2 KSP Residual norm 1.136083748449e-01 
      3 KSP Residual norm 3.495550503143e-02 
      4 KSP Residual norm 1.400045594653e-02 
  8 KSP Component U,V,W,P residual norm [ 1.543550495536e-09, 1.131917067202e-09, 2.223570469700e-09, 2.213185854149e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.756591849838e+00 
      1 KSP Residual norm 4.242390010753e-01 
      2 KSP Residual norm 1.279764345506e-01 
      3 KSP Residual norm 5.170902369315e-02 
      4 KSP Residual norm 2.249383949810e-02 
  9 KSP Component U,V,W,P residual norm [ 4.177976497491e-10, 8.886352013919e-10, 7.266342894637e-10, 1.237458594890e-08 ]
  Linear solve converged due to CONVERGED_RTOL iterations 9
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.92e-02 (sec)
  4 SNES Function norm 4.780560007257e-06 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.00e-02 (sec)
  0 KSP Component U,V,W,P residual norm [ 2.888961444287e-06, 3.591835300007e-06, 1.267368061636e-06, 1.237458594910e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.000770713462e+00 
      1 KSP Residual norm 2.158613361914e-01 
      2 KSP Residual norm 7.437397449985e-02 
      3 KSP Residual norm 2.350749148453e-02 
      4 KSP Residual norm 7.479638916087e-03 
  1 KSP Component U,V,W,P residual norm [ 7.960748546636e-08, 1.015003899468e-07, 3.627047688738e-08, 7.748246430575e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.938642801302e+00 
      1 KSP Residual norm 2.901880744173e-01 
      2 KSP Residual norm 6.794415582860e-02 
      3 KSP Residual norm 1.969796697144e-02 
      4 KSP Residual norm 6.426139451349e-03 
  2 KSP Component U,V,W,P residual norm [ 2.940373689577e-08, 3.706307811036e-08, 2.022219594791e-08, 4.329636542199e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.032862303048e+00 
      1 KSP Residual norm 4.908044379193e-01 
      2 KSP Residual norm 1.280742241616e-01 
      3 KSP Residual norm 3.252962260417e-02 
      4 KSP Residual norm 8.908240900665e-03 
  3 KSP Component U,V,W,P residual norm [ 7.795078107531e-09, 8.821090949012e-09, 5.102584641801e-09, 2.198384148096e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 3.067714470839e+00 
      1 KSP Residual norm 4.305892342276e-01 
      2 KSP Residual norm 1.263028496674e-01 
      3 KSP Residual norm 4.553892070592e-02 
      4 KSP Residual norm 1.612036351475e-02 
  4 KSP Component U,V,W,P residual norm [ 2.776043348937e-09, 4.973682414315e-09, 3.598210120893e-09, 1.156849949407e-07 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.719327151908e+00 
      1 KSP Residual norm 4.886851433303e-01 
      2 KSP Residual norm 1.545464072829e-01 
      3 KSP Residual norm 5.276431721090e-02 
      4 KSP Residual norm 1.866572893033e-02 
  5 KSP Component U,V,W,P residual norm [ 3.409683853410e-09, 3.755289810139e-09, 3.612982999572e-09, 5.557774078852e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.709325960400e+00 
      1 KSP Residual norm 4.454652189151e-01 
      2 KSP Residual norm 1.060356662725e-01 
      3 KSP Residual norm 3.201566064873e-02 
      4 KSP Residual norm 1.204318616481e-02 
  6 KSP Component U,V,W,P residual norm [ 7.453931648218e-10, 1.368541809533e-09, 8.813147443631e-10, 3.214745019841e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.710478254070e+00 
      1 KSP Residual norm 4.298562802429e-01 
      2 KSP Residual norm 1.327647829069e-01 
      3 KSP Residual norm 4.370705212513e-02 
      4 KSP Residual norm 1.957133536713e-02 
  7 KSP Component U,V,W,P residual norm [ 5.220859392759e-10, 1.167719461642e-09, 1.002929081798e-09, 1.663866055375e-08 ]
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 2.605667746002e+00 
      1 KSP Residual norm 4.783918096071e-01 
      2 KSP Residual norm 1.410185704058e-01 
      3 KSP Residual norm 3.914816841582e-02 
      4 KSP Residual norm 1.478496016680e-02 
  8 KSP Component U,V,W,P residual norm [ 5.268803751182e-10, 3.863059373787e-10, 5.373477247364e-10, 7.657923547424e-09 ]
  Linear solve converged due to CONVERGED_ATOL iterations 8
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.91e-02 (sec)
  5 SNES Function norm 1.899032790192e-06 
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 5
SNES Object: 4 MPI processes
  type: newtonls
  maximum iterations=5, maximum function evaluations=10000
  tolerances: relative=1e-05, absolute=1e-08, solution=1e-08
  total number of linear solver iterations=46
  total number of function evaluations=6
  norm schedule ALWAYS
  SNESLineSearch Object: 4 MPI processes
    type: basic
    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: 4 MPI processes
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=2000, initial guess is zero
    tolerances:  relative=0.001, absolute=1e-08, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 4 MPI processes
    type: fieldsplit
      FieldSplit with Schur preconditioner, factorization UPPER
      Preconditioner for the Schur complement formed from A11
      Split info:
      Split number 0 Defined by IS
      Split number 1 Defined by IS
      KSP solver for A00 block
        KSP Object: (fieldsplit_u_) 4 MPI processes
          type: fgmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            happy breakdown tolerance 1e-30
          maximum iterations=10, initial guess is zero
          tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
          right preconditioning
          using UNPRECONDITIONED norm type for convergence test
        PC Object: (fieldsplit_u_) 4 MPI processes
          type: mg
            type is MULTIPLICATIVE, levels=2 cycles=v
              Cycles per PCApply=1
              Not using Galerkin computed coarse grid matrices
          Coarse grid solver -- level -------------------------------
            KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=4, initial guess is zero
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using UNPRECONDITIONED norm type for convergence test
            PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
              type: bjacobi
                number of blocks = 4
                Local solve is same for all blocks, in the following KSP and PC objects:
              KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                        package used to perform factorization: petsc
                        total: nonzeros=126711, allocated nonzeros=126711
                        total number of mallocs used during MatSetValues calls =0
                          using I-node routines: found 225 nodes, limit used is 5
                linear system matrix = precond matrix:
                Mat Object: (Buu_) 1 MPI processes
                  type: seqaij
                  rows=675, cols=675, bs=3
                  total: nonzeros=126711, allocated nonzeros=126711
                  total number of mallocs used during MatSetValues calls =0
                    using I-node routines: found 225 nodes, limit used is 5
              linear system matrix = precond matrix:
              Mat Object: (Buu_) 4 MPI processes
                type: mpiaij
                rows=2187, cols=2187, bs=3
                total: nonzeros=533871, allocated nonzeros=533871
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
          Down solver (pre-smoother) on level 1 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: gmres
                restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                happy breakdown tolerance 1e-30
              maximum iterations=2, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              right preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: (fieldsplit_u_) 4 MPI processes
                type: shell
                rows=14739, cols=14739, bs=3
          Up solver (post-smoother) same as down solver (pre-smoother)
          linear system matrix = precond matrix:
          Mat Object: (fieldsplit_u_) 4 MPI processes
            type: shell
            rows=14739, cols=14739, bs=3
      KSP solver for S = A11 - A10 inv(A00) A01 
        KSP Object: (fieldsplit_p_) 4 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: (fieldsplit_p_) 4 MPI processes
          type: jacobi
          linear system matrix followed by preconditioner matrix:
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: schurcomplement
            rows=2048, cols=2048
              Schur complement A11 - A10 inv(A00) A01
              A11
                Mat Object: (fieldsplit_p_) 4 MPI processes
                  type: mpisbaij
                  rows=2048, cols=2048, bs=4
                  total: nonzeros=8192, allocated nonzeros=8192
                  total number of mallocs used during MatSetValues calls =0
              A10
                Mat Object: (Bpu_) 4 MPI processes
                  type: shell
                  rows=2048, cols=14739
              KSP of A00
                KSP Object: (fieldsplit_u_) 4 MPI processes
                  type: fgmres
                    restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                    happy breakdown tolerance 1e-30
                  maximum iterations=10, initial guess is zero
                  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
                  right preconditioning
                  using UNPRECONDITIONED norm type for convergence test
                PC Object: (fieldsplit_u_) 4 MPI processes
                  type: mg
                    type is MULTIPLICATIVE, levels=2 cycles=v
                      Cycles per PCApply=1
                      Not using Galerkin computed coarse grid matrices
                  Coarse grid solver -- level -------------------------------
                    KSP Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=4, initial guess is zero
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using UNPRECONDITIONED norm type for convergence test
                    PC Object: (fieldsplit_u_mg_coarse_) 4 MPI processes
                      type: bjacobi
                        number of blocks = 4
                        Local solve is same for all blocks, in the following KSP and PC objects:
                      KSP Object: (fieldsplit_u_mg_coarse_sub_) 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: (fieldsplit_u_mg_coarse_sub_) 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=675, cols=675, bs=3
                                package used to perform factorization: petsc
                                total: nonzeros=126711, allocated nonzeros=126711
                                total number of mallocs used during MatSetValues calls =0
                                  using I-node routines: found 225 nodes, limit used is 5
                        linear system matrix = precond matrix:
                        Mat Object: (Buu_) 1 MPI processes
                          type: seqaij
                          rows=675, cols=675, bs=3
                          total: nonzeros=126711, allocated nonzeros=126711
                          total number of mallocs used during MatSetValues calls =0
                            using I-node routines: found 225 nodes, limit used is 5
                      linear system matrix = precond matrix:
                      Mat Object: (Buu_) 4 MPI processes
                        type: mpiaij
                        rows=2187, cols=2187, bs=3
                        total: nonzeros=533871, allocated nonzeros=533871
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                  Down solver (pre-smoother) on level 1 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: gmres
                        restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                        happy breakdown tolerance 1e-30
                      maximum iterations=2, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      right preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_1_) 4 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: (fieldsplit_u_) 4 MPI processes
                        type: shell
                        rows=14739, cols=14739, bs=3
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  linear system matrix = precond matrix:
                  Mat Object: (fieldsplit_u_) 4 MPI processes
                    type: shell
                    rows=14739, cols=14739, bs=3
              A01
                Mat Object: (Bup_) 4 MPI processes
                  type: shell
                  rows=14739, cols=2048
          Mat Object: (fieldsplit_p_) 4 MPI processes
            type: mpisbaij
            rows=2048, cols=2048, bs=4
            total: nonzeros=8192, allocated nonzeros=8192
            total number of mallocs used during MatSetValues calls =0
    linear system matrix followed by preconditioner matrix:
    Mat Object: (stokes_Amf_) 4 MPI processes
      type: shell
      rows=16787, cols=16787
    Mat Object: 4 MPI processes
      type: nest
      rows=16787, cols=16787
        Matrix object: 
          type=nest, rows=2, cols=2 
          MatNest structure: 
          (0,0) : prefix="fieldsplit_u_", type=shell, rows=14739, cols=14739 
          (0,1) : prefix="Bup_", type=shell, rows=14739, cols=2048 
          (1,0) : prefix="Bpu_", type=shell, rows=2048, cols=14739 
          (1,1) : prefix="fieldsplit_p_", type=mpisbaij, rows=2048, cols=2048 
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 1.91e-02 (sec)
[[ModelOutput_SD3D]]
slab base range: gmin 6.9697e-01 4.9902e-01 -7.0740e-13
slab base range: gmax 7.5758e-01 4.9920e-01 3.7878e-01
slab edge range: gmax(x) 6.9715e-01
slab edge identified: pidx_w,rank_w 5180 1
[step 2] slab edge range: W,Dn 6.9715e-01 8.7815e-01
slab edge (back face) range: gmax(x) 6.9715e-01
slab edge (back face) identified: pidx_w,rank_w 5155 1
[step 2] slab edge (back face) range: W,Dn 6.9715e-01 8.7815e-01
slab edge (front face) range: gmin(z) 3.7847e-01
slab edge (front face) identified: pidx_w,rank_w 5103 1
[step 2] slab edge (front face) range: W,Dn 3.7847e-01 8.7822e-01
topo range: 1.0000e+00 1.0000e+00
vol,int v.v,vrms: 5.7392e-01 1.5117e-02 1.6230e-01
Update non-linearities (VPSTD) [mpoint]: (min,max)_eta 1.00e-03,2.67e+00; log10(max/min) 3.43e+00; npoints_yielded 0; cpu time 2.03e-02 (sec)
Phi: 0.0000e+00
  TimeStep control(StkCourant): | current = 1.0000e+32 : trial = 3.1164e-02 [accepted] | dt > dt_max : restricted to = 2.3891e-03 | ==>> dt used = 2.3891e-03 |
  TimeStep control(StkSurfaceCourant): | current = 2.3891e-03 : trial = 1.0000e+32 | ==>> dt used = 2.3891e-03 |
  timestep_stokes[2] dt_courant = 2.3891e-03 
pTatin3dMaxMemoryUsage = [0.0000e+00 , 0.0000e+00] (MB) 
[[ModelDestroy_Template]]
[[ModelDestroy_ViscousSinker]]
[[ModelDestroy_Gene3D]]
[[ModelDestroy_Indentor]]
[[ModelDestroy_Rift3D_T]]
[[ModelDestroy_AdvDiffExample]]
[[ModelDestroy_Delamination]]
[[ModelDestroy_Riftrh]]
[[ModelDestroy_Rift_oblique3d]]
[[ModelDestroy_MultilayerFolding]]
[[ModelDestroy_SubmarineLavaFlow]]
[[ModelDestroy_ExSubduction]]
[[ModelDestroy_iPLUS]]
[[ModelDestroy_Subduction_Initiation2d]]
[[ModelDestroy_Thermal_Convection2d]]
************************************************************************************************************************
***             WIDEN YOUR WINDOW TO 120 CHARACTERS.  Use 'enscript -r -fCourier9' to print this document            ***
************************************************************************************************************************

---------------------------------------------- PETSc Performance Summary: ----------------------------------------------

/Users/patrick/ptatin3d/arch-darwin-maint-double-opt/bin/ptatin_driver_nonlinear_ts.app on a arch-darwin-maint-double-opt named Windbauer.local with 4 processors, by patrick Thu Mar 15 08:18:03 2018
Using Petsc Release Version 3.8.3, unknown 

                         Max       Max/Min        Avg      Total 
Time (sec):           5.624e+01      1.00001   5.624e+01
Objects:              3.506e+03      1.00000   3.506e+03
Flop:                 1.634e+10      1.05802   1.586e+10  6.345e+10
Flop/sec:            2.906e+08      1.05802   2.820e+08  1.128e+09
MPI Messages:         4.479e+04      1.01822   4.437e+04  1.775e+05
MPI Message Lengths:  1.378e+08      1.12319   2.930e+03  5.201e+08
MPI Reductions:       4.460e+04      1.00000

Flop counting convention: 1 flop = 1 real number operation of type (multiply/divide/add/subtract)
                            e.g., VecAXPY() for real vectors of length N --> 2N flop
                            and VecAXPY() for complex vectors of length N --> 8N flop

Summary of Stages:   ----- Time ------  ----- Flop -----  --- Messages ---  -- Message Lengths --  -- Reductions --
                        Avg     %Total     Avg     %Total   counts   %Total     Avg         %Total   counts   %Total 
 0:      Main Stage: 1.4090e+01  25.1%  2.1248e+10  33.5%  4.210e+04  23.7%  7.491e+02       25.6%  7.531e+03  16.9% 
 1:        MG Apply: 4.0455e+01  71.9%  4.2206e+10  66.5%  1.354e+05  76.3%  2.181e+03       74.4%  1.207e+04  27.1% 

------------------------------------------------------------------------------------------------------------------------
See the 'Profiling' chapter of the users' manual for details on interpreting output.
Phase summary info:
   Count: number of times phase was executed
   Time and Flop: Max - maximum over all processors
                   Ratio - ratio of maximum to minimum over all processors
   Mess: number of messages sent
   Avg. len: average message length (bytes)
   Reduct: number of global reductions
   Global: entire computation
   Stage: stages of a computation. Set stages with PetscLogStagePush() and PetscLogStagePop().
      %T - percent time in this phase         %F - percent flop in this phase
      %M - percent messages in this phase     %L - percent message lengths in this phase
      %R - percent reductions in this phase
   Total Mflop/s: 10e-6 * (sum of flop over all processors)/(max time over all processors)
------------------------------------------------------------------------------------------------------------------------
Event                Count      Time (sec)     Flop                             --- Global ---  --- Stage ---   Total
                   Max Ratio  Max     Ratio   Max  Ratio  Mess   Avg len Reduct  %T %F %M %L %R  %T %F %M %L %R Mflop/s
------------------------------------------------------------------------------------------------------------------------

--- Event Stage 0: Main Stage

MatMultMFA11         670 1.0 5.7589e+00 1.0 1.48e+09 1.0 1.6e+04 3.4e+03 0.0e+00 10  9  9 11  0  41 28 38 42  0  1028
MatMultMFA           356 1.0 2.6307e+00 1.0 2.71e+09 1.0 8.5e+03 3.4e+03 0.0e+00  5 17  5  6  0  19 51 20 22  0  4127
MatMultMFA12         167 1.0 9.7993e-01 1.0 7.47e+08 1.0 2.0e+03 3.4e+03 0.0e+00  2  5  1  1  0   7 14  5  5  0  3049
DataExTopoSetup        1 1.0 2.6650e-03 1.1 0.00e+00 0.0 6.0e+01 5.6e+00 1.5e+01  0  0  0  0  0   0  0  0  0  0     0
DataExBegin            2 1.0 7.2956e-05 2.0 0.00e+00 0.0 1.2e+01 6.4e+03 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
DataExEnd              2 1.0 1.2929e-0315.4 0.00e+00 0.0 1.2e+01 6.4e+03 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelInit              1 1.0 5.3906e-04 3.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelInitSoln          1 1.0 0.0000e+00 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelInitMesh          1 1.0 2.1386e-04 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 2.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelInitMat           1 1.0 7.2055e-02 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelInitStkVar        1 1.0 0.0000e+00 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelBC                3 1.0 6.8581e-03 1.0 0.00e+00 0.0 5.8e+02 3.4e+03 0.0e+00  0  0  0  0  0   0  0  1  1  0     0
ModelBCMG              3 1.0 1.0354e-02 1.0 0.00e+00 0.0 1.2e+03 2.3e+03 0.0e+00  0  0  1  0  0   0  0  3  2  0     0
ModelMatBC             2 1.0 1.1921e-06 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelUpdateMesh        2 1.0 9.5367e-07 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
ModelOutput            5 1.0 2.9416e-01 1.0 3.33e+07 1.0 9.6e+02 2.2e+03 3.5e+02  1  0  1  0  1   2  1  2  2  5   453
MPAdvGCoord            2 1.0 1.1537e-02 1.9 0.00e+00 0.0 2.4e+01 3.4e+03 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MPAdvLCoord            2 1.0 2.1045e-02 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MPAdvComm              2 1.0 8.5580e-03 2.1 0.00e+00 0.0 4.8e+01 3.2e+03 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MPAdvRemove            2 1.0 3.0398e-0415.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MPPCInsert             2 1.0 1.0152e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 2.0e+00  0  0  0  0  0   0  0  0  0  0     0
VecMax                24 1.0 3.9694e-03 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 2.4e+01  0  0  0  0  0   0  0  0  0  0     0
VecMin                24 1.0 4.9632e-03 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 2.4e+01  0  0  0  0  0   0  0  0  0  0     0
VecMDot              846 1.0 2.1884e-01 1.8 2.34e+07 1.3 0.0e+00 0.0e+00 8.5e+02  0  0  0  0  2   1  0  0  0 11   384
VecNorm             1254 1.0 2.4684e-01 1.0 8.98e+06 1.3 0.0e+00 0.0e+00 1.3e+03  0  0  0  0  3   2  0  0  0 17   130
VecScale            1196 1.0 4.6833e-03 1.3 5.00e+06 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  3817
VecCopy              463 1.0 1.7769e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
VecSet              4584 1.0 1.1258e-02 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
VecAXPY              350 1.0 2.3184e-03 1.2 3.08e+06 1.2 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  4774
VecAYPX              183 1.0 3.2256e-03 1.3 8.50e+05 1.2 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   952
VecWAXPY              22 1.0 5.0688e-04 1.3 1.02e+05 1.2 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   729
VecMAXPY            1196 1.0 5.0111e-02 1.2 3.98e+07 1.2 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  1  0  0  0  2850
VecPointwiseMult     203 1.0 6.6447e-04 1.2 9.36e+04 1.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   554
VecScatterBegin     4695 1.0 1.5975e-01 1.1 0.00e+00 0.0 3.6e+04 3.1e+03 0.0e+00  0  0 20 22  0   1  0 86 85  0     0
VecScatterEnd       4695 1.0 7.0509e-01 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  1  0  0  0  0   4  0  0  0  0     0
VecReduceArith        48 1.0 3.2830e-04 1.2 4.46e+05 1.2 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  4908
VecReduceComm         16 1.0 3.1421e-03 1.8 0.00e+00 0.0 0.0e+00 0.0e+00 1.6e+01  0  0  0  0  0   0  0  0  0  0     0
VecNormalize           9 1.0 1.6849e-03 1.1 1.82e+04 1.6 0.0e+00 0.0e+00 9.0e+00  0  0  0  0  0   0  0  0  0  0    35
MatMult             1196 1.0 9.3658e+00 1.0 4.94e+09 1.0 2.7e+04 3.4e+03 0.0e+00 17 31 15 18  0  66 93 63 69  0  2110
MatMultTranspose      54 1.0 3.9477e-03 1.1 4.56e+05 1.2 2.7e+02 2.7e+02 0.0e+00  0  0  0  0  0   0  0  1  0  0   427
MatAssemblyBegin     114 1.0 1.6087e-01 2.3 0.00e+00 0.0 1.0e+03 1.7e+04 3.0e+02  0  0  1  3  1   1  0  2 13  4     0
MatAssemblyEnd       187 1.0 7.1136e-02 2.8 0.00e+00 0.0 1.1e+02 2.8e+02 7.0e+01  0  0  0  0  0   0  0  0  0  1     0
MatGetRow              1 1.0 6.1989e-06 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MatCreateSubMat      108 1.0 3.1543e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 1.5e+02  0  0  0  0  0   0  0  0  0  2     0
MatZeroEntries        35 1.0 2.2957e-03 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MatView               36 1.3 1.3487e-02 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 1.6e+01  0  0  0  0  0   0  0  0  0  0     0
DMCoarsen              1 1.0 4.7350e-03 1.0 0.00e+00 0.0 5.6e+01 4.3e+02 3.7e+01  0  0  0  0  0   0  0  0  0  0     0
DMCreateInterp         1 1.0 2.4972e-03 1.0 0.00e+00 0.0 2.0e+01 7.0e+01 2.1e+01  0  0  0  0  0   0  0  0  0  0     0
Warning -- total time of event greater than time of entire stage -- something is wrong with the timer
SNESSolve              4 1.0 5.3511e+01 1.0 1.63e+10 1.1 1.7e+05 2.9e+03 1.8e+04 95100 97 98 41 380297409382242  1181
SNESFunctionEval      24 1.0 1.0074e+00 1.0 1.60e+08 1.0 3.7e+03 1.9e+03 8.6e+02  2  1  2  1  2   7  3  9  5 11   635
SNESJacobianEval      16 1.0 1.7295e+00 1.0 6.04e+07 1.0 7.0e+03 3.1e+03 1.9e+03  3  0  4  4  4  12  1 17 17 25   140
SNESLineSearch        16 1.0 6.9398e-01 1.0 1.07e+08 1.0 2.5e+03 1.9e+03 5.9e+02  1  1  1  1  1   5  2  6  3  8   617
KSPGMRESOrthog       837 1.0 2.4484e-01 1.6 4.67e+07 1.3 0.0e+00 0.0e+00 8.4e+02  0  0  0  0  2   1  1  0  0 11   685
KSPSetUp              83 1.0 4.5462e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 2.4e+01  0  0  0  0  0   0  0  0  0  0     0
Warning -- total time of event greater than time of entire stage -- something is wrong with the timer
KSPSolve              16 1.0 5.0930e+01 1.0 1.61e+10 1.1 1.6e+05 3.0e+03 1.6e+04 91 98 91 92 35 361294385361207  1226
PCSetUp               51 1.0 4.3385e-01 1.0 4.33e+07 1.0 2.4e+02 3.9e+03 3.1e+02  1  0  0  0  1   3  1  1  1  4   399
Warning -- total time of event greater than time of entire stage -- something is wrong with the timer
PCApply              167 1.0 4.7915e+01 1.0 1.33e+10 1.1 1.5e+05 2.9e+03 1.4e+04 85 81 87 87 31 340242365338181  1074
Warning -- total time of event greater than time of entire stage -- something is wrong with the timer
KSPSolve_FS_0        167 1.0 4.6921e+01 1.0 1.26e+10 1.1 1.5e+05 2.9e+03 1.4e+04 83 76 85 85 30 333228360333181  1033
KSPSolve_FS_Schu     167 1.0 1.3466e-02 1.1 8.55e+04 1.0 0.0e+00 0.0e+00 6.0e+00  0  0  0  0  0   0  0  0  0  0    25
MGSetup Level 0       16 1.0 1.6086e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 6.0e+00  0  0  0  0  0   0  0  0  0  0     0
MGSetup Level 1       16 1.0 3.7569e-01 1.0 4.33e+07 1.0 1.6e+02 3.4e+03 6.0e+00  1  0  0  0  0   3  1  0  0  0   461

--- Event Stage 1: MG Apply

MatMultMFA11        4020 1.0 3.4692e+01 1.0 8.88e+09 1.0 9.6e+04 3.4e+03 0.0e+00 61 56 54 64  0  85 84 71 86  0  1024
VecMDot             5360 1.0 1.6606e+00 4.0 4.22e+07 1.3 0.0e+00 0.0e+00 5.4e+03  2  0  0  0 12   3  0  0  0 44    89
VecNorm             6700 1.0 1.5562e+00 1.3 3.22e+07 1.3 0.0e+00 0.0e+00 6.7e+03  2  0  0  0 15   3  0  0  0 56    73
VecScale            6775 1.0 1.5752e-02 1.2 1.64e+07 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  3670
VecCopy             6030 1.0 2.0100e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
VecSet             16756 1.0 4.0124e-02 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
VecAXPY             2680 1.0 1.0587e-02 1.2 1.75e+07 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  5873
VecAYPX              670 1.0 9.1031e-03 1.2 2.77e+06 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  1085
VecMAXPY            7370 1.0 1.0088e-01 1.3 6.80e+07 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  1  0  0  0  2364
VecPointwiseMult    4020 1.0 6.4389e-02 1.1 1.66e+07 1.3 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   920
VecScatterBegin    12063 1.0 5.4486e-01 1.1 0.00e+00 0.0 1.4e+05 2.9e+03 0.0e+00  1  0 76 74  0   1  0100100  0     0
VecScatterEnd      12063 1.0 1.8027e+00 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  3  0  0  0  0   4  0  0  0  0     0
VecNormalize        7370 1.0 1.5820e+00 1.3 4.86e+07 1.3 0.0e+00 0.0e+00 6.7e+03  2  0  0  0 15   3  0  0  0 56   108
MatMult             6700 1.0 3.6230e+01 1.0 9.79e+09 1.0 1.3e+05 3.0e+03 0.0e+00 64 60 72 74  0  89 91 95 99  0  1059
MatMultAdd           670 1.0 1.0687e-01 1.4 1.70e+07 1.2 3.4e+03 5.5e+02 0.0e+00  0  0  2  0  0   0  0  2  0  0   588
MatMultTranspose     670 1.0 9.5607e-02 1.5 1.70e+07 1.2 3.4e+03 5.5e+02 0.0e+00  0  0  2  0  0   0  0  2  0  0   657
MatSolve            3350 1.0 1.3614e+00 1.8 8.47e+08 1.8 0.0e+00 0.0e+00 0.0e+00  2  4  0  0  0   3  6  0  0  0  1876
MatLUFactorNum        16 1.0 2.7202e-01 2.0 1.49e+08 2.0 0.0e+00 0.0e+00 0.0e+00  0  1  0  0  0   0  1  0  0  0  1586
MatILUFactorSym        3 1.0 5.2209e-03 2.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MatResidual          670 1.0 5.8093e+00 1.0 1.48e+09 1.0 1.6e+04 3.4e+03 0.0e+00 10  9  9 11  0  14 14 12 14  0  1020
MatGetRowIJ            3 1.0 5.2452e-06 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
MatGetOrdering         3 1.0 2.6870e-04 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
KSPGMRESOrthog      5360 1.0 1.7324e+00 3.5 8.45e+07 1.3 0.0e+00 0.0e+00 5.4e+03  2  0  0  0 12   3  1  0  0 44   171
KSPSetUp              16 1.0 3.0041e-05 3.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0
KSPSolve            2010 1.0 3.4401e+01 1.0 9.51e+09 1.1 1.1e+05 2.9e+03 1.2e+04 61 57 63 63 27  85 86 83 85100  1051
PCSetUp               16 1.0 2.7783e-01 2.0 1.49e+08 2.0 0.0e+00 0.0e+00 0.0e+00  0  1  0  0  0   1  1  0  0  0  1553
PCSetUpOnBlocks      670 1.0 2.7896e-01 2.0 1.49e+08 2.0 0.0e+00 0.0e+00 0.0e+00  0  1  0  0  0   1  1  0  0  0  1547
PCApply             7370 1.0 1.7380e+00 1.6 8.73e+08 1.8 3.6e+01 3.4e+03 6.0e+00  2  4  0  0  0   3  6  0  0  0  1526
MGSmooth Level 0     670 1.0 4.1760e+00 1.0 1.93e+09 1.8 3.2e+04 1.6e+03 6.0e+03  7  9 18 10 14  10 14 24 13 50  1421
MGSmooth Level 1    1340 1.0 3.0384e+01 1.0 7.57e+09 1.0 8.0e+04 3.4e+03 6.0e+03 54 48 45 53 14  75 72 59 72 50   994
MGResid Level 1      670 1.0 5.8101e+00 1.0 1.48e+09 1.0 1.6e+04 3.4e+03 0.0e+00 10  9  9 11  0  14 14 12 14  0  1020
MGInterp Level 1    1340 1.0 1.9334e-01 1.1 3.40e+07 1.2 6.7e+03 5.5e+02 0.0e+00  0  0  4  1  0   0  0  5  1  0   650
------------------------------------------------------------------------------------------------------------------------

Memory usage is given in bytes:

Object Type          Creations   Destructions     Memory  Descendants' Mem.
Reports information only for process 0.

--- Event Stage 0: Main Stage

           Container     4              4         2272     0.
              Viewer    15             13        10816     0.
              Vector  1734           1749     37314312     0.
      Vector Scatter   201            201      1960168     0.
              Matrix    49             52     11102652     0.
   Matrix Null Space     3              3         2064     0.
    Distributed Mesh   189            189       974808     0.
           Index Set   466            475      1980420     0.
   IS L to G Mapping   185            185      1705368     0.
   Star Forest Graph   378            378       299632     0.
     Discrete System   189            189       169352     0.
                SNES     3              3         4092     0.
      SNESLineSearch     3              3         2976     0.
              DMSNES     3              3         1992     0.
       Krylov Solver    24             24       240504     0.
     DMKSP interface     3              3         1944     0.
      Preconditioner    24             24        23136     0.

--- Event Stage 1: MG Apply

              Vector    15              0            0     0.
              Matrix     3              0            0     0.
           Index Set    15              6        10104     0.
========================================================================================================================
Average time to get PetscTime(): 0.
Average time for MPI_Barrier(): 7.13825e-05
Average time for zero size MPI_Send(): 1.97291e-05
#PETSc Option Table entries:
-A11_operator_type 0,1
-checkpoint_disable
-dau_nlevels 2
-dt_max 0.0023890909091
-fieldsplit_p_ksp_type preonly
-fieldsplit_p_pc_type jacobi
-fieldsplit_u_ksp_max_it 10
-fieldsplit_u_ksp_monitor
-fieldsplit_u_ksp_rtol 1.0e-2
-fieldsplit_u_ksp_type fgmres
-fieldsplit_u_mg_coarse_ksp_max_it 4
-fieldsplit_u_mg_coarse_ksp_pc_side right
-fieldsplit_u_mg_coarse_ksp_type gmres
-fieldsplit_u_mg_coarse_pc_type bjacobi
-fieldsplit_u_mg_levels_ksp_max_it 2
-fieldsplit_u_mg_levels_ksp_pc_side right
-fieldsplit_u_mg_levels_ksp_type gmres
-fieldsplit_u_mg_levels_pc_type jacobi
-fieldsplit_u_pc_mg_levels 2
-fieldsplit_u_pc_mg_log
-fieldsplit_u_pc_type mg
-ksp_atol 1.0e-8
-ksp_converged_reason
-ksp_max_it 2000
-ksp_rtol 1.0e-3
-ksp_type fgmres
-lattice_layout_Nx 3
-lattice_layout_Ny 3
-lattice_layout_Nz 3
-lattice_layout_perturb 0.0
-log_view
-monitor_stages
-mp_popctrl_np_upper -1
-mx 8
-my 8
-mz 8
-nsteps 2
-options_left 1
-output_frequency 1
-output_path sd3d_test
-pc_fieldsplit_schur_fact_type upper
-pc_fieldsplit_type schur
-pc_type fieldsplit
-ptatin_model sd3d
-snes_atol 1.0e-8
-snes_converged_reason
-snes_linesearch_type basic
-snes_max_it 5
-snes_max_linear_solve_fail 100000
-snes_monitor
-snes_rtol 1.0e-5
-snes_type newtonls
-snes_view
-stk_velocity_dm_mat_type aij
#End of PETSc Option Table entries
Compiled without FORTRAN kernels
Compiled with full precision matrices (default)
sizeof(short) 2 sizeof(int) 4 sizeof(long) 8 sizeof(void*) 8 sizeof(PetscScalar) 8 sizeof(PetscInt) 4
Configure options: PETSC_DIR=/Users/patrick/petsc-maint PETSC_ARCH=arch-darwin-maint-double-opt --with-debugging=0 --with-precision=double --with-scalar-type=real --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich --COPTFLAGS="-g -O3 -march=native -mtune=native " --CXXOPTFLAGS="-g -O3 -march=native -mtune=native " --FOPTFLAGS="-g -O3 -march=native -mtune=native " --download-c2html --download-suitesparse --download-c2html --download-sowing
-----------------------------------------
Libraries compiled on Tue Mar 13 19:32:24 2018 on Windbauer.local 
Machine characteristics: Darwin-17.3.0-x86_64-i386-64bit
Using PETSc directory: /Users/patrick/petsc-maint
Using PETSc arch: arch-darwin-maint-double-opt
-----------------------------------------

Using C compiler: /Users/patrick/petsc-maint/arch-darwin-maint-double-opt/bin/mpicc    -Wall -Wwrite-strings -Wno-strict-aliasing -Wno-unknown-pragmas -fstack-protector -fvisibility=hidden  ${COPTFLAGS} ${CFLAGS}
Using Fortran compiler: /Users/patrick/petsc-maint/arch-darwin-maint-double-opt/bin/mpif90   -Wall -ffree-line-length-0 -Wno-unused-dummy-argument   ${FOPTFLAGS} ${FFLAGS} 
-----------------------------------------

Using include paths: -I/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/include -I/Users/patrick/petsc-maint/include -I/Users/patrick/petsc-maint/include -I/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/include -I/opt/X11/include
-----------------------------------------

Using C linker: /Users/patrick/petsc-maint/arch-darwin-maint-double-opt/bin/mpicc
Using Fortran linker: /Users/patrick/petsc-maint/arch-darwin-maint-double-opt/bin/mpif90
Using libraries: -Wl,-rpath,/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/lib -L/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/lib -lpetsc -Wl,-rpath,/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/lib -L/Users/patrick/petsc-maint/arch-darwin-maint-double-opt/lib -Wl,-rpath,/opt/X11/lib -L/opt/X11/lib -Wl,-rpath,/opt/local/lib/gcc6/gcc/x86_64-apple-darwin17/6.4.0 -L/opt/local/lib/gcc6/gcc/x86_64-apple-darwin17/6.4.0 -Wl,-rpath,/opt/local/lib/gcc6 -L/opt/local/lib/gcc6 -lumfpack -lklu -lcholmod -lbtf -lccolamd -lcolamd -lcamd -lamd -lsuitesparseconfig -llapack -lblas -lX11 -ldl -lmpifort -lgfortran -lquadmath -lm -lmpicxx -lstdc++ -lmpi -lpmpi -lSystem -lgcc_ext.10.5 -ldl
-----------------------------------------

#PETSc Option Table entries:
-A11_operator_type 0,1
-checkpoint_disable
-dau_nlevels 2
-dt_max 0.0023890909091
-fieldsplit_p_ksp_type preonly
-fieldsplit_p_pc_type jacobi
-fieldsplit_u_ksp_max_it 10
-fieldsplit_u_ksp_monitor
-fieldsplit_u_ksp_rtol 1.0e-2
-fieldsplit_u_ksp_type fgmres
-fieldsplit_u_mg_coarse_ksp_max_it 4
-fieldsplit_u_mg_coarse_ksp_pc_side right
-fieldsplit_u_mg_coarse_ksp_type gmres
-fieldsplit_u_mg_coarse_pc_type bjacobi
-fieldsplit_u_mg_levels_ksp_max_it 2
-fieldsplit_u_mg_levels_ksp_pc_side right
-fieldsplit_u_mg_levels_ksp_type gmres
-fieldsplit_u_mg_levels_pc_type jacobi
-fieldsplit_u_pc_mg_levels 2
-fieldsplit_u_pc_mg_log
-fieldsplit_u_pc_type mg
-ksp_atol 1.0e-8
-ksp_converged_reason
-ksp_max_it 2000
-ksp_rtol 1.0e-3
-ksp_type fgmres
-lattice_layout_Nx 3
-lattice_layout_Ny 3
-lattice_layout_Nz 3
-lattice_layout_perturb 0.0
-log_view
-monitor_stages
-mp_popctrl_np_upper -1
-mx 8
-my 8
-mz 8
-nsteps 2
-options_left 1
-output_frequency 1
-output_path sd3d_test
-pc_fieldsplit_schur_fact_type upper
-pc_fieldsplit_type schur
-pc_type fieldsplit
-ptatin_model sd3d
-snes_atol 1.0e-8
-snes_converged_reason
-snes_linesearch_type basic
-snes_max_it 5
-snes_max_linear_solve_fail 100000
-snes_monitor
-snes_rtol 1.0e-5
-snes_type newtonls
-snes_view
-stk_velocity_dm_mat_type aij
#End of PETSc Option Table entries
There are no unused options.
