1

Существует ли аналог алгоритма std::hypot() для вектора с 3+ элементами в стандартной библиотеке? Если нет, то можно предложить подобный по функционалу алгоритм?

6
  • А чем не нравится корень из суммы квадратов?
    – MBo
    23 янв 2019 в 17:15
  • @MBo переполнением
    – IIiiII
    23 янв 2019 в 17:17
  • @PavelGridin при больших значениях x и y сумма их квадратов теряет точность
    – IIiiII
    23 янв 2019 в 17:30
  • стёрлось сообщение, понятно, т. е. double не хватает 23 янв 2019 в 17:32
  • Может попробовать hypot из Intel IPP. Из описания: function flavor ippsHypot_64f_A53 guarantees 53 correctly rounded bits of significand, including the implied bit, with the maximum guaranteed error within 1 ulp 23 янв 2019 в 17:36

1 ответ 1

3

Раз могут возникнуть проблемы с переполнением, то можно использовать ту же технику, что и в hypot:

Выбирается максимальное по модулю значение из вектор

AMax = Max(Abs(V[i]))

Вычисляется сумма квадратов нормированных на этот максимум компонентов

for all:
   nv = V[i] / AMax
   sum += nv * nv

Извлекается корень и возвращается нужный порядок

  hyp = AMax * Sqrt(sum)

Можно также отсортировать компоненты и суммировать с минимальных по модулю либо использовать алгоритм Кэхэна.

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.