Peano 4
Loading...
Searching...
No Matches
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// #include "exahype/solvers/ADERDGSolver.h"
6#include "SolverInformation.h"
8// #include "kernels/GaussLegendreBasis.h"
9// #include "kernels/GaussLobattoBasis.h"
11
12template <int order>
14public:
16 SolverInformation(a_solver) {
17
18 std::copy_n(kernels::lobatto::nodes[order], order + 1, nodes);
19 for (int i = 0; i < order + 1; i++) {
20 std::copy_n(kernels::lobatto::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(exahype2::solvers::aderdg::ADERDGSolver *a_solver)