pdsyev ругается с вычислительной сеткой 1x4 (2x2 - в норме)

И scilab, и octave могут использовать lapack для вычислений.

распределённых вычислений?

для вычислений

в дебиане есть теперь упакованный openblas или надо именно mpi?

Такого ещё не встречал, спасибо за инфо («Unlike Atlas, OpenBLAS provides a multiple architecture library», но под amd64 им пока не удалось 0.1.0 собрать).

В контексте вопроса мне интересны именно параллельные варианты.

atlas умеет openmp и на blas level 3 выдает количество gflops почти в точном соответствии со спецификациями процессора :)

openblas из репов у меня собрался именно на амд, или я чего то не понимаю?

ну и как атлас?39? или как его там, именно в дебиане поиметь?

Тогда можно обрадовать дебианщиков, потому что у них 0.1.0 только для i386 собран.

Тем временем, обнаружил подозрительный комментарий в pdsyev.f

*  Alignment requirements
*  ======================
*
*  The distributed submatrices A(IA:*, JA:*) and Z(IZ:IZ+M-1,JZ:JZ+N-1)
*  must verify some alignment properties, namely the following
*  expressions should be true:
*
*  ( MB_A.EQ.NB_A.EQ.MB_Z .AND. IROFFA.EQ.IROFFZ .AND. IROFFA.EQ.0 .AND.
*    IAROW.EQ.IZROW )
*  where
*  IROFFA = MOD( IA-1, MB_A ) and ICOFFA = MOD( JA-1, NB_A ).

Как-то совсем негибко получается.

Помнится тут месяц-другой назад обсуждали кластеры в универах. Что они - эти кластеры - есть, а считать некому. Вспомнить бы ники, чтобы позвать этих обсуждающих сюда. Ну не может быть, чтобы никто ни разу не пробовал посчитать собственные векторы на кластере. Я уже молчу о том, что запустить orterun (или как раньше mpirun) -np 4 можно и на стареньком десктопе.

На parallel.ru тоже что-то не очень людно на форуме. Где же обитают те редкие кадры, работающие на кластерах?..

Они там вот только 0.1.1 собрали. Тоже исключительно для i386. Кстати, вспомнил: один из тех людей - это AIv.

Где же обитают те редкие кадры, работающие на кластерах?..

Ну мы работаем.

Так ScaLAPACK используется на практике для распределённых решений задач линейной алгебры? Или использовалась бы, но нет потребности (тогда в чём потребность)?

Нами не используется. Но мы и не решаем задачи линейной алгебры;-)

Ясно. А кто тут ещё пробегал с кластера?

Reset должен быть в теме,Evgueni, marat2345.

Спасибо, тогда попросим заглянуть сюда Evgueni, marat2345. И вот только что нашёлся ещё один товарищ morse.

Не, я от линейной алгебры бесконечно далек. Хотя диплом свой писал именно по лапаку. Я там сравнивал intel MKL, sun perflib, gotoblas, и amd'шную либу название которой уже не помню. Сравнивал на домашнем интелловском четырехядернике, не на кластере. Результат был такой: MKL всех уделал, gotoblas работал круто и параллельно только на некоторых функциях, а на некоторых - сливал по-черному, все остальные сливали MKL в два-четыре раза, даже в однопоточном режиме.