** ====================================================================================== 
**
**             ___________                          _______
**     _______/          /_____ ________ __ _   ___/       \ ____
**    /      /___    ___/      /__   __/  /  \ /  /\__ _   /     \
**   /  //  /   /   /  /  //  /  /  / /  /    /  /___/_   /  //  /
**  /  ___ /   /   /  /  _   /  /  / /  /  /    //       /  //  /
** /__/       /___/  /__//__/  /__/ /__/__/ \__//_______/______/
**
** 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                  = 200 
  buffer (max)       = 0 
  allocated          = 200 
  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                                                      = 8.00e-02 (MB) : <collective over 1 ranks>
[pTatin] Created output directory: pt3dout 
[pTatin] Created log file: pt3dout/ptatin.log-2018.03.14_19:33:17 
[pTatin] Created options file: pt3dout/ptatin.options-2018.03.14_19:33:17 
[pTatin] Created options file: pt3dout/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 "viscous_sinker" was detected
[[ModelInitialize_ViscousSinker]]
  MaterialPointsStokes: Using Q1 projection
[[Swarm initialization: 0.0003 (sec)]]
[[Swarm->coordinate assignment: 4096 points : 0.0014 (sec)]]
[[SwarmDMDA3dDataExchangerCreate: time = 1.6403e-04 (sec)]]
[[ModelApplyInitialMeshGeometry_ViscousSinker]]
RUNNING DEFORMED MESH EXAMPLE 
[[ViscousSinker_ApplyInitialMaterialGeometry_SingleInclusion]]
[[ModelApplyBoundaryCondition_ViscousSinker]]
Mesh size (8 x 8 x 8) : MG levels 3  
         level [ 0]: global Q2 elements (2 x 2 x 2) 
         level [ 1]: global Q2 elements (4 x 4 x 4) 
         level [ 2]: global Q2 elements (8 x 8 x 8) 
[r   0]: level [ 0]: local Q2 elements  (2 x 2 x 2) 
[r   0]: level [ 1]: local Q2 elements  (4 x 4 x 4) 
[r   0]: level [ 2]: local Q2 elements  (8 x 8 x 8) 
[r   0]: level [ 0]: element range [0 - 1] x [0 - 1] x [0 - 1] 
[r   0]: level [ 1]: element range [0 - 3] x [0 - 3] x [0 - 3] 
[r   0]: level [ 2]: element range [0 - 7] x [0 - 7] x [0 - 7] 
[[ModelApplyBoundaryConditionMG_ViscousSinker]]
[[ModelInitialCondition_ViscousSinker]]
*** Rheology update for RHEOLOGY_VISCOUS selected ***
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; cpu time 6.41e-05 (sec)
Level [2]: Coarse grid type :: Re-discretisation :: assembled operator 
Level [1]: Coarse grid type :: Galerkin :: assembled operator 
Level [0]: Coarse grid type :: Galerkin :: assembled operator 
[[ModelOutput_ViscousSinker]]
  writing pvdfilename pt3dout/timeseries_vp.pvd 
  writing pvdfilename pt3dout/timeseries_mpoints_std.pvd 
   [[ COMPUTING FLOW FIELD FOR STEP : 0 ]]
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; cpu time 5.10e-05 (sec)
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; cpu time 7.10e-05 (sec)
    0 KSP Residual norm 0.00328249 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1. 
      1 KSP Residual norm 0.361905 
      2 KSP Residual norm 0.00648797 
    1 KSP Residual norm 0.00328244 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.00913657 
      1 KSP Residual norm 0.00149661 
      2 KSP Residual norm 0.000286117 
      3 KSP Residual norm 7.51999e-06 
    2 KSP Residual norm 0.00328096 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.108788 
      1 KSP Residual norm 0.000844228 
    3 KSP Residual norm 0.00327805 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.379174 
      1 KSP Residual norm 0.00282463 
    4 KSP Residual norm 0.00327076 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.984649 
      1 KSP Residual norm 0.00958806 
    5 KSP Residual norm 0.00322628 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.28477 
      1 KSP Residual norm 0.0127001 
    6 KSP Residual norm 0.00314596 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 1.02788 
      1 KSP Residual norm 0.00728591 
    7 KSP Residual norm 0.00268119 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.871612 
      1 KSP Residual norm 0.00831313 
    8 KSP Residual norm 0.00180557 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.872015 
      1 KSP Residual norm 0.0140172 
      2 KSP Residual norm 0.00516406 
    9 KSP Residual norm 0.00107576 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.511082 
      1 KSP Residual norm 0.00533368 
      2 KSP Residual norm 0.00271548 
   10 KSP Residual norm 0.000532551 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.511023 
      1 KSP Residual norm 0.00865199 
      2 KSP Residual norm 0.00595491 
      3 KSP Residual norm 0.000105401 
   11 KSP Residual norm 0.00027117 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.501737 
      1 KSP Residual norm 0.00701739 
      2 KSP Residual norm 0.00398537 
   12 KSP Residual norm 0.00012517 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.487643 
      1 KSP Residual norm 0.00671948 
      2 KSP Residual norm 0.00464823 
   13 KSP Residual norm 4.58432e-05 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.365518 
      1 KSP Residual norm 0.00496607 
      2 KSP Residual norm 0.00356216 
   14 KSP Residual norm 1.93963e-05 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.444427 
      1 KSP Residual norm 0.00273843 
   15 KSP Residual norm 9.87314e-06 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.849056 
      1 KSP Residual norm 0.00961961 
      2 KSP Residual norm 0.00479929 
   16 KSP Residual norm 3.78996e-06 
      Residual norms for fieldsplit_u_ solve.
      0 KSP Residual norm 0.553595 
      1 KSP Residual norm 0.00208465 
   17 KSP Residual norm 2.6426e-06 
  Linear solve converged due to CONVERGED_RTOL iterations 17
