PlateHole

2D-serial

Two-dimensional

Serial

Validated using check script

./bin/mechanics-2d-g++ tests/PlateHole/input

2D-serial-coverage

Two-dimensional

Serial

Not validated

./bin/mechanics-2d-g++ tests/PlateHole/input  solver.fixed_iter="1"

3D-parallel

Three-dimensional

Parallel (4 procs)

Not validated

13.52s (beaker) 11.48s (statler) 7.05s (waldorf)

mpiexec -np 4 ./bin/mechanics-3d-g++ tests/PlateHole/input

3D-serial

Three-dimensional

Serial

Not validated

42.00s (beaker) 21.28s (waldorf)

./bin/mechanics-3d-g++ tests/PlateHole/input
Input file (../../tests/PlateHole/input)
#@  [2D-serial]
#@  exe    = mechanics
#@  dim    = 2 
#@  nprocs = 1
#@  check  = true
#@ 
#@  [2D-serial-coverage]
#@  exe    = mechanics
#@  dim    = 2 
#@  nprocs = 1
#@  check  = false
#@  coverage = true
#@  args=solver.fixed_iter = 1
#@
#@  [3D-parallel]
#@  exe    = mechanics
#@  dim    = 3
#@  nprocs = 4
#@  check  = false
#@  benchmark-beaker = 13.52
#@  benchmark-statler = 11.48
#@  benchmark-waldorf = 7.05
#@
#@  [3D-serial]
#@  exe    = mechanics
#@  dim    = 3
#@  nprocs = 1
#@  check  = false
#@  benchmark-beaker = 42.00
#@  benchmark-waldorf = 21.28




alamo.program = mechanics
alamo.program.mechanics.model = linear.isotropic

plot_file		    = tests/PlateHole/output

# this is not a time integration, so do
# exactly one timestep and then quit
timestep		    = 0.1
stop_time		    = 0.1
amr.plot_int		    = 1

amr.blocking_factor	    = 8
amr.regrid_int		    = -1
amr.grid_eff		    = 1.0
amr.node.all                = 1
amr.cell.all                = 1

# grid
amr.max_level		    = 5 #8
amr.n_cell		    = 32 32 32

# geometry
geometry.prob_lo	    = -16 -16 -16
geometry.prob_hi	    = 16 16 16
geometry.is_periodic	    = 0 0 0

# ellipse configuration
psi.ic.type        = ellipse
psi.ic.ellipse.a   = 1.0 1.0 1.0  # ellipse radius
psi.ic.ellipse.x0  = 0 0 0 # location of ellipse center
psi.ic.ellipse.eps = 0.2
psi.ic.ellipse.invert = 1

# elastic moduli
nmodels = 1
model1.E = 1.0
model1.nu = 0.3
model1.planestress = 1

solver.verbose = 3
solver.nriters = 1
solver.max_iter = 20
solver.bottom_solver = smoother
solver.tol_abs = 1E-16 # This is very important for near-singular problems!!!

print_residual = 1

bc.type = tensiontest
bc.tensiontest.type = uniaxial_stress
bc.tensiontest.disp=0.01

elasticop.small=1E-16