Poisson::FirstOrderSystem< Dim, BackgroundGeometry, DataType > Struct Template Reference

The Poisson equation formulated as a set of coupled first-order PDEs. More...

#include <FirstOrderSystem.hpp>

Public Types
using	primal_fields = implementation defined
using	primal_fluxes = implementation defined
using	background_fields = implementation defined
using	inv_metric_tag = implementation defined
using	fluxes_computer = Fluxes< Dim, BackgroundGeometry, DataType >
using	sources_computer = implementation defined
using	boundary_conditions_base = elliptic::BoundaryConditions::BoundaryCondition< Dim >
using	modify_boundary_data = void

Static Public Attributes
static constexpr size_t	volume_dim = Dim

Detailed Description

template<size_t Dim, Geometry BackgroundGeometry, typename DataType = DataVector>
struct Poisson::FirstOrderSystem< Dim, BackgroundGeometry, DataType >

The Poisson equation formulated as a set of coupled first-order PDEs.

Details

This system formulates the Poisson equation $-{\mathrm{\Delta }}_{\gamma }u(x)=f(x)$ on a background metric ${\gamma }_{ij}$ as the set of coupled first-order PDEs

$$\begin{array}{r}-{\partial }_i{v}^i-{\mathrm{\Gamma }}_{ij}^i{v}^j=f(x)\\ v^i={\gamma }^{ij}{\partial }_{j}u(x)\end{array}$$

where ${\mathrm{\Gamma }}_{jk}^i=\frac{1}{2}{\gamma }^{il}({\partial }_j{\gamma }_{kl}+{\partial }_k{\gamma }_{jl}-{\partial }_l{\gamma }_{jk})$ are the Christoffel symbols of the second kind of the background metric ${\gamma }_{ij}$. The background metric ${\gamma }_{ij}$ and the Christoffel symbols derived from it are assumed to be independent of the variables $u$, i.e. constant throughout an iterative elliptic solve.

The system can be formulated in terms of these fluxes and sources (see elliptic::protocols::FirstOrderSystem):

$$\begin{array}{rl}{F}^i& =v^i={\gamma }^{ij}{\partial }_{j}u\\ S& =-{\mathrm{\Gamma }}_{ij}^i{v}^j\\ f& =f(x)\text{.}\end{array}$$

The fluxes and sources simplify significantly when the background metric is flat and we employ Cartesian coordinates so ${\gamma }_{ij}={\delta }_{ij}$ and ${\mathrm{\Gamma }}_{jk}^i=0$. Set the template parameter BackgroundGeometry to Poisson::Geometry::FlatCartesian to specialize the system for this case. Set it to Poisson::Geometry::Curved for the general case.

Complex Poisson equation

This system can also be used to solve the complex Poisson equation where $u(x)$ and $f(x)$ are complex-valued, by setting the DataType template parameter to ComplexDataVector. Note that the real and imaginary sectors of the equations decouple, so they are essentially two independent Poisson equations. This is useful for testing the elliptic solver with complex-valued fields, and also as building blocks for other Poisson-like systems of equations that have additional complex-valued source terms.

---

The documentation for this struct was generated from the following file:

• src/Elliptic/Systems/Poisson/FirstOrderSystem.hpp