Add recoverable testing of complex grid functions

lindsayad · lindsayad · commit a1bca1b87be8 · 2026-04-09T19:11:47.000-07:00
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0000.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0000.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+0,0,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0010.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0010.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+1.7092102224128,1.7092102224128,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0020.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0020.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+1.4604882411775,1.4604882411775,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0030.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0030.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+1.3383797678165,1.3383797678165,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0040.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0040.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+1.2718650483004,1.2718650483004,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0050.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0050.csv
@@ -0,0 +1,2 @@
+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1
+1.2321741962316,1.2321741962316,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0000.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0000.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+0,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0010.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0010.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+1.7092102224128,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0020.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0020.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+1.4604882411775,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0030.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0030.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+1.3383797678165,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0040.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0040.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+1.2718650483004,0,0
diff --git a/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0050.csv b/test/tests/mfem/kernels/gold/OutputData/NLHeatConductionComplexAux_centre_temperature_0050.csv
@@ -0,0 +1,2 @@
+temperature_0,x_0,x_1
+1.2321741962316,0,0
diff --git a/test/tests/mfem/kernels/nl_heatconduction_complex_aux.i b/test/tests/mfem/kernels/nl_heatconduction_complex_aux.i
@@ -0,0 +1,127 @@
+# Transient nonlinear heat conduction with a complex auxiliary variable.
+# The complex aux variable exercises the cmplx_gridfunctions checkpoint/recover code path
+# in SolutionStateData.C without requiring a complex primary solve.
+
+kappa = 0.5
+alpha = 1e-2
+
+[Mesh]
+  type = MFEMMesh
+  file = ../mesh/star.mesh
+  uniform_refine = 1
+[]
+
+[Problem]
+  type = MFEMProblem
+[]
+
+[FESpaces]
+  [H1FESpace]
+    type = MFEMScalarFESpace
+    fec_type = H1
+    fec_order = SECOND
+  []
+[]
+
+[Variables]
+  [temperature]
+    type = MFEMVariable
+    fespace = H1FESpace
+  []
+[]
+
+[AuxVariables]
+  [complex_temperature]
+    type = MFEMComplexVariable
+    fespace = H1FESpace
+  []
+[]
+
+[AuxKernels]
+  [project_complex_temperature]
+    type = MFEMComplexScalarProjectionAux
+    variable = complex_temperature
+    coefficient_real = temperature
+    coefficient_imag = temperature
+    execute_on = TIMESTEP_END
+  []
+[]
+
+[Functions]
+  [initial]
+    type = ParsedFunction
+    expression = 'if((x*x + y*y > 0.251), 1.0, 2.0)'
+  []
+  [diffusivity_temperature_dependence]
+    type = MFEMParsedFunction
+    expression = 'alpha * temperature'
+    symbol_names = 'alpha temperature'
+    symbol_values = '${alpha} temperature'
+  []
+[]
+
+[ICs]
+  [diffused_ic]
+    type = MFEMScalarIC
+    coefficient = initial
+    variable = temperature
+  []
+[]
+
+[Kernels]
+  [nl_diffusion]
+    type = MFEMNLDiffusionKernel
+    variable = temperature
+    k_coefficient = diffusivity_temperature_dependence
+    dk_du_coefficient = ${alpha}
+  []
+  [linear_diffusion]
+    type = MFEMDiffusionKernel
+    variable = temperature
+    coefficient = ${kappa}
+  []
+  [dT_dt]
+    type = MFEMTimeDerivativeMassKernel
+    variable = temperature
+  []
+[]
+
+[Solver]
+  type = MFEMMUMPS
+  print_level = 0
+[]
+
+[Executioner]
+  type = MFEMTransient
+  device = cpu
+  assembly_level = legacy
+  dt = 1e-2
+  start_time = 0.0
+  end_time = 0.5
+
+  nl_max_its = 30
+  nl_abs_tol = 1.0e-5
+  nl_rel_tol = 1.0e-5
+  print_level = 1
+[]
+
+[VectorPostprocessors]
+  [centre_temperature]
+    type = MFEMPointValueSampler
+    variable = 'temperature'
+    points = '0.0 0.0 0.0'
+    execute_on = TIMESTEP_END
+  []
+  [centre_complex_temperature]
+    type = MFEMComplexPointValueSampler
+    variable = 'complex_temperature'
+    points = '0.0 0.0 0.0'
+    execute_on = TIMESTEP_END
+  []
+[]
+
+[Outputs]
+  file_base = OutputData/NLHeatConductionComplexAux
+  csv = true
+  time_step_interval = 10
+[]
diff --git a/test/tests/mfem/kernels/tests b/test/tests/mfem/kernels/tests
@@ -277,4 +277,26 @@
     compute_devices = 'cpu'
     restep = false
   []
+  [MFEMNonLinearHeatConductionComplexAux]
+    design = 'EquationSystem.md'
+    issues = '#30076'
+    requirement = 'The system shall correctly checkpoint and recover a transient MFEM problem that includes a complex auxiliary variable, exercising the cmplx_gridfunctions restore path.'
+    type = CSVDiff
+    input = nl_heatconduction_complex_aux.i
+    csvdiff = 'OutputData/NLHeatConductionComplexAux_centre_temperature_0000.csv
+                OutputData/NLHeatConductionComplexAux_centre_temperature_0010.csv
+                OutputData/NLHeatConductionComplexAux_centre_temperature_0020.csv
+                OutputData/NLHeatConductionComplexAux_centre_temperature_0030.csv
+                OutputData/NLHeatConductionComplexAux_centre_temperature_0040.csv
+                OutputData/NLHeatConductionComplexAux_centre_temperature_0050.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0000.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0010.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0020.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0030.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0040.csv
+                OutputData/NLHeatConductionComplexAux_centre_complex_temperature_0050.csv'
+    capabilities = 'mfem'
+    compute_devices = 'cpu'
+    restep = false
+  []
 []

Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+0,0,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+1.7092102224128,1.7092102224128,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+1.4604882411775,1.4604882411775,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+1.3383797678165,1.3383797678165,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+1.2718650483004,1.2718650483004,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+complex_temperature_imag_0,complex_temperature_real_0,x_0,x_1`
	`2`	`+1.2321741962316,1.2321741962316,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+temperature_0,x_0,x_1`
	`2`	`+1.7092102224128,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+temperature_0,x_0,x_1`
	`2`	`+1.4604882411775,0,0`
Original file line number	Diff line number	Diff line change
`@@ -0,0 +1,2 @@`
	`1`	`+temperature_0,x_0,x_1`
	`2`	`+1.3383797678165,0,0`