caxapa.ru :: Я очень извиняюсь, попав в толпу ПХПшников, а вот оборачивание матриц для МКЭ считается достойной задачей для серьёзных камней?

Крок (13.02.2017 21:14, просмотров: 332) ответил Evgeny_CD на В PHP все есть -> -->, а остальное - проблема идиотов, которые школоло на серьезные задачи поставило.

Я очень извиняюсь, попав в толпу ПХПшников, а вот оборачивание матриц для МКЭ считается достойной задачей для серьёзных камней?

Ответить

- Да, это типичная задача для FPU, которую можно считать серьезной, т.к. требует double64-арифметику из-за того, что данный процесс чувствителен к потере значащих цифр мантиссы. А с матрицами большой размерности задача становится еще труднее, т.к. Ксения(77 знак., 13.02.2017 21:39)
  - А что ж это за сигналы такие непредсказуемого размаха, что для них обязательна плавучка? Обычно высвобожденных битов порядка хватает на то, чтобы тупо создать в фикс. точке запас на все возможные случаи. - Николай Коровин(13.02.2017 23:00)
    - Нет. Обращение только в даблах. - Evgeny_CD(13.02.2017 23:05)
      - В принципе существуют еще и итерационные методы, когда обратная матрица вычисляется "халтурно" (с большой ошибкой), а затем ее умножают на исходную и смотрят на то, что получилось. Отклонение произведения от диагональной матрицы можно использовать Ксения(594 знак., 13.02.2017 23:45 - 23:50)
        
        Спасибо за пояснения. Но я все же про embedded, а там - недостаток ресурсов. Лобовой вариант с дабловым FPU выглядит быстрее всех, и по памяти разумно. Открытых исходников тоже много, можно брать и допиливать. - Evgeny_CD(14.02.2017 00:07)
        
        Кто бы спорил, но не я. До сих пор горюю из-за изгнания типа long double (float80) из компиляторов. Если уж заколодило из-за того, что 80 не является степенью двойки, то за 25 лет уж можно было сделать float128, благо что и стандарт такой есть. Но Ксения(92 знак., 14.02.2017 00:27)
        
        Купите машину Мир-2, там нет ограничений по разрядности - Крок(14.02.2017 08:51)
        
        Для меня, кстати, загадка, почему рядом с AVX, на худой конец AVX 512 не угнездился 128 битный float. Но, похоже, он нужен небольшому числу пользователей -> - Evgeny_CD(14.02.2017 00:44, ссылка)
        
        Потому что тяжело в реализации. Это AVX128-AVX256-AVX512 делаются простым клонированием AVX128 блока, т.к. "этажи" там работают независимо друг от друга. А умножение 128х128 даже выглядит страшно - где-то в квадрате должно быть сложнее, чем 64х64. Ксения(116 знак., 14.02.2017 01:39)
        
        5 лет назад => - Ксения(14.02.2017 01:43 - 01:53, ссылка)
        
        Как время быстро бежит :( - Evgeny_CD(14.02.2017 01:55)
      - А, обращение. - Николай Коровин(13.02.2017 23:28)
  - какие страсти Вы рассказываете! Люди попроще разбивают задачу на такие КЭ, чтоб матрица была более-менее обусловленной, а с недостатком точности борются вторым или третьим (крайне редко) прогоном. - Крок(13.02.2017 21:56)
    - Я не спец по КЭ, но если есть возможность минимизировать размер матриц и сделать их более обусловленными, то это большая удача. Однако в самом общем случае трудно оценить ранг матрицы, пока не начнешь ее обращать. - Ксения(13.02.2017 22:04)
      - Только в самом общем :-) - Крок(13.02.2017 22:08)
        
        Тогда вы сами себе ответили, что ваша задача несерьезна :) - Ксения(13.02.2017 22:25)
- У нас LPC3250 работал с кучкой матриц 120х120 в даблах. В том числе обращал их. Все получалось на пределе возможностей. - Evgeny_CD(13.02.2017 21:39)
- Зависит от числа КЭ. очевидно же.. :) - Хаос(13.02.2017 21:24, )

Средства и методы разработки