File Manager

X7ROOT File Manager

Current Path: /opt/golang/1.17.2/src/math

📁 ..
📄 abs.go(363 B)
📄 acos_s390x.s(3.73 KB)
📄 acosh.go(1.71 KB)
📄 acosh_s390x.s(4.32 KB)
📄 all_test.go(85.27 KB)
📄 arith_s390x.go(3.73 KB)
📄 arith_s390x_test.go(10.78 KB)
📄 asin.go(1.08 KB)
📄 asin_s390x.s(4.16 KB)
📄 asinh.go(1.92 KB)
📄 asinh_s390x.s(5.74 KB)
📄 atan.go(3.03 KB)
📄 atan2.go(1.52 KB)
📄 atan2_s390x.s(6.93 KB)
📄 atan_s390x.s(3.69 KB)
📄 atanh.go(1.99 KB)
📄 atanh_s390x.s(5.36 KB)
📁 big
📁 bits
📄 bits.go(1.87 KB)
📄 cbrt.go(2.31 KB)
📄 cbrt_s390x.s(4.89 KB)
📁 cmplx
📄 const.go(2.33 KB)
📄 const_test.go(1.29 KB)
📄 copysign.go(378 B)
📄 cosh_s390x.s(5.59 KB)
📄 dim.go(1.68 KB)
📄 dim_amd64.s(1.92 KB)
📄 dim_arm64.s(963 B)
📄 dim_asm.go(380 B)
📄 dim_noasm.go(450 B)
📄 dim_riscv64.s(1.16 KB)
📄 dim_s390x.s(1.97 KB)
📄 erf.go(11.5 KB)
📄 erf_s390x.s(8.5 KB)
📄 erfc_s390x.s(14.4 KB)
📄 erfinv.go(3.36 KB)
📄 example_test.go(3.66 KB)
📄 exp.go(5.37 KB)
📄 exp2_asm.go(268 B)
📄 exp2_noasm.go(301 B)
📄 exp_amd64.go(277 B)
📄 exp_amd64.s(4.24 KB)
📄 exp_arm64.s(5.36 KB)
📄 exp_asm.go(296 B)
📄 exp_noasm.go(333 B)
📄 exp_s390x.s(4.65 KB)
📄 expm1.go(7.9 KB)
📄 expm1_s390x.s(5.29 KB)
📄 export_s390x_test.go(732 B)
📄 export_test.go(357 B)
📄 floor.go(3.28 KB)
📄 floor_386.s(1.47 KB)
📄 floor_amd64.s(2 KB)
📄 floor_arm64.s(573 B)
📄 floor_asm.go(482 B)
📄 floor_noasm.go(589 B)
📄 floor_ppc64x.s(523 B)
📄 floor_s390x.s(579 B)
📄 floor_wasm.s(459 B)
📄 fma.go(4.46 KB)
📄 frexp.go(926 B)
📄 gamma.go(5.52 KB)
📄 huge_test.go(2.56 KB)
📄 hypot.go(845 B)
📄 hypot_386.s(1.81 KB)
📄 hypot_amd64.s(1.05 KB)
📄 hypot_asm.go(284 B)
📄 hypot_noasm.go(319 B)
📄 j0.go(13.6 KB)
📄 j1.go(13.3 KB)
📄 jn.go(7.17 KB)
📄 ldexp.go(1.05 KB)
📄 lgamma.go(11.02 KB)
📄 log.go(3.86 KB)
📄 log10.go(869 B)
📄 log10_s390x.s(4.73 KB)
📄 log1p.go(6.34 KB)
📄 log1p_s390x.s(5.15 KB)
📄 log_amd64.s(3.67 KB)
📄 log_asm.go(281 B)
📄 log_s390x.s(4.31 KB)
📄 log_stub.go(316 B)
📄 logb.go(1014 B)
📄 mod.go(900 B)
📄 modf.go(910 B)
📄 modf_arm64.s(447 B)
📄 modf_asm.go(322 B)
📄 modf_noasm.go(359 B)
📄 modf_ppc64x.s(440 B)
📄 nextafter.go(1.2 KB)
📄 pow.go(3.22 KB)
📄 pow10.go(1.24 KB)
📄 pow_s390x.s(16.27 KB)
📁 rand
📄 remainder.go(2.04 KB)
📄 signbit.go(302 B)
📄 sin.go(6.35 KB)
📄 sin_s390x.s(8.34 KB)
📄 sincos.go(1.75 KB)
📄 sinh.go(1.69 KB)
📄 sinh_s390x.s(5.98 KB)
📄 sqrt.go(4.9 KB)
📄 sqrt_386.s(304 B)
📄 sqrt_amd64.s(334 B)
📄 sqrt_arm.s(529 B)
📄 sqrt_arm64.s(310 B)
📄 sqrt_asm.go(416 B)
📄 sqrt_mipsx.s(409 B)
📄 sqrt_noasm.go(469 B)
📄 sqrt_ppc64x.s(362 B)
📄 sqrt_riscv64.s(308 B)
📄 sqrt_s390x.s(309 B)
📄 sqrt_wasm.s(273 B)
📄 stubs.go(2.59 KB)
📄 stubs_s390x.s(12.38 KB)
📄 tan.go(3.67 KB)
📄 tan_s390x.s(2.73 KB)
📄 tanh.go(2.65 KB)
📄 tanh_s390x.s(4.57 KB)
📄 trig_reduce.go(3.33 KB)
📄 unsafe.go(1.27 KB)

