Commit 788a283a authored by chinsp's avatar chinsp
Browse files

fixed test for Gauss Hermite rule

parent 7e1f499c
Pipeline #160039 failed with stages
in 11 minutes and 28 seconds
...@@ -17,6 +17,8 @@ protected: ...@@ -17,6 +17,8 @@ protected:
int outputs; int outputs;
//double alpha;
double weight = 0.0; double weight = 0.0;
RVector sample{}; RVector sample{};
......
...@@ -8,7 +8,7 @@ ...@@ -8,7 +8,7 @@
#include "SparseGridGenerator.hpp" #include "SparseGridGenerator.hpp"
#include "HybridFluxGenerator.hpp" #include "HybridFluxGenerator.hpp"
#include "CirculantEmbedding.hpp" #include "CirculantEmbedding.hpp"
#include "math.h"
class IStochasticProblem { class IStochasticProblem {
protected: protected:
...@@ -286,12 +286,13 @@ public: ...@@ -286,12 +286,13 @@ public:
class SparseGrid2DhermiteProblem : public SparseGridGeneratorProblem { class SparseGrid2DhermiteProblem : public SparseGridGeneratorProblem {
public: public:
explicit SparseGrid2DhermiteProblem(const Meshes &meshes) : explicit SparseGrid2DhermiteProblem(const Meshes &meshes) :
SparseGridGeneratorProblem(meshes, SparseGridGenerator(meshes, 1, 0, 5, TasGrid::type_qptotal, TasGrid::rule_gausshermite)) {} SparseGridGeneratorProblem(meshes, SparseGridGenerator(meshes, 1, 0, 8,
TasGrid::type_qptotal, TasGrid::rule_gausshermite)) {}
double FunctionEvaluation() override { double FunctionEvaluation() override {
RVector sample = this->generator.EvalSample(); RVector sample = this->generator.EvalSample();
//double y1 = sample[0]; //double y1 = sample[0];
return 1; return 1 / sqrt(M_PI);
} }
string Name() const override { string Name() const override {
......
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