mirror of
https://github.com/golang/go
synced 2024-10-04 16:31:22 -06:00
13c2e62966
Also: * document special cases for Frexp and Ldexp * handle ±Inf in Ldexp * correctly return -0 on underflow in Ldexp * test special cases for Ldexp * test boundary cases for Frexp, Ilogb, Ldexp, and Logb R=rsc CC=golang-dev https://golang.org/cl/3676041
55 lines
1.3 KiB
Go
55 lines
1.3 KiB
Go
// Copyright 2010 The Go Authors. All rights reserved.
|
|
// Use of this source code is governed by a BSD-style
|
|
// license that can be found in the LICENSE file.
|
|
|
|
package math
|
|
|
|
// Logb(x) returns the binary exponent of x.
|
|
//
|
|
// Special cases are:
|
|
// Logb(±Inf) = +Inf
|
|
// Logb(0) = -Inf
|
|
// Logb(NaN) = NaN
|
|
func Logb(x float64) float64 {
|
|
// TODO(rsc): Remove manual inlining of IsNaN, IsInf
|
|
// when compiler does it for us
|
|
// special cases
|
|
switch {
|
|
case x == 0:
|
|
return Inf(-1)
|
|
case x < -MaxFloat64 || x > MaxFloat64: // IsInf(x, 0):
|
|
return Inf(1)
|
|
case x != x: // IsNaN(x):
|
|
return x
|
|
}
|
|
return float64(ilogb(x))
|
|
}
|
|
|
|
// Ilogb(x) returns the binary exponent of x as an integer.
|
|
//
|
|
// Special cases are:
|
|
// Ilogb(±Inf) = MaxInt32
|
|
// Ilogb(0) = MinInt32
|
|
// Ilogb(NaN) = MaxInt32
|
|
func Ilogb(x float64) int {
|
|
// TODO(rsc): Remove manual inlining of IsNaN, IsInf
|
|
// when compiler does it for us
|
|
// special cases
|
|
switch {
|
|
case x == 0:
|
|
return MinInt32
|
|
case x != x: // IsNaN(x):
|
|
return MaxInt32
|
|
case x < -MaxFloat64 || x > MaxFloat64: // IsInf(x, 0):
|
|
return MaxInt32
|
|
}
|
|
return ilogb(x)
|
|
}
|
|
|
|
// logb returns the binary exponent of x. It assumes x is finite and
|
|
// non-zero.
|
|
func ilogb(x float64) int {
|
|
x, exp := normalize(x)
|
|
return int((Float64bits(x)>>shift)&mask) - bias + exp
|
|
}
|