Ошибка сегментации (ядро сброшено) во время устранения Гаусса - PullRequest
Новая
       21

Ошибка сегментации (ядро сброшено) во время устранения Гаусса

0 голосов
/ 08 марта 2011

Во время работы этой программы я получаю ошибку сегментации.Я пытаюсь решить линейное уравнение, с числами, полученными из файла .dat.Я могу добраться до части прямо перед функцией гаусса, но потом она сбрасывает.Любая помощь?Благодарю.

    void prob2(void)
{
        FILE *matrix;

        matrix=fopen("matrix.dat", "r");
        double a1, a2, a3, a4, **A, *a, *rhs, length;
        int k, row;
        A=(double **)malloc(4*sizeof(double *));
        printf("\n");
        for(k=1;k<=4; k++) 
        {
                a=(double *)malloc(4*sizeof(double));
                fscanf(matrix, "%lf %lf %lf %lf", &a1, &a2, &a3, &a4);
                a[0]=a1;
                a[1]=a2;
                a[2]=a3;
                a[3]=a4;
                printf(" a[%d][] = %5.2f %5.2f %5.2f %5.2f\n", k, a[0], a[1], a[2], a[3]);
        }
        rhs=(double *)malloc(4*sizeof(double));
        fscanf(matrix, "%lf %lf %lf %lf", &a1, &a2, &a3, &a4);
        rhs[0]=a1;
        rhs[1]=a2;
        rhs[2]=a3;
        rhs[3]=a4;
        printf("\nb[]={  %.3f,   %.3f,  %.3f,   %.3f  }\n", rhs[0],rhs[1],rhs[2],rhs[3]);
        printf("hiii");
        gauss(4, A, rhs);
        /* print the solution x[] stored in rhs[]:*/
        printf("x[]={%7.3f, %7.3f, %7.3f, %7.3f  }\n", 
        rhs[0], rhs[1], rhs[2], rhs[3]);
        length=sqrt(pow(rhs[0],2)+pow(rhs[1],2)+pow(rhs[2],2)+pow(rhs[3],2));
        printf("The length of x[] is  %.6lf", length);
        /* free memory */
        for(row=0; row<4; row++) free(A[row]);
        free(A); 
        free(rhs);
        fclose(matrix);

}

void gauss(int n, double **A, double *rhs)
{
/* By Gauss elimination, solve a system of equations:
   A[][]*x[] = rhs[] where A[][] (n x n)
   and rhs[] (n x 1) are input */

/* x[] is stored in rhs[] */
printf("hiii");
double one = 1.0, zero=0.0;
double b, c, d;
int nm, row, col, krow;
nm = n - 1;
if(n == 1)
{
rhs[0] /= A[0][0]; A[0][0] = one;
return;
}
/* forward reduction */
for(row=0; row<nm; row++)
{
b = A[row][row]; A[row][row] = one;
for(col=row+1; col<n; col++) A[row][col] /= b;
rhs[row] /= b;
/* sweep rows of A[row+1][] to A[nm][] */
for(krow=row+1; krow<n; krow++)
{
c = A[krow][row]; A[krow][row] = zero;
for(col=row+1; col<n; col++) A[krow][col] -= c * A[row][col];
rhs[krow] -= c * rhs[row];
}
}
/* back substitution */

rhs[nm] /=A[nm][nm];
for(row=nm-1; row>=0; row--)
{
for(col=row+1; col<n; col++) rhs[row] -= A[row][col]*rhs[col];
}
return;
}

Ответы [ 3 ]

2 голосов
/ 08 марта 2011

Похоже, матрица A не назначена.Я думаю, что вам нужно это после выделения a: 

A[k-1] = a;
0 голосов
/ 08 марта 2011

Вы никогда ничего не храните в A. Вы используете malloc в цикле k = 1 ... 4, но вы никогда не сохраняете то, что вы используете malloc.

0 голосов
/ 08 марта 2011

Обычно это происходит, когда алгоритм встречает 0 на диагонали. Вы должны использовать поворот, чтобы избежать этого. Подробнее см. статью в Википедии об исключении Гаусса.

...