2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Тестовые задачи для солвера
Сообщение26.05.2008, 23:20 


26/11/06
76
Наконец таки дописал свой решатель систем линейных алгебраических уравнений. В нем есть как прямые так и итерационные методы. В основном он направлен на решение больших разреженных СЛАУ. Но есть одно ограничение, накладываемое на класс решаемых задач. Матрицы должны быть симметричные и положительно определены. Теперь бы хотелось протестировать его. Точнее я уже протестировал на множестве матриц, скачанных в интернете. Всё работает. Но всё таки хотелось бы взять реальные задачи, имеено из области строительной механики, из какого - нибудь коммерческого комплекса (SCAD,Lira, ES Stark и др. ) , занимающегося КЭ - анализом строительных конструкций. Только вот, существует две проблемы:
1. Нет примеров больших детальных 3D моделей зданий и сооружений.
2. Даже если бы такие примеры и были, не один из перечисленных коммерческих продуктов не позволяет импортировать матрицу жесткости и вектор правых частей.
Пробывал связаться с разработчиками но они на отрез отказались предоставить мне такие тестовые данные.
Может кто - нибудь из вас сталкивался с такой проблемой или знает где взять такие тесты?

 Профиль  
                  
 
 Re: Тестовые задачи для солвера
Сообщение27.05.2008, 00:12 


05/08/07

194
В пакет LAPACK включены довольно жесткие тесты.

 Профиль  
                  
 
 
Сообщение27.05.2008, 12:28 


26/11/06
76
Цитата:
В пакет LAPACK включены довольно жесткие тесты.


Ну во 1 - ых LAPACK предназначен больше для работы с плотными матрицами а у меня матрицы разреженные. И "жестких " тестов я там не нашел. Во основным не большие тестовые примеры для проверки работы каждой функции. Во 2-ых мне нужны задачи именно из области строительной... а не какой - либо другой.

 Профиль  
                  
 
 
Сообщение27.05.2008, 19:05 


05/08/07

194
vitaly333 писал(а):
Цитата:
В пакет LAPACK включены довольно жесткие тесты.


Ну во 1 - ых LAPACK предназначен больше для работы с плотными матрицами а у меня матрицы разреженные. И "жестких " тестов я там не нашел. Во основным не большие тестовые примеры для проверки работы каждой функции. Во 2-ых мне нужны задачи именно из области строительной... а не какой - либо другой.

Вы писали: "В основном он направлен на решение больших разреженных СЛАУ". Из этой фразы я сделал вывод, что Вы имеете дело и с плотными матрицами. Далее Вы писали: "Матрицы должны быть симметричные и положительно определены". Во многих задачах такие матрицы очень плохо обусловлены и с ними необходимо работать с ювелирной точностью. Поэтому я и написал про жесткие тесты. По поводу: "жестких " тестов я там не нашел". Плохо искали. Посмотрите мою страницу, посвященную диагонализации, конкретней - пункт 8 (для процессора P4). Тестовая программа генерируется с помощью Visual Studio (необходимо только реализовать "обертку" для своих функций). Т.к. на этом форуме не разрешают оставлять ссылки на свою страницу, то в Google задайте для поиска "Ю.Ф. Сиголаев".

 Профиль  
                  
 
 
Сообщение27.05.2008, 19:30 


26/11/06
76
Цитата:
Из этой фразы я сделал вывод, что Вы имеете дело и с плотными матрицами.


Нет.. только с разреженными.

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

Именно для таких вещей итерационные методы применяютя с предобуславливателями. На прямое решение плохая обусловленность не влияет.

Цитата:
Поэтому я и написал про жесткие тесты. По поводу: "жестких " тестов я там не нашел". Плохо искали. Посмотрите мою страницу, посвященную диагонализации, конкретней - пункт 8 (для процессора P4). Тестовая программа генерируется с помощью Visual Studio (необходимо только реализовать "обертку" для своих функций). Т.к. на этом форуме не разрешают оставлять ссылки на свою страницу, то в Google задайте для поиска "Ю.Ф. Сиголаев".

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

 Профиль  
                  
 
 Re: Тестовые задачи для солвера
Сообщение28.05.2008, 14:45 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
vitaly333 писал(а):
Но всё таки хотелось бы взять реальные задачи

Каков размер необходимой Вам матрицы?
Как Вы сравните свое решение с решением других пакетов?
vitaly333 писал(а):
существует две проблемы:
1. Нет примеров больших детальных 3D моделей зданий и сооружений.
2. Даже если бы такие примеры и были, не один из перечисленных коммерческих продуктов не позволяет импортировать матрицу жесткости и вектор правых частей.


Найдите 3D модель хотя бы одного здания или сооружения. Это тоже проблема.
Больших матриц в строительных МКЭ расчетах нет, т.к. используются линейные и оболочечные элементы и нет особой нужды в очень больших матрицах.

