0

Cでスパース行列の割り当てを完了しようとしています。値と座標のリストとしてスパース行列を保持し、それをエール形式に変換しています。

私は、これまで誰も見たことがないような奇妙なメモリ割り当ての問題に遭遇しました。私のコードは次のとおりです。

yale* convertMatrix(matrix_list* input){
int matrix_elements = input->elements;
int matrix_rows = input->m;

yale* yale = (struct y*)calloc(1, sizeof(yale));

int* A = (int*)calloc(matrix_elements, sizeof(int));
int* IA = (int*)calloc(matrix_rows + 1, sizeof(int));    
int* JA = (int*)calloc(matrix_elements, sizeof(int));

printf("%d elements\n",matrix_elements);

yale->A = A;      // Value
yale->IA = IA;          // Row (X)
yale->JA = JA;     // Column (Y)
yale->elements = matrix_elements;
yale->m = matrix_rows;
yale->n = input->n;

list* tmp_list = input->first;

for(int i = 0, j = 0, tmp_y = 0; i < matrix_elements && tmp_list!=NULL; i++){
    printf("Input Value: %d \n",tmp_list->point.value);
    A[i] = tmp_list->point.value;
    // Initialise the first row
    if(i == 0) IA[0] = tmp_list->point.x;
    else{
        // Add a new row index
        if(tmp_y != tmp_list->point.x){
            j++;
            IA[j] = i;
            tmp_y = tmp_list->point.x;
        }
    }
    JA[i] = tmp_list->point.y;
    tmp_list = tmp_list->next;
}

for(int i = 0; i < matrix_elements; i++)
    printf("%d,",yale->A[i]);
printf("\n");
for(int i = 0; i < matrix_rows + 1; i++)
    printf("%d,",yale->IA[i]);
printf("\n");
for(int i = 0; i < matrix_elements; i++)
    printf("%d,",yale->JA[i]);

return yale;
}

そして、ここにエールの構造体があります：

typedef struct y{
int n;
int m;
int elements;
int *IA;
int *JA;
int *A;
} yale;

ただし、プログラムは、ループの最初の反復で最初の関連するprintfでsegfaultを実行します。

      printf("%d,",yale->A[i]);

私は肯定的です：matrix_elementsは整数です（私のテストケースでは9）matrix_rowsは整数ですA / IA / JAはすべて正しい値で満たされています（printfでyale-> AをAに交換すると正常に動作します）。配列を構造体ポインターに直接呼び出しても、結果には影響しません。型キャストではなく、マロッキング、キャロシング、すべて効果がありません。

Xcodeとgdbのおかげで、セグメンテーション違反の時点でもそれを確認できます。構造体ポインタが配列を指していないようです

4

2 に答える 2

0

コンパイルするには、これをスニペットの前に追加する必要があります。

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

typedef struct y{
        int n;
        int m;
        int elements;
        int *IA;
        int *JA;
        int *A;
        } yale;

typedef struct list {
        struct list *next;
        struct point { int x,y,value; } point;
        } list;

typedef struct matrix_list {
        int elements;
        int m;
        int n;
        struct list *first;
        int *point;
        } matrix_list;

更新：私はプログラムをより読みやすいものに変換しました（IMHO）。IAとJAが何をするのか、私にはわかりませんが、以下のフラグメントはOPと同等であるはずです。

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

struct y {
        unsigned int n;
        unsigned int m;
        unsigned int elements;
        unsigned int *IA;
        unsigned int *JA;
        int *A;
        } ;

struct list {
        struct list *next;
        struct point { unsigned int x,y; int value; } point;
        } ;

struct matrix_list {
        unsigned int elements;
        unsigned int m;
        unsigned int n;
        struct list *first;
        } ;

struct y *convertMatrix(struct matrix_list* input)
{
unsigned int matrix_elements = input->elements;
unsigned int matrix_rows = input->m;
unsigned int ii,jj,tmp_y;

struct y *yale ;
struct list *tmp_list ;

yale = calloc(1, sizeof *yale);
assert (yale != NULL);

printf("%u elements\n",matrix_elements);

yale->A = calloc(matrix_elements, sizeof *yale->A);
assert (yale->A != NULL);
yale->IA = calloc(matrix_rows + 1, sizeof *yale->IA);
assert (yale->IA != NULL);
yale->JA = calloc(matrix_elements, sizeof *yale->JA);
assert (yale->JA != NULL);

yale->elements = matrix_elements;
yale->m = matrix_rows;
yale->n = input->n;

    // Initialise the first row, set start condition
        // FIXME: this ignores the empty list or size=0 cases
yale->IA[0] = tmp_y = input->first->point.x;
ii = jj = 0;
for(tmp_list = input->first ;tmp_list; tmp_list = tmp_list->next) {
    printf("Input Value: %d \n",tmp_list->point.value);
    yale->A[ii] = tmp_list->point.value;
        // Add a new row index
    if(tmp_y != tmp_list->point.x){
        jj++;
        yale->IA[jj] = ii;
        tmp_y = tmp_list->point.x;
    }
    yale->JA[ii] = tmp_list->point.y;
    if (++ii >= matrix_elements ) break;
}

for(int i = 0; i < matrix_elements; i++)
    printf("%d,",yale->A[i]);
printf("\n");
for(int i = 0; i < matrix_rows + 1; i++)
    printf("%u,",yale->IA[i]);
printf("\n");
for(int i = 0; i < matrix_elements; i++)
    printf("%u,",yale->JA[i]);

return yale;
}

注：（ii == 0）{}条件をループの外に移動し、1文字のインデックスを2文字の同等のインデックスに置き換えました。また、すべてのインデックスは署名されていません（本来あるべきです）

于 2012-01-25T13:55:10.597 に答える
0

Valgrindでコードを実行することをお勧めします。これにより、バッファオーバーフローエラーが報告されます。（バッファオーバーフローは、配列の終わりを超えて書き込む場所です）。

また、コードの単体テストを作成することをお勧めします。それらはバグの検出に非常に役立ちます。特に、すべての位置に値を持つ3x3の入力行列を使用してテストを作成することをお勧めします。あなたが得る値があなたが期待するものであることを確認してください。

于 2012-01-29T00:52:53.310 に答える