2010-03-09 13:49:24 -07:00
|
|
|
// 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.
|
|
|
|
|
|
|
|
#include "runtime.h"
|
|
|
|
|
2010-06-18 16:46:00 -06:00
|
|
|
typedef struct Complex128 Complex128;
|
|
|
|
|
2010-03-09 13:49:24 -07:00
|
|
|
void
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
runtime·complex128div(Complex128 n, Complex128 d, Complex128 q)
|
2010-03-09 13:49:24 -07:00
|
|
|
{
|
2010-06-18 16:46:00 -06:00
|
|
|
int32 ninf, dinf, nnan, dnan;
|
2010-03-09 13:49:24 -07:00
|
|
|
float64 a, b, ratio, denom;
|
|
|
|
|
2010-06-18 16:46:00 -06:00
|
|
|
// Special cases as in C99.
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
ninf = runtime·isInf(n.real, 0) || runtime·isInf(n.imag, 0);
|
|
|
|
dinf = runtime·isInf(d.real, 0) || runtime·isInf(d.imag, 0);
|
2010-06-18 16:46:00 -06:00
|
|
|
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
nnan = !ninf && (runtime·isNaN(n.real) || runtime·isNaN(n.imag));
|
|
|
|
dnan = !dinf && (runtime·isNaN(d.real) || runtime·isNaN(d.imag));
|
2010-06-18 16:46:00 -06:00
|
|
|
|
|
|
|
if(nnan || dnan) {
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
q.real = runtime·NaN();
|
|
|
|
q.imag = runtime·NaN();
|
2010-06-18 16:46:00 -06:00
|
|
|
} else if(ninf && !dinf && !dnan) {
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
q.real = runtime·Inf(0);
|
|
|
|
q.imag = runtime·Inf(0);
|
2010-06-18 16:46:00 -06:00
|
|
|
} else if(!ninf && !nnan && dinf) {
|
|
|
|
q.real = 0;
|
|
|
|
q.imag = 0;
|
|
|
|
} else if(d.real == 0 && d.imag == 0) {
|
|
|
|
if(n.real == 0 && n.imag == 0) {
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
q.real = runtime·NaN();
|
|
|
|
q.imag = runtime·NaN();
|
2010-06-18 16:46:00 -06:00
|
|
|
} else {
|
runtime: ,s/[a-zA-Z0-9_]+/runtime·&/g, almost
Prefix all external symbols in runtime by runtime·,
to avoid conflicts with possible symbols of the same
name in linked-in C libraries. The obvious conflicts
are printf, malloc, and free, but hide everything to
avoid future pain.
The symbols left alone are:
** known to cgo **
_cgo_free
_cgo_malloc
libcgo_thread_start
initcgo
ncgocall
** known to linker **
_rt0_$GOARCH
_rt0_$GOARCH_$GOOS
text
etext
data
end
pclntab
epclntab
symtab
esymtab
** known to C compiler **
_divv
_modv
_div64by32
etc (arch specific)
Tested on darwin/386, darwin/amd64, linux/386, linux/amd64.
Built (but not tested) for freebsd/386, freebsd/amd64, linux/arm, windows/386.
R=r, PeterGo
CC=golang-dev
https://golang.org/cl/2899041
2010-11-04 12:00:19 -06:00
|
|
|
q.real = runtime·Inf(0);
|
|
|
|
q.imag = runtime·Inf(0);
|
2010-06-18 16:46:00 -06:00
|
|
|
}
|
2010-03-09 13:49:24 -07:00
|
|
|
} else {
|
2010-06-18 16:46:00 -06:00
|
|
|
// Standard complex arithmetic, factored to avoid unnecessary overflow.
|
|
|
|
a = d.real;
|
|
|
|
if(a < 0)
|
|
|
|
a = -a;
|
|
|
|
b = d.imag;
|
|
|
|
if(b < 0)
|
|
|
|
b = -b;
|
|
|
|
if(a <= b) {
|
|
|
|
ratio = d.real/d.imag;
|
|
|
|
denom = d.real*ratio + d.imag;
|
|
|
|
q.real = (n.real*ratio + n.imag) / denom;
|
|
|
|
q.imag = (n.imag*ratio - n.real) / denom;
|
|
|
|
} else {
|
|
|
|
ratio = d.imag/d.real;
|
|
|
|
denom = d.imag*ratio + d.real;
|
|
|
|
q.real = (n.imag*ratio + n.real) / denom;
|
|
|
|
q.imag = (n.imag - n.real*ratio) / denom;
|
|
|
|
}
|
2010-03-09 13:49:24 -07:00
|
|
|
}
|
2010-06-18 16:46:00 -06:00
|
|
|
FLUSH(&q);
|
2010-03-09 13:49:24 -07:00
|
|
|
}
|