Commit 2ccb9383 authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

back to default problem

parent e2d2ad37
......@@ -263,22 +263,22 @@ private:
public:
GaussHat2D(Meshes &meshes) :
mu(VectorField(0.75, 0.75, 0.0)),
mu(VectorField(0.5, 0.5, 0.0)),
sigma(Tensor()),
vectorField(VectorField(1 / sqrt(2), 1 / sqrt(2), 0)),
IStochasticTransportProblem(GeneratorNames{"DummyGenerator"},
meshes) {
sigma[0][0] = 0.1;
sigma[0][0] = 0.01;
sigma[0][1] = 0.0;
sigma[1][0] = 0.0;
sigma[1][1] = 0.1;
sigma[1][1] = 0.01;
}
Scalar Solution(double t, const Point &x) const override {
double factor = 1 / sqrt(pow(2 * Pi, 2) * sigma.det());
VectorField diff = (x - mu);
VectorField temp = Invert(sigma) * diff;
Scalar exponent = - 0.5 * diff * temp;
Scalar exponent = -0.5 * diff * temp;
return factor * exp(exponent);
}
......
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