Determinant matice – C / C++ – Fórum – Programujte.com

Dobré ráno,
pachtím se s algoritmem pro výpočet determinantu čtvercové matice, ale za boha mi to nechce šlapat a bůhví jestli je to vůbec dobře navrhnuto :/ pokud by měl někdo nějaký kus kódu, který řeší tento problém, byl bych mu velice vděčný.
Děkuji

int i, vysl = 1, diag=0, poddiag, mult = 1, add;

for (i = 0; i<m->rows; i++)

{

for (int j = 0; j<m->cols; j++)

{

diag = m->matrix[i][i];

if (m->matrix[i][j] != 0)

m->matrix[i][j] = diag/m->matrix[i][j];

}

mult = mult * diag;

for (int a = i; a<m->rows; a++)

{

poddiag = m->matrix[a][i];

for (int b = 0; b < m->cols; b++)

{

add = (-poddiag) * m->matrix[i];

m->matrix[a] = add+m->matrix[a];

}

}

}

for (i=0; i<m->rows;i++)

{

vysl *= m->matrix[i][i];

}

vysl *= mult;