** ====================================================================================== 
**
**             ___________                          _______
**     _______/          /_____ ________ __ _   ___/       \ ____
**    /      /___    ___/      /__   __/  /  \ /  /\__ _   /     \
**   /  //  /   /   /  /  //  /  /  / /  /    /  /___/_   /  //  /
**  /  ___ /   /   /  /  _   /  /  / /  /  /    //       /  //  /
** /__/       /___/  /__//__/  /__/ /__/__/ \__//_______/______/
**
** 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 = -DTATIN_HAVE_CUDA -DTATIN_HAVE_OPENCL -DTATIN_HAVE_AVX
**                                                                       
** ====================================================================================== 
RheologyConstantsInitialise: global viscosity cut-off, min= 1.000000e-100, max = 1.000000e+100  
[pTatin] Created output directory: pt3dout 
[pTatin] Created log file: pt3dout/ptatin.log-2017.07.14_09:36:51 
[pTatin] Writing output to existing directory: pt3dout 
[pTatin] Created options file: pt3dout/ptatin.options-2017.07.14_09:36:51 
[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) = 1.69e-01,2.14e-01 (MB) 
VolumeQuadratureCreate_GaussLegendreStokes:
	Using 3x3 pnt Gauss Legendre quadrature
DataBucketView(MPI): ("GaussLegendre StokesCoefficients")
  L                  = 55296 
  buffer (max)       = 12 
  allocated          = 55308 
  nfields registered = 1 
    [  0]: field name  ==>>              QPntVolCoefStokes : Mem. usage = 2.07e-01 (MB) : rank0
  Total mem. usage                                                      = 2.65e+00 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 0]")
  L                  = 1152 
  buffer (max)       = 12 
  allocated          = 1170 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 2.24e-04 (MB) : rank0
  Total mem. usage                                                      = 1.31e-01 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 1]")
  L                  = 1152 
  buffer (max)       = 12 
  allocated          = 1170 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 2.03e-02 (MB) : rank0
  Total mem. usage                                                      = 1.31e-01 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 2]")
  L                  = 1152 
  buffer (max)       = 12 
  allocated          = 1172 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 2.24e-04 (MB) : rank0
  Total mem. usage                                                      = 1.31e-01 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 3]")
  L                  = 1152 
  buffer (max)       = 12 
  allocated          = 1172 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 3.24e-02 (MB) : rank0
  Total mem. usage                                                      = 1.31e-01 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 4]")
  L                  = 2304 
  buffer (max)       = 12 
  allocated          = 2322 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 2.24e-04 (MB) : rank0
  Total mem. usage                                                      = 2.60e-01 (MB) : collective
DataBucketView(MPI): ("SurfaceGaussLegendre StokesCoefficients[face 5]")
  L                  = 2304 
  buffer (max)       = 12 
  allocated          = 2322 
  nfields registered = 1 
    [  0]: field name  ==>>             QPntSurfCoefStokes : Mem. usage = 4.04e-02 (MB) : rank0
  Total mem. usage                                                      = 2.60e-01 (MB) : collective
  MaterialPointsStokes: Using Q1 projection
[[Swarm initialization: 0.0002 (sec)]]
SwarmMPntStd_AssignUniquePointIdentifiers : max_pid = 0 
[[Swarm->coordinate assignment: 1280 points : 0.0004 (sec)]]
************************** Starting _DataExCompleteCommunicationMap ************************** 
max_nnz = 11 
Mat Object: 12 MPI processes
  type: mpiaij
  rows=12, cols=12
  total: nonzeros=100, allocated nonzeros=112
  total number of mallocs used during MatSetValues calls =0
    not using I-node (on process 0) routines
