public override IEnumerable <INumericalResult> Solve(int resultsCount)
{
List <INumericalResult> results = new List <INumericalResult>();
if (_mesh.IsMeshGenerated)
{
List <int> indeciesToDelete;
Matrix GeneralMatrix = GetGeneralMatrix(out indeciesToDelete);
Vector init = InitValue();
CauchyProblemResult cauchyProblemResult = CauchyProblemSolver.HeunMethodSolve((t, v) => GeneralMatrix * v, init, _maxTime, _intervalsTime);
//CauchyProblemResult cauchyProblemResult = CauchyProblemSolver.AdamsBashforthMethodsSolve((t, v) => GeneralMatrix*v, init, _maxTime, _intervalsTime);
//CauchyProblemResult cauchyProblemResult = CauchyProblemSolver.GirMethodsSolve((t, v) => GeneralMatrix*v, init, _maxTime, _intervalsTime);
SemidiscreteVibrationsNumericalResult result = new SemidiscreteVibrationsNumericalResult(_mesh.Elements, cauchyProblemResult.DeltaTime, _maxTime);
foreach (Vector v in cauchyProblemResult.Results)
{
AddToResult(result, v, indeciesToDelete);
}
results.Add(result);
}
return(results);
}
public override IEnumerable <INumericalResult> Solve(int resultsCount)
{
List <INumericalResult> results = new List <INumericalResult>();
if (_mesh.IsMeshGenerated)
{
GetConstantMatrix();
_indeciesToDelete = getIndeciesWithStaticBoundaryConditions();
Matrix StiffnessMatrix = AsumeStaticBoundaryConditions(GetStiffnessMatrix(), _indeciesToDelete);
Matrix MassMatrix = AsumeStaticBoundaryConditions(GetMassMatrix(), _indeciesToDelete);
_init = AsumeStaticBoundaryConditionsToVector(_init, _indeciesToDelete);
CauchyProblemResult newmarkResult = NewmarkSolver(MassMatrix, StiffnessMatrix, _init, 0.5, 0.5, _maxTime, _intervalsTime);
SemidiscreteVibrationsNumericalResult result = new SemidiscreteVibrationsNumericalResult(_mesh.Elements, newmarkResult.DeltaTime, _maxTime);
foreach (Vector v in newmarkResult.Results)
{
Vector resultWithStaticPoints = addStaticPoints(v, _indeciesToDelete);
result.AddResult(resultWithStaticPoints);
}
results.Add(result);
}
return(results);
}
private void AddToResult(SemidiscreteVibrationsNumericalResult result, Vector V, List <int> indeciesToDelete)
{
Vector u = GetU(V);
addStaticPoints(u, indeciesToDelete);
result.AddResult(u);
}
