caxapa.ru :: Да, это типичная задача для FPU, которую можно считать серьезной, т.к. требует double64-арифметику из-за того, что данный процесс чувствителен к потере значащих цифр мантиссы. А с матрицами большой размерности задача становится еще труднее, т.к.

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

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

Ответить

- А что ж это за сигналы такие непредсказуемого размаха, что для них обязательна плавучка? Обычно высвобожденных битов порядка хватает на то, чтобы тупо создать в фикс. точке запас на все возможные случаи. - Николай Коровин(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)

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