TestPDESolver.cpp

#include "assemble/LagrangeEllipticAssemble.hpp"
#include "assemble/MixedEllipticAssemble.hpp"
#include "assemble/HybridEllipticAssemble.hpp"
#include "assemble/DGEllipticAssemble.hpp"
#include "assemble/DGTransportAssemble.hpp"
#include "assemble/DGReactionAssemble.hpp"
#include "assemble/PGReactionAssemble.hpp"
#include "montecarlo/PDESolver.hpp"
#include "TestEnvironment.hpp"


class TestPDESolver : public Test {
protected:
    PDESolver *pdeSolver;

    Meshes *meshes;

    MatrixGraphs *mGraphs;

    SampleID dummyID;

    TestPDESolver() {
        meshes = createTestMeshes("Square_DirichletBC", 0, 0);
        dummyID.level = Level(1 ,0);
        dummyID.coarse = false;
        dummyID.number = -1;
    }

    void TearDown() override {
        delete pdeSolver;
        delete meshes;
    }
};

class TestEllipticPDESolver : public TestPDESolver {
protected:
    TestEllipticPDESolver() : TestPDESolver() {
        pdeSolver = new EllipticPDESolver(
            new LagrangeEllipticAssemble(
                new LagrangeDiscretization(*meshes, 1),
                CreateStochasticEllipticProblem("StochasticLaplace2D", *meshes)
            )
        );

        mGraphs = pdeSolver->CreateSolutionMatrixGraph(*meshes);
    }
};

class TestTransportPDESolver : public TestPDESolver {
    void SetUp() override {
        // Todo
    }

    void TearDown() override {
        // Todo
    }
};

TEST_F(TestEllipticPDESolver, TestLagrangeDiscretization) {
    SampleSolution solution(*mGraphs, dummyID.level.mGfine);
    pdeSolver->Run(solution);
}

int main(int argc, char **argv) {
    MppTest mppTest = MppTestBuilder(argc, argv).WithPPM().WithScreenLogging();
    return mppTest.RUN_ALL_MPP_TESTS();
}