31 const ext::shared_ptr<FdmMesher>& mesher,
32 Real valueOnBoundary,
Size direction,
35 valueOnBoundary_(valueOnBoundary),
37 direction, side).getIndices()) {
40 xExtreme_ = mesher->locations(direction)[0];
44 ->locations(direction)[mesher->layout()->dim()[direction]-1];
55 for (
unsigned long indice :
indices_) {
1-D array used in linear algebra.
void applyAfterSolving(array_type &) const override
BoundaryCondition< FdmLinearOp >::Side Side
void applyAfterApplying(array_type &) const override
FdmDirichletBoundary(const ext::shared_ptr< FdmMesher > &mesher, Real valueOnBoundary, Size direction, Side side)
const Real valueOnBoundary_
void applyBeforeApplying(operator_type &) const override
void applyBeforeSolving(operator_type &, array_type &) const override
const std::vector< Size > indices_
#define QL_FAIL(message)
throw an error (possibly with file and line information)
Dirichlet boundary conditions for differential operators.
helper class to extract the indices on a boundary
memory layout of a fdm linear operator
std::size_t Size
size of a container