Quelles sont les meilleures bibliothèques pour une grande itérative d'itératif (gradient de conjugué, minres, GMRES, etc.) de résolution de système algébra linéaire? J'ai souvent codé mes propres routines, mais je suis intéressé de savoir quels paquets "hors-la-étagère" préfèrent. J'ai entendu parler de Petsc, Tâcs, IML ++ et quelques autres. Je me demande comment ces empilement et quoi d'autre est là-bas. Ma préférence est la facilité d'utilisation et le logiciel librement disponible. P>
5 Réponses :
Je vérifierais la Fondation de Microsoft Solver. C'est gratuit pour bon marché pour même de très gros problèmes. La version illimitée est la force industrielle et est basée sur GURROBI et bien sûr n'est pas bon marché. P>
De votre lien, il semble que Microsoft Solver Foundation est une bibliothèque d'optimisation linéaire et non linéaire, connectée à une algèbre linéaire numérique, mais n'est pas le bon outil de résolution du système linéaire AX = B code> pour le vecteur
x code> donné une matrice de taille (code> A code> et le vecteur
b code>.
Vue d'ensemble des forfaits solvants de système linéaire itératif serait probablement un bon endroit pour commencer. P>
C'est génial, bien qu'un peu obsolète.
Vous voudrez peut-être aussi regarder des trilinos http://trilinos.sandia.gov/ p>
Il est conçu par un excellent artisan logiciel, en utilisant moderne Techniques de conception. P>
En outre, à partir de trilinos, vous pouvez appeler petsc si vous le souhaitez. P>
NIST a un peu de logiciel algébrique linéaire rare que vous pouvez télécharger Ici: http://math.nist.gov/sparselib++/ et ici: http://math.nist.gov/spblas/ p>
Je n'ai pas utilisé ces paquets moi-même, mais j'ai entendu de bonnes choses à leur sujet. P>
SPARSELIB ++ n'est qu'une bibliothèque de classe matricielle clairsemée. La bibliothèque de solveur itérative pertinente est iml ++ ( math.nist.gov/iml++ ), qui a été mentionné dans la question.
http: //www.cise.ufl.edu/research/sparse/umfpack/ p>
Umfpack est un ensemble de routines pour Résoudre peu sympatétrique linéaire Systèmes, AX = B, à l'aide du non -memetrique Méthode multifontalale. Écrit en ANSI / ISO C, avec un matlab (version 6.0 et ultérieure) interface. Apparaît comme un Routine intégrée (pour Lu, Backslash, et avancer slash) à Matlab. Inclut une interface MATLAB, un C-appelable interface et un appelant-appelable interface. Notez que "UMFPACK" est prononcé dans deux syllabes, "Umph Pack ". Ce n'est pas" VOUS EM EF PACK ". P> blockQuote>
Je l'utilise pour le code FEM. P>
Merci, bien que je pense que Umfpack est un solutionneur direct (pas un solveur itératif).