save progress on SF selection pairs in compressible FHD

Author: isriva

src_compressible/compressible_functions.cpp

@@ -13,6 +13,12 @@ AMREX_GPU_MANAGED bool compressible::do_reservoir = false;
 AMREX_GPU_MANAGED amrex::Real compressible::zeta_ratio = -1.0;
 AMREX_GPU_MANAGED amrex::Real compressible::p_init = -1.0;
 AMREX_GPU_MANAGED amrex::GpuArray<amrex::Real, MAX_SPECIES> compressible::Xk_init;
+AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 2 + MAX_SPECIES> cons_SF_pairA;
+AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 2 + MAX_SPECIES> cons_SF_pairB;
+AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 3 + 2*MAX_SPECIES> prim_SF_pairA;
+AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 3 + 2*MAX_SPECIES> prim_SF_pairB;
+AMREX_GPU_MANAGED int do_SF_pair_cons = -1;
+AMREX_GPU_MANAGED int do_SF_pair_prim = -1;
 
 
 void InitializeCompressibleNamespace()
@@ -108,6 +114,16 @@ void InitializeCompressibleNamespace()
         amrex::Print() << "\n";
     }
 
+    // Get structure factor pairs
+    for (int i=0; i< AMREX_SPACEDIM + 2 + MAX_SPECIES; ++i) {
+      cons_SF_pairA[i] = -1;
+      cons_SF_pairB[i] = -1;
+    }
+    for (int i=0; i< AMREX_SPACEDIM + 3 + 2*MAX_SPECIES; ++i) {
+      prim_SF_pairA[i] = -1;
+      prim_SF_pairB[i] = -1;
+    }
+
     return;
 }
 

src_compressible/compressible_namespace.H

@@ -12,6 +12,11 @@ namespace compressible {
     extern AMREX_GPU_MANAGED amrex::Real zeta_ratio;
     extern AMREX_GPU_MANAGED amrex::Real p_init;
     extern AMREX_GPU_MANAGED amrex::GpuArray<amrex::Real, MAX_SPECIES> Xk_init;
-
+    extern AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 2 + MAX_SPECIES> cons_SF_pairA;
+    extern AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 2 + MAX_SPECIES> cons_SF_pairB;
+    extern AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 3 + 2*MAX_SPECIES> prim_SF_pairA;
+    extern AMREX_GPU_MANAGED amrex::GpuArray<int, AMREX_SPACEDIM + 3 + 2*MAX_SPECIES> prim_SF_pairB;
+    extern AMREX_GPU_MANAGED int do_SF_pair_cons;
+    extern AMREX_GPU_MANAGED int do_SF_pair_prim;
 }