Segmentation fault (core dumped)というエラーが出ます

Question

Cのプログラムについて質問です。
以下のように行列の行数、列数、そして行列そのものを読み込んでそれらをそのままprintfするプログラムをコンパイルして実行してみたところ、n(行数)とm(列数)はちゃんとprintfしてくれるのですが、その直後に行列をprintfする段階でSegmentation fault (core dumped)というエラーが出ます。

間違えている箇所がありましたら教えてください。
宜しくお願いします。

#include <stdio.h>
#include <stdlib.h>

#define NEW(p,n){p=malloc((n)*sizeof(p[0]));}

typedef double* row_vector;
typedef row_vector* matrix;

row_vector create_row_vector(int m)
{
    int i;
    row_vector V;
    NEW(V, m);
    for (i=0; i<m; i++) V[i]=0.0;
    return V;
}

matrix create_matrix(int n, int m)
{
    int i;
    matrix A;
    NEW(A, n);
    for (i=0; i<n; i++) A[i]=create_row_vector(m);
    return A;
}

void read_matrix(int n, int m, matrix A)
{
    int i, j;
    A=create_matrix(n, m);
    for (i=0; i<n; i++) {
        for(j=0; j<m; j++){
            scanf("%lf", &A[i][j]);
        }
    }
}

void print_matrix(int n, int m, matrix A)
{
    int i, j;
    printf("%d %d\n", n, m);
    for (i=0; i<n; i++) {
        for(j=0; j<m; j++){
            printf("%lf", A[i][j]);
        }
        printf("\n");
    }
}

int main()
{
    int n, m;
    matrix A;

    scanf("%d", &n);
    scanf("%d", &m);

    read_matrix(n, m, A);
    print_matrix(n, m, A);

    free(A);
    return 0;
}

Answer 1

read_matrix() 中で、確保された matrix Aですが、 ここで確保された領域は、 この関数の外では参照できません。
その結果、print_matrix()で、参照しようとした時、エラーとなります。
少なくとも、
read_matrix(int n, int m, matrix *A) とし、他もそれに合わせて、修正すべきです。

Comments

[aardvark]

ありがとうございます。

Answer 2

void read_matrix(int n, int m, matrix A)

Aはポインタの値渡しなので、この関数内でAに設定したアドレスは呼び出し元には返りません。

Comments

[aardvark]

ありがとうございます。
解決しました。

Answer 3

read_matrix(n, m, A);

print_matrix(n, m, A);
free(A);

ローカル変数Aが全く初期化されないままprint_matrixされ、freeされてます

Comments

[aardvark]

ありがとうございます。