public Vector Solve(Slae <ProfileMatrix> SLAE, Vector Initial, int maxiter, double eps)
{
int iterNum;
double pr, pr1, alpha, betta, residual;
Vector Az;
Vector r = new Vector(SLAE.RightPart.Differ(SLAE.Matrix.Multiply(Initial)));
Vector result = new Vector(Initial);
residual = r.Norm() / SLAE.RightPart.Norm();
Vector z = new Vector(r);
Vector s = new Vector(r);
Vector p = new Vector(r);
pr1 = p.Scalar(r);
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
Az = SLAE.Matrix.Multiply(z);
alpha = pr1 / Az.Scalar(s);
result = result.Sum(z.Mult(alpha));
r = r.Differ(Az.Mult(alpha));
p = p.Differ(SLAE.Matrix.TMultiply(s).Mult(alpha));
pr = p.Scalar(r);
betta = pr / pr1;
pr1 = pr;
z = r.Sum(z.Mult(betta));
s = p.Sum(s.Mult(betta));
residual = r.Norm() / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
return(result);
}
public Vector Solve(Slae <DiagonalMatrix> SLAE, Vector Initial, int maxiter, double eps)
{
int iterNum, vec_index;
double residual, diff, ML, MU;
Vector result = new Vector(Initial);
residual = SLAE.RightPart.Differ(SLAE.Matrix.Multiply(result)).Norm() / SLAE.RightPart.Norm();
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
for (vec_index = 0; vec_index < result.size; vec_index++)
{
ML = SLAE.Matrix.MultiplyL(vec_index, result);
MU = SLAE.Matrix.MultiplyU(vec_index, result);
diff = SLAE.RightPart.values[vec_index] - ML - MU;
result.values[vec_index] += diff / SLAE.Matrix.get_di(vec_index);
}
residual = SLAE.RightPart.Differ(SLAE.Matrix.Multiply(result)).Norm() / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
return(result);
}
public Vector Solve(Slae <ProfileMatrix> SLAE, Vector Initial, int maxiter, double eps)
{
int iterNum;
double iner_r, iner_mr, alpha, betta, residual;
Vector Az, Mr;
Vector r = new Vector(SLAE.RightPart.Differ(SLAE.Matrix.Multiply(Initial)));
Vector result = new Vector(Initial);
residual = r.Norm() / SLAE.RightPart.Norm();
if (Data.preconditioner == 0)//нет предобуславливания
{
Vector z = new Vector(r);
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
iner_r = r.Scalar(r);
Az = SLAE.Matrix.Multiply(z);
alpha = r.Scalar(r) / Az.Scalar(z);
result = result.Sum(z.Mult(alpha));
r = r.Differ(Az.Mult(alpha));
betta = r.Scalar(r) / iner_r;
z = r.Sum(z.Mult(betta));
residual = r.Norm() / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
else
if (Data.preconditioner == 1 || Data.preconditioner == 2 || Data.preconditioner == 3 || Data.preconditioner == 4) //диагональное предобуславливание или Неполное разложение Холесского
{
Vector z = new Vector(SLAE.PMatrix.ReverseProgress(SLAE.PMatrix.DirectProgress(r)));
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
Mr = SLAE.PMatrix.ReverseProgress(SLAE.PMatrix.DirectProgress(r)); //M^-1*r(k-1)
Az = SLAE.Matrix.Multiply(z); //A*z(k-1)
alpha = Mr.Scalar(r) / Az.Scalar(z);
result = result.Sum(z.Mult(alpha));
iner_mr = Mr.Scalar(r); //(M^-1*r(k-1),r(k-1))
r = r.Differ(Az.Mult(alpha));
Mr = SLAE.PMatrix.ReverseProgress(SLAE.PMatrix.DirectProgress(r));
betta = Mr.Scalar(r) / iner_mr;
z = Mr.Sum(z.Mult(betta));
residual = r.Norm() / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
return(result);
}
public Vector Solve(Slae <DisperseMatrix> SLAE, Vector Initial, int maxiter, double eps)
{
int iterNum;
double iner_p, alpha, betta, residual;
Vector Ar;
Vector result = new Vector(Initial);
if (Data.preconditioner == 0)//нет предобуславливания
{
Vector r = new Vector(SLAE.RightPart.Differ(SLAE.Matrix.Multiply(Initial)));
Vector z = new Vector(r);
Vector p = new Vector(SLAE.Matrix.Multiply(z));
residual = Math.Sqrt(r.Scalar(r)) / SLAE.RightPart.Norm();
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
iner_p = p.Scalar(p);
alpha = p.Scalar(r) / iner_p;
result = result.Sum(z.Mult(alpha));
r = r.Differ(p.Mult(alpha));
Ar = SLAE.Matrix.Multiply(r);
betta = -p.Scalar(Ar) / iner_p;
z = r.Sum(z.Mult(betta));
p = Ar.Sum(p.Mult(betta));
residual = Math.Sqrt(r.Scalar(r)) / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
else
if (Data.preconditioner == 1 || Data.preconditioner == 2 || Data.preconditioner == 3 || Data.preconditioner == 4) //LU или LU(sq)-предобуславливание
{
Vector fAx = new Vector(SLAE.RightPart.Differ(SLAE.Matrix.Multiply(result)));
Vector r = new Vector(SLAE.PMatrix.DirectProgress(fAx));
Vector z = new Vector(SLAE.PMatrix.ReverseProgress(r));
Vector p = new Vector(SLAE.PMatrix.DirectProgress(SLAE.Matrix.Multiply(z)));
if ((r == null) || (z == null) || (p == null))
{
return(null);
}
residual = Math.Sqrt(r.Scalar(r)) / SLAE.RightPart.Norm();
for (iterNum = 0; iterNum < maxiter && residual >= eps; iterNum++)
{
iner_p = p.Scalar(p);
alpha = p.Scalar(r) / iner_p;
result = result.Sum(z.Mult(alpha));
r = r.Differ(p.Mult(alpha));
Ar = SLAE.PMatrix.DirectProgress(SLAE.Matrix.Multiply(SLAE.PMatrix.ReverseProgress(r)));//L^-1*A*U^-1*r(k)
betta = -p.Scalar(Ar) / iner_p;
z = SLAE.PMatrix.ReverseProgress(r).Sum(z.Mult(betta));
p = Ar.Sum(p.Mult(betta));
residual = Math.Sqrt(r.Scalar(r)) / SLAE.RightPart.Norm();
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum, residual, maxiter, false))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
if (!autotest)
{
if (!InputOutput.OutputIterationToForm(iterNum - 1, residual, maxiter, true))
{
MessageBox.Show("Ошибка при выводе данных на форму.", "Опаньки...", MessageBoxButtons.OK);
}
}
}
return(result);
}