Hallo,

ich habe hier zur Berechnung der Determinante einer quadratischen Matrix zwei Moeglichkeiten:

Einmal Berechnung ueber das Gauss-Jordan Verfahren und einmal ueber LU Decomposition.

Welches von beiden Verfahren ist denn schneller ? Die Matritzen sind so in der Groessenordnung von 1000 x 1000 und beinhalten nur natuerliche Zahlen.

Danke und Gruss HCN

Hab ich vor ca. 20 Jahren gemacht. Damals war LU-Dekomposition wesentlich schneller. Da sich das Verfahren nicht, sondern nur die CPUs in der Zwisc henzeit geändert haben, denke ich, Du müsstest mit LU-Dekomp hinkommen.

Das Gauß-Verfahren ist zudem wegen Rundungsfehlern viel ungenauer. Es könnte sein, das auch die LU-Dekomposition an seine Grenzen stößt bzw. Du mit höherer Genauigkeit rechnen musst.

In den Numerical Recipes (Kapitel 2.3) gibt es noch folgenden Hinweis:

Zitat:

For a matrix of any substantial size, it is quite likely that the determinant will overflow or underflow your computer’s floating-point dynamic range. In this case you can modify the loop of the above fragment and (e.g.) divide by powers of ten, to keep track of the scale separately, or (e.g.) accumulate the sum of logarithms of the absolute values of the factors and the sign separately.

(Die benutzen auch LU Zerlegung für die Berechnung der Determinante)

Darf ich fragen, wofür du den Wert brauchst?

Die Matrize beinhaltet HU Units aus der med. Bildgebung. Mehr darf ich leider nicht sagen....