Задача вроде не особо большая. Вы уверены, что без FMM не обойтись? Интересно было-бы посчитать напрямую, если бы была доступна сама модель.

дык Alex Z и писал, что "забил" на MLFMM и считал стандартным MOM.

меня вот больше другое удивляет, почему на одну итерацию тратится уйма времени.
из выходного файла видно, что оперативки для хранения матриц хватает.

2Alex Z: а в "основной" задаче сколько памяти затрачивается?
и вообще вы писали что у вас 4ГБ ОЗУ....а винда 32битная разве позволяет ее полностью использовать?
такое ощущение, что памяти не хватает и прога на винт ачинает сохранять....тогда точно долго считаться будет.

ОС 32 -битная, а памяти извиняюсь 3 Гига (при загрузке комп показывает около 3.2 Гига) комп Р3000.

Прога на винт может по-разному сохранять - бывает что просто система дает возможность аллокировать больше памяти чем есть и уходит в swapping. Программа работает как-будто с нормальной памятью, а на самом деле с винтом. Тогда все реально очень тормозит. В другом случае использаются специальные солверы ориентированные на out-of-core решение - сами аккуратно записывают что надо на винт и работают только с реальной памятью. В ЕМС Studiо это просто переключается автоматически (или если очень надо вручную), а в ФЕКО как я помню можно указать лимит памяти в pre файле через переменную #maxallocm = 1024 или сколько нужно (может сейчас уже есть и интерфесная поддержка). Нормальная производительность достигается - вполне терпимая. Начиная с версии 5.3 ФЕКО может аллокировать как 32-битная программа около 1.8 Гб, предыдущие версии - не более ~1.4 Gb.