Peano 4
Loading...
Searching...
No Matches
old__SolverInformationADERDG.h
Go to the documentation of this file.
1#ifndef EXASEIS_SOLVERADERDG_INFORMATION_HEADER
2#define EXASEIS_SOLVERADERDG_INFORMATION_HEADER
3
4#include <algorithm>
5
6#include "SolverInformation.h"
7#include "kernels/GaussLegendreBasis.h"
8#include "exahype/solvers/ADERDGSolver.h"
9#include "../../PML_LOH1_Marc_test/AbstractElasticWaveSolver.h"
10
11template <int order>
13public:
14 // SolverInformationADERDG(exahype::solvers::ADERDGSolver* a_solver): SolverInformation(a_solver){
15 SolverInformationADERDG(Elastic::AbstractElasticWaveSolver* a_solver):
16 SolverInformation(a_solver) {
17
18 std::copy_n(a_solver->nodes[order], order + 1, nodes);
19 for (int i = 0; i < order + 1; i++) {
20 std::copy_n(a_solver->dudx[order][i], order + 1, dudx[i]);
21 }
22
23#ifdef Asserts
24 for (int i = 0; i < order + 1; i++) {
25 assertion2(std::isfinite(nodes[i]), nodes[i], i);
26 }
27
28 for (int i = 0; i < order + 1; i++) {
29 for (int j = 0; j < order + 1; j++) {
30 assertion3(std::isfinite(dudx[i][j]), dudx[i][j], i, j);
31 }
32 }
33#endif
34 };
35
36 double getNodes(int i) { return nodes[i]; }
37
38 double getDuDx(int i, int j) { return dudx[i][j]; }
39
40 bool isDG() override { return true; }
41
42 double nodes[order + 1];
43 double dudx[order + 1][order + 1];
44};
45#endif
#define assertion2(expr, param0, param1)
#define assertion3(expr, param0, param1, param2)
And from this we can write down f$ nabla phi_i nabla phi_i dx but since we are constructing matrix let s investigate the f$ j
SolverInformationADERDG(Elastic::AbstractElasticWaveSolver *a_solver)