libgo/go/math/hypot.go

   1 // Copyright 2010 The Go Authors. All rights reserved.
   2 // Use of this source code is governed by a BSD-style
   3 // license that can be found in the LICENSE file.
   4
   5 package math
   6
   7 /*
   8         Hypot -- sqrt(p*p + q*q), but overflows only if the result does.
   9 */
  10
  11 // Hypot returns Sqrt(p*p + q*q), taking care to avoid
  12 // unnecessary overflow and underflow.
  13 //
  14 // Special cases are:
  15 //      Hypot(±Inf, q) = +Inf
  16 //      Hypot(p, ±Inf) = +Inf
  17 //      Hypot(NaN, q) = NaN
  18 //      Hypot(p, NaN) = NaN
  19 func Hypot(p, q float64) float64 {
  20         return hypot(p, q)
  21 }
  22
  23 func hypot(p, q float64) float64 {
  24         // special cases
  25         switch {
  26         case IsInf(p, 0) || IsInf(q, 0):
  27                 return Inf(1)
  28         case IsNaN(p) || IsNaN(q):
  29                 return NaN()
  30         }
  31         p, q = Abs(p), Abs(q)
  32         if p < q {
  33                 p, q = q, p
  34         }
  35         if p == 0 {
  36                 return 0
  37         }
  38         q = q / p
  39         return p * Sqrt(1+q*q)
  40 }