[cig-commits] [commit] master: remove accidentally added broken test (d14051b)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Sun Jun 8 07:15:13 PDT 2014


Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/277ad233ad48d1c5db9e23cd1e74d6914cb94919...d14051bb8dc3b2a63654bfac6174f2c4371133f0

>---------------------------------------------------------------

commit d14051bb8dc3b2a63654bfac6174f2c4371133f0
Author: Timo Heister <timo.heister at gmail.com>
Date:   Sun Jun 8 10:15:06 2014 -0400

    remove accidentally added broken test


>---------------------------------------------------------------

d14051bb8dc3b2a63654bfac6174f2c4371133f0
 tests/petsc_direct_solver.prm           | 164 --------------------------------
 tests/petsc_direct_solver/screen-output |  98 -------------------
 tests/petsc_direct_solver/statistics    |  18 ----
 3 files changed, 280 deletions(-)

diff --git a/tests/petsc_direct_solver.prm b/tests/petsc_direct_solver.prm
deleted file mode 100644
index 9aab906..0000000
--- a/tests/petsc_direct_solver.prm
+++ /dev/null
@@ -1,164 +0,0 @@
-# test that the direct solver works in parallel with PETSc
-
-# LINEAR ALGEBRA: PETSC
-# MPI: 3
-
-set Dimension = 2
-
-set Use direct solver for Stokes system = true
-
-set Start time                             = 0
-set End time                               = 6e15
-set Use years in output instead of seconds = false
-
-set Output directory                       = output
-set Pressure normalization                 = surface
-
-
-subsection Geometry model
-  set Model name = box
-
-  subsection Box
-    set X extent = 1000000
-    set Y extent = 1000000
-  end
-end
-
-
-subsection Gravity model
-  set Model name = vertical
-  subsection Vertical
-    set Magnitude = 10.0
-  end
-end
-
-
-subsection Model settings
-
-  # As we only want to look at the effects of latent heating, we disable all
-  # the other heating terms. 
-  set Include adiabatic heating               = false
-  set Include latent heat                     = true
-  set Include shear heating                   = false
-
-  # We only fix the temperature at the upper boundary, the other boundaries
-  # are isolating. To guarantuee a steady downward flow, we fix the velocity
-  # at the top and bottom, and set it to free slip on the sides. 
-  set Fixed temperature boundary indicators   = 3
-  set Prescribed velocity boundary indicators = 2:function, 3:function
-  set Tangential velocity boundary indicators = 0, 1
-end
-
-
-############### Boundary conditions
-# We set the top temperature to T1=1000K. 
-subsection Boundary temperature model
-  set Model name = box
-  subsection Box
-    set Top temperature = 1000
-  end
-end
-
-# We prescribe a constant downward flow.
-subsection Boundary velocity model
-  subsection Function
-    set Function expression = 0;-2.1422e-11
-    set Variable names      = x,y
-  end
-end
-
-subsection Initial conditions
-  set Model name = function
-  subsection Function
-    set Function expression = 1000.0
-    set Variable names      = x,y
-  end
-end
-
-
-subsection Material model
-  set Model name = latent heat
-  subsection Latent heat
-
-    # The change of density across the phase transition. Together with the
-    # Clapeyron slope, this is what determines the entropy change.
-    set Phase transition density jumps                 = 115.6
-    set Corresponding phase for density jump           = 0
-
-    # If the temperature is equal to the phase transition temperature, the 
-    # phase transition will occur at the phase transition depth. However, 
-    # if the temperature deviates from this value, the Clapeyron slope 
-    # determines how much the pressure (and depth) of the phase boundary
-    # changes. Here, the phase transition will be in the middle of the box
-    # for T=T1. 
-    set Phase transition depths                        = 500000
-    set Phase transition temperatures                  = 1000
-    set Phase transition Clapeyron slopes              = 1e7
-
-    # We set the width of the phase transition to 5 km. You may want to 
-    # change this parameter to see how latent heating depends on the width
-    # of the phase transition. 
-    set Phase transition widths                        = 5000
-
-    set Reference density                              = 3400
-    set Reference specific heat                        = 1000
-    set Reference temperature                          = 1000
-    set Thermal conductivity                           = 2.38
-
-    # We set the thermal expansion amd the compressibility to zero, so that 
-    # all temperature (and density) changes are caused by advection, diffusion 
-    # and latent heating. 
-    set Thermal expansion coefficient                  = 0.0
-    set Compressibility                                = 0.0
-
-    # Viscosity is constant. 
-    set Thermal viscosity exponent                     = 0.0
-    set Viscosity                                      = 8.44e21
-    set Viscosity prefactors                           = 1.0, 1.0
-    set Composition viscosity prefactor                = 1.0
-    set Activation enthalpies                          = 3.9473e-3, 3.9473e-3 
-  end
-end
-
-
-subsection Mesh refinement
-  set Initial adaptive refinement        = 0 
-  set Initial global refinement          = 5
-  set Time steps between mesh refinement = 0
-
-end
-
-
-subsection Discretization
-  subsection Stabilization parameters
-    # The exponent $\alpha$ in the entropy viscosity stabilization. Units:
-    # None.
-    set alpha = 2
-
-    # The $\beta$ factor in the artificial viscosity stabilization. An
-    # appropriate value for 2d is 0.052 and 0.078 for 3d. Units: None.
-    set beta  = 0.078
-
-    # The $c_R$ factor in the entropy viscosity stabilization. Units: None.
-    set cR    = 0.5   # default: 0.11
-  end
-end
-
-
-subsection Postprocess
-
-  set List of postprocessors = visualization
-
-  subsection Visualization
-
-    set Number of grouped files       = 0
-    set Output format                 = vtu
-
-    # We are only interested in the last timestep (when the system hast reached
-    # a steady state). For following the development of the system or checking
-    # if the solution already reached steady state, this parameter can be set 
-    # to a smaller value.
-    set Time between graphical output = 0
-    set List of output variables      = density
-  end
-end
diff --git a/tests/petsc_direct_solver/screen-output b/tests/petsc_direct_solver/screen-output
deleted file mode 100644
index 0ca5c9b..0000000
--- a/tests/petsc_direct_solver/screen-output
+++ /dev/null
@@ -1,98 +0,0 @@
------------------------------------------------------------------------------
--- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
---     . running in DEBUG mode
---     . running with 3 MPI processes
---     . using Trilinos
------------------------------------------------------------------------------
-
-Number of active cells: 1,024 (on 6 levels)
-Number of degrees of freedom: 13,764 (9,539+4,225)
-
-*** Timestep 0:  t=0 seconds
-   Solving temperature system... 0 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00000
-
-*** Timestep 1:  t=7.26222e+14 seconds
-   Solving temperature system... 189 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00001
-
-*** Timestep 2:  t=1.4246e+15 seconds
-   Solving temperature system... 22 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00002
-
-*** Timestep 3:  t=2.11579e+15 seconds
-   Solving temperature system... 20 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00003
-
-*** Timestep 4:  t=2.79859e+15 seconds
-   Solving temperature system... 21 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00004
-
-*** Timestep 5:  t=3.48103e+15 seconds
-   Solving temperature system... 20 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00005
-
-*** Timestep 6:  t=4.16434e+15 seconds
-   Solving temperature system... 20 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00006
-
-*** Timestep 7:  t=4.85169e+15 seconds
-   Solving temperature system... 19 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00007
-
-*** Timestep 8:  t=5.54591e+15 seconds
-   Solving temperature system... 19 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00008
-
-*** Timestep 9:  t=6e+15 seconds
-   Solving temperature system... 16 iterations.
-   Solving Stokes system... done.
-
-   Postprocessing:
-     Writing graphical output: output-direct_solver_2/solution-00009
-
-Termination requested by criterion: end time
-
-
-+---------------------------------------------+------------+------------+
-| Total wallclock time elapsed since start    |      3.75s |            |
-|                                             |            |            |
-| Section                         | no. calls |  wall time | % of total |
-+---------------------------------+-----------+------------+------------+
-| Assemble Stokes system          |        10 |      0.65s |        17% |
-| Assemble temperature system     |        10 |     0.497s |        13% |
-| Build temperature preconditioner|        10 |    0.0327s |      0.87% |
-| Solve Stokes system             |        10 |      1.88s |        50% |
-| Solve temperature system        |        10 |    0.0984s |       2.6% |
-| Initialization                  |         2 |    0.0393s |         1% |
-| Postprocessing                  |        10 |      0.24s |       6.4% |
-| Setup dof systems               |         1 |    0.0692s |       1.8% |
-+---------------------------------+-----------+------------+------------+
-
diff --git a/tests/petsc_direct_solver/statistics b/tests/petsc_direct_solver/statistics
deleted file mode 100644
index 3735ebc..0000000
--- a/tests/petsc_direct_solver/statistics
+++ /dev/null
@@ -1,18 +0,0 @@
-# 1: Time step number
-# 2: Time (seconds)
-# 3: Number of mesh cells
-# 4: Number of Stokes degrees of freedom
-# 5: Number of temperature degrees of freedom
-# 6: Iterations for temperature solver
-# 7: Time step size (seconds)
-# 8: Visualization file name
-0 0.0000e+00 1024 9539 4225   0 7.2622e+14 output-direct_solver_2/solution-00000 
-1 7.2622e+14 1024 9539 4225 189 6.9838e+14 output-direct_solver_2/solution-00001 
-2 1.4246e+15 1024 9539 4225  22 6.9119e+14 output-direct_solver_2/solution-00002 
-3 2.1158e+15 1024 9539 4225  20 6.8280e+14 output-direct_solver_2/solution-00003 
-4 2.7986e+15 1024 9539 4225  21 6.8244e+14 output-direct_solver_2/solution-00004 
-5 3.4810e+15 1024 9539 4225  20 6.8331e+14 output-direct_solver_2/solution-00005 
-6 4.1643e+15 1024 9539 4225  20 6.8734e+14 output-direct_solver_2/solution-00006 
-7 4.8517e+15 1024 9539 4225  19 6.9422e+14 output-direct_solver_2/solution-00007 
-8 5.5459e+15 1024 9539 4225  19 4.5409e+14 output-direct_solver_2/solution-00008 
-9 6.0000e+15 1024 9539 4225  16 7.1111e+14 output-direct_solver_2/solution-00009 



More information about the CIG-COMMITS mailing list