Skip to content

Commit

Permalink
Panzer :: Start reproducer
Browse files Browse the repository at this point in the history
Signed-off-by: Bryan Reuter <bwreute@sandia.gov>
  • Loading branch information
reuterb committed Feb 21, 2025
1 parent ad0de88 commit a8d12b2
Show file tree
Hide file tree
Showing 4 changed files with 354 additions and 0 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -36,6 +36,9 @@ TRIBITS_COPY_FILES_TO_BINARY_DIR(main_driver_files
periodic_wedge.xml
periodic_wedge.pam
energy-periodic-32bit-limit.xml
energy-ss-tp-multiblock-ic-bc-issue.xml
test.jou
make-mesh.sh
EXEDEPS main_driver
)

Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,310 @@
<ParameterList name="Drekar">

<ParameterList name="Mesh">

<Parameter name="Source" type="string" value="Exodus File" />

<ParameterList name="Exodus File">
<Parameter name="File Name" type="string" value="tri3-combined-centers.exo" />
</ParameterList>

</ParameterList>

<ParameterList name="Initial Conditions">
<ParameterList name="left">
<ParameterList name="TEMPERATURE">
<Parameter name="Value" type="double" value="0.0"/>
</ParameterList>
</ParameterList>
<ParameterList name="right">
<ParameterList name="TEMPERATURE">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
</ParameterList>
</ParameterList>

<ParameterList name="Responses">
<ParameterList name="Volume Integral">
<Parameter name="Type" type="string" value="Functional"/>
<Parameter name="Field Name" type="string" value="Volume_Integral"/>
<Parameter name="Element Blocks" type="string" value="left"/>
<Parameter name="Evaluation Types" type="string" value="Residual"/>
</ParameterList>
<ParameterList name="Volume Integral2">
<Parameter name="Type" type="string" value="Functional"/>
<Parameter name="Field Name" type="string" value="Volume_Integral"/>
<Parameter name="Element Blocks" type="string" value="left"/>
<Parameter name="Evaluation Types" type="string" value="Residual"/>
</ParameterList>
</ParameterList>

<ParameterList name="Block ID to Physics ID Mapping">
<Parameter name="left" type="string" value="solid"/>
<Parameter name="right" type="string" value="solid"/>
</ParameterList>

<ParameterList name="Assembly">
<Parameter name="Workset Size" type="int" value="99"/>
<Parameter name="Use Tpetra" type="bool" value="true"/>
</ParameterList>

<ParameterList name="Physics Blocks">

<ParameterList name="solid">

<ParameterList name="EQ 0">
<Parameter name="Type" type="string" value="Energy"/>
<Parameter name="Basis Type" type="string" value="HGrad"/>
<Parameter name="Basis Order" type="int" value="1"/>
<Parameter name="Integration Order" type="int" value="2"/>
<Parameter name="Model ID" type="string" value="fluid model"/>
<Parameter name="Prefix" type="string" value=""/>
</ParameterList>

</ParameterList>

</ParameterList>

<ParameterList name="Closure Models">

<ParameterList name="fluid model">

<ParameterList name="Volume Integral">
</ParameterList>

<ParameterList name="SOURCE_TEMPERATURE">
<Parameter name="Type" type="string" value="Parameter"/>
</ParameterList>
<ParameterList name="Heat Capacity">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
<ParameterList name="Thermal Conductivity">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
<ParameterList name="DENSITY">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
<ParameterList name="HEAT_CAPACITY">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
<ParameterList name="Global Statistics">
<Parameter name="Value" type="string" value="TEMPERATURE,DENSITY"/>
</ParameterList>

</ParameterList>

</ParameterList>

<ParameterList name="User Data">

<ParameterList name="function data one">

</ParameterList>

<ParameterList name="IP Coordinates">
<Parameter name="Integration Order" type="int" value="2"/>
</ParameterList>


</ParameterList>

<ParameterList name="Boundary Conditions">

<ParameterList>
<Parameter name="Type" type="string" value="Dirichlet"/>
<Parameter name="Sideset ID" type="string" value="ymax"/>
<Parameter name="Element Block ID" type="string" value="left"/>
<Parameter name="Equation Set Name" type="string" value="TEMPERATURE"/>
<Parameter name="Strategy" type="string" value="Constant"/>
<ParameterList name="Data">
<Parameter name="Value" type="double" value="0.0"/>
</ParameterList>
</ParameterList>

<ParameterList>
<Parameter name="Type" type="string" value="Dirichlet"/>
<Parameter name="Sideset ID" type="string" value="ymin"/>
<Parameter name="Element Block ID" type="string" value="left"/>
<Parameter name="Equation Set Name" type="string" value="TEMPERATURE"/>
<Parameter name="Strategy" type="string" value="Constant"/>
<ParameterList name="Data">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
</ParameterList>