Editing: pow.go

// Copyright 2009 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

func isOddInt(x float64) bool {
	xi, xf := Modf(x)
	return xf == 0 && int64(xi)&1 == 1
}

// Special cases taken from FreeBSD's /usr/src/lib/msun/src/e_pow.c
// updated by IEEE Std. 754-2008 "Section 9.2.1 Special values".

// Pow returns x**y, the base-x exponential of y.
//
// Special cases are (in order):
//	Pow(x, ±0) = 1 for any x
//	Pow(1, y) = 1 for any y
//	Pow(x, 1) = x for any x
//	Pow(NaN, y) = NaN
//	Pow(x, NaN) = NaN
//	Pow(±0, y) = ±Inf for y an odd integer < 0
//	Pow(±0, -Inf) = +Inf
//	Pow(±0, +Inf) = +0
//	Pow(±0, y) = +Inf for finite y < 0 and not an odd integer
//	Pow(±0, y) = ±0 for y an odd integer > 0
//	Pow(±0, y) = +0 for finite y > 0 and not an odd integer
//	Pow(-1, ±Inf) = 1
//	Pow(x, +Inf) = +Inf for |x| > 1
//	Pow(x, -Inf) = +0 for |x| > 1
//	Pow(x, +Inf) = +0 for |x| < 1
//	Pow(x, -Inf) = +Inf for |x| < 1
//	Pow(+Inf, y) = +Inf for y > 0
//	Pow(+Inf, y) = +0 for y < 0
//	Pow(-Inf, y) = Pow(-0, -y)
//	Pow(x, y) = NaN for finite x < 0 and finite non-integer y
func Pow(x, y float64) float64 {
	if haveArchPow {
		return archPow(x, y)
	}
	return pow(x, y)
}

func pow(x, y float64) float64 {
	switch {
	case y == 0 || x == 1:
		return 1
	case y == 1:
		return x
	case IsNaN(x) || IsNaN(y):
		return NaN()
	case x == 0:
		switch {
		case y < 0:
			if isOddInt(y) {
				return Copysign(Inf(1), x)
			}
			return Inf(1)
		case y > 0:
			if isOddInt(y) {
				return x
			}
			return 0
		}
	case IsInf(y, 0):
		switch {
		case x == -1:
			return 1
		case (Abs(x) < 1) == IsInf(y, 1):
			return 0
		default:
			return Inf(1)
		}
	case IsInf(x, 0):
		if IsInf(x, -1) {
			return Pow(1/x, -y) // Pow(-0, -y)
		}
		switch {
		case y < 0:
			return 0
		case y > 0:
			return Inf(1)
		}
	case y == 0.5:
		return Sqrt(x)
	case y == -0.5:
		return 1 / Sqrt(x)
	}

yi, yf := Modf(Abs(y))
	if yf != 0 && x < 0 {
		return NaN()
	}
	if yi >= 1<<63 {
		// yi is a large even int that will lead to overflow (or underflow to 0)
		// for all x except -1 (x == 1 was handled earlier)
		switch {
		case x == -1:
			return 1
		case (Abs(x) < 1) == (y > 0):
			return 0
		default:
			return Inf(1)
		}
	}

// ans = a1 * 2**ae (= 1 for now).
	a1 := 1.0
	ae := 0

// ans *= x**yf
	if yf != 0 {
		if yf > 0.5 {
			yf--
			yi++
		}
		a1 = Exp(yf * Log(x))
	}

// ans *= x**yi
	// by multiplying in successive squarings
	// of x according to bits of yi.
	// accumulate powers of two into exp.
	x1, xe := Frexp(x)
	for i := int64(yi); i != 0; i >>= 1 {
		if xe < -1<<12 || 1<<12 < xe {
			// catch xe before it overflows the left shift below
			// Since i !=0 it has at least one bit still set, so ae will accumulate xe
			// on at least one more iteration, ae += xe is a lower bound on ae
			// the lower bound on ae exceeds the size of a float64 exp
			// so the final call to Ldexp will produce under/overflow (0/Inf)
			ae += xe
			break
		}
		if i&1 == 1 {
			a1 *= x1
			ae += xe
		}
		x1 *= x1
		xe <<= 1
		if x1 < .5 {
			x1 += x1
			xe--
		}
	}

// ans = a1*2**ae
	// if y < 0 { ans = 1 / ans }
	// but in the opposite order
	if y < 0 {
		a1 = 1 / a1
		ae = -ae
	}
	return Ldexp(a1, ae)
}

X7ROOT File Manager

Editing: pow.go

Upload File

Create Folder