Speed up reg tests and add unit tests for moments/vlasov/pkpm

moments/creg/rt_10m_burch.c

@@ -126,8 +126,8 @@ create_ctx(void)
     + n1 * (Ti1 + Te1) - n2 * Ti2) / n2; // Magnetosheath electron temperature (so that the system is in force balance).
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Nx = 16; // Cell count (x-direction).
+  int Ny = 16; // Cell count (y-direction).
   double Lx = 40.96 * di; // Domain size (x-direction).
   double Ly = 20.48 * di; // Domain size (y-direction).
   double k0_elc = 1.0; // Closure parameter for electrons.

moments/creg/rt_10m_burch_grad_closure.c

@@ -126,8 +126,8 @@ create_ctx(void)
     + n1 * (Ti1 + Te1) - n2 * Ti2) / n2; // Magnetosheath electron temperature (so that the system is in force balance).
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Nx = 64; // Cell count (x-direction).
+  int Ny = 64; // Cell count (y-direction).
   double Lx = 40.96 * di; // Domain size (x-direction).
   double Ly = 20.48 * di; // Domain size (y-direction).
   double k0_elc = 10.0; // Closure parameter for electrons.

moments/creg/rt_10m_checkerboard.c

@@ -83,7 +83,7 @@ create_ctx(void)
   double cfl_frac = 1.0; // CFL coefficient.
 
   double t_end = 1000.0; // Final simulation time.
-  int num_frames = 100; // Number of output frames.
+  int num_frames = 1; // Number of output frames.
   int field_energy_calcs = INT_MAX; // Number of times to calculate field energy.
   int integrated_mom_calcs = INT_MAX; // Number of times to calculate integrated moments.
   double dt_failure_tol = 1.0e-4; // Minimum allowable fraction of initial time-step.

moments/creg/rt_10m_expanding_axi_sodshock.c

@@ -80,7 +80,7 @@ create_ctx(void)
   double pr = 1.0; // Right/outer fluid pressure.
 
   // Simulation parameters.
-  int Nr = 128; // Cell count (radial direction).
+  int Nr = 16; // Cell count (radial direction).
   int Ntheta = 128 * 6; // Cell count (angular direction).
   double Lr = 1.0; // Domain size (radial direction).
   double Ltheta = 2.0 * pi; // Domain size (angular direction).

moments/creg/rt_10m_gem.c

@@ -108,7 +108,7 @@ create_ctx(void)
   double omega_ci = fabs(charge_ion * B0 / mass_ion); // Ion cyclotron frequency.
 
   // Simulation parameters.
-  int Nx = 128; // Cell count (x-direction).
+  int Nx = 64; // Cell count (x-direction).
   int Ny = 64; // Cell count (y-direction).
   double Lx = 25.6 * di; // Domain size (x-direction).
   double Ly = 12.8 * di; // Domain size (y-direction).

moments/creg/rt_10m_gem_grad_closure.c

@@ -108,7 +108,7 @@ create_ctx(void)
   double omega_ci = fabs(charge_ion * B0 / mass_ion); // Ion cyclotron frequency.
 
   // Simulation parameters.
-  int Nx = 128; // Cell count (x-drection).
+  int Nx = 64; // Cell count (x-drection).
   int Ny = 64; // Cell count (y-direction).
   double Lx = 25.6 * di; // Domain size (x-direction).
   double Ly = 12.8 * di; // Domain size (y-direction).

moments/creg/rt_10m_lhdi.c

@@ -142,7 +142,7 @@ create_ctx(void)
 
   // Simulation parameters.
   int Nx = 64; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Ny = 64; // Cell count (y-direction).
   double Lx = 6.4 * l; // Domain size (x-direction).
   double Ly = 12.8 * l; // Domain size (y-direction).
   double k0_elc = 1.0; // Closure parameter for electrons.

moments/creg/rt_10m_lhdi_grad_closure.c

@@ -141,8 +141,8 @@ create_ctx(void)
   double iy = 1.0; // Current (y-direction).
 
   // Simulation parameters.
-  int Nx = 64; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Nx = 32; // Cell count (x-direction).
+  int Ny = 64; // Cell count (y-direction).
   double Lx = 6.4 * l; // Domain size (x-direction).
   double Ly = 12.8 * l; // Domain size (y-direction).
   double k0_elc = 1.0; // Closure parameter for electrons.

moments/creg/rt_10m_ot_grad_closure.c

@@ -102,8 +102,8 @@ create_ctx(void)
   double delta_u0 = 0.2 * vAi; // Reference fluid velocity perturbation.
 
   // Simulation parameters.
-  int Nx = 128; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Nx = 64; // Cell count (x-direction).
+  int Ny = 64; // Cell count (y-direction).
   double Lx = 20.48 * d_i; // Domain size (x-direction).
   double Ly = 20.48 * d_i; // Domain size (y-direction).
   double k0 = 5.0; // Closure parameter.

