PARALLEL.RU

Дискуссионный клуб по параллельным вычислениям
Текущее время: 18 авг 19 4:21

Часовой пояс: UTC + 4 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 9 ] 
Автор Сообщение
 Заголовок сообщения: Параллельные решатели СЛАУ
СообщениеДобавлено: 15 окт 06 19:29 
Не в сети

Зарегистрирован: 22 май 06 1:34
Сообщения: 4
Откуда: г Саров
Здравствуйте уважаемые параллельщики!
Занимаюсь сабжем вот уже порядка пяти лет в математическом отделении РФЯЦ-ВНИИЭФ ( Российский Федеральный Ядерный Центр - Всероссийский Научно Исследовательский Институт Экспериментальной Физики).

До сих пор занимались созданием собственного продукта, и по ряду причин практически не оглядывались по сторонам. В качестве базы взяли неполное LU разложение + BiCGStab.
Сейчас настало время сравнить результат работы с какими либо внешними разработками, чтобы понять чего мы достигли и куда имеет смысл двигаться дальше.
Вот хотелось бы поинтересоваться у общественности, чем в мире люди считают СЛАУ, возникающие при решении задач газодинамики и теплопроводности на нерегулярных и неортогональных сетках, в параллельном режиме. Если можно ссылки на описания , обзоры. Для свободных продуктов буду благодарен за ссылки на исходники!


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 15 окт 06 22:50 
Не в сети

Зарегистрирован: 30 янв 04 2:32
Сообщения: 27
Откуда: Donetsk, Ukraine
PetSC, PBLAS.


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 16 окт 06 7:34 
Не в сети

Зарегистрирован: 12 янв 06 11:26
Сообщения: 98
Откуда: Хабаровск, ВЦ ДВО РАН
Интересно, что это за причины, по которым вы не огладывались по сторонам? Обычно сначала смотрят, что уже сделано в данной области, а уж потом начинают заниматься велосипедоизобретательством.


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 16 окт 06 12:53 
Не в сети

Зарегистрирован: 30 ноя 05 16:09
Сообщения: 130
Откуда: Ростов-на-Дону
Мы используем библиотеку Aztec.
Правда, в ней кроме того же BICGSTAB и других
методов Крылова ничего нет, но эти методы
реализованы весьма неплохо.
http://www.cs.sandia.gov/CRF/aztec1.html


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 16 окт 06 17:07 
Не в сети

Зарегистрирован: 27 апр 03 18:47
Сообщения: 41
Откуда: НЦ ПЭ, Петербург
PLAPACK использую для задач электродинамики. Однако назрел переход на более подходящий инструмент, потому что там общие методы, а у нас матрицы получаются разреженные.

http://www.cs.utexas.edu/~plapack/

_________________
Don't shoot the pianist. He's doing his best...


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 16 окт 06 18:24 
Не в сети

Зарегистрирован: 22 май 06 1:34
Сообщения: 4
Откуда: г Саров
"Интересно, что это за причины, по которым вы не огладывались по сторонам? Обычно сначала смотрят, что уже сделано в данной области, а уж потом начинают заниматься велосипедоизобретательством"

++ Я немного не корректно выразился. Имелось ввиду, что к моему приходу, уже было протестировано несколько библиотек, и они видимо были признаны неудовлетворяющими требованиям. Сравнение я затеял в частном порядке, поскольку видимо буду делать диссер и хочется не только расширить кругозор в теоретическом плане, но и пощщупать руками готовые творения других коллективов, на конкретных задачах, решаемых у нас в настоящее время. А до сего дня лично я этим не занимался просто за недостатком времени.

"Мы используем библиотеку Aztec.
Правда, в ней кроме того же BICGSTAB и других
методов Крылова ничего нет, но эти методы
реализованы весьма неплохо. "

Пакет я себе скачал, сейчас пытаюсь как раз откомпилировать.

"PLAPACK использую для задач электродинамики. Однако назрел переход на более подходящий инструмент, потому что там общие методы, а у нас матрицы получаются разреженные."

С такими же имеем дело мы.


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 16 окт 06 18:48 
Не в сети

Зарегистрирован: 30 ноя 05 16:09
Сообщения: 130
Откуда: Ростов-на-Дону
По-поводу библиотеки Aztec я уже писал
на этом форуме, что она содержит небольшие
ошибки. По крайней мере версия 2.0, с которой мы
работаем. По-моему, они не устранены и в более поздней
версии. Ошибка заключается в том, что внутри функций
библиотеки при каждом вызове выполняется аллокатирование
и освобождение памяти Так вот, не во всех функциях
это делается корректно и при каждом вызове программа
откусывет небольшой кусочек оперативной памяти.
Если итерационый процесс по временному параметру
длится достаточно долго ~ 500000 итераций, то с некоторого
момента происходит исчерпание оперативной памяти,
программа начинает сильно тормозить, а затем
вообще вылетает. Мы эту ошибку устранили.
Где-то на форуме я уже помещал сообщение об этом.


Вернуться к началу
 Профиль  
 
 Заголовок сообщения:
СообщениеДобавлено: 19 окт 06 14:18 
Не в сети

Зарегистрирован: 27 апр 03 18:47
Сообщения: 41
Откуда: НЦ ПЭ, Петербург
Когда-то было желание попробовать вот это:

Watson Sparse Matrix Package ------- http://www.alphaworks.ibm.com/eval/wsmp ... aw-eml-evl

даже скачал пробную версию, но руки не дошли.

_________________
Don't shoot the pianist. He's doing his best...


Вернуться к началу
 Профиль  
 
 Заголовок сообщения: Re: Параллельные решатели СЛАУ
СообщениеДобавлено: 12 мар 07 2:04 
Не в сети

Зарегистрирован: 21 дек 06 15:16
Сообщения: 121
Откуда: Нижний Новгород
stas9ih писал(а):
Здравствуйте уважаемые параллельщики!
Занимаюсь сабжем вот уже порядка пяти лет в математическом отделении РФЯЦ-ВНИИЭФ ( Российский Федеральный Ядерный Центр - Всероссийский Научно Исследовательский Институт Экспериментальной Физики). ...

Здоровеньки булы, сосед! :D :D
:arrow: Приглашаю включиться в работу сетевого совещания-семинара СОВРЕМЕННАЯ МАТЕМАТИЧЕСКАЯ ФИЗИКА И ЕЕ ПРИЛОЖЕНИЯ

:idea: А также приглашаю принимать участие и в наших очных заседаниях в Нижнем Новгороде --- следующее в середине мая.

Я думаю, два атомщика найдут, о чем поговорить и в виртуале, и в реале.
Атомная энергетика у нас здесь. Раздел ведет парень из Обнинска --- возможно, Вы его лично знаете.

ShapovalovTS писал(а):
Интересно, что это за причины, по которым вы не огладывались по сторонам? Обычно сначала смотрят, что уже сделано в данной области, а уж потом начинают заниматься велосипедоизобретательством.

:shock: Не надо давить нижегородцам на психику, глубокоуважаемый ShapovalovTS! Все правильно он делал --- сначала реши сам, потом посмотри, как делали другие. В физике это называется методологией Ландау. :lol:

_________________
Если тебе дадут линованную бумагу --- пиши поперек


Вернуться к началу
 Профиль  
 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 9 ] 

Часовой пояс: UTC + 4 часа [ Летнее время ]


Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 2


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  
cron
Создано на основе phpBB® Forum Software © phpBB Group
Русская поддержка phpBB