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

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

    AMax = Max(Abs(V[i]))
Вычисляется сумма квадратов нормированных на этот максимум компонентов
  
    for all:
       nv = V[i] / AMax
       sum += nv * nv
Извлекается корень и возвращается нужный порядок

      hyp = AMax * Sqrt(sum)

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