public void BuildMatrices()
{
if (timeIncrement < BDForder - 1)
{
var coeffs = new ImplicitIntegrationCoefficients
{
Mass = 1 / timeStep / timeStep,
Damping = -1,
Stiffness = 1 / timeStep
};
foreach (ILinearSystem linearSystem in linearSystems.Values)
{
linearSystem.Matrix = provider.LinearCombinationOfMatricesIntoStiffness(coeffs, linearSystem.Subdomain);
}
}
else
{
var coeffs = new ImplicitIntegrationCoefficients
{
Mass = BDFcoeff[BDForder - 1, BDForder] / timeStep / timeStep,
Damping = -BDFcoeff[BDForder - 1, BDForder + 1],
Stiffness = BDFcoeff[BDForder - 1, BDForder + 1] / timeStep
};
foreach (ILinearSystem linearSystem in linearSystems.Values)
{
linearSystem.Matrix = provider.LinearCombinationOfMatricesIntoStiffness(coeffs, linearSystem.Subdomain);
}
}
}
public void BuildMatrices()
{
var coeffs = new ImplicitIntegrationCoefficients
{
Mass = 1,
Stiffness = 0
};
foreach (ILinearSystem linearSystem in linearSystems.Values)
{
linearSystem.Matrix = provider.LinearCombinationOfMatricesIntoStiffness(coeffs, linearSystem.Subdomain);
}
}