ThunderEgg  1.0.0
MatrixHelper.h
Go to the documentation of this file.
1 /***************************************************************************
2  * ThunderEgg, a library for solvers on adaptively refined block-structured
3  * Cartesian grids.
4  *
5  * Copyright (c) 2017-2021 Scott Aiton
6  *
7  * This program is free software: you can redistribute it and/or modify
8  * it under the terms of the GNU General Public License as published by
9  * the Free Software Foundation, either version 3 of the License, or
10  * (at your option) any later version.
11  *
12  * This program is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15  * GNU General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License
18  * along with this program. If not, see <https://www.gnu.org/licenses/>.
19  ***************************************************************************/
20 
21 #ifndef THUNDEREGG_POISSON_MATRIXHELPER_H
22 #define THUNDEREGG_POISSON_MATRIXHELPER_H
23 
29 #include <ThunderEgg/Domain.h>
30 #include <petscmat.h>
31 
32 namespace ThunderEgg::Poisson {
39 {
40 private:
44  Domain<3> domain;
48  std::bitset<6> neumann;
49 
50 public:
57  explicit MatrixHelper(const Domain<3>& domain, std::bitset<6> neumann);
58 
64  Mat formCRSMatrix();
65 };
66 } // namespace ThunderEgg::Poisson
67 #endif
ThunderEgg::Poisson::MatrixHelper
Create a matrix for the 3D second-order Laplacian operator.
Definition: MatrixHelper.h:38
ThunderEgg::Domain< 3 >
ThunderEgg::Poisson
Classes specific to the Poisson equation.
Definition: DFTPatchSolver.h:40
ThunderEgg::Poisson::MatrixHelper::formCRSMatrix
Mat formCRSMatrix()
Form the matrix for the domain.
ThunderEgg::Poisson::MatrixHelper::MatrixHelper
MatrixHelper(const Domain< 3 > &domain, std::bitset< 6 > neumann)
Create a MatrixHelper for a given 3D domain.
Domain.h
Domain class.