2009-06-07 4 views
0

J'utilise LAPACK pour diagonaliser des matrices hermitiennes complexes. Je peux choisir entre ZHEEV et ZHEEVD. Laquelle de ces routines est la plus précise pour les matrices de taille 40 et une gamme de valeurs propres de 1E-2 à 1E1?Précision de ZHEEV et ZHEEVD

Répondre

2

ZHEEVD utilise une méthode de division et de conquête pour calculer des valeurs propres.

Si vos matrices sont 40 x 40 et les valeurs propres sont dans l'intervalle [1E-2, 1e1] alors vous devriez avoir absolument aucun problème numérique. Vous pouvez utiliser n'importe quelle routine.

+0

Merci, c'est ce que je voulais entendre. –

0

Je ne connais pas la réponse mais,

  • Il dépend probablement de quelle bibliothèque LAPACK que vous utilisez. Il y en a plusieurs, optimisés pour diverses plateformes. Utilisez-vous Netlib, MKL, ACML, ??
  • Pourquoi prendriez-vous un mot d'un inconnu pour cela quand vous pouvez le mesurer vous-même?
+0

Je suppose que la précision varie beaucoup moins entre les implémentations que les performances. W/r à la mesure, je n'ai pas l'expérience pour préparer un ensemble correct de matrices de test. –

Questions connexes