Fixed resonator example for multiple MPI processes.

share/examples/resonator/params_resonator_gmshfem_comparison.lua

@@ -16,13 +16,20 @@ postpro["E"].silo_mode = "zonal"
 postpro["H"].silo_mode = "none"
 postpro["J"].silo_mode = "none"
 
-local epsr = 1
-local mur = 1
+local epsr = 1.0
+local mur = 1.0
 
-materials[1] = {}
-materials[1].epsilon = epsr
-materials[1].mu = mur
-materials[1].sigma = 0
+materials.epsilon = function(tag, x, y, z)
+    return epsr;
+end
+
+materials.mu = function(tag, x, y, z)
+    return mur;
+end
+
+materials.sigma = function(tag, x, y, z)
+    return 0.0;
+end
 
 function electric_initial_condition(x, y, z)
     local Ex = 0
@@ -57,8 +64,8 @@ local resonance_f = omega/(2*math.pi)
 local resonance_MHz = resonance_f/1e6
 local cycle_timesteps = 1/(resonance_f*sim.dt)
 -- Compute impedance
-local eps = materials[1].epsilon * const.eps0
-local mu = materials[1].mu * const.mu0
+local eps = epsr * const.eps0
+local mu = mur * const.mu0
 local k_sq = (omega*mu)*(omega*eps)
 local kc_sq = u*u + v*v
 local beta = math.sqrt(k_sq - kc_sq)

share/examples/resonator/params_test_maxwell_resonator.lua

@@ -6,7 +6,7 @@ sim.name = "test_maxwell_resonator"             -- simulation name
 sim.dt = 1e-12                                  -- timestep size
 sim.timesteps = 501                           -- num of iterations
 sim.gmsh_model = "resonator.geo"                -- gmsh model filename
-sim.use_gpu = 1                                 -- 0: cpu, 1: gpu
+sim.use_gpu = 0                                 -- 0: cpu, 1: gpu
 sim.approx_order = 3                            -- approximation order
 sim.time_integrator = "leapfrog"
 postpro.silo_output_rate = 100
@@ -16,13 +16,20 @@ postpro["E"].silo_mode = "nodal"
 postpro["H"].silo_mode = "nodal"
 postpro["J"].silo_mode = "none"
 
-local epsr = 1
-local mur = 1
+local epsr = 1.0
+local mur = 1.0
 
-materials[1] = {}
-materials[1].epsilon = epsr
-materials[1].mu = mur
-materials[1].sigma = 0
+materials.epsilon = function(tag, x, y, z)
+    return epsr;
+end
+
+materials.mu = function(tag, x, y, z)
+    return mur;
+end
+
+materials.sigma = function(tag, x, y, z)
+    return 0.0;
+end
 
 function electric_initial_condition(x, y, z)
     local Ex = 0
@@ -57,8 +64,8 @@ local resonance_f = omega/(2*math.pi)
 local resonance_MHz = resonance_f/1e6
 local cycle_timesteps = 1/(resonance_f*sim.dt)
 -- Compute impedance
-local eps = materials[1].epsilon * const.eps0
-local mu = materials[1].mu * const.mu0
+local eps = epsr * const.eps0
+local mu = mur * const.mu0
 local k_sq = (omega*mu)*(omega*eps)
 local kc_sq = u*u + v*v
 local beta = math.sqrt(k_sq - kc_sq)