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

вопрос по sparse solver


Гость suvorov.as

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

Вы опять ничего не поняли.

Плохо обусловленные матрицы быстрее решать ПРЯМЫМ методом.

Ладно, поехали по новой:

прямой метод продувает с ростом размерности итерационному методу.

И также непонятна его эффективность на системах с распределеннной памятью.

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


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

В МГСУ была недавно конференция, где люди из Абакуса рассказывали про свои инструменты решателя.

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

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

Утром приходишь - невязки улетели в космос на пустом месте. Начали копаться и вскрыли в Ланцоше некоторые неприятные особенности. После разборок ввели ограничения некоторые на итерирование в отдельных местах. Но все пляски с бубном посчитали неправильным, и в итоге пришли к созданию SOFGMRES, где если сходится, то сходится.

Мне вот интересно, к задачам акустики в сильно неоднородных средах впринципе итерационные методы применимы?

Выхлоп в трубе сильно неоднородная среда?

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

Отрыл на просторах гугла помоему очень достойый современный курс лекций французов по СЛАУ

<noindex>http://perso.ens-lyon.fr/~bucar/CR07/index.html</noindex>

Может кому интересно будет

Кстати, может кто знает, как на русский корректно, применительно к матрицам, переводятся термины

elimination tree (Гугл говорит, что это - уничтожение деревьев :rolleyes: )

и

supernode

Спасибо!

elimination tree -- дерево исключения (<noindex>здесь</noindex> и <noindex>тут</noindex>)

supernode -- в какой лекции он впервые встречается?

PS Посмотрел сам внимательнее свои ссылки.. Походу, вся наука в Нижнем)

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

Supernode <noindex>тоже в лоб перевели</noindex>-- суперузел.

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

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

Тоже выкладываю одну главу книженции.

2 Solution of Linear Algebraic Equations

2.0 Introduction 32

2.1 Gauss-Jordan Elimination 36

2.2 Gaussian Elimination with Backsubstitution 41

2.3 LU Decomposition and Its Applications 43

2.4 Tridiagonal and Band Diagonal Systems of Equations 50

2.5 Iterative Improvement of a Solution to Linear Equations 55

2.6 Singular Value Decomposition 59

2.7 Sparse Linear Systems 71

2.8 Vandermonde Matrices and Toeplitz Matrices 90

2.9 Cholesky Decomposition 96

2.10 QR Decomposition 98

2.11 Is Matrix Inversion an $N^3$ Process? 102

1.rar

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

Вуаля

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

suvorov.as, для umfpack параллельной версии вроде не было. Или нашли?
Ссылка на сообщение
Поделиться на других сайтах

Гораздо хуже, что у россиян таких библиотек нет.

А если они есть, но их в интернете найти нельзя, то это совсем хреново.

Если верить докладу ВНИИЭФ на ПаВТ'12, то через пару годиков они выложат на суд общественности отечественную библиотеку решателей для СЛАУ всех мастей, которуя сейчас собирают из разработок различных авторских коллективов по стране и собственных. Туда войдут и итерационные схемы ТЕСИСа.

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

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

Фиг знает.

Вроде много лет уже считают этим методом течения в каналах.

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

У нас, просто, не принято фуфло толкать заказчику.

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

Нет, спасибо. Не интересует. В любом случае у них расчет с экспериментом не сошелся. Уж незнаю кто там прав, расчетчики или экспериментаторы.

Пусть возятся.

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

Продавец компании ТЕСИС возможно поможет и завтра, а возможно и через полгода, а возможно и через два. Ибо у ВНИИЭФ свои дела, а у ТЕСИСа свои. А возможно, что и вообще не поможет - ведь чудес на свете не бывает. Просто на форуме деловые вопросы не обсждаются. Взрослеть надо, ребята. Если нет формализованных отношений, то будет все течь само собой, как сейчас.

А касаемо данных Интела по распределенному Пардизо, то чуда не случилось.

post-1864-1333454500.png

Много лет назад Abaqus/Explicit показывал схожее ускорение на нелинейной задаче с контактом на неявном решателе.

Свежих данных под рукой нет.

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

Ссылка на сообщение
Поделиться на других сайтах
suvorov.as , подробности будут? Уже столько времени следим за сериалом. Они, наконец, поженились. Но интересно же как он делал предложение, как она думала, как родители были против.
Ссылка на сообщение
Поделиться на других сайтах
  • 2 недели спустя...

А мне стало интересно. Допустим все то же самое, стержень под 45, но сила все время под 30 градусов действует. ИСПА такую задачу решит?

30гр. к оси стержня ? Как бы вы стали решать эту динамику в АНСИС ?
Ссылка на сообщение
Поделиться на других сайтах

В нелинейной постановке интересно решит? Я такие не пробовал решать.

Естественное же решение, что стержень должен повернуться вдоль действующей силы

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

А без этого кому нужны миллиарды степеней.

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

За вас это делает Ансис. Ансис подкрепляет тысячи степеней свободы автоматически.

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

Подумайте сами. Вы только меняете систему координат и решатель отказывается решать корректную задачу.

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

Что в этом криминального?

*** WARNING *** CP = 49.531 TIME= 13:30:56

Small equation solver pivot term encountered at UY DOF of node 1.

ANSYS automatically constrained this DOF. Check for an insufficiently

constrained model.

----------

Factor Small Pivots Constrained List factor number= 3

Constrained zero pivot User node 1 DOF UY

Constrained zero pivot User node 2 DOF UX

Constrained zero pivot User node 2 DOF UY

ИСПА, зачем вы удаляете свои сообщения ?

suvorov.as

всегда отвечайте с цитированием сообщений от ИСПА, а то потом непонятно будет.

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

Борман какое сообщение я сегодня удалил?

Да я про ту тему, где обсуждали эту задачу.

В глобальной системе координат

[1 1; 0 0]{ux;uy}={f;0}

В локальной системе координат

[1 0; 0 0]{u;v}={f1;0}

А чего сила то не поворачивается ? А глобальной СК есть fx и fy.
Ссылка на сообщение
Поделиться на других сайтах

ИСПА, вы наверное яиц объелись. Я никого ничем не поливал. А если кто то вас и поливал, то я то тут каким боком? Это целиком ваша заслуга.

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

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




  • Сообщения

    • gudstartup
      а вы хоть станок проверяли по программе на изделии на точность прежде чем товарищей этих выгнать? если нет то грешите на самих себя! система в наших краях еще не распространенная поэтому и тем тут нет надо в поднебесную писать
    • AlexArt
      Ну допустим, ты и на другом ресурсе это опубликовал. А не коммуниздил. Но вот продвигать воровство от государства, ворующее из Вики, это верх мерзости.
    • maxx2000
      Ах, да. Фильтры выбора добавили. Теперь можно выбрать только то что видно на первом плане, а не вместе с тем что с обратной стороны детали. В общем надо обновляться. Как раз работёнка на прессформу нарисовалась 
    • maxx2000
      Причина того - Кроилово. Кроилово всегда приводит к попадалову. Месяц простоял сколько мильонов деревянных потеряли? Вопрос риторический. И ещё будет стоять. Как памятник человеческой глупости и жадности.
    • AlexKaz
      "9 июля 1968 года на мышах был проведен самый знаменитый эксперимент американского ученого-этолога Джона Кэлхуна «Вселенная-25». Суть опыта заключалась в создании идеальных условий, где мыши могли бы жить и размножаться, не ведая никаких забот, вдали от хищников и в отсутствие эпидемий и заболеваний. Для этих целей ученый построил специальный загон, куда были помещены четыре пары белых мышей (самцов и самок). В распоряжении мышей всегда была чистая вода и еда в изобилии, специальные гнезда, где можно обустроить себе жилище ― гнезд в загоне хватало для проживания нескольких тысяч мышей. Температура в загоне в среднем составляла около 20 ℃ и была комфортной для мышей. Животные не подвергались никаким влияниям извне и жили в идеальных условиях в свое удовольствие. А дальше началось самое интересное. На первом этапе эксперимента мыши хорошо размножались, вели активный образ жизни, охотно играли. На следующей фазе эксперимента мыши стали есть меньше, перестали наедаться до отвала. На третьей фазе эксперимента, когда в загоне были уже сотни мышей, произошло распределение социальных ролей, стала ярко выраженной иерархия, клановость. Появились так называемые отверженные ― молодые особи, которых другие, взрослые мыши сгоняли в центр загона, не давали им вести нормальный образ жизни, причиняли физический вред. В природе такое, наверное, было бы невозможно, ведь эти мыши-агрессоры просто не дожили бы до старости: их бы съели хищники. Но в загоне Кэлхуна хищников не было, и взрослые мыши начали попросту издеваться над молодняком. Образовались две большие группировки: самцы-одиночки и самки-одиночки. При этом самки-одиночки отказывались спариваться <с менее статусными многочисленными молодыми самцами и с оставшимися старыми статусными> и отвергали ухаживания самцов. У мышей стал проявляться тотальный индивидуализм, мыши не стремились создать семью. На последней, четвертой стадии мышиная популяция стала сокращаться. Появились самцы, которых сам Кэлхун назвал «красивыми» (англ. beautiful ones), из-за отсутствия ран и рубцов. <В оригинале: They never engaged in sexual approaches toward females, and they never engaged in fighting, and so they had no wound or scar tissue. Thus their pelage remained in excellent condition. - Дословный перевод: Они никогда не прибегали к сексуальным подходам к самкам, и они никогда не участвовали в боях, и поэтому у них не было ран или рубцовой ткани. Таким образом, их шерсть сохранилась в отличном состоянии.> Эти мыши не вступали в борьбу за самок и территорию, не проявляли активности к размножению и только питались, спали и чистили шёрстку. У мышей стали проявляться различные формы девиантного поведения, вспышки агрессии. Самки стали проявлять агрессию, защищать себя сами, стали умерщвлять своих детенышей, а затем окончательно отказались размножаться. На пике эксперимента в загоне одновременно проживало чуть более двух тыс. мышей. Еды и гнезд было достаточно для дальнейшего роста популяции, но через четыре года после начала эксперимента Кэлхун остановил свой опыт, потому что в загоне осталось чуть более сотни мышей, и все они уже вышли из репродуктивного возраста. По итогам эксперимента Кэлхун пришел к выводу, что достижение определенной плотности населения и заполнение социальных ролей в популяции приводит к распаду общества" https://physicsoflife.pl/dict/pic/calhoun/calhoun.. https://scientificrussia.ru/articles/utopiya-dlya-mys.. https://ru.wikipedia.org/wiki/Кэлхун,_Джон_(этолог)
    • gudstartup
      @Koels вот в чем дело пока ds609 это предупреждение поэтому F может и не появится если sv601 это значит ошибка. возможно при нагреве радиатора серво определяет это как предупреждение или ваш вентилятор крутиться медленнее чем оригинальный и серва думает что он встал хотяпри этом обычно на экране в строке состояния FAN.мигает больше у меня вариантов нет....  
    • ДОБРЯК
      Решите любым алгоритмом. Тогда будет конструктивный разговор. :=)
    • Fedor
      https://en.wikipedia.org/wiki/List_of_numerical_analysis_topics#Eigenvalue_algorithms     :) 
    • Юрий К.Ф.
      Добрый день. Не нашёл тут тему по стойке Китайско Китайской)) Lynuc N3ME. Видать мне так повезло с её наличием)) Приобрели 5-ти осевой Китаец. В б/у состоянии после удара по оси Z. Отремонтировали по механике, заменили батареи на драйверах, выставили лимиты. Всё Ок. Пригласили со стороны людей которые бы разобрались по операторской части. Те два выходных ковырялись, после сказали покажут расскажут, но за огромные деньги. Не сошлись. После месяц станок простоял, когда включили перестал реагировать на регулировку скорости шпинделя. То есть в режиме Jog, включаем обороты, которые стандартно 2140-2149 об/мин. При регулировке процетности не меняются (сама процентность показывает на мониторе). Так же при включении оборотов через команду M03S300 или другое значение, скорость так же показывает 2140-2149 об/мин. Грешить на тех товарищей с которыми не сошлись по деньгам для обучения, как то не хочется. Поковырялся в настройках шпинделя, вроде всё в норме. Проводку на шпинделе прозвонил, целая. В чём причина, не понятна. Кто нибудь сталкивался с подобным, или с подобной стойкой? Может подсказать варианты причины подобного?
    • ДОБРЯК
×
×
  • Создать...