Skip to content
Snippets Groups Projects
Commit 1a9c7b5c authored by Tom Fischer's avatar Tom Fischer
Browse files

[T] XDMF/HDF5 gml output test

parent 3565937c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 270 additions and 0 deletions
ProcessLib/SteadyStateDiffusion/Tests.cmake
+ 13
− 0
View file @ 1a9c7b5c
......@@ -963,6 +963,19 @@ AddTest(
cube_1x1x1_hex_1e2_inner_bar_domain_ts_1_t_1.000000.vtu cube_1x1x1_hex_1e2_inner_bar_domain_3d_submesh_output_test_ts_1_t_1.000000.vtu v v 1e-13 1e-13
)
AddTest(
NAME SteadyStateDiffusion_square_1x1_1e2_GMRES_GML_output_xdmf-hdf5
PATH Elliptic/square_1x1_SteadyStateDiffusion
EXECUTABLE ogs
EXECUTABLE_ARGS square_1e2_GMRES_GML_output_xdmf-hdf5.prj
WRAPPER time
TESTER xdmfdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
square_1x1_quad_1e2_GMRES_GML_output_square_1x1_quad_1e2.xdmf square_1x1_quad_1e2_GMRES_GML_output_square_1x1_quad_1e2.xdmf pressure pressure 1e-14 1e-14
square_1x1_quad_1e2_GMRES_GML_output_square_1x1_geometry_left.xdmf square_1x1_quad_1e2_GMRES_GML_output_square_1x1_geometry_left.xdmf pressure pressure 1e-14 1e-14
)
if(OGS_USE_MPI AND OGS_USE_PIP)
NotebookTest(NOTEBOOKFILE Notebooks/SimplePETSc.ipynb RUNTIME 10)
else()
......
Tests/Data/Elliptic/square_1x1_SteadyStateDiffusion/square_1e2_GMRES_GML_output_xdmf-hdf5.prj 0 → 100644
+ 137
− 0
View file @ 1a9c7b5c
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh>square_1x1_quad_1e2.vtu</mesh>
<geometry>square_1x1.gml</geometry>
<processes>
<process>
<name>SteadyStateDiffusion</name>
<type>STEADY_STATE_DIFFUSION</type>
<integration_order>2</integration_order>
<process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="v"/>
</secondary_variables>
</process>
</processes>
<media>
<medium id="0">
<phases/>
<properties>
<property>
<name>diffusion</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="SteadyStateDiffusion">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<abstol>1.e-6</abstol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>SingleStep</type>
</time_stepping>
</process>
</processes>
<output>
<type>XDMF</type>
<prefix>{:meshname}_GMRES_GML_output</prefix>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
<meshes>
<mesh>square_1x1_quad_1e2</mesh>
</meshes>
<geometrical_sets>
<geometrical_set>
<name>square_1x1_geometry</name>
<geometry>left</geometry>
</geometrical_set>
</geometrical_sets>
<variables>
<variable> pressure </variable>
<variable> v </variable>
</variables>
</output>
</time_loop>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>10</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
<eigen>
<solver_type>GMRES</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-16</error_tolerance>
<restart>30</restart>
</eigen>
<petsc>
<prefix>gw</prefix>
<parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
<parameters>
<parameter>
<name>p0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p_Dirichlet_left</name>
<type>Constant</type>
<value>1</value>
</parameter>
<parameter>
<name>p_Dirichlet_right</name>
<type>Constant</type>
<value>-1</value>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>p0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>square_1x1_geometry</geometrical_set>
<geometry>left</geometry>
<type>Dirichlet</type>
<parameter>p_Dirichlet_left</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>square_1x1_geometry</geometrical_set>
<geometry>right</geometry>
<type>Dirichlet</type>
<parameter>p_Dirichlet_right</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
</OpenGeoSysProject>
Tests/Data/Elliptic/square_1x1_SteadyStateDiffusion/square_1x1_quad_1e2_GMRES_GML_output.h5 0 → 100644
+ 0
− 0
View file @ 1a9c7b5c
File added
Tests/Data/Elliptic/square_1x1_SteadyStateDiffusion/square_1x1_quad_1e2_GMRES_GML_output_square_1x1_geometry_left.xdmf 0 → 100644
+ 58
− 0
View file @ 1a9c7b5c
<?xml version="1.0" encoding="utf-8"?>
<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
<Domain>
<Information Name="OGS_VERSION" Value="6.4.3-761-g38d23b54.dirty"/>
<Grid CollectionType="Temporal" GridType="Collection" Name="Collection">
<Grid Name="Grid" GridType="Uniform">
<Time Value="0"/>
<Geometry Origin="" Type="XYZ">
<DataItem DataType="Float" Dimensions="11 3" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/geometry|0 0 0:1 1 1:1 11 3:2 11 3</DataItem>
</Geometry>
<Topology Dimensions="40" Type="Mixed">
<DataItem DataType="Int" Dimensions="40" Format="HDF" Precision="4">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/topology|0 0:1 1:1 40:2 40</DataItem>
</Topology>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="D1_left_bottom_N1_right" Type="None">
<DataItem DataType="Float" Dimensions="11" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/D1_left_bottom_N1_right|0 0:1 1:1 11:2 11</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="Linear_1_to_minus1" Type="None">
<DataItem DataType="Float" Dimensions="11" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/Linear_1_to_minus1|0 0:1 1:1 11:2 11</DataItem>
</Attribute>
<Attribute Center="Cell" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="bulk_element_ids" Type="None">
<DataItem DataType="UInt" Dimensions="10" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/bulk_element_ids|0 0:1 1:1 10:2 10</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="bulk_node_ids" Type="None">
<DataItem DataType="UInt" Dimensions="11" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/bulk_node_ids|0 0:1 1:1 11:2 11</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
<DataItem DataType="Float" Dimensions="11" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/pressure|0 0:1 1:1 11:2 11</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
<DataItem DataType="Float" Dimensions="11 2" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/v|0 0 0:1 1 1:1 11 2:2 11 2</DataItem>
</Attribute>
</Grid>
<Grid Name="Grid" GridType="Uniform">
<Time Value="1"/>
<xi:include xpointer="element(/1/1/2/1/2)"/>
<xi:include xpointer="element(/1/1/2/1/3)"/>
<xi:include xpointer="element(/1/1/2/1/4)"/>
<xi:include xpointer="element(/1/1/2/1/5)"/>
<xi:include xpointer="element(/1/1/2/1/6)"/>
<xi:include xpointer="element(/1/1/2/1/7)"/>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
<DataItem DataType="Float" Dimensions="11" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/pressure|1 0:1 1:1 11:2 11</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
<DataItem DataType="Float" Dimensions="11 2" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_geometry_left/v|1 0 0:1 1 1:1 11 2:2 11 2</DataItem>
</Attribute>
</Grid>
</Grid>
</Domain>
</Xdmf>
\ No newline at end of file
Tests/Data/Elliptic/square_1x1_SteadyStateDiffusion/square_1x1_quad_1e2_GMRES_GML_output_square_1x1_quad_1e2.xdmf 0 → 100644
+ 62
− 0
View file @ 1a9c7b5c
<?xml version="1.0" encoding="utf-8"?>
<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
<Domain>
<Information Name="OGS_VERSION" Value="6.4.3-761-g38d23b54.dirty"/>
<Grid CollectionType="Temporal" GridType="Collection" Name="Collection">
<Grid Name="Grid" GridType="Uniform">
<Time Value="0"/>
<Geometry Origin="" Type="XYZ">
<DataItem DataType="Float" Dimensions="121 3" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/geometry|0 0 0:1 1 1:1 121 3:2 121 3</DataItem>
</Geometry>
<Topology Dimensions="500" Type="Mixed">
<DataItem DataType="Int" Dimensions="500" Format="HDF" Precision="4">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/topology|0 0:1 1:1 500:2 500</DataItem>
</Topology>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="D1_left_bottom_N1_right" Type="None">
<DataItem DataType="Float" Dimensions="121" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/D1_left_bottom_N1_right|0 0:1 1:1 121:2 121</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="Linear_1_to_minus1" Type="None">
<DataItem DataType="Float" Dimensions="121" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/Linear_1_to_minus1|0 0:1 1:1 121:2 121</DataItem>
</Attribute>
<Attribute Center="Cell" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="MaterialIDs" Type="None">
<DataItem DataType="Int" Dimensions="100" Format="HDF" Precision="4">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/MaterialIDs|0 0:1 1:1 100:2 100</DataItem>
</Attribute>
<Attribute Center="Cell" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="bulk_element_ids" Type="None">
<DataItem DataType="UInt" Dimensions="100" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/bulk_element_ids|0 0:1 1:1 100:2 100</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="bulk_node_ids" Type="None">
<DataItem DataType="UInt" Dimensions="121" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/bulk_node_ids|0 0:1 1:1 121:2 121</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
<DataItem DataType="Float" Dimensions="121" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/pressure|0 0:1 1:1 121:2 121</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
<DataItem DataType="Float" Dimensions="121 2" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/v|0 0 0:1 1 1:1 121 2:2 121 2</DataItem>
</Attribute>
</Grid>
<Grid Name="Grid" GridType="Uniform">
<Time Value="1"/>
<xi:include xpointer="element(/1/1/2/1/2)"/>
<xi:include xpointer="element(/1/1/2/1/3)"/>
<xi:include xpointer="element(/1/1/2/1/4)"/>
<xi:include xpointer="element(/1/1/2/1/5)"/>
<xi:include xpointer="element(/1/1/2/1/6)"/>
<xi:include xpointer="element(/1/1/2/1/7)"/>
<xi:include xpointer="element(/1/1/2/1/8)"/>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
<DataItem DataType="Float" Dimensions="121" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/pressure|1 0:1 1:1 121:2 121</DataItem>
</Attribute>
<Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
<DataItem DataType="Float" Dimensions="121 2" Format="HDF" Precision="8">square_1x1_quad_1e2_GMRES_GML_output.h5:/meshes/square_1x1_quad_1e2/v|1 0 0:1 1 1:1 121 2:2 121 2</DataItem>
</Attribute>
</Grid>
</Grid>
</Domain>
</Xdmf>
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment