Перейти к публикации

Суперкомпьютер своими руками для расчетов в ansys и 3d приложениях


Рекомендованные сообщения

Думаю, что METIS и создает А'. И именно из-за этого преобразования данный метод не решает матрицы с комплексными числами.

Там рядышком ссылка на эту книженцию..

George, Alan and Liu, Joseph W-H, Computer Solution of Large Sparse Positive Definite Systems, Prentice-Hall, Inc. (1981).

Товарищ майор, тебе пакет -

Ссылка на сообщение
Поделиться на других сайтах


Так оно же вроде гарантирует сходимость метода сопр. градиентов. А для прямых методов вроде без разницы.

Ссылка на сообщение
Поделиться на других сайтах

"Почему матрица должна быть положительно определенной" - корень квадратный даст комплексное значение в противном случае. Матрицы мкэ полуопределенные и достаточно часто такое выскакивает если влоб решать. Можно слегка возмутить диагональ например умножив все ее члены на 1.00000001 или что-то подобное и тогда проблем не возникает, а лучше по моему просто пользоваться LDLt разложением вместо LLt и тогда вообще проблем не возникает. Интересно еще поэкспериментировать с несколько диагональными D, которые легко обращаются. У Икрамова было описание трехдиагональной D в тоненькой книжке. Там же кстати было описание работы с матричными уравнениями, что позволило элементарно строить базисные функции для любых элементов c наперед заданными свойствами :)

Ну теперь-то они будут выдавать ваши правильные мысли за свои никчемные.

Это вряд ли. Просто гуманитарня всегда выступает тормозом прогресса. Когда-то Карамзин на 50 лет затормозил развитие капитализма в России, убедив Александра 1 отказаться от отмены крепостного права. Сейчас они цепляются за архаичную систему организации науки и делают то же самое...

Надо разделить на естественные и неестественные в организационном плане, иначе вечно будем следовать методу авторитетов, а не логики, а это не прагматично, как Пирс указал стрелкой Пирса и и заложив Основы прагматизма, когда знание полезных следствий и есть знание о предмете :)

Ссылка на сообщение
Поделиться на других сайтах

Это вряд ли. Просто гуманитарня всегда выступает тормозом прогресса. Когда-то Карамзин на 50 лет затормозил развитие капитализма в России, убедив Александра 1 отказаться от отмены крепостного права. Сейчас они цепляются за архаичную систему организации науки и делают то же самое...

У нас в кабинете истории всегда висела надпись, на которую я смотрел несколько лет...

"Можно не знать тысячи наук и быть образованным человеком, но не любить историю может только человек совершенно не развитый умственно". / Н.Г. Чернышевский

С тех пор я не люблю Чернышевского...

Ссылка на сообщение
Поделиться на других сайтах

Целиться стоит на парочку таких <noindex> таких </noindex> наверное, а для начала и одной хватит, то есть прикинуть программу максимум и минимум и брать железо по минимуму с прицелом на максимум развития :unsure:

Ссылка на сообщение
Поделиться на других сайтах
  • 2 недели спустя...

Хотел потестить режим GPU на 13-ой версии и вот что получилось:

While 1 GPU devices have been detected on this machine, no Tesla-series

GPU devices have been detected. Only Tesla-series GPU devices are

allowed. Please verify that a supported GPU device is installed on

this machine or rerun without using the -acc command line option.

Это что, только Tesla-series поддерживаются? А у меня GT220, 48 процессоров, должно же работать??! :confused_1: Кто-нибудь пробовал подключать видеокарту к расчетам?

PS: На fea.ru обещают скоро выложить результаты тестовых расчетов в ansys с использованием Tesla.

Ссылка на сообщение
Поделиться на других сайтах

Ну теперь-то они будут выдавать ваши правильные мысли за свои никчемные.

Научил редакцию ЛГ конституцию уважать в которой цензура в любой форме запрещена. Не захотели быть преступниками :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

Обычно конечноэлементные матрицы полученные грамотно, являются положительно обусловлеными.

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

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

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

Даю ссылочку на ребят, которые вывели, сформулировали и реализовали, собственную КЭ прогу для механики пористых деформируемых сред. Так вот - получили положительно определённую матрицу. Проводили расчёты на gpu и cpu. Максимальное увеличение быстродействия - в 11 раз - <noindex>http://daim.idi.ntnu.no/masteroppgave?id=4773</noindex>

Ссылка на сообщение
Поделиться на других сайтах

" являются положительно обусловлеными" - из производных их в принципе не сделать, так как с точностью до константы. Посмотрите как строятся метрики <noindex> Соболева </noindex> и матрица Грама в энергетическом пространстве :)

Грамотная запись - положительная определенность для обусловленности :unsure:

Видел даже какую-то нелепую докторскую в которой изобретались ограничения на форму элемента чтобы L Lt проходило, вместо того, чтобы не парясь использовать L D Lt и перейти на двойную точность :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

Видел даже какую-то нелепую докторскую в которой изобретались ограничения на форму элемента чтобы L Lt проходило, вместо того, чтобы не парясь использовать L D Lt и перейти на двойную точность :rolleyes:

Нам тоже препод по "Истории математики" (и такое было) расказывал... Типа прибегает к нему мужик и кричит - "Я решил <noindex>неразрешимую задачу трисекции угла</noindex>! Вы только поглядите какая точность!!!!"
Ссылка на сообщение
Поделиться на других сайтах
Ио́сиф Виссарио́нович

Думаю вопрос будет по теме.

Как заставить ANSYS вообще не использовать файл подкачки? Чтоб он торчал полностью в оперативке?

Вот поставлю себе, к примеру, 128 Гб ддэра, так ведь всёравно при расчёте в сисе будет расти файл подкачки.....Как его полностью обэлектронить и загнать в быструю память?

Может в этих 128 гигибайтах создать RAM диск и туда установить ansys ...? но так ведь после перезагрузки он оттуда сотрётся :(

Изменено пользователем Ио́сиф Виссарио́нович
Ссылка на сообщение
Поделиться на других сайтах

Так может просто на нем размещать рабочий каталог. Он и базу там держать будет, а программу то уж пусть с диска таскает и дллельки, это дело быстрое.

Да и файл подкачки можно на нем разместить :)

Нам тоже препод по "Истории математики" (и такое было) расказывал... Типа прибегает к нему мужик и кричит - "Я решил <noindex>неразрешимую задачу трисекции угла</noindex>! Вы только поглядите какая точность!!!!"

Мне массажист в бане рассказывал, что придумал способ рисовать звезду с помощью циркуля и линейки. Я принес ему

<noindex> листки </noindex> и показал ему, что он не Евклид, так как ошибся в рассуждениях. Теперь опасаюсь у него массаж делать :)

А ведь топологически элементарно "Правильный пятиугольник можно получить, завязав узлом полоску бумаги" не сложнее чем восемьсот разделить на три без остатка :rolleyes:

Ссылка на сообщение
Поделиться на других сайтах

" новый ПРЯМОЙ метод решения уравнений" - народился <noindex> новый Гаусс или Холецкий </noindex> ? :)

Ссылка на сообщение
Поделиться на других сайтах

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

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

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

А вторую половину - и вы одобряете, поэтому не цитировали.

Заметьте, я же не использую слов значения которых не понимаю, а вот прочие участники дискуссии - бывает:)

Я приветствую ваш предыдущий комментарий Фёдор.

Ссылка на сообщение
Поделиться на других сайтах

" метрики и матрицы никакого смысла не имеют применительно к моим рассуждениям" - это и печально :)

"требуется получать грамотно - то есть формулировать или выводить" - да это в любом букваре по мкэ есть. Квадратичный функционал записал от деформаций, продифференцировал и все дела :)

Попробуйте разложить любую элементную матрицу и убедитесь, что ничего не получится. Краевые условия ее улучшат, но ввиду ошибок округления легко получить слегка возмущенную и она не будет разлагаться в действительных числах, только комплексных, но это дело канительное. А норма комплексной части будет на уровне ошибки, поэтому ее можно без ущерба отбросить. Ошибки связаны с обусловленностью элементных матриц, а они зависят от искажений относительно порождающей. Чего тут переписчиков читать, взяли да поэкспериментировали в какой-нибудь математической программе. Делов то :rolleyes:

"Положительная обусловленность очень хорошее свойство матриц - ускоряет в несколько раз решение, даёт уверенность в том, что небольшие ошибки округления не извратят решение" - ускорение дает учет симметрии, а мкэшные матрицы легко вылетают при чистом LLt разложении. Раньше многие в это впиливались при реальном программировании.

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

Ссылка на сообщение
Поделиться на других сайтах

Ну так и напишите про новую программую реализацию такого-то метода, а не про новый метод :rolleyes:

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

Ссылка на сообщение
Поделиться на других сайтах

Всё верно, я потому и написала что положительная определённость избавляет нас с вами от влияния ошибок округления и от много того - что вы Фёдор написали. А кроме того, когда вы разогнались на грамотной записи, мне показалось что вы как-то сформулировали кое-чего, а матрицы-то оказались не очень-то хороши - сожалею. Мне конечно везло побольше - чаще всего я решаю свою слау по дедовски - оцениваю число обусловленности, работаю с cg или гауссом, а потом исправляю решение применяя итерационное уточнение. Если оказывается, что слау не совсем корректна:) используя регуляризацию тихонова - если получается здорово, если же не - переформулирую - то бишь перевывожу функционал в другую форму или применяю что-нибудь иное для дискретизации. Я с этим впервые стокнулась для случая когда i3=0, так что вы меня поймёте

Я с вами согласно по тому пункту, что

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

Потому, так важно чтобы получались именно положительные матрицы.
Ссылка на сообщение
Поделиться на других сайтах

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

"как-то сформулировали кое-чего, а матрицы-то оказались не очень-то хороши - сожалею" - я просто использовал LDLt и никаких проблем не испытывал, это же всего одну строчку кода немного иначе записать. Те кто используют чужое, так поступить не могут и мучаются, пока не сообразят умножить диагональные на 1.000000001 или что-нибудь подобное. Обычно в мкэ сумма внедиагональных равна диагональному, поэтому их и не считают, а просто заполняют диагональную матрицу и слегка при этом возмущают, чтобы не было проблем при разложении.

Удобно флажок в программе организовать, чтобы получать или LLt или LDLt. n проверок логического условия при вычислениях и под микроскопом не разглядеть на общем фоне.

" используя регуляризацию тихонова" - увеличение слегка диагонального примерно то же самое и никаких заморочек c дополнительными условиями :)

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

Есть единственный надежный способ убедится в положительной определенности матрицы - попробовать разложить ее в LLt :)

Увы, этот мир не совершенный, а грязная математика его отражение :unsure:

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах
  • MFS открепил тему
Гость
Эта тема закрыта для публикации сообщений.
  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • streamdown
      Коллеги приветствую! IPS 8. Подскажите пожалуйста, кто какое серверное железо использует? Интересуют параметры при одновременной работе, ну например, 400 пользователей онлайн
    • gudstartup
      так он так и позиционируется по m19 pmc выдает g70.6 а чпу отвечает f45.7 но ориентацию и смещение в 4077 он отрабатывает нормально шпиндель встает ровно и смена происходит хорошо. вопрос почему после ввода команды управления скоростью он все еще продолжает контролировать число импульсов между нуль метками хотя в принципе уже должен отменить позиционный контроль и просто считать обороты по 0 метке как он это делает без М19? это все понятно но почему оно продолжает проверять это после завершения ориентации мне непонятно
    • Александр1979
      SP9047 SSPA:47 ILLEGAL SIGNAL OF POSITION CODER "The relationship between the A/B phase and 1-rotation signal is incorrect (Pulse interval mismatch)." "Неправильное значение счетчика импульсов сигнала на энкодере ALPHAi. На фазах A и B энкодера за один оборот шпинделя насчитывается 4096 импульсов обратной связи. Программное обеспечение по управлению шпинделем проверяет количество импульсов на фазах A и B, соответствующее энкодеру, при каждой генерации сигнала одного оборота. Данный аварийный сигнал срабатывает, если регистрируется число импульсов, нарушающее заданный диапазон."
    • vs3dpro
      Добрый день! У нас на есть SLA принтер 600х600х400мм. Можно напечатать мастер- модели, и можно приехать посмотреть. mail@iges.space
    • gudstartup
      после ввода 789 в edit mod надо перейти в режим ручное упр. нажать mod а потом помощь и откроется wzw.hlp  с возможностью выполнения команд управления atc. если не появляется содержимое wzw.hlp то возможно что у вас другой код хотя 789 'это стандарт для dmg в этом случае надо смотреть ваш plc раздел у меня пароль прописан в файле wzw_rett.mod
    • Anat2015
      SERVO TURRET AMP ALARM - переводить не пробовали, учитывая, что AMP - это и есть сервопривод? Что тут не ясного, язык не знаком?
    • Viktor2004
      я бы попробовал позиционировать шпиндель сигналом G70.6 Ответный сигнал, шпиндель сориентирован F45.7
    • gudstartup
      Здравствуйте уважаемые участники форума попрошу вас поделится мнениями. Имеется фрезерный шпиндель с позиционным кодером (зубчатое колесо и считывающая головка). Ошибку 9047 выдает после того отрабатывает M19 c последующим включением оборотов Sxxx M3 причем обороты могут быть любыми от 100-макс. Если задавать обороты без предварительной ориентации то 9047 не возникает и шпиндель нормально вращается на любой скорости при этом обороты не гуляют точность поддержания +-1 об\мин. Кабель заменен заменить сам датчик нет возможности. Мой вопрос в следующем почему шпиндель после M19 и включении управления скоростью Sxxx M3 остается в режиме позиционирования а не сразу переходит врежим управления скоростью и даже явная отмена ориентации кодом M20 ничего не меняет ведь режим позиционного управления должен быть отменен или он остается и отменяется не сразу после подачи команды M3 и кода скорости? Станок достаточно старый 2009 года плату управления шпинделем меняли на соседний ничего не изменилось. Это может быть логика PMC или скорее всего это датчик?   PMC-SB.001
    • Killerchik
      А мы-то думали, что пропорциональный гидравлический.....
    • Grig75
×
×
  • Создать...