Существует ли аналог алгоритма std::hypot()
для вектора с 3+ элементами в стандартной библиотеке? Если нет, то можно предложить подобный по функционалу алгоритм?
-
А чем не нравится корень из суммы квадратов?– MBo23 янв 2019 в 17:15
-
@MBo переполнением– IIiiII23 янв 2019 в 17:17
-
@PavelGridin при больших значениях x и y сумма их квадратов теряет точность– IIiiII23 янв 2019 в 17:30
-
стёрлось сообщение, понятно, т. е. double не хватает– Pavel Gridin23 янв 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– Pavel Gridin23 янв 2019 в 17:36
1 ответ
Раз могут возникнуть проблемы с переполнением, то можно использовать ту же технику, что и в hypot
:
Выбирается максимальное по модулю значение из вектор
AMax = Max(Abs(V[i]))
Вычисляется сумма квадратов нормированных на этот максимум компонентов
for all:
nv = V[i] / AMax
sum += nv * nv
Извлекается корень и возвращается нужный порядок
hyp = AMax * Sqrt(sum)
Можно также отсортировать компоненты и суммировать с минимальных по модулю либо использовать алгоритм Кэхэна.