update the test results

Author: YiminJin

source/utilities.cc

@@ -3306,6 +3306,7 @@ namespace aspect
       }
 
 
+
       template <>
       const Tensor<3,3> &levi_civita<3>()
       {
@@ -3327,6 +3328,7 @@ namespace aspect
       }
 
 
+
       template <int dim>
       SymmetricTensor<2,dim> consistent_deviator(const SymmetricTensor<2,dim> &input)
       {
@@ -3340,6 +3342,7 @@ namespace aspect
       }
 
 
+
       template <int dim>
       double consistent_second_invariant_of_deviatoric_tensor(const SymmetricTensor<2,dim> &t)
       {

tests/Picard_compressible/screen-output

@@ -9,13 +9,13 @@ Number of degrees of freedom: 53,952 (25,344+3,264+12,672+12,672)
       Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.0979426
 
    Solving Stokes system (GMG)... 70+0 iterations.
-      Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.000868458
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.000868383
 
    Solving Stokes system (GMG)... 65+0 iterations.
-      Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 5.63416e-05
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 5.63357e-05
 
    Solving Stokes system (GMG)... 51+0 iterations.
-      Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 7.73424e-06
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 7.73425e-06
 
 
    Postprocessing:

tests/boundary_traction_function_cartesian_free_surface/screen-output

@@ -108,9 +108,9 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_xx_old system ... 1 iterations.
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
-   Solving Stokes system (GMG)... 108+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.46859e-16, 8.48765e-08, 8.48765e-08, 1.46642e-07, 0, 0, 0, 1.3304e-05
-      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.3304e-05
+   Solving Stokes system (GMG)... 109+0 iterations.
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.46859e-16, 8.73173e-08, 8.73173e-08, 1.43995e-07, 0, 0, 0, 1.35446e-05
+      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.35446e-05
 
    Solving temperature system... 0 iterations.
    Solving ve_stress_xx system ... 1 iterations.
@@ -119,13 +119,13 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_xx_old system ... 1 iterations.
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
-   Solving Stokes system (GMG)... 24+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.59591e-16, 5.86833e-06, 5.86833e-06, 3.49505e-06, 0, 0, 0, 4.27217e-11
-      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.86833e-06
+   Solving Stokes system (GMG)... 25+0 iterations.
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.55204e-16, 5.84216e-06, 5.84216e-06, 3.43134e-06, 0, 0, 0, 4.32824e-11
+      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.84216e-06
 
 
    Postprocessing:
-     Compositions min/max/mass: -3.512e+04/3.512e+04/-1.351e+08 // -3.512e+04/3.512e+04/1.351e+08 // -2.271e+04/3.336e+04/3.431e+12 // -3.512e+04/3.512e+04/-1.351e+08 // -3.512e+04/3.512e+04/1.351e+08 // -2.271e+04/3.336e+04/3.431e+12
+     Compositions min/max/mass: -3.612e+04/4.354e+04/-1.335e+08 // -3.612e+04/4.354e+04/1.367e+08 // -2.271e+04/3.336e+04/3.431e+12 // -3.612e+04/4.354e+04/-1.335e+08 // -3.612e+04/4.354e+04/1.367e+08 // -2.271e+04/3.336e+04/3.431e+12
      Temperature min/avg/max:   293 K, 293 K, 293 K
      Topography min/max:        -7.451e-09 m, 9.313e-10 m
      RMS, max velocity:         0.0109 m/year, 0.0377 m/year
@@ -144,12 +144,12 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
    Solving Stokes system (GMG)... 113+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.30109e-16, 4.08776e-07, 4.08776e-07, 3.87956e-07, 0, 0, 0, 7.76729e-06
-      Relative nonlinear residual (total system) after nonlinear iteration 1: 7.76729e-06
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.33736e-16, 4.05723e-07, 4.05723e-07, 3.7722e-07, 0, 0, 0, 8.08466e-06
+      Relative nonlinear residual (total system) after nonlinear iteration 1: 8.08466e-06
 
 
    Postprocessing:
-     Compositions min/max/mass: -1.06e+05/1.06e+05/-7.371e+07 // -1.06e+05/1.06e+05/7.371e+07 // -8.863e+04/5.522e+04/-4.078e+14 // -1.06e+05/1.06e+05/-7.371e+07 // -1.06e+05/1.06e+05/7.371e+07 // -8.863e+04/5.522e+04/-4.078e+14
+     Compositions min/max/mass: -1.066e+05/1.066e+05/-1.526e+08 // -1.066e+05/1.066e+05/-1.45e+07 // -8.863e+04/5.229e+04/-4.078e+14 // -1.066e+05/1.066e+05/-1.526e+08 // -1.066e+05/1.066e+05/-1.45e+07 // -8.863e+04/5.229e+04/-4.078e+14
      Temperature min/avg/max:   293 K, 293 K, 293 K
      Topography min/max:        -0.09426 m, 0.02027 m
      RMS, max velocity:         0.011 m/year, 0.0383 m/year

tests/boundary_traction_function_cartesian_free_surface/statistics

@@ -42,6 +42,6 @@
 # 42: Maximum topography (m)
 # 43: RMS velocity (m/year)
 # 44: Max. velocity (m/year)
-0 0.000000000000e+00 0.000000000000e+00 4380 41098 4644 236520 2 0 0 0 0 0 0 0 216 219 219  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00 0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00 0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.40488011e-15 -9.31322575e-10 9.31322575e-10 3.44414753e-05 5.26615309e-04 
-1 2.500000000000e+00 2.500000000000e+00 4380 41098 4644 236520 2 0 2 2 2 2 2 2 132 135 135 -3.51213277e+04 3.51213271e+04 -1.35090965e+08 -3.51213271e+04 3.51213277e+04 1.35090965e+08 -2.27140438e+04 3.33570536e+04  3.43054006e+12 -3.51213281e+04 3.51213275e+04 -1.35090998e+08 -3.51213275e+04 3.51213281e+04 1.35090998e+08 -2.27140442e+04 3.33570533e+04  3.43053620e+12 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.61580189e-15 -7.45058060e-09 9.31322575e-10 1.08910272e-02 3.76945054e-02 
-2 5.000000000000e+00 2.500000000000e+00 4380 41098 4644 236520 1 0 1 1 1 1 1 1 113 115 115 -1.06045937e+05 1.06045937e+05 -7.37124683e+07 -1.06045937e+05 1.06045937e+05 7.37124683e+07 -8.86308305e+04 5.52207115e+04 -4.07785493e+14 -1.06045914e+05 1.06045914e+05 -7.37125503e+07 -1.06045914e+05 1.06045914e+05 7.37125503e+07 -8.86308305e+04 5.52207121e+04 -4.07787409e+14 2.93000000e+02 2.93000000e+02 2.93000000e+02 -1.35592571e-16 -9.42646312e-02 2.02718014e-02 1.10291575e-02 3.82968158e-02 
+0 0.000000000000e+00 0.000000000000e+00 4380 41098 4644 236520 2 0 0 0 0 0 0 0 216 219 219  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.40488011e-15 -9.31322575e-10 9.31322575e-10 3.44414753e-05 5.26615309e-04 
+1 2.500000000000e+00 2.500000000000e+00 4380 41098 4644 236520 2 0 2 2 2 2 2 2 134 137 137 -3.61236472e+04 4.35399634e+04 -1.33521411e+08 -3.61236462e+04 4.35399638e+04  1.36660518e+08 -2.27140439e+04 3.33570525e+04  3.43053703e+12 -3.61236472e+04 4.35399636e+04 -1.33521445e+08 -3.61236462e+04 4.35399641e+04  1.36660551e+08 -2.27140442e+04 3.33570523e+04  3.43053412e+12 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.61580189e-15 -7.45058060e-09 9.31322575e-10 1.08910239e-02 3.76945053e-02 
+2 5.000000000000e+00 2.500000000000e+00 4380 41098 4644 236520 1 0 1 1 1 1 1 1 113 115 115 -1.06640941e+05 1.06573682e+05 -1.52557581e+08 -1.06640941e+05 1.06573682e+05 -1.45034263e+07 -8.86307931e+04 5.22853238e+04 -4.07787457e+14 -1.06640918e+05 1.06573644e+05 -1.52557651e+08 -1.06640918e+05 1.06573644e+05 -1.45033571e+07 -8.86307931e+04 5.22853244e+04 -4.07789373e+14 2.93000000e+02 2.93000000e+02 2.93000000e+02 -1.50658412e-16 -9.42646312e-02 2.02718014e-02 1.10291555e-02 3.82968120e-02 

tests/boundary_traction_function_cartesian_free_surface/topography.00002

@@ -302,8 +302,8 @@
 4.30419e+06 4.69719e+06 -0.00387885
 4.29735e+06 4.70345e+06 -0.00254236
 4.29735e+06 4.70345e+06 -0.00254236
-4.2905e+06 4.70969e+06 -0.00127911
-4.2905e+06 4.70969e+06 -0.00127911
+4.2905e+06 4.70969e+06 -0.0012791
+4.2905e+06 4.70969e+06 -0.0012791
 4.28365e+06 4.71593e+06 -8.42009e-05
 4.28365e+06 4.71593e+06 -8.42009e-05
 4.27679e+06 4.72216e+06 0.00104675
@@ -320,8 +320,8 @@
 4.24233e+06 4.75313e+06 0.00586973
 4.23541e+06 4.7593e+06 0.00668936
 4.23541e+06 4.7593e+06 0.00668936
-4.22849e+06 4.76545e+06 0.0074667
-4.22849e+06 4.76545e+06 0.0074667
+4.22849e+06 4.76545e+06 0.00746671
+4.22849e+06 4.76545e+06 0.00746671
 4.22155e+06 4.7716e+06 0.008204
 4.22155e+06 4.7716e+06 0.008204
 4.21461e+06 4.77773e+06 0.0089033

tests/boundary_traction_function_spherical_free_surface/screen-output

@@ -108,9 +108,9 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_xx_old system ... 1 iterations.
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
-   Solving Stokes system (GMG)... 108+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.46859e-16, 8.48765e-08, 8.48765e-08, 1.46642e-07, 0, 0, 0, 1.3304e-05
-      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.3304e-05
+   Solving Stokes system (GMG)... 109+0 iterations.
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.46859e-16, 8.73173e-08, 8.73173e-08, 1.43995e-07, 0, 0, 0, 1.35446e-05
+      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.35446e-05
 
    Solving temperature system... 0 iterations.
    Solving ve_stress_xx system ... 1 iterations.
@@ -119,13 +119,13 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_xx_old system ... 1 iterations.
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
-   Solving Stokes system (GMG)... 24+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.59591e-16, 5.86833e-06, 5.86833e-06, 3.49505e-06, 0, 0, 0, 4.27217e-11
-      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.86833e-06
+   Solving Stokes system (GMG)... 25+0 iterations.
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.55204e-16, 5.84216e-06, 5.84216e-06, 3.43134e-06, 0, 0, 0, 4.32824e-11
+      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.84216e-06
 
 
    Postprocessing:
-     Compositions min/max/mass: -3.512e+04/3.512e+04/-1.351e+08 // -3.512e+04/3.512e+04/1.351e+08 // -2.271e+04/3.336e+04/3.431e+12 // -3.512e+04/3.512e+04/-1.351e+08 // -3.512e+04/3.512e+04/1.351e+08 // -2.271e+04/3.336e+04/3.431e+12
+     Compositions min/max/mass: -3.612e+04/4.354e+04/-1.335e+08 // -3.612e+04/4.354e+04/1.367e+08 // -2.271e+04/3.336e+04/3.431e+12 // -3.612e+04/4.354e+04/-1.335e+08 // -3.612e+04/4.354e+04/1.367e+08 // -2.271e+04/3.336e+04/3.431e+12
      Temperature min/avg/max:   293 K, 293 K, 293 K
      Topography min/max:        -7.451e-09 m, 9.313e-10 m
      RMS, max velocity:         0.0109 m/year, 0.0377 m/year
@@ -144,12 +144,12 @@ Number of mesh deformation degrees of freedom: 9,288
    Solving ve_stress_yy_old system ... 1 iterations.
    Solving ve_stress_xy_old system ... 1 iterations.
    Solving Stokes system (GMG)... 113+0 iterations.
-      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.30109e-16, 4.08776e-07, 4.08776e-07, 3.87956e-07, 0, 0, 0, 7.76729e-06
-      Relative nonlinear residual (total system) after nonlinear iteration 1: 7.76729e-06
+      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.33736e-16, 4.05723e-07, 4.05723e-07, 3.7722e-07, 0, 0, 0, 8.08466e-06
+      Relative nonlinear residual (total system) after nonlinear iteration 1: 8.08466e-06
 
 
    Postprocessing:
-     Compositions min/max/mass: -1.06e+05/1.06e+05/-7.371e+07 // -1.06e+05/1.06e+05/7.371e+07 // -8.863e+04/5.522e+04/-4.078e+14 // -1.06e+05/1.06e+05/-7.371e+07 // -1.06e+05/1.06e+05/7.371e+07 // -8.863e+04/5.522e+04/-4.078e+14
+     Compositions min/max/mass: -1.066e+05/1.066e+05/-1.526e+08 // -1.066e+05/1.066e+05/-1.45e+07 // -8.863e+04/5.229e+04/-4.078e+14 // -1.066e+05/1.066e+05/-1.526e+08 // -1.066e+05/1.066e+05/-1.45e+07 // -8.863e+04/5.229e+04/-4.078e+14
      Temperature min/avg/max:   293 K, 293 K, 293 K
      Topography min/max:        -0.09426 m, 0.02027 m
      RMS, max velocity:         0.011 m/year, 0.0383 m/year

tests/boundary_traction_function_spherical_free_surface/statistics

@@ -42,6 +42,6 @@
 # 42: Maximum topography (m)
 # 43: RMS velocity (m/year)
 # 44: Max. velocity (m/year)
-0 0.000000000000e+00 0.000000000000e+00 4380 41098 4644 236520 2 0 0 0 0 0 0 0 216 219 219  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00 0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00 0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.40488011e-15 -9.31322575e-10 9.31322575e-10 3.44414753e-05 5.26615309e-04 
-1 2.500000000000e+00 2.500000000000e+00 4380 41098 4644 236520 2 0 2 2 2 2 2 2 132 135 135 -3.51213277e+04 3.51213271e+04 -1.35090965e+08 -3.51213271e+04 3.51213277e+04 1.35090965e+08 -2.27140438e+04 3.33570536e+04  3.43054006e+12 -3.51213281e+04 3.51213275e+04 -1.35090998e+08 -3.51213275e+04 3.51213281e+04 1.35090998e+08 -2.27140442e+04 3.33570533e+04  3.43053620e+12 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.61580189e-15 -7.45058060e-09 9.31322575e-10 1.08910272e-02 3.76945054e-02 
-2 5.000000000000e+00 2.500000000000e+00 4380 41098 4644 236520 1 0 1 1 1 1 1 1 113 115 115 -1.06045937e+05 1.06045937e+05 -7.37124683e+07 -1.06045937e+05 1.06045937e+05 7.37124683e+07 -8.86308305e+04 5.52207115e+04 -4.07785493e+14 -1.06045914e+05 1.06045914e+05 -7.37125503e+07 -1.06045914e+05 1.06045914e+05 7.37125503e+07 -8.86308305e+04 5.52207121e+04 -4.07787409e+14 2.93000000e+02 2.93000000e+02 2.93000000e+02 -1.35592571e-16 -9.42646312e-02 2.02718014e-02 1.10291575e-02 3.82968158e-02 
+0 0.000000000000e+00 0.000000000000e+00 4380 41098 4644 236520 2 0 0 0 0 0 0 0 216 219 219  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00  0.00000000e+00 0.00000000e+00  0.00000000e+00 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.40488011e-15 -9.31322575e-10 9.31322575e-10 3.44414753e-05 5.26615309e-04 
+1 2.500000000000e+00 2.500000000000e+00 4380 41098 4644 236520 2 0 2 2 2 2 2 2 134 137 137 -3.61236472e+04 4.35399634e+04 -1.33521411e+08 -3.61236462e+04 4.35399638e+04  1.36660518e+08 -2.27140439e+04 3.33570525e+04  3.43053703e+12 -3.61236472e+04 4.35399636e+04 -1.33521445e+08 -3.61236462e+04 4.35399641e+04  1.36660551e+08 -2.27140442e+04 3.33570523e+04  3.43053412e+12 2.93000000e+02 2.93000000e+02 2.93000000e+02  3.61580189e-15 -7.45058060e-09 9.31322575e-10 1.08910239e-02 3.76945053e-02 
+2 5.000000000000e+00 2.500000000000e+00 4380 41098 4644 236520 1 0 1 1 1 1 1 1 113 115 115 -1.06640941e+05 1.06573682e+05 -1.52557581e+08 -1.06640941e+05 1.06573682e+05 -1.45034263e+07 -8.86307931e+04 5.22853238e+04 -4.07787457e+14 -1.06640918e+05 1.06573644e+05 -1.52557651e+08 -1.06640918e+05 1.06573644e+05 -1.45033571e+07 -8.86307931e+04 5.22853244e+04 -4.07789373e+14 2.93000000e+02 2.93000000e+02 2.93000000e+02 -1.50658412e-16 -9.42646312e-02 2.02718014e-02 1.10291555e-02 3.82968120e-02 

tests/boundary_traction_function_spherical_free_surface/topography.00002

@@ -302,8 +302,8 @@
 4.30419e+06 4.69719e+06 -0.00387885
 4.29735e+06 4.70345e+06 -0.00254236
 4.29735e+06 4.70345e+06 -0.00254236
-4.2905e+06 4.70969e+06 -0.00127911
-4.2905e+06 4.70969e+06 -0.00127911
+4.2905e+06 4.70969e+06 -0.0012791
+4.2905e+06 4.70969e+06 -0.0012791
 4.28365e+06 4.71593e+06 -8.42009e-05
 4.28365e+06 4.71593e+06 -8.42009e-05
 4.27679e+06 4.72216e+06 0.00104675
@@ -320,8 +320,8 @@
 4.24233e+06 4.75313e+06 0.00586973
 4.23541e+06 4.7593e+06 0.00668936
 4.23541e+06 4.7593e+06 0.00668936
-4.22849e+06 4.76545e+06 0.0074667
-4.22849e+06 4.76545e+06 0.0074667
+4.22849e+06 4.76545e+06 0.00746671
+4.22849e+06 4.76545e+06 0.00746671
 4.22155e+06 4.7716e+06 0.008204
 4.22155e+06 4.7716e+06 0.008204
 4.21461e+06 4.77773e+06 0.0089033

tests/continental_extension/screen-output

@@ -36,12 +36,12 @@ Number of mesh deformation degrees of freedom: 462
    Solving upper_crust system ... 9 iterations.
    Solving lower_crust system ... 9 iterations.
    Solving lithospheric_mantle system ... 9 iterations.
-   Initial Newton Stokes residual = 8.15787e+16, v = 8.15787e+16, p = 1.24975e+13
+   Initial Newton Stokes residual = 8.15787e+16, v = 8.15787e+16, p = 1.24042e+13
 
    Rebuilding Stokes preconditioner...
    Solving Stokes system (AMG)... 0+18 iterations.
-      Newton system information: Norm of the rhs: 2.11532e+13
-      Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.000259297
+      Newton system information: Norm of the rhs: 2.11455e+13
+      Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.000259204
 
 
    WARNING: The nonlinear solver in the current timestep failed to converge.