Hallo,
Fuer ein Projekt muss ich Gleichungen der Form A*X = B loesen, wobei A, X und B spaerlich besetzte Matrizen mit komplexen Zahlen und ohne bestimmte Symmetrien sind.
Was fuer eine Bibliothek kann ich dafuer benutzen?
Mein Programm ist in C++ geschrieben, alles was C- oder C++ bindings hat waere geeignet.
Ich habe bisher Eigen2 benutzt, das unter der Haube SuperLU benutzt, aber das kommt anscheinend nur fuer dicht besetzte X und B zurecht.
Meschach habe ich mir auch angeschaut, aber das sie fast so aus, als ob die Matrizen entweder Komplex oder Sparse sein koennten, aber nicht beides.
Habt ihr noch irgend welche Tipps fuer mich?
(Open source waere toll, kostenlos fuer akademische Zwecke waere auch akzeptabel).
Gruesse,
Moritz
Hallo Moritz
Habt ihr noch irgend welche Tipps fuer mich?
Sparse BLAS => http://math.nist.gov/spblas/
[C++]
...
The distribution is quite small (less than 14 Kbytes)
and it is meant as a starting for developing optimized
and architecture-dependent version. (C++ was used,
rather than C, as it has support for complex arithmetic
and templates to facilitate to creation of various
precision codes.) The library includes support for all
four precision types (single, double precision real
and complex) and Level 1, 2, and 3 operations.
...
(Open source waere toll, kostenlos fuer akademische Zwecke
waere auch akzeptabel).
Ich habe sie noch nicht verwendet (bei mir steht
wenigstens überall was in den Matrizen drin 
,
aber die Beschreibung liest sich nicht schlecht.
Grüße
CMb
Hallo Semjon,
Habt ihr noch irgend welche Tipps fuer mich?
Sparse BLAS => http://math.nist.gov/spblas/
9) B
Bei mir ist leider B sparse, nicht dense. Mit Algorithmen fuer dense B wird mein Programm leider zu langsam. Fuer dense B Bibliotheken zu finden ist relativ einfach, fuer sparse nicht - weshalb ich hier gefragt habe.
Trotzdem vielen Dank fuer deinen Tipp.
Gruesse,
Moritz