Российские САЕ

[Fedor 1 613]

Если у человека матрицы идут от сеток то там обычно ленты и стрелки. То есть структура фиксирована и достаточно эффективна. Быстрее прогонки Вам все равно не сосчитать на таких матрицах

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

[Eugeen 6]

Начните с полной факторизации. Переставьте строки и столбцы для уменьшения заполнения матрицы при решении. И решайте разреженным решателем.

Почитать можно здесь

<noindex>http://www.ispa-soft.ru/statxi/statxq15.htm</noindex>

Почитайте, в нашем докладе все так и сделано. Из этой процедуры мы высосали все что модно. Ускорились в 125 раз.

Но надо двигаться дальше. Нужны алгоритмы распараллеливания для работы уже с оптимизированной разреженной матрицей.

По крайней мере ленточные вроде остаются ленточными с заполнением внутри ленты. Насколько понял у человека речь не об ускорении, а о возможности вообще что-нибудь решить. Все-таки это разные задачи...

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

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

Ребята, мы истребитель давно освоили, нам нужен фотонный околосветовой звездолет. Давайте не предлагать тривиальные решения, мы это давно прошли.

Если у человека матрицы идут от сеток то там обычно ленты и стрелки. То есть структура фиксирована и достаточно эффективна. Быстрее прогонки Вам все равно не сосчитать на таких матрицах

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

Нет у нас никаких лент Я же привел для скачивания матрицу. Не матрица, а крякозабр какой-то. Это потому, что одновременно решается связная задача - гидродинамика - температурные поля в конструкциях - нейтронная кинетика. Ну откуда там взятся ленте, это больше похоже на архипелаг экзотических островов (Таити).

[Fedor 1 613]

Наверное паритесь в поисках новизны да актуальности. Флаг в руки cо <noindex> штучной уникалкой </noindex>

Можете попробовать лимон баксов с Боинга срубить за фундаментальное решение уравнения Навье-Стокса

"можно ли найти некую другую матрицу, более простого вида (3-х или 5-ти диагональную, например) для использования ее в итерациях вместо исходной?" у Икрамова было разложение типа L D Lt где D - трехдиагональная. Можно поискать и для 5 и вообще для любого числа диагоналей. Взять Mathematica и поискать алгоритм например для размерности 5-7 а потом распространить по подобию на n . Для Ваших целей можно будет L загрублять до единичной разумно полагаясь, что основная информация группируется в окрестности узла, а более дальними взаимодействиями можно пренебречь ...

[Eugeen 6]

Наверное паритесь в поисках новизны да актуальности. Флаг в руки cо <noindex> штучной уникалкой </noindex>

Можете попробовать лимон баксов с Боинга срубить за фундаментальное решение уравнения Навье-Стокса

Нет, цель чисто утилитарная - проведение многовариантных расчетов для обоснования безопасности АЭС.

[Fedor 1 613]

Домики и конструкция не проблема, ну а с жидкостями и газами трубки потолще сделайте и все дела. Куда они денутся. :)

"для обоснования безопасности АЭС" - <noindex> Сначала прикиньте </noindex> насколько надежность чувствительна к повышению точности расчетов

[Eugeen 6]

Домики и конструкция не проблема, ну а с жидкостями и газами трубки потолще сделайте и все дела. Куда они денутся. :)

С толстыми трубами Росатом не впарит наш проект даже в северную Корею. Это как покупать ИЛ-86 вместо Боинга-787.

[Fedor 1 613]

Повышение толщины трубки в сметах под микроскопом не разглядеть, а люди увидят, что все толстое и крепкое, наоборот можно будет надежность приподнять, и цену, соответственно. Витрувий учил что главное в строительстве прочность - польза - красота, нескупое расходование материалов и опираться на материк :)

[Eugeen 6]

Вам нужно определитьcя. Для начала написать техническое задание. А то задача меняется каждый день.

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

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

Домики и конструкция не проблема, ну а с жидкостями и газами трубки потолще сделайте и все дела. Куда они денутся. :)

