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); }