1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
// compile
// Copyright 2016 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.
// This tickles (a version of) the PPC64 back end to
// emit a BVS instruction.
package foo
type Flag int
const (
Identity Flag = iota - 2 // H is the identity matrix; no rotation is needed.
Rescaling // H specifies rescaling.
)
type DrotmParams struct {
Flag
}
func Drotmg(d1, d2, x1, y1 float64) (p DrotmParams, rd1, rd2, rx1 float64) {
const (
gam = 4.0
gamsq = 16.0
rgamsq = 5e-8
)
if d1 < 0 {
p.Flag = Rescaling
return
}
for rd1 <= rgamsq || rd1 >= gamsq {
if rd1 <= rgamsq {
rd1 *= gam * gam
rx1 /= gam
} else {
rd1 /= gam * gam
rx1 *= gam
}
}
return
}
|
