public DoubleDenseMatrix GetRightCorrectionFunc()
{
if (_cfRight == null)
{
_cfRight = new DoubleDenseMatrix(_polynomialOrder + 1, 1);
for (var i = 0; i < _polynomialOrder + 1; i++)
{
var deriv = Polynomials.Derivation(Polynomials.LeftRadau, _polynomialOrder + 1, _coordinates[i, 0]);
_cfRight[i, 0] = deriv;
}
}
return(_cfRight);
}
private DoubleDenseMatrix CreateDifferentiationMatrix()
{
var Dx = Polynomials.DiffLagrange(_coordinates);
return(Dx);
}