Question

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

Était-ce utile?

La solution

ZHEEVD utilise un procédé de division pour régner pour calculer des valeurs propres.

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

Autres conseils

Je ne connais pas la réponse mais,

  • Il dépend probablement de quelle bibliothèque LAPACK que vous utilisez. Il y a un certain nombre d'entre eux là-bas, optimisé pour différentes plates-formes. Utilisez-vous Netlib, MKL, ACML, ??
  • Pourquoi voudriez-vous prendre prendre la parole d'un inconnu pour cela quand vous pouvez mesurer vous-même?
Licencié sous: CC-BY-SA avec attribution
Non affilié à StackOverflow
scroll top