moments/creg/rt_10m_ot_nn_closure.c

@@ -102,8 +102,8 @@ create_ctx(void)
   double delta_u0 = 0.2 * vAi; // Reference fluid velocity perturbation.
 
   // Simulation parameters.
-  int Nx = 128; // Cell count (x-direction).
-  int Ny = 128; // Cell count (y-direction).
+  int Nx = 8; // Cell count (x-direction).
+  int Ny = 16; // Cell count (y-direction).
   double Lx = 20.48 * d_i; // Domain size (x-direction).
   double Ly = 20.48 * d_i; // Domain size (y-direction).
   double k0 = 5000.0; // Closure parameter.

moments/creg/rt_10m_riem_nn_closure_p1.c

@@ -98,7 +98,7 @@ create_ctx(void)
   double rhor_elc = rhor_ion * mass_elc / mass_ion; // Right electron mass density.
 
   // Simulation parameters.
-  int Nx = 1024; // Cell count (x-direction).
+  int Nx = 64; // Cell count (x-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double k0 = 500.0; // Closure parameter.
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_10m_riem_nn_closure_p2.c

@@ -98,7 +98,7 @@ create_ctx(void)
   double rhor_elc = rhor_ion * mass_elc / mass_ion; // Right electron mass density.
 
   // Simulation parameters.
-  int Nx = 1024; // Cell count (x-direction).
+  int Nx = 64; // Cell count (x-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double k0 = 500.0; // Closure parameter.
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_5m_burch.c

@@ -119,7 +119,7 @@ create_ctx(void)
     + n1 * (Ti1 + Te1) - n2 * Ti2) / n2; // Magnetosheath electron temperature (so that the system is in force balance).
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 64; // Cell count (x-direction).
   int Ny = 128; // Cell count (y-direction).
   double Lx = 40.96 * di; // Domain size (x-direction).
   double Ly = 20.48 * di; // Domain size (y-direction).

moments/creg/rt_5m_expanding_axi_sodshock.c

@@ -72,7 +72,7 @@ create_ctx(void)
   double pr = 1.0; // Right/outer fluid pressure.
 
   // Simulation parameters.
-  int Nr = 128; // Cell count (radial direction).
+  int Nr = 32; // Cell count (radial direction).
   int Ntheta = 128 * 6; // Cell count (angular direction).
   double Lr = 1.0; // Domain size (radial direction).
   double Ltheta = 2.0 * pi; // Domain size (angular direction).

moments/creg/rt_5m_rt.c

@@ -98,7 +98,7 @@ create_ctx(void)
   double grav = g_hat * omega_ci * vAi; // Gravitational acceleration.
 
   // Simulation parameters.
-  int Nx = 64; // Cell count (x-direction).
+  int Nx = 32; // Cell count (x-direction).
   int Ny = 64; // Cell count (y-direction).
   double Lx = 3.0; // Domain size (x-direction).
   double Ly = 3.75; // Domain size (y-direction).

moments/creg/rt_euler_embedded_surface.c

@@ -59,8 +59,8 @@ create_ctx(void)
   double p1 = 1.0;
 
   // Simulation parameters.
-  int Nx = 300; // Cell count (x-direction).
-  int Ny = 300; // Cell count (y-direction).
+  int Nx = 150; // Cell count (x-direction).
+  int Ny = 150; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.

moments/creg/rt_euler_embedded_surface_mapc2p.c

@@ -59,8 +59,8 @@ create_ctx(void)
   double p1 = 1.0;
 
   // Simulation parameters.
-  int Nx = 300; // Cell count (x-direction).
-  int Ny = 300; // Cell count (y-direction).
+  int Nx = 150; // Cell count (x-direction).
+  int Ny = 150; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.

moments/creg/rt_euler_mixture_shock_bubble.c

@@ -85,7 +85,7 @@ create_ctx(void)
   double p_bub = 1.0 / gas_gamma1; // Bubble fluid pressure.
 
   // Simulation parameters.
-  int Nx = 325; // Cell count (x-direction).
+  int Nx = 89; // Cell count (x-direction).
   int Ny = 89; // Cell count (y-direction).
   double Lx = 0.325; // Domain size (x-direction).
   double Ly = 0.089; // Domain size (y-direction).

moments/creg/rt_euler_rgfm_shock_bubble.c

@@ -86,7 +86,7 @@ create_ctx(void)
   double p_bub = 1.0 / gas_gamma1; // Bubble fluid pressure.
 
   // Simulation parameters.