************************** Ending _DataExCompleteCommunicationMap [setup time: 2.8751e-03 (sec)] ************************** 
[[SwarmDMDA3dDataExchangerCreate: time = 2.9981e-03 (sec)]]
[[ModelApplyInitialMeshGeometry_ViscousSinker]]
RUNNING DEFORMED MESH EXAMPLE 
[[ViscousSinker_ApplyInitialMaterialGeometry_SingleInclusion]]
[[ModelApplyBoundaryCondition_ViscousSinker]]



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

+  Test [compare_mf_A11]: Mesh 16 x 16 x 8 
MatMultA11(MF):      time 1.2195e-02 (sec): ratio 9.9764e+01%: min/max 1.2177e-02 1.2206e-02 (sec)
MatAssemblyA11(ASM): time 6.8427e-02 (sec): ratio 9.9660e+01%: min/max 6.8276e-02 6.8509e-02 (sec)
MatMultA11(ASM):     time 3.4571e-03 (sec): ratio 9.2113e+01%: min/max 3.2609e-03 3.5400e-03 (sec)
  y.y    = +1.31651556e+05 [mfo]
  y2.y2  = +1.31651556e+05 [asm]
  min[A11_mfo.x-A11_asm.x]  = -2.26767390e-14 
  max[A11_mfo.x-A11_asm.x]  = +2.25052182e-14 

+  Test [compare_mf_A21]: Mesh 16 x 16 x 8 
StokesA21Preallocation_basic: using nnz = 81 and onnz = 30 
  y.y    = +9.43076063e-01 [mfo]
  y2.y2  = +9.43076063e-01 [asm]
  min[A21_mfo.x-A21_asm.x]  = -4.09828421e-17 
  max[A21_mfo.x-A21_asm.x]  = +7.00394603e-17 

+  Test [compare_mf_A12]: Mesh 16 x 16 x 8 
StokesA12Preallocation_basic: using nnz = 32 and onnz = 16 
  y.y    = +1.75148072e+00 [mfo]
  y2.y2  = +1.75148072e+00 [asm]
  min[A12_mfo.x-A12_asm.x]  = -1.76941795e-16 
  max[A12_mfo.x-A12_asm.x]  = +1.35308431e-16 

+  Test [compare_mf_A]: Mesh 16 x 16 x 8 
  y.y    = +4.98077749e+03 [mfo]
  y2.y2  = +4.98077749e+03 [asm]
  min[A_mfo.x-A_asm.x]  = -8.88178420e-16 
  max[A_mfo.x-A_asm.x]  = +1.33226763e-15 

+  Test [compare_mf_diagA11]: Mesh 16 x 16 x 8 
  y.y    = +3.82278542e+04 [mfo]
  y2.y2  = +3.82278542e+04 [asm]
  min[diagA11_mfo-diagA11_asm]  = -3.55271368e-15 
  max[diagA11_mfo-diagA11_asm]  = +4.44089210e-15 



====================================================================
[[ModelDestroy_Template]]
[[ModelDestroy_ViscousSinker]]
[[ModelDestroy_Gene3D]]
[[ModelDestroy_Gene3DNueve]]
[[ModelDestroy_Indentor]]
[[ModelDestroy_Rift3D]]
[[ModelDestroy_Rift3D_T]]
[[ModelDestroy_Sierra]]
[[ModelDestroy_AdvDiffExample]]
[[ModelDestroy_Delamination]]
[[ModelDestroy_Riftrh]]
[[ModelDestroy_MultilayerFolding]]
[[ModelDestroy_SubmarineLavaFlow]]
[[ModelDestroy_ExSubduction]]
[[ModelDestroy_iPLUS]]
[[ModelDestroy_Subduction_Initiation2d]]
[[ModelDestroy_Thermal_Convection2d]]
#PETSc Option Table entries:
-a11_op subrepart
-compare_operators
-lattice_layout_perturb 0.0
-model_viscous_sinker_eta1 10
-mx 16
-my 16
-mz 8
-options_left
-output_path pt3dout
-ptatin_model viscous_sinker
#End of PETSc Option Table entries
There are no unused options.
