internal override void Solve_BdEqualsa(double[] a)
{
// We have LB=UE or B =L(-1)UE. We need to solve L(-1)UEd=a, or UEd=La
//calculating La
foreach (Matrix m in LP)
{
a = m * a;
}
//solving Ud=r
U.SolveRightSystem(a);
//solving Ed=d
for (int i = 0; i < this.etaList.Count; i++)
{
etaList[i].SolveRightSystem(a);
}
}
