Cette question est une extension de deux discussions qui ont récemment été abordées dans les réponses à " C ++ vs Fortran pour HPC ". Et c'est un peu plus un défi qu'une question ... L'un des arguments les plus souvent entendus en faveur de Fortran est que les compilateurs sont simplement...
J'ai étudié les bibliothèques d'algèbre linéaire C ++ pour un projet sur lequel je travaillais. Quelque chose que je n'ai pas encore saisi est la connexion de BLAS et LAPACK à d'autres bibliothèques d'algèbre linéaire. En parcourant cet article sur les bibliothèques d'algèbre linéaire, j'ai trouvé...
Je voudrais savoir s'il existe un moyen rapide de calculer la distance euclidienne de deux vecteurs en octave. Il semble qu'il n'y ait pas de fonction spéciale pour cela, alors devrais-je simplement utiliser la formule avec
Étant donné deux implémentations BLAS différentes, pouvons-nous nous attendre à ce qu'elles effectuent exactement les mêmes calculs en virgule flottante et renvoient les mêmes résultats? Ou peut-il arriver, par exemple, que l'on calcule un produit scalaire comme et un comme donnant ainsi...
Désolé pour le long post, mais je voulais inclure tout ce que je pensais pertinent au premier coup. Ce que je veux J'implémente une version parallèle des méthodes de sous-espace de Krylov pour les matrices denses. Principalement GMRES, QMR et CG. J'ai réalisé (après profilage) que ma routine DGEMV...
Partout où j'ai vu, le didacticiel PETSc / documents, etc., dit qu'il est utile pour l'algèbre linéaire et spécifie généralement que les systèmes clairsemés en bénéficieront. Et les matrices denses? Je veux résoudreA x = bUNEX=bAx=b pour dense UNEUNEA. J'ai écrit mon propre code pour CG et QMR à...
nVidia, par exemple, a CUBLAS, qui promet une accélération de 7 à 14 fois. Naïvement, cela est loin du débit théorique de l'une des cartes GPU de nVidia. Quels sont les défis liés à l'accélération de l'algèbre linéaire sur les GPU, et existe-t-il déjà des routages d'algèbre linéaire plus...
L'implémentation de netlib BLAS est une excellente référence, étant pour la plupart non optimisée et bien documentée (par exemple zgemm ). Cependant, c'est dans Fortran 77, ce qui le rend quelque peu inaccessible à ceux qui ont une formation en programmation plus moderne. Existe-t-il une...
Le produit intérieur complexe a deux définitions différentes décidées par des conventions: ou . Dans BLAS, j'ai trouvé les routines cdotu, zdotu et cdotc, zdotc. Les deux anciennes routines calculent en fait (un faux produit intérieur!) Et les deux dernières routines conjuguent le premier vecteur...
La compilation de PETSc avec une bibliothèque BLAS / LAPACK externe affecte-t-elle de manière significative les performances sur des matrices clairsemées, ou utilise-t-elle uniquement ces bibliothèques pour les mathématiques matricielles