public override void SetSup(TridiagOperator1D tridiaOp, double dx)
{
tridiaOp.SetSupBoundaryCondition((dx - 2.0) / (dx - 1.0), 1.0 / (dx - 1.0), 0.0);
}
public override void SetInf(TridiagOperator1D tridiaOp, double dx)
{
tridiaOp.SetInfBoundaryCondition((2.0 + dx) / (1.0 + dx), -1.0 / (1.0 + dx), 0.0);
}
public override void SetInf(TridiagOperator1D tridiaOp, double dx)
{
// V[-1] = 2 * V[0] - V[1]
tridiaOp.SetInfBoundaryCondition(2.0, -1.0, 0.0);
}
public override void SetSup(TridiagOperator1D tridiaOp, double dx)
{
// V[N+1] = 2 * V[N] - V[N-1]
tridiaOp.SetSupBoundaryCondition(2.0, -1.0, 0.0);
}
public abstract void SetSup(TridiagOperator1D tridiaOp, double dx);