internal PcpgAlgorithm(IMaxIterationsProvider maxIterationsProvider, double residualTolerance,
IFetiPcgConvergence convergence)
{
this.residualTolerance = residualTolerance;
this.convergence = convergence;
this.maxIterationsProvider = maxIterationsProvider;
}
public Vector CalcLagrangeMultipliers(Feti1FlexibilityMatrix flexibility, IFetiPreconditioner preconditioner,
Feti1Projection projection, Vector disconnectedDisplacements, Vector rigidBodyModesWork, double globalForcesNorm,
SolverLogger logger)
{
// PCPG starts from the particular lagrange multipliers: λ0 = Q * G * inv(G^T * Q * G) * e
Vector lagranges = projection.CalcParticularLagrangeMultipliers(rigidBodyModesWork);
IFetiPcgConvergence pcpgConvergenceStrategy =
pcgConvergenceStrategyFactory.CreateConvergenceStrategy(globalForcesNorm);
var pcpg = new PcpgAlgorithm(maxIterationsProvider, pcpgConvergenceTolerance, pcpgConvergenceStrategy);
IterativeStatistics stats =
pcpg.Solve(flexibility, preconditioner, projection, disconnectedDisplacements, lagranges);
if (!stats.HasConverged)
{
throw new IterativeSolverNotConvergedException(Feti1Solver.name + " did not converge to a solution. PCPG"
+ $" algorithm run for {stats.NumIterationsRequired} iterations and the residual norm ratio was"
+ $" {stats.ResidualNormRatioEstimation}");
}
// Log statistics about PCG execution
logger.LogIterativeAlgorithm(stats.NumIterationsRequired, stats.ResidualNormRatioEstimation);
return(lagranges);
}