В коммерческих продуктах матрица жесткости и правая часть записываются в двоичном формате. При желании их можно считать другими программами.

 Профиль  
                  
 
 
Сообщение28.05.2008, 17:08 


05/08/07

194
vitaly333 писал(а):
abc_qmost писал(а):

Из этой фразы я сделал вывод, что Вы имеете дело и с плотными матрицами.



Нет.. только с разреженными.


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

 Профиль  
                  
 
 
Сообщение28.05.2008, 17:57 


26/11/06
76
Цитата:
Каков размер необходимой Вам матрицы?

От 10.000 до нескольких миллионов уравнений.

Цитата:
Как Вы сравните свое решение с решением других пакетов?

Импортирую ту же матрицу в другие пакеты (тот же TAUCS,MUMPS, INTEL MKL и др.). Если вы имели ввиду САПР то можно просто посчитать задачку в них (в том из которого и была вытащена матрица) и сравнить со своим.

Цитата:
Найдите 3D модель хотя бы одного здания или сооружения. Это тоже проблема.

Пока не могу нигде найти....

Цитата:
Больших матриц в строительных МКЭ расчетах нет, т.к. используются линейные и оболочечные элементы и нет особой нужды в очень больших матрицах.

Особой нужды нет, но всё таки большие матрицы там возникают, например при изучении совместной работы зданий и основания, при КЭ анализе сложных многоэтажных пространственных конструкций и др....в этих случаях порядок систем уравнений которые приходится решать достигает 600.000 и более. У меня даже есть несколько примеров (не самих моделей).

Цитата:
В коммерческих продуктах матрица жесткости и правая часть записываются в двоичном формате. При желании их можно считать другими программами.

Интересно что это за программы и с каких именно КЭ комплексов они могут произвести считывание, в какой формат пишут... и.т.д.

 Профиль  
                  
 
 
Сообщение29.05.2008, 09:16 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
vitaly333 писал(а):
Интересно что это за программы


Практически в каждом пакете есть утилиты работы с бинарными файлами. Например:
http://www.padtinc.com/epubs/focus/comm ... ticle1.htm
Цитата:
If you want matrices from ANSYS, use the new HBWRITE command in ANSYS 8.0 to create a file in Harwell-Boeing format.

 Профиль  
                  
 
 
Сообщение29.05.2008, 21:29 


26/11/06
76
Цитата:
If you want matrices from ANSYS, use the new HBWRITE command in ANSYS 8.0 to create a file in Harwell-Boeing format.


Да, но это же только в Ансисе. А для Скада, Лиры, МикроФЭ такие программы есть? То что из Ансиса матрицу вытащить можно это я знал. А вот найти модели зданий , созданных в Ансисе ещё сложнее чем найти аналогичные для Скада, Лиры.....

 Профиль  
                  
 
 
Сообщение30.05.2008, 05:58 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
vitaly333 писал(а):
найти модели зданий , созданных в Ансисе ещё сложнее чем найти аналогичные для Скада, Лиры.....


Модели из других пакетов можно перевести в модели ANSYSа.

 Профиль  
                  
 
 
Сообщение30.05.2008, 13:04 


26/11/06
76
Цитата:
Модели из других пакетов можно перевести в модели ANSYSа.


Это каким же образом?

 Профиль  
                  
 
 
Сообщение30.05.2008, 14:06 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
В коммерческих пакетах предусматривается архивация моделей в ASCII формате (данные по узлам, элементам, материалам и граничным условиям). Необходимы описания форматов и программа трансляции данных. Один пакет создает архивную модель. Программа обрабатывает и создает архивную модель для другого пакета. Другой пакет считывает модель.

 Профиль  
                  
 
 
Сообщение30.05.2008, 20:52 


26/11/06
76
Цитата:
В коммерческих пакетах предусматривается архивация моделей в ASCII формате (данные по узлам, элементам, материалам и граничным условиям). Необходимы описания форматов и программа трансляции данных. Один пакет создает архивную модель. Программа обрабатывает и создает архивную модель для другого пакета. Другой пакет считывает модель.


И что есть такие например SCAD->Ansys или Stark -> Ansys? Если есть то где их взять?

 Профиль  
                  
 
 
Сообщение02.06.2008, 12:00 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
Готовых программ SCAD->Ansys или Stark -> Ansys Вы вряд ли обнаружите. SCAD,Lira, ES Stark это российские программы. Возможно что ANSYS CivilFEM (http://www.civilfem.com) имеет транслятор в Ansys. У SCADа текстовые данные моделей *.txt в директории C:\WinSCAD\Examples\ можно перевести в стандартный формат команд Ansys, но для этого нужно знать формат данных SCADа.

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

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



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

Сейчас этот форум просматривают: DariaRychenkova


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

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group