public Vector MultiplyTransposedFIrc(Vector vector)
{
Preconditions.CheckMultiplicationDimensions(Order, vector.Length);
// FIrc[s]^T * x = sum_over_s( Lc[s]^T * (Krc[s]^T * (inv(Krr[s]) * (Br[s]^T * x))) )
var result = Vector.CreateZero(numCornerDofs);
foreach (int s in Br.Keys)
{
IFetiDPSubdomainMatrixManager matrices = matrixManagers[s];
Vector temp = Br[s].Multiply(vector, true);
temp = matrices.MultiplyInverseKrrTimes(temp);
temp = matrices.MultiplyKcrTimes(temp);
temp = Lc[s].Multiply(temp, true);
result.AddIntoThis(temp);
}
return(result);
}
public Vector CreateCoarseProblemRhs(FetiDPDofSeparator dofSeparator,
Dictionary <int, IFetiDPSubdomainMatrixManager> matrixManagers,
Dictionary <int, Vector> fr, Dictionary <int, Vector> fbc)
{
// Static condensation for the force vectors
var globalFcStar = Vector.CreateZero(dofSeparator.NumGlobalCornerDofs);
for (int s = 0; s < subdomains.Count; ++s)
{
IFetiDPSubdomainMatrixManager matrices = matrixManagers[s];
// fcStar[s] = fbc[s] - Krc[s]^T * inv(Krr[s]) * fr[s]
// globalFcStar = sum_over_s(Lc[s]^T * fcStar[s])
UnsignedBooleanMatrix Lc = dofSeparator.CornerBooleanMatrices[s];
Vector temp = matrices.MultiplyInverseKrrTimes(fr[s]);
temp = matrices.MultiplyKcrTimes(temp);
Vector fcStar = fbc[s] - temp;
globalFcStar.AddIntoThis(Lc.Multiply(fcStar, true));
}
return(globalFcStar);
}