2009-06-07 8 views
0

Ich verwende LAPACK, um komplexe hermitesche Matrizen zu diagonalisieren. Ich kann zwischen ZHEEV und ZHEEVD wählen. Welche dieser Routinen ist genauer für Matrizen der Größe 40 und eine Reihe von Eigenwerten von 1E-2 bis 1E1?Genauigkeit von ZHEEV und ZHEEVD

Antwort

2

ZHEEVD verwendet eine Divide-and-Conquer-Methode zur Berechnung von Eigenwerten.

Wenn Ihre Matrizen 40 x 40 sind und die Eigenwerte innerhalb des Bereichs [1e-2, 1e1] liegen, dann sollten Sie absolut keine numerischen Probleme haben, . Sie können jede Routine verwenden.

+0

Danke, das wollte ich hören. –

0

Ich weiß nicht, die Antwort aber,

  • Es hängt wahrscheinlich von dem LAPACK-Bibliothek Sie verwenden. Es gibt eine Reihe von ihnen, optimiert für verschiedene Plattformen. Verwenden Sie Netlib, MKL, ACML, ??
  • Warum sollten Sie das Wort eines Fremden nehmen, wenn Sie es selbst messen können?
+0

Ich würde vermuten, dass die Genauigkeit zwischen Implementierungen viel weniger als die Leistung variieren wird. W/r zum Messen, ich habe nicht die Erfahrung, um einen korrekten Satz von Testmatrizen vorzubereiten. –