<ParameterList>
<Parameter name="Type" type="string" value="Dirichlet"/>
<Parameter name="Sideset ID" type="string" value="ymax"/>
<Parameter name="Element Block ID" type="string" value="right"/>
<Parameter name="Equation Set Name" type="string" value="TEMPERATURE"/>
<Parameter name="Strategy" type="string" value="Constant"/>
<ParameterList name="Data">
<Parameter name="Value" type="double" value="0.0"/>
</ParameterList>
</ParameterList>

<ParameterList>
<Parameter name="Type" type="string" value="Dirichlet"/>
<Parameter name="Sideset ID" type="string" value="ymin"/>
<Parameter name="Element Block ID" type="string" value="right"/>
<Parameter name="Equation Set Name" type="string" value="TEMPERATURE"/>
<Parameter name="Strategy" type="string" value="Constant"/>
<ParameterList name="Data">
<Parameter name="Value" type="double" value="1.0"/>
</ParameterList>
</ParameterList>

</ParameterList>

<ParameterList name="Output">
<Parameter name="File Name" type="string" value="energy-ss-tp.exo"/>
<!-- This is inefficient: writes to file every residual evaluation: Used for debug, added to exercise capability -->
<ParameterList name="Cell Average Vectors">
<Parameter name="left" type="string" value="GRAD_TEMPERATURE"/>
</ParameterList>
<!-- -->
</ParameterList>

<ParameterList name="Options">
<Parameter name="Write Volume Assembly Graphs" type="bool" value="true"/>
<Parameter name="Volume Assembly Graph Prefix" type="string" value="energy-ss-tp"/>
</ParameterList>

<ParameterList name="Active Parameters">
<Parameter name="Number of Parameter Vectors" type="int" value="1"/>
<ParameterList name="Parameter Vector 0">
<Parameter name="Number" type="int" value="1"/>
<Parameter name="Parameter 0" type="string" value="SOURCE_TEMPERATURE"/>
<Parameter name="Initial Value 0" type="double" value="1.0"/>
</ParameterList>
</ParameterList>

