[Fix] The three uncapitalized files

Author: lukbrb

init/hydro/Blast.h

@@ -0,0 +1,43 @@
+#pragma once
+
+#include "../../SimInfo.h"
+#include "../InitData.h"
+
+namespace fv2d
+{
+
+/**
+ * @brief Sedov blast initial conditions
+ */
+KOKKOS_INLINE_FUNCTION
+void initBlast(Array Q, int i, int j, const DeviceParams &params, const InitData &init_data)
+{
+  real_t xmid = 0.5 * (params.xmin + params.xmax);
+  real_t ymid = 0.5 * (params.ymin + params.ymax);
+
+  Pos pos  = getPos(params, i, j);
+  real_t x = pos[IX];
+  real_t y = pos[IY];
+
+  real_t xr = xmid - x;
+  real_t yr = ymid - y;
+  real_t r  = sqrt(xr * xr + yr * yr);
+
+  if (r < 0.2)
+  {
+    Q(j, i, IR) = 1.0;
+    Q(j, i, IU) = 0.0;
+    Q(j, i, IV) = 0.0;
+    Q(j, i, IP) = 10.0;
+  }
+  else
+  {
+    Q(j, i, IR) = 1.2;
+    Q(j, i, IU) = 0.0;
+    Q(j, i, IV) = 0.0;
+    Q(j, i, IP) = 0.1;
+  }
+}
+
+REGISTER_PROBLEM(initBlast, "blast")
+} // namespace fv2d

init/hydro/Diffusion.h

@@ -0,0 +1,37 @@
+#pragma once
+
+#include "../../SimInfo.h"
+#include "../InitData.h"
+
+namespace fv2d
+{
+
+/**
+ * @brief Simple diffusion test with a structure being advected on the grid
+ */
+KOKKOS_INLINE_FUNCTION
+void initDiffusion(Array Q, int i, int j, const DeviceParams &params, const InitData &init_data)
+{
+  real_t xmid = 0.5 * (params.xmin + params.xmax);
+  real_t ymid = 0.5 * (params.ymin + params.ymax);
+
+  Pos pos = getPos(params, i, j);
+
+  real_t x0 = (pos[IX] - xmid);
+  real_t y0 = (pos[IY] - ymid);
+
+  real_t r = sqrt(x0 * x0 + y0 * y0);
+
+  if (r < 0.2)
+    Q(j, i, IR) = 1.0;
+  else
+    Q(j, i, IR) = 0.1;
+
+  Q(j, i, IP) = 1.0;
+  Q(j, i, IU) = 1.0;
+  Q(j, i, IV) = 1.0;
+}
+
+REGISTER_PROBLEM(initDiffusion, "diffusion")
+
+} // namespace fv2d

init/hydro/Sod.h

@@ -0,0 +1,53 @@
+#pragma once
+
+#include "../../SimInfo.h"
+#include "../InitData.h"
+
+namespace fv2d
+{
+
+/**
+ * @brief Sod Shock tube aligned along the X axis
+ */
+KOKKOS_INLINE_FUNCTION
+void initSodX(Array Q, int i, int j, const DeviceParams &params, const InitData &init_data)
+{
+  if (getPos(params, i, j)[IX] <= 0.5)
+  {
+    Q(j, i, IR) = 1.0;
+    Q(j, i, IP) = 1.0;
+    Q(j, i, IU) = 0.0;
+  }
+  else
+  {
+    Q(j, i, IR) = 0.125;
+    Q(j, i, IP) = 0.1;
+    Q(j, i, IU) = 0.0;
+  }
+}
+
+REGISTER_PROBLEM(initSodX, "sod_x")
+
+/**
+ * @brief Sod Shock tube aligned along the Y axis
+ */
+KOKKOS_INLINE_FUNCTION
+void initSodY(Array Q, int i, int j, const DeviceParams &params, const InitData &init_data)
+{
+  if (getPos(params, i, j)[IY] <= 0.5)
+  {
+    Q(j, i, IR) = 1.0;
+    Q(j, i, IP) = 1.0;
+    Q(j, i, IU) = 0.0;
+  }
+  else
+  {
+    Q(j, i, IR) = 0.125;
+    Q(j, i, IP) = 0.1;
+    Q(j, i, IU) = 0.0;
+  }
+}
+
+REGISTER_PROBLEM(initSodY, "sod_y")
+
+} // namespace fv2d