public double EstimateResidualNormRatio(PcgAlgorithmBase pcg)
{
var lagrangesBar = Vector.CreateZero(pcg.Solution.Length);
lagrangesBar.CopyFrom(pcg.Solution);
Vector lagranges = interfaceProblemSolver.CombineLagrangeMultipliers(lagrangesParticular, lagrangesBar, projection);
return(CalcExactResidualNorm(lagranges));
}
public void Initialize(PcgAlgorithmBase pcg)
{
} // Do nothing
public double EstimateResidualNormRatio(PcgAlgorithmBase pcg)
{
Debug.Assert(globalForcesNorm != 0.0, "norm2(globalForces) must be set first");
return(pcg.PrecondResidual.Norm2() / globalForcesNorm);
}