SNES Object: 1 MPI processes
  type: ksponly
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=17
  total number of function evaluations=1
  norm schedule ALWAYS
  SNESLineSearch Object: 1 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=1
  KSP Object: 1 MPI processes
    type: fgmres
      restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
      happy breakdown tolerance 1e-30
    maximum iterations=10000, nonzero initial guess
    tolerances:  relative=0.001, absolute=1e-50, divergence=10000.
    right preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 1 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_) 1 MPI processes
          type: fgmres
            restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
            happy breakdown tolerance 1e-30
          maximum iterations=10000, initial guess is zero
          tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
          right preconditioning
          using UNPRECONDITIONED norm type for convergence test
        PC Object: (fieldsplit_u_) 1 MPI processes
          type: mg
            type is MULTIPLICATIVE, levels=3 cycles=v
              Cycles per PCApply=1
              Not using Galerkin computed coarse grid matrices
          Coarse grid solver -- level -------------------------------
            KSP Object: (fieldsplit_u_mg_coarse_) 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_) 1 MPI processes
              type: lu
                out-of-place factorization
                tolerance for zero pivot 2.22045e-14
                using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
                matrix ordering: nd
                factor fill ratio given 5., needed 1.50824
                  Factored matrix follows:
                    Mat Object: 1 MPI processes
                      type: seqaij
                      rows=375, cols=375, bs=3
                      package used to perform factorization: petsc
                      total: nonzeros=93105, allocated nonzeros=93105
                      total number of mallocs used during MatSetValues calls =0
                        using I-node routines: found 118 nodes, limit used is 5
              linear system matrix = precond matrix:
              Mat Object: 1 MPI processes
                type: seqaij
                rows=375, cols=375, bs=3
                total: nonzeros=61731, allocated nonzeros=0
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
                  using I-node routines: found 125 nodes, limit used is 5
          Down solver (pre-smoother) on level 1 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_1_) 1 MPI processes
              type: chebyshev
                eigenvalue estimates used:  min = 0.525085, max = 2.88796
                eigenvalues estimate via gmres min 0.164191, max 2.62542
                eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
                KSP Object: (fieldsplit_u_mg_levels_1_esteig_) 1 MPI processes
                  type: gmres
                    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=1e-12, absolute=1e-50, divergence=10000.
                  left preconditioning
                  using NONE norm type for convergence test
                estimating eigenvalues using noisy right hand side
              maximum iterations=10, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              left preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_1_) 1 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: 1 MPI processes
                type: seqaij
                rows=2187, cols=2187, bs=3
                total: nonzeros=533871, allocated nonzeros=0
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
                  using I-node routines: found 729 nodes, limit used is 5
          Up solver (post-smoother) same as down solver (pre-smoother)
          Down solver (pre-smoother) on level 2 -------------------------------
            KSP Object: (fieldsplit_u_mg_levels_2_) 1 MPI processes
              type: chebyshev
                eigenvalue estimates used:  min = 0.521622, max = 2.86892
                eigenvalues estimate via gmres min 0.0738372, max 2.60811
                eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
                KSP Object: (fieldsplit_u_mg_levels_2_esteig_) 1 MPI processes
                  type: gmres
                    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=1e-12, absolute=1e-50, divergence=10000.
                  left preconditioning
                  using NONE norm type for convergence test
                estimating eigenvalues using noisy right hand side
              maximum iterations=10, nonzero initial guess
              tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
              left preconditioning
              using NONE norm type for convergence test
            PC Object: (fieldsplit_u_mg_levels_2_) 1 MPI processes
              type: jacobi
              linear system matrix = precond matrix:
              Mat Object: (fieldsplit_u_) 1 MPI processes
                type: seqaij
                rows=14739, cols=14739, bs=3
                total: nonzeros=4437351, allocated nonzeros=4437351
                total number of mallocs used during MatSetValues calls =0
                  has attached near null space
                  using I-node routines: found 4913 nodes, limit used is 5
          Up solver (post-smoother) same as down solver (pre-smoother)
          linear system matrix = precond matrix:
          Mat Object: (fieldsplit_u_) 1 MPI processes
            type: seqaij
            rows=14739, cols=14739, bs=3
            total: nonzeros=4437351, allocated nonzeros=4437351
            total number of mallocs used during MatSetValues calls =0
              has attached near null space
              using I-node routines: found 4913 nodes, limit used is 5
      KSP solver for S = A11 - A10 inv(A00) A01 
        KSP Object: (fieldsplit_p_) 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_p_) 1 MPI processes
          type: jacobi
          linear system matrix followed by preconditioner matrix:
          Mat Object: (fieldsplit_p_) 1 MPI processes
            type: schurcomplement
            rows=2048, cols=2048
              Schur complement A11 - A10 inv(A00) A01
              A11
                Mat Object: (fieldsplit_p_) 1 MPI processes
                  type: seqsbaij
                  rows=2048, cols=2048, bs=4
                  total: nonzeros=8192, allocated nonzeros=8192
                  total number of mallocs used during MatSetValues calls =0
                      block size is 4
              A10
                Mat Object: (Bpu_) 1 MPI processes
                  type: shell
                  rows=2048, cols=14739
              KSP of A00
                KSP Object: (fieldsplit_u_) 1 MPI processes
                  type: fgmres
                    restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
                    happy breakdown tolerance 1e-30
                  maximum iterations=10000, initial guess is zero
                  tolerances:  relative=0.01, absolute=1e-50, divergence=10000.
                  right preconditioning
                  using UNPRECONDITIONED norm type for convergence test
                PC Object: (fieldsplit_u_) 1 MPI processes
                  type: mg
                    type is MULTIPLICATIVE, levels=3 cycles=v
                      Cycles per PCApply=1
                      Not using Galerkin computed coarse grid matrices
                  Coarse grid solver -- level -------------------------------
                    KSP Object: (fieldsplit_u_mg_coarse_) 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_) 1 MPI processes
                      type: lu
                        out-of-place factorization
                        tolerance for zero pivot 2.22045e-14
                        using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
                        matrix ordering: nd
                        factor fill ratio given 5., needed 1.50824
                          Factored matrix follows:
                            Mat Object: 1 MPI processes
                              type: seqaij
                              rows=375, cols=375, bs=3
                              package used to perform factorization: petsc
                              total: nonzeros=93105, allocated nonzeros=93105
                              total number of mallocs used during MatSetValues calls =0
                                using I-node routines: found 118 nodes, limit used is 5
                      linear system matrix = precond matrix:
                      Mat Object: 1 MPI processes
                        type: seqaij
                        rows=375, cols=375, bs=3
                        total: nonzeros=61731, allocated nonzeros=0
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                          using I-node routines: found 125 nodes, limit used is 5
                  Down solver (pre-smoother) on level 1 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_1_) 1 MPI processes
                      type: chebyshev
                        eigenvalue estimates used:  min = 0.525085, max = 2.88796
                        eigenvalues estimate via gmres min 0.164191, max 2.62542
                        eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
                        KSP Object: (fieldsplit_u_mg_levels_1_esteig_) 1 MPI processes
                          type: gmres
                            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=1e-12, absolute=1e-50, divergence=10000.
                          left preconditioning
                          using NONE norm type for convergence test
                        estimating eigenvalues using noisy right hand side
                      maximum iterations=10, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      left preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_1_) 1 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: 1 MPI processes
                        type: seqaij
                        rows=2187, cols=2187, bs=3
                        total: nonzeros=533871, allocated nonzeros=0
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                          using I-node routines: found 729 nodes, limit used is 5
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  Down solver (pre-smoother) on level 2 -------------------------------
                    KSP Object: (fieldsplit_u_mg_levels_2_) 1 MPI processes
                      type: chebyshev
                        eigenvalue estimates used:  min = 0.521622, max = 2.86892
                        eigenvalues estimate via gmres min 0.0738372, max 2.60811
                        eigenvalues estimated using gmres with translations  [0. 0.2; 0. 1.1]
                        KSP Object: (fieldsplit_u_mg_levels_2_esteig_) 1 MPI processes
                          type: gmres
                            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=1e-12, absolute=1e-50, divergence=10000.
                          left preconditioning
                          using NONE norm type for convergence test
                        estimating eigenvalues using noisy right hand side
                      maximum iterations=10, nonzero initial guess
                      tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
                      left preconditioning
                      using NONE norm type for convergence test
                    PC Object: (fieldsplit_u_mg_levels_2_) 1 MPI processes
                      type: jacobi
                      linear system matrix = precond matrix:
                      Mat Object: (fieldsplit_u_) 1 MPI processes
                        type: seqaij
                        rows=14739, cols=14739, bs=3
                        total: nonzeros=4437351, allocated nonzeros=4437351
                        total number of mallocs used during MatSetValues calls =0
                          has attached near null space
                          using I-node routines: found 4913 nodes, limit used is 5
                  Up solver (post-smoother) same as down solver (pre-smoother)
                  linear system matrix = precond matrix:
                  Mat Object: (fieldsplit_u_) 1 MPI processes
                    type: seqaij
                    rows=14739, cols=14739, bs=3
                    total: nonzeros=4437351, allocated nonzeros=4437351
                    total number of mallocs used during MatSetValues calls =0
                      has attached near null space
                      using I-node routines: found 4913 nodes, limit used is 5
              A01
                Mat Object: (Bup_) 1 MPI processes
                  type: shell
                  rows=14739, cols=2048
          Mat Object: (fieldsplit_p_) 1 MPI processes
            type: seqsbaij
            rows=2048, cols=2048, bs=4
            total: nonzeros=8192, allocated nonzeros=8192
            total number of mallocs used during MatSetValues calls =0
                block size is 4
    linear system matrix followed by preconditioner matrix:
    Mat Object: (stokes_Amf_) 1 MPI processes
      type: shell
      rows=16787, cols=16787
    Mat Object: 1 MPI processes
      type: nest
      rows=16787, cols=16787
        Matrix object: 
          type=nest, rows=2, cols=2 
          MatNest structure: 
          (0,0) : prefix="fieldsplit_u_", type=seqaij, 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=seqsbaij, rows=2048, cols=2048 
Update rheology (viscous) [mpoint]: (min,max)_eta 1.00e-03,1.00e+00; log10(max/min) 3.00e+00; cpu time 5.82e-05 (sec)
[[ModelOutput_ViscousSinker]]
  TimeStep control(StkCourant): | current = 1.0000e+30 : trial = 8.0328e-02 [accepted] | ==>> dt used = 8.0328e-02 |
  TimeStep control(StkSurfaceCourant): | current = 8.0328e-02 : trial = 1.0000e+32 | ==>> dt used = 8.0328e-02 |
  timestep[] dt_courant = 8.0328e-02 
[[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]]