"для обоснования безопасности АЭС" - <noindex> Сначала прикиньте </noindex> насколько надежность чувствительна к повышению точности расчетов

В Атомной энергетике это не прокатит. Там есть скурпулезное перечисление процедур, норм и правил при обосновании безопасности. И будь Вы даже Президентом Всемирной Академии Наук, Вас пошлют даааалеко при непрохождении даже одной процедуры, какой бы дурацкой она Вам ни казалась.

С другой стороны Вам скажут, что при всем соблюдении норм и правил, оборудование АЭС должно быть дешевым и экономичным. Иначе Арева, Вестингауз, Дженерал Электрик займут Ваше место на рынке. У них все тип-топ с правилами, а трубы тоньше (дешевле). Вот здесь то и нужны расчеты, модели, точность, быстродействие и т.п. и т.д.

[Fedor 1 613]

А потом как в Японии разгребают :)

[a_schelyaev 367]

Так выложили уже матрицу - качайте.

[a_schelyaev 367]

Дык, у них может быть там ветки кода в принципе не параллятся.

По уму все надо изначально под нее родную писать, в том числе под MPI.

[Eugeen 6]

Так нужен код, который нужно распараллелить.

Выкладывайте. Поможем.

Вот этот код мы используем:

<noindex>http://www.cise.ufl.edu/research/sparse/umfpack/</noindex>

[Eugeen 6]

Если вам действительно нужен "звездолет" для решения больших задач, то вам нужно менять оптимизатор перестановки строк и столбцов.

Для матриц 1000000х1000000 и большего размера алгоритм минимальной степени не подходит. Скорость решения для таких матриц в основном зависит от оптимизатора.

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

Бесплатно вы найдете исходный код только для маленьких матриц.

У кого купить?

[Eugeen 6]

Так как вы работаете на PC то посмотрите пакет от Интел. Они утверждают, что на 8-ми ядрах считают в 7 раз быстрее.

Или в Интернете поищите.

Мое мнение, что на ускорение двнного кода вы зря потратите и время и деньги. Но я его не навязываю.

а НАЗВАНИЕ ПАКЕТА?

[Fedor 1 613]

<noindex>http://software.intel.com/ru-ru/articles/intel-mkl/</noindex> наверное речь об этом?

[Борман 2 389]

a_schelyaev

Побывал тут на днях в окресностях Сарова, в результате чего узрел кроме всего прочего продукт Л0Г0C-Adaptive. Я так понимаю, это FV в чистом виде. Даже код тот же. Есть ли участие ТЕТИСа в прочих продуктах Л0Г0C-CFD ?

[a_schelyaev 367]

Сгружаем кроме всего прочего решалку СЛАУ, которую они могут прикручивать где угодно - нам не сообщают.

:)

[Apologist 0]

a_schelyaev, а у них в пресс-релизах про вас и решалку где-нибудь упоминается?

Нет, я не ёрничаю, просто интересно :)

А код вы им в исходниках отдаёте? А вдруг там "закладка"? ВНИИЭФ же, не шутки, секретность и всё такое....

[a_schelyaev 367]

Пресс-релизы ихние я не отслеживал, а вот в своих докладах про нас упоминают.

Касаемо закладок, то все сгружаем в соответствии с их требованиями.

[Apologist 0]

a_schelyaev, нагуглил сходу пару презентаций... Интересно, по цифрам круто и... неожиданно с точки зрения секретности. Обычно деталей от ВНИИЭФ не узнать, а тут вы чтоль на них влияете? Довольно подробно расписан функционал (техническая часть) и план усовершенствований.

Хотя может не планировали афишировать, а про вездесущий гугл забыли :)

Или вот финансовый вопрос.

a_schelyaev, это в какие деньги обходится этот проект... Вы же не дураки и понимаете, что растите конкурента, который уделает вас по всем статьям на российском рынке. Оборонку и госпредприятия принудят покупать Логос, а это же ваш сектор (он мог бы быть и нашим частично если б не экспорт-контроль) и огромные деньги.

Или они в ход утюг пустили? :)