Add definition of Cjr_half

src/language/modules/SchemeModule.cpp

@@ -472,7 +472,7 @@ SchemeModule::SchemeModule()
                                                                      explicit_zone_list,
                                                                      dt}
                                   .solver()
-                                  .apply(dt, rho, u, E);
+                                  .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
                               }
 
                               ));
@@ -502,7 +502,7 @@ SchemeModule::SchemeModule()
                                              std::vector<std::shared_ptr<const IZoneDescriptor>>{},
                                              dt}
           .solver()
-          .apply(dt, rho, u, E);
+          .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
       }
 
       ));
@@ -540,7 +540,7 @@ SchemeModule::SchemeModule()
                                                                      explicit_zone_list,
                                                                      dt}
                                   .solver()
-                                  .apply(dt, rho, u, E);
+                                  .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
                               }
 
                               ));
@@ -570,7 +570,7 @@ SchemeModule::SchemeModule()
                                              std::vector<std::shared_ptr<const IZoneDescriptor>>{},
                                              dt}
           .solver()
-          .apply(dt, rho, u, E);
+          .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
       }
 
       ));
@@ -608,7 +608,7 @@ SchemeModule::SchemeModule()
                                                                    explicit_zone_list,
                                                                    dt}
                                 .solver()
-                                .apply(dt, rho, u, E);
+                                .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
                             }
 
                             ));
@@ -638,7 +638,7 @@ SchemeModule::SchemeModule()
                                              std::vector<std::shared_ptr<const IZoneDescriptor>>{},
                                              dt}
           .solver()
-          .apply(dt, rho, u, E);
+          .apply(dt, rho, u, E, c, p, pi, gamma, Cv, entropy);
       }
 
       ));

src/scheme/ImplicitAcousticSolver.hpp

@@ -50,7 +50,13 @@ class ImplicitAcousticSolverHandler
     apply(const double& dt,
           const std::shared_ptr<const IDiscreteFunction>& rho,
           const std::shared_ptr<const IDiscreteFunction>& u,
-          const std::shared_ptr<const IDiscreteFunction>& E) const = 0;
+          const std::shared_ptr<const IDiscreteFunction>& E,
+          const std::shared_ptr<const IDiscreteFunction>& c,
+          const std::shared_ptr<const IDiscreteFunction>& p,
+          const std::shared_ptr<const IDiscreteFunction>& pi,
+          const std::shared_ptr<const IDiscreteFunction>& gamma,
+          const std::shared_ptr<const IDiscreteFunction>& Cv,
+          const std::shared_ptr<const IDiscreteFunction>& entropy) = 0;
 
     IImplicitAcousticSolver()                          = default;
     IImplicitAcousticSolver(IImplicitAcousticSolver&&) = default;
@@ -65,8 +71,8 @@ class ImplicitAcousticSolverHandler
   std::unique_ptr<IImplicitAcousticSolver> m_implicit_acoustic_solver;
 
  public:
-  const IImplicitAcousticSolver&
-  solver() const
+  IImplicitAcousticSolver&
+  solver()
   {
     return *m_implicit_acoustic_solver;
   }