forked from jcosborn/qhmc
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ratfunsqrt.lua
83 lines (81 loc) · 1.7 KB
/
ratfunsqrt.lua
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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
-- approximate x^1/4 on [1e-6,100]
GRcoeffs = {
5.3299169928573225e+00,
-2.0882452575579878e-08,
-2.8850595278430710e-07,
-3.6166144644815880e-06,
-4.4935689953819157e-05,
-5.5771630080104194e-04,
-6.9230623110420357e-03,
-8.6124589723631892e-02,
-1.0891980155601380e+00,
-1.5604401397480311e+01,
-6.2553104306358671e+02
}
GRshifts = {
false,
1.0375540499266855e-06,
1.0655725087980704e-05,
8.3053651903913846e-05,
6.2597840618800252e-04,
4.6972757114629908e-03,
3.5231680272542594e-02,
2.6449044533148985e-01,
2.0001629302403892e+00,
1.5997607306850448e+01,
2.0936246640846309e+02
}
-- approximate x^-1/4 on [1e-6,100]
FAcoeffs = {
1.8762018270455441e-01,
1.4270895504485702e-05,
6.0972929291128314e-05,
2.7225515021822935e-04,
1.2311393174732779e-03,
5.5773967472347082e-03,
2.5276750652566956e-02,
1.1467604887817914e-01,
5.2430591709875252e-01,
2.5409070274230348e+00,
1.9398138448352206e+01
}
FAshifts = {
false,
4.7764053278251641e-07,
6.2509347855524823e-06,
4.9995927075791560e-05,
3.7808549142358830e-04,
2.8383545498377469e-03,
2.1288935575138838e-02,
1.5974991950435852e-01,
1.2040409748109777e+00,
9.3846264964921655e+00,
9.6380521098699489e+01
}
-- approximate x^-1/2 on [1e-6,100]
MDcoeffs = {
3.2123938669006216e-02,
7.2602402436934858e-04,
1.5304033363450992e-03,
4.0303538496748387e-03,
1.0995649647171764e-02,
3.0143710622588880e-02,
8.2697454896435818e-02,
2.2705444225309490e-01,
6.2671751975410206e-01,
1.8004316874890463e+00,
7.0567520057763060e+00,
}
MDshifts = {
false,
2.7698093990089493e-07,
4.6740638660503145e-06,
3.8352164081297480e-05,
2.9190241457477520e-04,
2.2002344497914885e-03,
1.6564127735668888e-02,
1.2473765234158118e-01,
9.4265100858885620e-01,
7.3175377111423652e+00,
7.0681491790812657e+01,
}