============================================================================== Test 1: With multiple elastic moduli but only prefactor has multiple values ============================================================================== ==> snippet from prm: subsection Visco Plastic set Prefactors for dislocation creep = 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-24 set Stress exponents for dislocation creep = 1.0 set Activation energies for dislocation creep = 0.0 set Activation volumes for dislocation creep = 0.0 set Elastic shear moduli = 1.0e11, 1.0e11, 1.0e11, 1.0e11, 1.0e10 ===> This runs without any problem ============================================================================== Test 2: With only one elastic moduli, but multiple pre-factors ============================================================================== ==> snippet from prm: subsection Visco Plastic set Prefactors for dislocation creep = 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-24 set Stress exponents for dislocation creep = 1.0 set Activation energies for dislocation creep = 0.0 set Activation volumes for dislocation creep = 0.0 set Elastic shear moduli = 1.0e11 #, 1.0e11, 1.0e11, 1.0e11, 1.0e10 ===> This run aborts in Stokes solve of first time step...here's the output.... (base) /Users/billen/Box-Sync/Aspect-Projects/BendingBeamTests % aspect-debug vep_testnew.prm ----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.5.0-pre (main, 029548bd7) -- . using deal.II 9.4.0 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.0.0 -- . using p4est 2.8.0 -- . running in DEBUG mode -- . running with 1 MPI process ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.5.0-pre&dg=1&sha=029548bd7&src=code ----------------------------------------------------------------------------- Number of active cells: 3,750 (on 1 levels) Number of degrees of freedom: 173,254 (30,502+3,876+3,876+33,750+33,750+33,750+33,750) *** Timestep 0: t=0 years, dt=0 years Solving temperature system... 0 iterations. Skipping ve_stress_xx composition solve because RHS is zero. Skipping ve_stress_yy composition solve because RHS is zero. Skipping ve_stress_xy composition solve because RHS is zero. Solving beam system ... 0 iterations. Rebuilding Stokes preconditioner... Solving Stokes system... zsh: abort /Users/billen/Box-Sync/Aspect-Src/aspect/build-debug/aspect vep_testnew.prm ============================================================================== Test 3: With multiple viscosity inputs for ALL stress components and compositions -------------------------------------------------------------------------------- ==> snippet from prm: subsection Visco Plastic set Prefactors for dislocation creep = 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-18, 0.5e-24 set Stress exponents for dislocation creep = 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 set Activation energies for dislocation creep = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 set Activation volumes for dislocation creep = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 set Elastic shear moduli = 1.0e11, 1.0e11, 1.0e11, 1.0e11, 1.0e10 ==> Aborts before running the first time step... run output... (base) /Users/billen/Box-Sync/Aspect-Projects/BendingBeamTests % aspect-debug vep_testnew.prm ----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.5.0-pre (main, 029548bd7) -- . using deal.II 9.4.0 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.0.0 -- . using p4est 2.8.0 -- . running in DEBUG mode -- . running with 1 MPI process ----------------------------------------------------------------------------- zsh: abort /Users/billen/Box-Sync/Aspect-Src/aspect/build-debug/aspect vep_testnew.prm