СЛАУ методом Зейделем

[Nadya13]

Всем добрый вечер. Надо решить СЛАУ методом Зейделя. Есть код, помогите, пожалуйста, найти ошибку. Ошибка в функции converge, когда эту функцию вызываю в цикле do...while функции iteration, происходит зацикливание. Подскажите, пожалуйста, как исправить ошибку.

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
double eps=0.001;
using namespace std;

bool InputFileMatr(double **A, int n)
{
    int i, j;
    FILE*f=fopen("matrix.txt", "r");
    if(!f)
        return false;
    for(i=0; i<n; i++)
    {
        for(j=0; j<n; j++)
        {
            fscanf(f, "%lf", &A[i][j]);
        }
    }
    fclose(f);
    return true;
}
bool InputFileVect(double *b, int n)
{
    int i;
    FILE *f=fopen("vector.txt", "r");
    if(!f)

        return false;
    for(i=0; i<n; i++)
    {
        fscanf(f, "%lf", &b[i]);
    }
    fclose(f);
    return true;
}
void OutputDescMatrix(double **A, int n)
{
    int i,j;
    printf("\nThe output matrix\n");
    for(i=0; i<n; i++)
    {
        for(j=0; j<n; j++)
        {
            printf("%lf ", A[i][j]);
        }
        printf("\n");
    }
    printf("\n");
}
void OutputDescVect(double *b, int n)
{
    int i;
    printf("\nThe output vector\n");
    for(i=0; i<n; i++)
    {
        printf("%lf ", b[i]);
        printf("\n");
    }
    printf("\n");
}
bool converge(double *xk, double *xkp, int n)
{
    double norm = 0;
    for (int i = 0; i < n; i++)
    {
        norm += (xk[i] - xkp[i])*(xk[i] - xkp[i]);
    }
    if (sqrt(norm) >= eps)
        return false;
    return true;
}
void iteration(double **a, double *b, double *x, int n)
{
 double p[n];
 for(int i=0; i<n; i++)
      x[i]=0;
do
{
    for (int i = 0; i < n; i++)
        p[i] = x[i];

    for (int i = 0; i < n; i++)
    {
        double var = 0.;
        for (int j = 0; j < i; j++)
            var += (a[i][j] * x[j]);
        for (int j = i + 1; j < n; j++)
            var += (a[i][j] * p[j]);
        x[i] = (b[i] - var) / a[i][i];
    }
}
while (!converge(x, p, n));
}

int main()
{
 int n;
    printf("Enter the size of matrix and vector\n");
    scanf("%d",&n);

    double **a=new double *[n];
    for(int i=0; i<=n-1; i++)
    {
        a[i]=new double [n];
    }

    double *b=new double [n];
    double *x=new double [n];

    InputFileMatr(a, n);
    InputFileVect(b, n);
    OutputDescMatrix(a, n);
    OutputDescVect(b, n);
    iteration(a, b, x, n);
    OutputDescVect(x, n);

    for(int j=0; j<=n-1; j++)
    {
        delete [] a[j];
    }
    delete [] a; a=NULL;
    delete [] b; b=NULL;
    delete [] x; x=NULL;
 return 0;
}

_____
Код программы нужно выделять (форматировать) тегами [CODE] (читать FAQ)
Модератор

[Sciv]

Это потому что у Вас converge всегда истинна. А истинна, потому что else забыли:

    
if (sqrt(norm) >= eps){
  return false;
}
else {
  return true;
}

а вообще - можно написать так и не заморачиваться:

return (sqrt(norm) >= eps);

[gantzboss]

Сообщение от Sciv

Это потому что у Вас converge всегда истинна. А истинна, потому что else забыли

Разве при return он не выходит, то есть если >= то он бы вышел и вернул false не идя дальше. Также там цикл while отрицательный если было бы тру он бы не циклился

[Sciv]

Сообщение от gantzboss

Разве при return он не выходит, то есть если >= то он бы вышел и вернул false не идя дальше. Также там цикл while отрицательный если было бы тру он бы не циклился

Ай чёрт, Ваша правда . Старею, наверное, азбуку забывать начал