24.09., 9:00 - 11:00: Due to updates GitLab will be unavailable for some minutes between 09:00 and 11:00.

Commit 9ffeff65 authored by Jean-Matthieu Gallard's avatar Jean-Matthieu Gallard

Mixed Euler - update test case to use predictor recompute + vect pde - p2

parent 954bc81d
...@@ -30,8 +30,6 @@ private: ...@@ -30,8 +30,6 @@ private:
* See also chapter 7.13.2 in "I do like CFD, VOL.1" by Katate Masatsuka. * See also chapter 7.13.2 in "I do like CFD, VOL.1" by Katate Masatsuka.
*/ */
static void entropyWave(const double* const x,double t, double* const Q); static void entropyWave(const double* const x,double t, double* const Q);
// traveling sin waves in x and y directions
static void sinWave(const double* const x,double t, double* const Q);
/** /**
* Log device * Log device
...@@ -52,7 +50,7 @@ public: ...@@ -52,7 +50,7 @@ public:
* calls ::entropyWave * calls ::entropyWave
* ErrorWriter, ErrorPlotter write errors of numerical solution. * ErrorWriter, ErrorPlotter write errors of numerical solution.
*/ */
static void referenceSolution(const double* const x, const double t, double* const Q); inline static void referenceSolution(const double* const x, const double t, double* const Q) {entropyWave(x,t,Q);}
/** /**
* Adjust the conserved variables and parameters (together: Q) at a given time t at the (quadrature) point x. * Adjust the conserved variables and parameters (together: Q) at a given time t at the (quadrature) point x.
...@@ -76,8 +74,12 @@ public: ...@@ -76,8 +74,12 @@ public:
* \param[inout] F the fluxes at that point as C array (already allocated). * \param[inout] F the fluxes at that point as C array (already allocated).
*/ */
virtual void flux(const double* const Q,double** const F); virtual void flux(const double* const Q,double** const F);
inline void flux2(const double* const Q,const double* const P,double** const F){flux(Q,F);} //add P array
void flux_vect(const double* const Q,const double* const P,double** const F); //add P array
virtual void nonConservativeProduct(const double* const Q, const double* const gradQ, double* const BgradQ); virtual void nonConservativeProduct(const double* const Q, const double* const gradQ, double* const BgradQ);
void nonConservativeProduct2(const double* const Q, const double* const P, const double* const * const gradQ, double* const BgradQ); // add P array and gradQ like F
void nonConservativeProduct_vect(const double* const Q, const double* const P, const double* const * const gradQ, double* const BgradQ); // add P array and gradQ like F
/** /**
* Compute the eigenvalues of the flux tensor per coordinate direction \p d. * Compute the eigenvalues of the flux tensor per coordinate direction \p d.
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment