Comme le titre le dit, j'ai besoin de former une décomposition de LDL cholesky pour ma matrice définie positive A (Comme normal cholesky, mais il y a ceux une diagonale de L, et D est matrice diagonale). J'ai trouvé une seule fonction dans Lapack qui fait cela, mais elle dit que la matrice A doit être tridiagonale. Y at-il un peu de fonction qui fait cela dans certaines librairies de sous-routines libres comme lapack?Existe-t-il un sous-programme Fortran dans LAPACK/BLAS ou ailleurs pour calculer la décomposition des LDL?
1
A
Répondre
2
Vérifiez la fonction SSPTRF de LAPACK:
Objet
=======
SSPTRF calcule la factorisation d'une matrice symétrique réelle A stocké en format compressé avec le Bunch- procédé de pivotement diagonale Kaufman:
a = U * D * U ** T ou a = L * D * L ** T
où U (ou L) est un produit de permutat ion et unité supérieure (inférieure) matrices triangulaires, et D est symétrique et bloc diagonale avec blocs 1-by-1 et 2-by-2 diagonale.
En outre, le livre "Matrice Computations" Golub et Van Loan donne de un algorithme pour faire la décomposition. Dans ma troisième édition, c'est à la page 138, section 4.1.2 «Symétrie et factorisation LDL».
Questions connexes
- 1. mtrace pour un programme fortran
- 2. Décomposition des applications Adobe AIR
- 3. bibliothèque ou algorithme pour calculer des satellites GPS visibles
- 4. Décomposition de la réponse HTTP
- 5. Comment identifier la norme fortran - '77, '90 ou '95?
- 6. Pour calculer jours dans un mois
- 7. décomposition WPF - DataGridTemplateColumn
- 8. Qu'est-ce que la décomposition fonctionnelle?
- 9. Télécharger un fichier dans Django, puis l'envoyer ailleurs via REST?
- 10. Comment utiliser les compilateurs Fortran pour analyser ou imprimer des fichiers sources?
- 11. Déclarations de type dans Fortran
- 12. coclass dans l'interface d'importation .idl définie ailleurs?
- 13. Motif de conception ou odeur de code, données dénormalisées résultant de la décomposition fonctionnelle
- 14. Comment puis-je supprimer la sortie due à un SIGSEGV ou à un SIGFPE dans un programme Fortran 90?
- 15. FORTRAN WRITE()
- 16. Réflexion sur la propriété pour obtenir des attributs. Comment faire quand ils sont définis ailleurs? Je
- 17. Rationaliser Javascript pour calculer des gammes
- 18. Calculer le hachage ou la somme de contrôle pour une table dans SQL Server
- 19. Débogueur Fortran pour Mac OS X
- 20. Calculer la TVA dans SQL
- 21. Fonctionnalité de type Autodocumentation pour Fortran?
- 22. Calculer la complexité cyclomatique pour Javascript
- 23. Recommandations pour une bibliothèque XML Fortran
- 24. Convertir la valeur dans un format convivial URL - décomposition Unicode ähhh
- 25. Décomposition de la valeur singulière (SVD) en PHP
- 26. MySQL Requête pour calculer la semaine précédente
- 27. Difficulté avec la variable variable fortran
- 28. Date Différence dans MySQL pour calculer l'âge
- 29. Faire un double-défaut par défaut dans g77, compilateur Fortran
- 30. Problème avec la fonction de décomposition Eigenvalue de Jama