tests work

lanl · Dec 16, 2024 · 3d18cae · 3d18cae
1 parent 4b4a2fa
commit 3d18cae
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 6 deletions.
diff --git a/singularity-eos/eos/modifiers/zsplit_eos.hpp b/singularity-eos/eos/modifiers/zsplit_eos.hpp
@@ -238,22 +238,22 @@ class ZSplit : public EosBase<ZSplit<ztype, T>> {
 
  private:
   static constexpr const int NL = 1;
-  template <typename Indexer_t>
-  PORTABLE_FORCEINLINE_FUNCTION Real GetIonizationState_(const Indexer_t &&lambda) const {
+  template <typename Indexer_t = Real *>
+  PORTABLE_FORCEINLINE_FUNCTION Real GetIonizationState_(Indexer_t &&lambda) const {
     return std::max(0.0, lambda[t_.nlambda()]);
   }
   // TODO(JMM): Runtime?
-  template <typename Indexer_t>
-  PORTABLE_FORCEINLINE_FUNCTION Real GetScale_(const Indexer_t &&lambda) const {
+  template <typename Indexer_t = Real *>
+  PORTABLE_FORCEINLINE_FUNCTION Real GetScale_(Indexer_t &&lambda) const {
     Real Z = GetIonizationState_(lambda);
     if constexpr (ztype == ZSplitComponent::Electrons) {
       return robust::ratio(Z, Z + 1);
     } else {
       return robust::ratio(1.0, Z + 1);
     }
   }
-  template <typename Indexer_t>
-  PORTABLE_FORCEINLINE_FUNCTION Real GetInvScale_(const Indexer_t &&lambda) const {
+  template <typename Indexer_t = Real *>
+  PORTABLE_FORCEINLINE_FUNCTION Real GetInvScale_(Indexer_t &&lambda) const {
     Real Z = GetIonizationState_(lambda);
     if constexpr (ztype == ZSplitComponent::Electrons) {
       return robust::ratio(Z + 1, Z);

diff --git a/test/CMakeLists.txt b/test/CMakeLists.txt
@@ -24,6 +24,7 @@ add_executable(
   test_eos_vinet.cpp
   test_eos_mgusup.cpp
   test_eos_powermg.cpp
+  test_eos_zsplit.cpp
   )
 
 add_executable(