고수님들 부탁드려요 ㅠmmult (희소행렬 곱셈)
파라나
sparse matrix의 곱셈을 구현하는 함수에요
[ C로 쓴 자료구조론 (Fundamentals of Data Structures in C번역본) ]
p.85 두 희소 행렬의 곱셈꼭좀 부탁드릴게요 내일이 시험인데
도저히 이 소스가 이해가 안가요 ㅠㅠㅠ
앞에 내용 그러니까 메트릭스를 sparse matrix로 바꾸는것까진 이해가 갔는데...
갑자기 곱이 나오는데 도통 뭔말인지 ㅠ 원리좀 알려주세요..
----------------------------------------------------------------------------------------------------main함수 생략.
void mmult (term a[], term b[], term d[])
{//두 희소 행렬을 곱한다.
int i, j, column, totalB=b[0].value, totalD=0;
int rowsA=a[0].row, colsA=a[0].col;
totalA=a[0].value; int colsB=b[0].col;
int rowBegin=1, row=a[1].row, sum=0;
int newB[MAX_TERMS][3];
if(colsA!=b[0].row) {
fprintf(stderr, Incompatible matrices\n);
exit(1);
}
fastTranspose(b,newB); //b행렬을 전치시킴. Transpose시키는 함수에요.
/* 경계 조건 설정*/ //이게 뭔말이에요? 경계조건이라뇽??;
a[TotalA+1].row=rowsA;
newB[totalB+1].row=colsB;
newB[totalB+1].col=0;
for(i=1; i=totalA; ) {
column=newB[1].row;
for(j=1; jtotalB+1;) {
/* a의 행과 b의 열ㅇ르 곱한다. */
if(a[i].row!= row) {
storeSum(d, &totalD,row,column,&sum);
i=rowBegin;
for(; newB[j].row==column; j++);
column=nowB[j].row;
}
else if(newB[j].row!=column) {
storeSum(d, &totalD, row, column, &sum);
i=rowBegin;
column=newB[j].row;
}
else switch(COMPARE(a[i].col, newB[j].col)) {
case -1 ://a의 다음 항으로 이동
i++; break;
case 0 ://항을 더하고, a와b를 다음 항으로 이동
sum+=(a[i++].value*newB[j++].value);
break;
case 1 : //b의 다음 항으로 이동
j++;
} //for j =totalB+1 end
for(; a[i].row==row; i++);
rowBegin=i; row=api].row;
} //for i=totalA end
d[0].row=rowsA;
d[0].col=colsB;
d[0].value=totalD;
}
아 쓰다보니 tab쓰니까 이상하게 되네요...
고수님들 부탁해요.. ㅜㅠㅠㅠ
