Я имею дело с матрицами до N = 10 ^ 7 x N = 10 ^ 7; количество ненулевых элементов составляет примерно 6 x N. (Эти элементы сгруппированы по диагонали.) Моя RAM имеет размер 16 Gbt; поэтому мне явно нужен решатель разреженных матриц. Я запускаю Ubuntu LINUX и использую fortran90 (gfortran), точнее говоря, ratfor90.
У меня есть LAPACK, но он, похоже, не поддерживает решение разреженных матриц. (я ошибаюсь?) MATLAB должен быть хорош, но я не хочу тратить много времени на знакомство с ним; время поджимает. У меня установлен old / gold slatec, и я использую его для спецификаций. функции; есть ли у него редкие матричные процедуры?
Я слышал об ARPACK, но можно ли его использовать как простой решатель? можно ли его вызвать из гфортрана?
Любое другое предложение?
Спасибо, - Алекс