수다닷컴

질문 제목 : 희소행렬의 곱셈 프로그램입니다질문 요약 :1.실행 후 오류를 분석하라는데 어디가 틀린지 감도 안오네요
어디가 틀린지 가르쳐주시구 수정 좀해주세요(오타가 있을수도있는데 오타수정이 중요한게아니라 이 프로그램이
왜 실행이 안되는지를 알아야합니다)질문 내용 :
void mmult(term a[], term b[], term d[])
{// a, b: 입력, d: 출력, d = a * b
int i, j, column, totald = 0;
int rows_a = a[0].row, cols_a = a[0].col, totala = a[0].value;
int cols_b = b[0].col, totalb = b[0].value;
int row_begin = 1, row = a[1].row, sum = 0; // row: a의 현재 행
term new_b[max_terms];
if (cols_a != b[0].row) {// a의 열의 수와 b의 행의 수는 동일
fpritf(stderr, “incompatible matricses\n”); exit(1);
}
fast_transpose(b, new_b);
// 경계 조건을 설정
a[totala+1].row = rows_a;
new_b[totalb+1].row = cols_b;
new_b[totalb+1].col = 0;
for (i = 1; i = totala; ) {
column = new_b[1].row;// b의 현재 열
for (j = 1; j = totalb + 1; ) {
// a의 현재 행과 b의 현재 열에 대해 곱셈 수행
if (a[i].row != row) { // a의 현재 행을 벗어남.
storesum(d, &totald, row, column, &sum);
i = row_begin;// b는 다음 열로. a는 원 위치로.
for (; new_b[j].row == cloumn; j++) ;
column = new_b[j].row;
}
else if (new_b[j].row != column) { // b의 현재 열을 벗어남.
storesum(d, &totald, row, column, &sum);
i = row_begin;// a는 원 위치
column = new_b[j].row // b는 다음 열로.
}
else switch (compare(a[i].col, new_b[j].col)){
case -1: // a[i].col new_b[j].col. a 증가
i++;break;
case 0: // 계산 후, a와 b를 모두 진행
sum += (a[i++].value * new_b[j++].value); break;
case 1: // a[i].col new_b[j].col. b 증가
j++;
}
}
for (; a[i].row == row; i++) ;// b의 모든 원소를 처리한 후,
row_begin = i;row = a[i].row;// a의 현재 행을 다음 행으로.
} /* end of for i = totala */
d[0].row = rows_a;
d[0].col = cols_b; d[0].value = totald;
}

void storesum(term d[], int *totald, int row, int column, int *sum)
{
/* sum이 0이 아니면, d 배열의 *totald+1 위치에 row, column 값과 함께 저장 */
if (*sum)
if (*totald max_terms) {
d[++*totald].row = row;
d[*totald].col = column;
d[*totald].value = *sum;
*sum = 0;
}
else {
fprintf(stderr, “numbers of terms in product
exceeds %d\n”, max_terms);
exit(1);
}
}

신청하기