Commit 329d79f0 authored by Houssein Zaghdane's avatar Houssein Zaghdane

solution in fv

parent 1f660718
......@@ -23,7 +23,7 @@ void SWE::MySWESolver_FV::adjustSolution(const double* const x,const double t,co
if (tarch::la::equals(t,0.0)) {
//std::cout << "Warning: using FV in initial step" << std::endl;
/*MySWESolver_FV::Variables vars(Q);
MySWESolver_FV::Variables vars(Q);
const double d = 0.15;
const double H = 0.3 * d;
......@@ -47,8 +47,17 @@ void SWE::MySWESolver_FV::adjustSolution(const double* const x,const double t,co
vars.b() = 0;
}
vars.hu() = -eta * std::sqrt(grav_FV/d) * vars.h();
vars.hv() = 0.0;*/
vars.hv() = 0.0;
std::vector<std::vector<double>> probe_point = { {1.5,0.5}, {2,0.5}, {2.5,0.5} //, {3, 0.5}, {4, 0.5}
};
for (int i = 0; i< probe_point.size(); i++){
if (abs(x[0] - probe_point[i][0]) < 0.00001 && abs(x[1] - probe_point[i][1]) < 0.00001){
muq::solution[i] = vars.h() + vars.b();
}
}
}
else{
......
......@@ -20,6 +20,12 @@
*/
#include "tarch/logging/Log.h"
namespace muq{
extern std::vector<double> solution;
extern std::vector<double> param;
}
namespace SWE{
class MySWESolver_FV;
......
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