<ParameterList name="Solution Control">
<Parameter name="Piro Solver" type="string" value="NOX"/>
<Parameter name="Compute Sensitivities" type="bool" value="0"/>
<Parameter name="Jacobian Operator" type="string" value="Have Jacobian"/>
<ParameterList name="LOCA">
<ParameterList name="Bifurcation"/>
<ParameterList name="Constraints"/>
<ParameterList name="Predictor">
<Parameter name="Method" type="string" value="Constant"/>
</ParameterList>
<ParameterList name="Stepper">
<Parameter name="Continuation Method" type="string" value="Natural"/>
<Parameter name="Initial Value" type="double" value="1.0"/>
<Parameter name="Continuation Parameter" type="string" value="Parameter 0"/>
<Parameter name="Max Steps" type="int" value="6"/>
<Parameter name="Max Value" type="double" value="12.25"/>
<Parameter name="Min Value" type="double" value="0.5"/>
<Parameter name="Compute Eigenvalues" type="bool" value="1"/>
<ParameterList name="Eigensolver">
<Parameter name="Method" type="string" value="Anasazi"/>
<Parameter name="Operator" type="string" value="Shift-Invert"/>
<Parameter name="Num Blocks" type="int" value="3"/>
<Parameter name="Num Eigenvalues" type="int" value="1"/>
<Parameter name="Block Size" type="int" value="1"/>
<Parameter name="Maximum Restarts" type="int" value="0"/>
<Parameter name="Shift" type="double" value="1.0"/>
</ParameterList>
</ParameterList>
<ParameterList name="Step Size">
<Parameter name="Initial Step Size" type="double" value="0.5"/>
<Parameter name="Aggressiveness" type="double" value="2.0"/>
</ParameterList>
</ParameterList>
<ParameterList name="NOX">
<ParameterList name="Direction">
<Parameter name="Method" type="string" value="Newton"/>
<ParameterList name="Newton">
<Parameter name="Forcing Term Method" type="string" value="Constant"/>
<Parameter name="Rescue Bad Newton Solve" type="bool" value="1"/>
<ParameterList name="Stratimikos Linear Solver">
<ParameterList name="NOX Stratimikos Options">
</ParameterList>
<ParameterList name="Stratimikos">
<Parameter name="Linear Solver Type" type="string" value="Belos"/>
<Parameter name="Preconditioner Type" type="string" value="Ifpack2"/>
<ParameterList name="Linear Solver Types">
<ParameterList name="Belos">
<Parameter name="Solver Type" type="string" value="Pseudo Block GMRES"/>
<ParameterList name="Solver Types">
<ParameterList name="Pseudo Block GMRES">
<Parameter name="Convergence Tolerance" type="double" value="1e-5"/>
<Parameter name="Output Frequency" type="int" value="10"/>
<Parameter name="Output Style" type="int" value="1"/>
<Parameter name="Verbosity" type="int" value="33"/>
<Parameter name="Maximum Iterations" type="int" value="100"/>
<Parameter name="Block Size" type="int" value="1"/>
<Parameter name="Num Blocks" type="int" value="20"/>
<Parameter name="Flexible Gmres" type="bool" value="0"/>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>
<ParameterList name="Preconditioner Types">
<ParameterList name="Ifpack2">
<Parameter name="Overlap" type="int" value="1"/>
<Parameter name="Prec Type" type="string" value="ILUT"/>
<ParameterList name="Ifpack2 Settings">
<Parameter name="fact: drop tolerance" type="double" value="0"/>
<Parameter name="fact: ilut level-of-fill" type="double" value="1"/>
<Parameter name="fact: level-of-fill" type="int" value="1"/>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>
<ParameterList name="Line Search">
<ParameterList name="Full Step">
<Parameter name="Full Step" type="double" value="1"/>
</ParameterList>
<Parameter name="Method" type="string" value="Full Step"/>
</ParameterList>
<Parameter name="Nonlinear Solver" type="string" value="Line Search Based"/>
<ParameterList name="Printing">
<Parameter name="Output Precision" type="int" value="3"/>
<Parameter name="Output Processor" type="int" value="0"/>
<ParameterList name="Output Information">
<Parameter name="Error" type="bool" value="1"/>
<Parameter name="Warning" type="bool" value="1"/>
<Parameter name="Outer Iteration" type="bool" value="1"/>
<Parameter name="Parameters" type="bool" value="1"/>
<Parameter name="Details" type="bool" value="1"/>
<Parameter name="Linear Solver Details" type="bool" value="1"/>
<Parameter name="Stepper Iteration" type="bool" value="1"/>
<Parameter name="Stepper Details" type="bool" value="1"/>
<Parameter name="Stepper Parameters" type="bool" value="1"/>
</ParameterList>
</ParameterList>
<ParameterList name="Solver Options">
<Parameter name="Status Test Check Type" type="string" value="Minimal"/>
</ParameterList>
<ParameterList name="Status Tests">
<Parameter name="Test Type" type="string" value="Combo"/>
<Parameter name="Combo Type" type="string" value="OR"/>
<Parameter name="Number of Tests" type="int" value="2"/>
<ParameterList name="Test 0">
<Parameter name="Test Type" type="string" value="Combo"/>
<Parameter name="Combo Type" type="string" value="AND"/>
<Parameter name="Number of Tests" type="int" value="2"/>
<ParameterList name="Test 0">
<Parameter name="Test Type" type="string" value="NormF"/>
<Parameter name="Tolerance" type="double" value="1.0e-8"/>
</ParameterList>
<ParameterList name="Test 1">
<Parameter name="Test Type" type="string" value="RelativeNormF"/>
<Parameter name="Tolerance" type="double" value="1.0e-4"/>
</ParameterList>
</ParameterList>
<ParameterList name="Test 1">
<Parameter name="Test Type" type="string" value="MaxIters"/>
<Parameter name="Maximum Iterations" type="int" value="10"/>
</ParameterList>
</ParameterList>
</ParameterList>
</ParameterList>


</ParameterList>
2 changes: 2 additions & 0 deletions packages/panzer/adapters-stk/example/main_driver/make-mesh.sh
Original file line number Diff line number Diff line change
@@ -0,0 +1,2 @@
cubit -nogui -nographics -nojournal -batch test.jou
decomp --processors 3 --rcb tri3-combined-centers.exo
39 changes: 39 additions & 0 deletions packages/panzer/adapters-stk/example/main_driver/test.jou
Original file line number Diff line number Diff line change
@@ -0,0 +1,39 @@
#{ctype = "tri3"}
#{scheme = "triadvance"}

create surface rectangle width 2 height 1 zplane
split surface 1 across location position -0.05 -2 0 location position -0.05 2 0
merge surface 3 2

set duplicate block elements off
block 1 add surface 3
block 1 name "left"
block 1 element type {ctype}

set duplicate block elements off
block 2 add surface 2
block 2 name "right"
block 2 element type {ctype}

sideset 1 add curve 2 wrt surface 3
sideset 1 name "xmin"
sideset 2 add curve 4 wrt surface 2
sideset 2 name "xmax"
sideset 3 add curve 9 wrt surface 3
sideset 3 add curve 6 wrt surface 2
sideset 3 name "ymin"
sideset 4 add curve 8 wrt surface 3
sideset 4 add curve 7 wrt surface 2
sideset 4 name "ymax"
sideset 5 add curve 5 wrt surface all
sideset 5 name "center"

surface all interval 7
surface all scheme {scheme}
mesh surface all

#{tolower(ctype)}

set exodus netcdf4 on

export mesh "{ctype}-combined-centers.exo" dimension 2 overwrite

0 comments on commit a8d12b2

Please sign in to comment.