private void solve_button_Click(object sender, EventArgs e)
{
//Запуск решения
bool success;
if (checkinput())
{
success = true;
richTextBox1.Clear();
Data.richtextbox = richTextBox1;
Data.ZedGraph.GraphPane.CurveList.Clear();
InputOutput.points.Clear();
Vector Initial;
ISolver solver;
double eps;
eps = 1.0;
for (int i = 0; i < eps_numericUpDown.Value; i++)
{
eps /= 10.0;
}
switch (Data.matrixformat)
{
case 0: //Плотная
{
Slae<DenseMatrix> SLAE = new Slae<DenseMatrix>();
SLAE.Matrix = new DenseMatrix(Data.matrixPath);
SLAE.RightPart = new Vector(Data.rightpartPath);
if (initial_checkBox.Checked) Initial = new Vector(SLAE.Matrix.getMatrix().N);
else Initial = new Vector(Data.initialPath);
if ((SLAE.Matrix.getMatrix().N == 0) || (SLAE.RightPart.SIZE == 0) || (Initial.SIZE == 0))
{
MessageBox.Show("Размерность матрицы/вектора правой части/вектора приближения равна нулю.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
if ((SLAE.Matrix.getMatrix().N != SLAE.RightPart.SIZE) || (SLAE.Matrix.getMatrix().N != Initial.SIZE))
{
MessageBox.Show("Размерности матрицы/вектора правой части/вектора приближения не совпадают.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
switch (Data.preconditioner)
{
case 0: //Нет предобуславливателя
{
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2: //БСГ стаб
{
solver = new BSGstab();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 3: //Гаусс-Зейдель
{
solver = new GaussZeidel();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
}
break;
}
case 1: //Диагональный
{
IPreconditioner<DenseMatrix> preconditioner = new DensePreconditioner();
if (!preconditioner.createDiag(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 2: //LLT
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного предобуславливателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
IPreconditioner<DenseMatrix> preconditioner = new DensePreconditioner();
if (!preconditioner.createLLT(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 3: //LU
{
IPreconditioner<DenseMatrix> preconditioner = new DensePreconditioner();
if (!preconditioner.createLU(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 4: //LUsq
{
IPreconditioner<DenseMatrix> preconditioner = new DensePreconditioner();
if (!preconditioner.createLUsq(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
}
break;
}
case 1: //Профильная
{
Slae<ProfileMatrix> SLAE = new Slae<ProfileMatrix>();
SLAE.Matrix = new ProfileMatrix(Data.matrixPath);
SLAE.RightPart = new Vector(Data.rightpartPath);
if (initial_checkBox.Checked) Initial = new Vector(SLAE.Matrix.getMatrix().N);
else Initial = new Vector(Data.initialPath);
if ((SLAE.Matrix.getMatrix().N == 0) || (SLAE.RightPart.SIZE == 0) || (Initial.SIZE == 0))
{
MessageBox.Show("Размерность матрицы/вектора правой части/вектора приближения равна нулю.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
if ((SLAE.Matrix.getMatrix().N != SLAE.RightPart.SIZE) || (SLAE.Matrix.getMatrix().N != Initial.SIZE))
{
MessageBox.Show("Размерности матрицы/вектора правой части/вектора приближения не совпадают.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
switch (Data.preconditioner)
{
case 0: //Нет предобуславливателя
{
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2: //БСГ стаб
{
solver = new BSGstab();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 3: //Гаусс-Зейдель
{
solver = new GaussZeidel();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
}
break;
}
case 1: //Диагональный
{
IPreconditioner<ProfileMatrix> preconditioner = new ProfilePreconditioner();
if (!preconditioner.createDiag(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 2: //LLT
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного предобуславливателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
IPreconditioner<ProfileMatrix> preconditioner = new ProfilePreconditioner();
if (!preconditioner.createLLT(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
//if (MessageBox.Show("Ваша матрица должна быть положительно определённой.\nЭто так?\nЕсли не уверены, исользуйте решатель МСГ.", "Обратите внимание!", MessageBoxButtons.YesNo, MessageBoxIcon.Information) == DialogResult.No) break;
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 3: //LU
{
IPreconditioner<ProfileMatrix> preconditioner = new ProfilePreconditioner();
if (!preconditioner.createLU(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 4: //LUsq
{
IPreconditioner<ProfileMatrix> preconditioner = new ProfilePreconditioner();
if (!preconditioner.createLUsq(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
}
break;
}
case 2: //Диагональная
{
Slae<DiagonalMatrix> SLAE = new Slae<DiagonalMatrix>();
SLAE.Matrix = new DiagonalMatrix(Data.matrixPath);
SLAE.RightPart = new Vector(Data.rightpartPath);
if (initial_checkBox.Checked) Initial = new Vector(SLAE.Matrix.getMatrix().N);
else Initial = new Vector(Data.initialPath);
if ((SLAE.Matrix.getMatrix().N == 0) || (SLAE.RightPart.SIZE == 0) || (Initial.SIZE == 0))
{
MessageBox.Show("Размерность матрицы/вектора правой части/вектора приближения равна нулю.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
if ((SLAE.Matrix.getMatrix().N != SLAE.RightPart.SIZE) || (SLAE.Matrix.getMatrix().N != Initial.SIZE))
{
MessageBox.Show("Размерности матрицы/вектора правой части/вектора приближения не совпадают.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
switch (Data.preconditioner)
{
case 0: //Нет предобуславливателя
{
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2: //БСГ стаб
{
solver = new BSGstab();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 3: //Гаусс-Зейдель
{
solver = new GaussZeidel();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
}
break;
}
case 1: //Диагональный
{
IPreconditioner<DiagonalMatrix> preconditioner = new DiagonalPreconditioner();
if (!preconditioner.createDiag(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 2: //LLT
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного предобуславливателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
IPreconditioner<DiagonalMatrix> preconditioner = new DiagonalPreconditioner();
if (!preconditioner.createLLT(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 3: //LU
{
IPreconditioner<DiagonalMatrix> preconditioner = new DiagonalPreconditioner();
if (!preconditioner.createLU(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 4: //LUsq
{
IPreconditioner<DiagonalMatrix> preconditioner = new DiagonalPreconditioner();
if (!preconditioner.createLUsq(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
}
break;
}
case 3: //Разреженая
{
Slae<DisperseMatrix> SLAE = new Slae<DisperseMatrix>();
SLAE.Matrix = new DisperseMatrix(Data.matrixPath);
SLAE.RightPart = new Vector(Data.rightpartPath);
if (initial_checkBox.Checked) Initial = new Vector(SLAE.Matrix.getMatrix().N);
else Initial = new Vector(Data.initialPath);
if ((SLAE.Matrix.getMatrix().N == 0) || (SLAE.RightPart.SIZE == 0) || (Initial.SIZE == 0))
{
MessageBox.Show("Размерность матрицы/вектора правой части/вектора приближения равна нулю.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
if ((SLAE.Matrix.getMatrix().N != SLAE.RightPart.SIZE) || (SLAE.Matrix.getMatrix().N != Initial.SIZE))
{
MessageBox.Show("Размерности матрицы/вектора правой части/вектора приближения не совпадают.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
switch (Data.preconditioner)
{
case 0: //Нет предобуславливателя
{
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2: //БСГ стаб
{
solver = new BSGstab();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 3: //Гаусс-Зейдель
{
solver = new GaussZeidel();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
}
break;
}
case 1: //Диагональный
{
IPreconditioner<DisperseMatrix> preconditioner = new DispersePreconditioner();
if (!preconditioner.createDiag(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 2: //LLT
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного предобуславливателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
IPreconditioner<DisperseMatrix> preconditioner = new DispersePreconditioner();
if (!preconditioner.createLLT(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 3: //LU
{
IPreconditioner<DisperseMatrix> preconditioner = new DispersePreconditioner();
if (!preconditioner.createLU(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
case 4: //LUsq
{
IPreconditioner<DisperseMatrix> preconditioner = new DispersePreconditioner();
if (!preconditioner.createLUsq(SLAE))
{
success = false;
break;
}
switch (Data.solver)
{
case 0: //МСГ
{
if (!SLAE.Matrix.CheckSymmetry())
{
MessageBox.Show("Для выбранного решателя ваша матрица должна быть симметричной.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
success = false;
break;
}
solver = new MSG();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 1: //ЛОС
{
solver = new LOS();
Data.result = solver.Solve(SLAE, Initial, (int)maxiter_numericUpDown.Value, eps);
break;
}
case 2:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
case 3:
{
MessageBox.Show("Для выбранного решателя предобуславливание недоступно.", "Такие дела.", MessageBoxButtons.OK, MessageBoxIcon.Information);
success = false;
break;
}
}
break;
}
}
break;
}
}
}
else success = false;
if (success)
MessageBox.Show("Сохраните ответ в файл на вкладке \"Вывод\".", "Решение завершено!", MessageBoxButtons.OK, MessageBoxIcon.Information);
else
MessageBox.Show("В процессе решения возникла ошибка.", "Решение не завершено!", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
private void Form1_KeyDown(object sender, KeyEventArgs e)
{
if (e.KeyCode == Keys.F2)
{
int i, n, iter;
Vector Initial;
ISolver solver;
IPreconditioner<DisperseMatrix> preconditioner = new DispersePreconditioner();
double eps;
if (System.IO.File.Exists("./auto/info.txt"))
{
string[] lines = System.IO.File.ReadAllLines("./auto/info.txt");
n = Int32.Parse(lines[0]);
for (i = 1; i <= n; i++)
{
eps = 1e-16;
Data.matrixformat = 3;
Data.matrixPath = "./auto/m" + i + ".txt";
Data.rightpartPath = "./auto/r" + i + ".txt";
System.IO.Directory.CreateDirectory("./auto/ans");
preconditioner.set_autotest(true);
if ((System.IO.File.Exists(Data.matrixPath)) && (System.IO.File.Exists(Data.rightpartPath)))
{
Slae<DisperseMatrix> SLAE = new Slae<DisperseMatrix>();
SLAE.Matrix = new DisperseMatrix(Data.matrixPath);
SLAE.RightPart = new Vector(Data.rightpartPath);
Initial = new Vector(SLAE.Matrix.getMatrix().N);
if ((SLAE.Matrix.getMatrix().N == 0) || (SLAE.RightPart.SIZE == 0) || (Initial.SIZE == 0))
{
MessageBox.Show("Размерность матрицы/вектора правой части/вектора приближения равна нулю.\nОшибка ввода матрицы.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
continue;
}
if ((SLAE.Matrix.getMatrix().N != SLAE.RightPart.SIZE) || (SLAE.Matrix.getMatrix().N != Initial.SIZE))
{
MessageBox.Show("Размерности матрицы/вектора правой части/вектора приближения не совпадают.", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Error);
continue;
}
if (SLAE.Matrix.getMatrix().N <= 1000) iter = 20000;
else iter = 2000;
if (SLAE.Matrix.CheckSymmetry())
{
Data.preconditioner = 0;
solver = new MSG();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
InputOutput.OutputVector("./auto/ans/x" + i + " MSG.txt", Data.result);
Data.preconditioner = 1;
if (preconditioner.createDiag(SLAE))
{
solver = new MSG();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " MSG_DIAG.txt", Data.result);
}
Data.preconditioner = 2;
if (preconditioner.createLLT(SLAE))
{
solver = new MSG();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " MSG_LLT.txt", Data.result);
}
Data.preconditioner = 3;
if (preconditioner.createLU(SLAE))
{
solver = new MSG();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " MSG_LU.txt", Data.result);
}
Data.preconditioner = 4;
if (preconditioner.createLUsq(SLAE))
{
solver = new MSG();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " MSG_LUsq.txt", Data.result);
}
}
Data.preconditioner = 0;
solver = new LOS();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " LOS.txt", Data.result);
Data.preconditioner = 1;
if (preconditioner.createDiag(SLAE))
{
solver = new LOS();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " LOS_DIAG.txt", Data.result);
}
if (SLAE.Matrix.CheckSymmetry())
{
Data.preconditioner = 2;
if (preconditioner.createLLT(SLAE))
{
solver = new LOS();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " LOS_LLT.txt", Data.result);
}
}
Data.preconditioner = 3;
if (preconditioner.createLU(SLAE))
{
solver = new LOS();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " LOS_LU.txt", Data.result);
}
Data.preconditioner = 4;
if (preconditioner.createLUsq(SLAE))
{
solver = new LOS();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " LOS_LUsq.txt", Data.result);
}
Data.preconditioner = 0;
solver = new BSGstab();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " BSGStab.txt", Data.result);
if (SLAE.Matrix.getMatrix().N < 50)
{
solver = new GaussZeidel();
solver.set_autotest(true);
Data.result = solver.Solve(SLAE, Initial, iter, eps);
solver.set_autotest(false);
if (Data.result != null) InputOutput.OutputVector("./auto/ans/x" + i + " GaussZeidel.txt", Data.result);
}
preconditioner.set_autotest(false);
if (i % (n / 10) == 0) MessageBox.Show("Матрица " + i + " протестирована.", "OK.", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
else
{
MessageBox.Show("Матрица " + i + " потеряна.", "OK.", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
MessageBox.Show("Автотестирование завершено.", "OK.", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
else
{
MessageBox.Show("Потерялся файл /auto/info.txt\n", "Ошибка.", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
}