gcc/libgo/go/runtime/complex.go
Ian Lance Taylor f8d9fa9e80 libgo, compiler: Upgrade libgo to Go 1.4, except for runtime.
This upgrades all of libgo other than the runtime package to
the Go 1.4 release.  In Go 1.4 much of the runtime was
rewritten into Go.  Merging that code will take more time and
will not change the API, so I'm putting it off for now.

There are a few runtime changes anyhow, to accomodate other
packages that rely on minor modifications to the runtime
support.

The compiler changes slightly to add a one-bit flag to each
type descriptor kind that is stored directly in an interface,
which for gccgo is currently only pointer types.  Another
one-bit flag (gcprog) is reserved because it is used by the gc
compiler, but gccgo does not currently use it.

There is another error check in the compiler since I ran
across it during testing.

gotools/:
	* Makefile.am (go_cmd_go_files): Sort entries.  Add generate.go.
	* Makefile.in: Rebuild.

From-SVN: r219627
2015-01-15 00:27:56 +00:00

53 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 runtime
func complex128div(n complex128, d complex128) complex128 {
// Special cases as in C99.
ninf := real(n) == posinf || real(n) == neginf ||
imag(n) == posinf || imag(n) == neginf
dinf := real(d) == posinf || real(d) == neginf ||
imag(d) == posinf || imag(d) == neginf
nnan := !ninf && (real(n) != real(n) || imag(n) != imag(n))
dnan := !dinf && (real(d) != real(d) || imag(d) != imag(d))
switch {
case nnan || dnan:
return complex(nan, nan)
case ninf && !dinf:
return complex(posinf, posinf)
case !ninf && dinf:
return complex(0, 0)
case real(d) == 0 && imag(d) == 0:
if real(n) == 0 && imag(n) == 0 {
return complex(nan, nan)
} else {
return complex(posinf, posinf)
}
default:
// Standard complex arithmetic, factored to avoid unnecessary overflow.
a := real(d)
if a < 0 {
a = -a
}
b := imag(d)
if b < 0 {
b = -b
}
if a <= b {
ratio := real(d) / imag(d)
denom := real(d)*ratio + imag(d)
return complex((real(n)*ratio+imag(n))/denom,
(imag(n)*ratio-real(n))/denom)
} else {
ratio := imag(d) / real(d)
denom := imag(d)*ratio + real(d)
return complex((imag(n)*ratio+real(n))/denom,
(imag(n)-real(n)*ratio)/denom)
}
}
}