АЛУ и большая разрядность данных Опции

[bva_rzn]

Доброго времени суток.

ПЛИС Cyclone IV E.

Решаю задачу решения СЛАУ. Ax=B. Система 5 на 5.
В матрице A все числа целые (разрядность 52 бита).
Взял метод Холецкого без вычисления квадратного корня (разложение LDL^t). Получилось около 20 делений.
Если числа переводить во floaе, то кажное умножение займёт 33 такта.
Причём в проекте необходимо решить три таких системы.

И здесь возник вопрос: при такой разрядности данных (52 бита) есть ли смысл переходить на фиксированную точку?

Спасибо.

[iiv]

Решаю задачу решения СЛАУ. Ax=B. Система 5 на 5.

У Вас есть больше информации о матрице - например, если Вы взяли Холецкого, значит - положительно определенная, но вот об обусловленности, диагональном преобладании, структуре Вы ничего не сказали?

Как Вариант, взять метод сопряженного градиента, параллелится легко, сойдется гарантированно за 5 шагов, не вывалится как Холецкий если будет почти сингулярная, и если сильно постараться, то будет только одно деление и квадратный корень. Без последнего тоже можно обойтись, но нужно явно прорисовать весь алгоритм и тащить эту нормировку с шага на шаг, а в последнем - искать его итеративно (ИМХО, проще все-таки с квадратным корнем).

EDIT и самое главное, сопряженный градиент тривиально пишется в целочисленной арифметике - все векторы у Вас нормированные, то есть, при экспоненте 2**52 Вы можете не задумываясь соскалировать каждый вектор на эту же величину и взять целое. Есть куча более элегантных решений и для CG и без оного, но они обычно на 5*5 матрицах не применяются