mnożenie macierzy - język C

Question

Czy ktoś mógłby wytłumaczyć mi ten kod mnożenia macierzy? Jeśli chciałabym pomnożyć macierz przez wektor jak powinnam go zmienić?

void dgemm_naive (int n, double* A, double* B, double* C)
{
int i,j,k;
double cij;
for (i=0; i<n; ++i)
 for (j=0; j<n; ++j)
 {
cij = C[i+j*n]; // cij = C[i][j]
for(k=0; k<n; ++k )
 cij += A[i+k*n] * B[k+j*n]; // cij += A[i][k]*B[k][j]
C[i+j*n] = cij; // C[i][j] = cij
 }
}

 

Comments

Nauczyć się i napisać samemu, skoro to twoje zadanie.
Zapraszam do lektury https://eduinf.waw.pl/inf/utils/002_roz/2008_20.php

---

@natalia2002.,

1. Ten kod to taka "akumulacja" - nie tyle wykonuje samo mnożenie co też dodaje wynik do aktualnej zawartości macierzy wynikowej C (linia 8).

2. Operuje on na macierzach kwadratowych - jest tylko jeden parametr określający wielkość macierzy. Najbardziej elegancko byłoby dodać parametry, by można było mnożyć też macierze nie-kwadratowe. Wtedy macierz o jednym wymiarzy równym 1 byłaby wektorem i efekt osiągnięty. Można co prawda też okroić ten kod tylko do podanego przypadku mnożenia przez wektor - w sumie mnożenie macierzy to mnożenie pierwszej macierzy przez każdą kolumnę drugiej a kolumna to wektor.