** ====================================================================================== 
**
**             ___________                          _______
**     _______/          /_____ ________ __ _   ___/       \ ____
**    /      /___    ___/      /__   __/  /  \ /  /\__ _   /     \
**   /  //  /   /   /  /  //  /  /  / /  /    /  /___/_   /  //  /
**  /  ___ /   /   /  /  _   /  /  / /  /  /    //       /  //  /
** /__/       /___/  /__//__/  /__/ /__/__/ \__//_______/______/
**
** 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 -g -Wall -Wno-unused-variable
**                                                                       
** 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 
**                                                                       
** ====================================================================================== 
RheologyConstantsInitialise: global viscosity cut-off, min= 1.000000e-100, max = 1.000000e+100  
[pTatin] Writing output to existing directory: pt3dout 
[pTatin] Created log file: pt3dout/ptatin.log-2015.07.15_15:35:06 
[pTatin] Writing output to existing directory: pt3dout 
[pTatin] Created options file: pt3dout/ptatin.options-2015.07.15_15:35:06 
[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 "Gene3DNueve"
  [pTatinModel]: Registering model [4] with name "indentor"
  [pTatinModel]: Registering model [5] with name "rift3D"
  [pTatinModel]: Registering model [6] with name "rift3D_T"
  [pTatinModel]: Registering model [7] with name "sierra"
  [pTatinModel]: Registering model [8] with name "advdiff_example"
  [pTatinModel]: Registering model [9] with name "delamination"
  [pTatinModel]: Registering model [10] with name "Riftrh"
  [pTatinModel]: Registering model [11] with name "geomod2008"
  [pTatinModel]: Registering model [12] with name "multilayer_folding"
  [pTatinModel]: Registering model [13] with name "submarinelavaflow"
  [pTatinModel]: Registering model [14] with name "ex_subduction"
  [pTatinModel]: Registering model [15] with name "iplus"
  [pTatinModel]: Registering model [16] with name "subduction_initiation2d"
  [pTatinModel]: Registering model [17] with name "convection2d"
  [pTatinModel]: Registering model [18] with name "thermal_sb"
  [pTatinModel]: Registering model [19] with name "sd3d"
  [pTatinModel]: Registering model [20] with name "pas"
  [pTatinModel]: Registering model [21] with name "pd"
  [pTatinModel]: -ptatin_model "viscous_sinker" was detected
[[ModelInitialize_ViscousSinker]]
BCList: Mem. usage (min,max) = 4.72e-01,4.72e-01 (MB) 
VolumeQuadratureCreate_GaussLegendreStokes:
	Using 3x3 pnt Gauss Legendre quadrature
DataBucketView(SEQ): ("GaussLegendre StokesCoefficients")
  L                  = 13824 
  buffer             = 1 
  allocated          = 13825 
  nfields registered = 1 
    [  0]: field name  ==>>              QPntVolCoefStokes : Mem. usage = 6.64e-01 (MB) 
  Total mem. usage                                                      = 6.64e-01 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 0]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 1]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 2]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 3]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 4]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
DataBucketView(SEQ): ("SurfaceGaussLegendre StokesCoefficients[face 5]")
  L                  = 576 
  buffer             = 1 
  allocated          = 577 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 6.46e-02 (MB) 
  Total mem. usage                                                      = 6.46e-02 (MB) 
  MaterialPointsStokes: Using Q1 projection
[[Swarm initialization: 0.0002 (sec)]]
SwarmMPntStd_AssignUniquePointIdentifiers : max_pid = 0 
[[Swarm->coordinate assignment: 4096 points : 0.0014 (sec)]]
************************** Starting _DataExCompleteCommunicationMap ************************** 
max_nnz = 0 
Mat Object: 1 MPI processes
  type: seqaij
  rows=1, cols=1
  total: nonzeros=0, allocated nonzeros=0
  total number of mallocs used during MatSetValues calls =0
    not using I-node routines
************************** Ending _DataExCompleteCommunicationMap [setup time: 6.9757e-04 (sec)] ************************** 
[[SwarmDMDA3dDataExchangerCreate: time = 7.2878e-04 (sec)]]
[[ModelApplyInitialMeshGeometry_ViscousSinker]]
RUNNING DEFORMED MESH EXAMPLE 
[[ViscousSinker_ApplyInitialMaterialGeometry_SingleInclusion]]
[[ModelApplyBoundaryCondition_ViscousSinker]]



====================================================================

+  Test [perform_viscous_solve]: Mesh 8 x 8 x 8 
MatAssemblyA11(ASM):                time 1.3019e+00 (sec): ratio 1.0000e+02%: min/max 1.3019e+00 1.3019e+00 (sec)
KSPSetUpA11:                        time 1.0144e+00 (sec): ratio 1.0000e+02%: min/max 1.0144e+00 1.0144e+00 (sec)
  0 KSP Residual norm 61086.4 
  1 KSP Residual norm 26276.9 
  2 KSP Residual norm 20345.9 
  3 KSP Residual norm 11666.8 
  4 KSP Residual norm 5094.93 
  5 KSP Residual norm 2397.25 
  6 KSP Residual norm 917.31 
  7 KSP Residual norm 286.629 
  8 KSP Residual norm 95.5837 
  9 KSP Residual norm 48.189 
 10 KSP Residual norm 27.1225 
 11 KSP Residual norm 15.1705 
 12 KSP Residual norm 5.51475 
 13 KSP Residual norm 1.18395 
 14 KSP Residual norm 0.383348 
 15 KSP Residual norm 0.176362 
 16 KSP Residual norm 0.0846772 
 17 KSP Residual norm 0.032297 
 18 KSP Residual norm 0.0151644 
 19 KSP Residual norm 0.00834386 
 20 KSP Residual norm 0.0045717 
 21 KSP Residual norm 0.00177446 
 22 KSP Residual norm 0.000696058 
 23 KSP Residual norm 0.000243403 
 24 KSP Residual norm 8.84972e-05 
 25 KSP Residual norm 2.39719e-05 
 26 KSP Residual norm 5.18094e-06 
 27 KSP Residual norm 1.44638e-06 
 28 KSP Residual norm 3.7809e-07 
 29 KSP Residual norm 1.14989e-07 
 30 KSP Residual norm 2.59811e-08 
KSPSolveA11(its = 30,cycles = 5)   time 3.6433e+00 (sec): ratio 1.0000e+02%: min/max 3.6433e+00 3.6433e+00 (sec)
KSPSolveA11: average                time 7.2865e-01 (sec): ratio 1.0000e+02%: min/max 7.2865e-01 7.2865e-01 (sec)
KSP Object: 1 MPI processes
  type: gmres
    GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
    GMRES: happy breakdown tolerance 1e-30
  maximum iterations=30, initial guess is zero
  tolerances:  relative=1e-20, absolute=1e-50, divergence=10000
  left preconditioning
  using PRECONDITIONED norm type for convergence test
PC Object: 1 MPI processes
  type: ilu
    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=14739, cols=14739, bs=3
          package used to perform factorization: petsc
          total: nonzeros=4.43735e+06, allocated nonzeros=4.43735e+06
          total number of mallocs used during MatSetValues calls =0
            using I-node routines: found 4913 nodes, limit used is 5
  linear system matrix = precond matrix:
  Mat Object:   1 MPI processes
    type: seqaij
    rows=14739, cols=14739, bs=3
    total: nonzeros=4.43735e+06, allocated nonzeros=4.43735e+06
    total number of mallocs used during MatSetValues calls =0
      using I-node routines: found 4913 nodes, limit used is 5



====================================================================