-  int Nx = 325; // Cell count (x-direction).
+  int Nx = 89; // Cell count (x-direction).
   int Ny = 89; // Cell count (y-direction).
   double Lx = 0.325; // Domain size (x-direction).
   double Ly = 0.089; // Domain size (y-direction).

moments/creg/rt_euler_riem_2d_hll.c

@@ -92,7 +92,7 @@ create_ctx(void)
   double p_lr = 0.3; // Lower-right fluid pressure.
 
   // Simulation parameters.
-  int Nx = 200; // Cell count (x-direction).
+  int Nx = 100; // Cell count (x-direction).
   int Ny = 200; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).

moments/creg/rt_euler_riem_2d_hllc.c

@@ -92,7 +92,7 @@ create_ctx(void)
   double p_lr = 0.3; // Lower-right fluid pressure.
 
   // Simulation parameters.
-  int Nx = 200; // Cell count (x-direction).
+  int Nx = 100; // Cell count (x-direction).
   int Ny = 200; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).

moments/creg/rt_euler_riem_2d_lax.c

@@ -92,7 +92,7 @@ create_ctx(void)
   double p_lr = 0.3; // Lower-right fluid pressure.
 
   // Simulation parameters.
-  int Nx = 200; // Cell count (x-direction).
+  int Nx = 100; // Cell count (x-direction).
   int Ny = 200; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).

moments/creg/rt_euler_riem_2d_roe.c

@@ -92,7 +92,7 @@ create_ctx(void)
   double p_lr = 0.3; // Lower-right fluid pressure.
 
   // Simulation parameters.
-  int Nx = 200; // Cell count (x-direction).
+  int Nx = 100; // Cell count (x-direction).
   int Ny = 200; // Cell count (y-direction).
   double Lx = 1.0; // Domain size (x-direction).
   double Ly = 1.0; // Domain size (y-direction).

moments/creg/rt_euler_riem_3d.c

@@ -67,7 +67,7 @@ create_ctx(void)
   double pr = 1.0; // Right/outer fluid pressure.
 
   // Simulation parameters.
-  int Nx = 37; // Cell count (x-direction).
+  int Nx = 18; // Cell count (x-direction).
   int Ny = 37; // Cell count (y-direction).
   int Nz = 25; // Cell count (z-direction).
   double Lx = 1.5; // Domain size (x-direction).

moments/creg/rt_gr_bhl_spinning.c

@@ -96,7 +96,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_bhl_spinning_neutronstar.c

@@ -120,8 +120,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_neutronstar_new(false, mass, spin, mass_quadrupole, spin_octupole, mass_hexadecapole, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_bhl_spinning_tetrad.c

@@ -96,7 +96,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_bhl_static.c

@@ -96,7 +96,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_bhl_static_neutronstar.c

@@ -120,8 +120,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_neutronstar_new(false, mass, spin, mass_quadrupole, spin_octupole, mass_hexadecapole, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_bhl_static_tetrad.c

@@ -96,7 +96,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_blackhole_spinning.c

@@ -107,7 +107,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_blackhole_static.c

@@ -107,7 +107,7 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
+  int Nx = 128; // Cell count (x-direction).
   int Ny = 256; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).

moments/creg/rt_gr_bz_monopole_fast.c

@@ -92,8 +92,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 10.0; // Domain size (x-direction).
   double Ly = 10.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_bz_monopole_fast_tetrad.c

@@ -92,8 +92,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 10.0; // Domain size (x-direction).
   double Ly = 10.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_bz_monopole_slow.c

@@ -92,8 +92,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 10.0; // Domain size (x-direction).
   double Ly = 10.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_bz_monopole_slow_tetrad.c

@@ -92,8 +92,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 10.0; // Domain size (x-direction).
   double Ly = 10.0; // Domain size (y-direction).
   double cfl_frac = 0.95; // CFL coefficient.

moments/creg/rt_gr_mhd_bhl_spinning.c

@@ -99,8 +99,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.

moments/creg/rt_gr_mhd_bhl_spinning_mhd.c

@@ -99,8 +99,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.

moments/creg/rt_gr_mhd_bhl_spinning_tetrad.c

@@ -99,8 +99,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.

moments/creg/rt_gr_mhd_bhl_static.c

@@ -99,8 +99,8 @@ create_ctx(void)
   struct gkyl_gr_spacetime *spacetime = gkyl_gr_blackhole_new(false, mass, spin, pos_x, pos_y, pos_z);
 
   // Simulation parameters.
-  int Nx = 256; // Cell count (x-direction).
-  int Ny = 256; // Cell count (y-direction).
+  int Nx = 128; // Cell count (x-direction).
+  int Ny = 128; // Cell count (y-direction).
   double Lx = 5.0; // Domain size (x-direction).
   double Ly = 5.0; // Domain size (y-direction).
   double cfl_frac = 0.9; // CFL coefficient.