gsl4r 0.0.1 → 0.0.2
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- data/README +18 -7
- data/Rakefile +7 -3
- data/changelog +17 -0
- data/lib/gsl4r/block.rb +59 -0
- data/lib/gsl4r/complex.rb +16 -54
- data/lib/gsl4r/const.rb +14 -0
- data/lib/gsl4r/const_cgs.rb +102 -0
- data/lib/gsl4r/const_cgsm.rb +108 -0
- data/lib/gsl4r/const_mks.rb +112 -0
- data/lib/gsl4r/const_mksa.rb +112 -0
- data/lib/gsl4r/const_num.rb +29 -0
- data/lib/gsl4r/platform.rb +3 -2
- data/lib/gsl4r/util.rb +70 -4
- data/lib/gsl4r/vector.rb +17 -12
- data/lib/t.rb +10 -0
- data/test/README +8 -0
- data/test/complex_test.rb +257 -257
- data/test/gsl_complex_tests_gen +0 -0
- data/test/gsl_complex_tests_gen.c +476 -476
- metadata +13 -4
- data/lib/gsl4r/f +0 -174
- data/test/complex_tests.rb +0 -28
@@ -5,101 +5,39 @@
|
|
5
5
|
int main( int argc, char **argv )
|
6
6
|
{
|
7
7
|
puts("$: << File.join('..','lib')\nrequire 'test/unit'\nrequire 'test/unit/autorunner'\nrequire 'gsl4r/complex'\ninclude GSL4r::Complex\nclass ComplexTests < Test::Unit::TestCase\n EPSILON = 5.0e-15");
|
8
|
-
/*
|
9
|
-
gsl_complex
|
10
|
-
GSL_SET_COMPLEX(&
|
11
|
-
|
12
|
-
|
13
|
-
|
14
|
-
puts("
|
15
|
-
puts("
|
16
|
-
puts("
|
17
|
-
puts(" r_r2 = GSL_Complex.new");
|
18
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r63),GSL_IMAG(c_r63));
|
19
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
20
|
-
puts("end");
|
21
|
-
/* c_test_gsl_complex_pow_real */
|
22
|
-
gsl_complex v37;
|
23
|
-
GSL_SET_COMPLEX(&v37, 2.0, 2.0);
|
24
|
-
double v38 = (double)2.0;
|
25
|
-
gsl_complex c_r36;
|
26
|
-
c_r36 = gsl_complex_pow_real(v37,v38);
|
27
|
-
puts("def test_gsl_complex_pow_real()");
|
28
|
-
puts(" v37 = GSL_Complex.create");
|
29
|
-
puts(" v37.set(2.0,2.0)");
|
8
|
+
/* c_test_gsl_complex_add_real */
|
9
|
+
gsl_complex v12;
|
10
|
+
GSL_SET_COMPLEX(&v12, 2.0, 2.0);
|
11
|
+
double v13 = (double)2.0;
|
12
|
+
gsl_complex c_r11;
|
13
|
+
c_r11 = gsl_complex_add_real(v12,v13);
|
14
|
+
puts("def test_gsl_complex_add_real()");
|
15
|
+
puts(" v12 = GSL_Complex.create");
|
16
|
+
puts(" v12.set(2.0,2.0)");
|
30
17
|
puts("");
|
31
|
-
puts("
|
32
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
33
|
-
puts(" r_r2 = GSL_Complex.new");
|
34
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r36),GSL_IMAG(c_r36));
|
35
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
36
|
-
puts("end");
|
37
|
-
/* c_test_gsl_complex_sinh */
|
38
|
-
gsl_complex v60;
|
39
|
-
GSL_SET_COMPLEX(&v60, 2.0, 2.0);
|
40
|
-
gsl_complex c_r59;
|
41
|
-
c_r59 = gsl_complex_sinh(v60);
|
42
|
-
puts("def test_gsl_complex_sinh()");
|
43
|
-
puts(" v60 = GSL_Complex.create");
|
44
|
-
puts(" v60.set(2.0,2.0)");
|
45
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sinh(v60)");
|
46
|
-
puts(" r_r2 = GSL_Complex.new");
|
47
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r59),GSL_IMAG(c_r59));
|
48
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
49
|
-
puts("end");
|
50
|
-
/* c_test_gsl_complex_sqrt */
|
51
|
-
gsl_complex v31;
|
52
|
-
GSL_SET_COMPLEX(&v31, 2.0, 2.0);
|
53
|
-
gsl_complex c_r30;
|
54
|
-
c_r30 = gsl_complex_sqrt(v31);
|
55
|
-
puts("def test_gsl_complex_sqrt()");
|
56
|
-
puts(" v31 = GSL_Complex.create");
|
57
|
-
puts(" v31.set(2.0,2.0)");
|
58
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sqrt(v31)");
|
59
|
-
puts(" r_r2 = GSL_Complex.new");
|
60
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r30),GSL_IMAG(c_r30));
|
61
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
62
|
-
puts("end");
|
63
|
-
/* c_test_gsl_complex_arccsc */
|
64
|
-
gsl_complex v56;
|
65
|
-
GSL_SET_COMPLEX(&v56, 2.0, 2.0);
|
66
|
-
gsl_complex c_r55;
|
67
|
-
c_r55 = gsl_complex_arccsc(v56);
|
68
|
-
puts("def test_gsl_complex_arccsc()");
|
69
|
-
puts(" v56 = GSL_Complex.create");
|
70
|
-
puts(" v56.set(2.0,2.0)");
|
71
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsc(v56)");
|
18
|
+
puts(" v13 = 2.0");
|
19
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_add_real(v12,v13)");
|
72
20
|
puts(" r_r2 = GSL_Complex.new");
|
73
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
21
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r11),GSL_IMAG(c_r11));
|
74
22
|
printf(" assert r_r1.equals(r_r2)\n");
|
75
23
|
puts("end");
|
76
|
-
/*
|
77
|
-
gsl_complex
|
78
|
-
GSL_SET_COMPLEX(&
|
79
|
-
|
80
|
-
|
81
|
-
|
82
|
-
|
83
|
-
puts("
|
84
|
-
puts("
|
85
|
-
puts("");
|
86
|
-
puts("
|
87
|
-
puts("
|
24
|
+
/* c_test_gsl_complex_pow */
|
25
|
+
gsl_complex v35;
|
26
|
+
GSL_SET_COMPLEX(&v35, 2.0, 2.0);
|
27
|
+
gsl_complex v36;
|
28
|
+
GSL_SET_COMPLEX(&v36, 2.0, 2.0);
|
29
|
+
gsl_complex c_r34;
|
30
|
+
c_r34 = gsl_complex_pow(v35,v36);
|
31
|
+
puts("def test_gsl_complex_pow()");
|
32
|
+
puts(" v35 = GSL_Complex.create");
|
33
|
+
puts(" v35.set(2.0,2.0)");
|
34
|
+
puts(" v36 = GSL_Complex.create");
|
35
|
+
puts(" v36.set(2.0,2.0)");
|
36
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_pow(v35,v36)");
|
88
37
|
puts(" r_r2 = GSL_Complex.new");
|
89
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
38
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r34),GSL_IMAG(c_r34));
|
90
39
|
printf(" assert r_r1.equals(r_r2)\n");
|
91
40
|
puts("end");
|
92
|
-
/* c_test_gsl_complex_abs2 */
|
93
|
-
gsl_complex v3;
|
94
|
-
GSL_SET_COMPLEX(&v3, 2.0, 2.0);
|
95
|
-
double c_r2;
|
96
|
-
c_r2 = gsl_complex_abs2(v3);
|
97
|
-
puts("def test_gsl_complex_abs2()");
|
98
|
-
puts(" v3 = GSL_Complex.create");
|
99
|
-
puts(" v3.set(2.0,2.0)");
|
100
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs2(v3)");
|
101
|
-
printf(" assert_in_delta r_r1, %.15g, EPSILON\n", c_r2);
|
102
|
-
puts("end");
|
103
41
|
/* c_test_gsl_complex_arccos */
|
104
42
|
gsl_complex v52;
|
105
43
|
GSL_SET_COMPLEX(&v52, 2.0, 2.0);
|
@@ -113,17 +51,17 @@ double v27 = (double)2.0;
|
|
113
51
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r51),GSL_IMAG(c_r51));
|
114
52
|
printf(" assert r_r1.equals(r_r2)\n");
|
115
53
|
puts("end");
|
116
|
-
/*
|
117
|
-
gsl_complex
|
118
|
-
GSL_SET_COMPLEX(&
|
119
|
-
gsl_complex
|
120
|
-
|
121
|
-
puts("def
|
122
|
-
puts("
|
123
|
-
puts("
|
124
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
54
|
+
/* c_test_gsl_complex_arcsinh */
|
55
|
+
gsl_complex v66;
|
56
|
+
GSL_SET_COMPLEX(&v66, 2.0, 2.0);
|
57
|
+
gsl_complex c_r65;
|
58
|
+
c_r65 = gsl_complex_arcsinh(v66);
|
59
|
+
puts("def test_gsl_complex_arcsinh()");
|
60
|
+
puts(" v66 = GSL_Complex.create");
|
61
|
+
puts(" v66.set(2.0,2.0)");
|
62
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsinh(v66)");
|
125
63
|
puts(" r_r2 = GSL_Complex.new");
|
126
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
64
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r65),GSL_IMAG(c_r65));
|
127
65
|
printf(" assert r_r1.equals(r_r2)\n");
|
128
66
|
puts("end");
|
129
67
|
/* c_test_gsl_complex_div_real */
|
@@ -142,17 +80,29 @@ double v19 = (double)2.0;
|
|
142
80
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r17),GSL_IMAG(c_r17));
|
143
81
|
printf(" assert r_r1.equals(r_r2)\n");
|
144
82
|
puts("end");
|
145
|
-
/*
|
146
|
-
gsl_complex
|
147
|
-
GSL_SET_COMPLEX(&
|
148
|
-
gsl_complex
|
149
|
-
|
150
|
-
puts("def
|
151
|
-
puts("
|
152
|
-
puts("
|
153
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
83
|
+
/* c_test_gsl_complex_log */
|
84
|
+
gsl_complex v40;
|
85
|
+
GSL_SET_COMPLEX(&v40, 2.0, 2.0);
|
86
|
+
gsl_complex c_r39;
|
87
|
+
c_r39 = gsl_complex_log(v40);
|
88
|
+
puts("def test_gsl_complex_log()");
|
89
|
+
puts(" v40 = GSL_Complex.create");
|
90
|
+
puts(" v40.set(2.0,2.0)");
|
91
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log(v40)");
|
154
92
|
puts(" r_r2 = GSL_Complex.new");
|
155
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
93
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r39),GSL_IMAG(c_r39));
|
94
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
95
|
+
puts("end");
|
96
|
+
/* c_test_gsl_complex_arcsec_real */
|
97
|
+
double v55 = (double)2.0;
|
98
|
+
gsl_complex c_r54;
|
99
|
+
c_r54 = gsl_complex_arcsec_real(v55);
|
100
|
+
puts("def test_gsl_complex_arcsec_real()");
|
101
|
+
puts("");
|
102
|
+
puts(" v55 = 2.0");
|
103
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsec_real(v55)");
|
104
|
+
puts(" r_r2 = GSL_Complex.new");
|
105
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r54),GSL_IMAG(c_r54));
|
156
106
|
printf(" assert r_r1.equals(r_r2)\n");
|
157
107
|
puts("end");
|
158
108
|
/* c_test_gsl_complex_arcsech */
|
@@ -168,21 +118,20 @@ double v19 = (double)2.0;
|
|
168
118
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r68),GSL_IMAG(c_r68));
|
169
119
|
printf(" assert r_r1.equals(r_r2)\n");
|
170
120
|
puts("end");
|
171
|
-
/*
|
172
|
-
gsl_complex
|
173
|
-
GSL_SET_COMPLEX(&
|
174
|
-
|
175
|
-
|
176
|
-
|
177
|
-
|
178
|
-
puts("
|
179
|
-
puts("
|
180
|
-
puts("
|
181
|
-
puts("
|
182
|
-
puts("
|
183
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div(v10,v11)");
|
121
|
+
/* c_test_gsl_complex_mul_imag */
|
122
|
+
gsl_complex v24;
|
123
|
+
GSL_SET_COMPLEX(&v24, 2.0, 2.0);
|
124
|
+
double v25 = (double)2.0;
|
125
|
+
gsl_complex c_r23;
|
126
|
+
c_r23 = gsl_complex_mul_imag(v24,v25);
|
127
|
+
puts("def test_gsl_complex_mul_imag()");
|
128
|
+
puts(" v24 = GSL_Complex.create");
|
129
|
+
puts(" v24.set(2.0,2.0)");
|
130
|
+
puts("");
|
131
|
+
puts(" v25 = 2.0");
|
132
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_imag(v24,v25)");
|
184
133
|
puts(" r_r2 = GSL_Complex.new");
|
185
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
134
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r23),GSL_IMAG(c_r23));
|
186
135
|
printf(" assert r_r1.equals(r_r2)\n");
|
187
136
|
puts("end");
|
188
137
|
/* c_test_gsl_complex_sin */
|
@@ -198,49 +147,102 @@ double v19 = (double)2.0;
|
|
198
147
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r43),GSL_IMAG(c_r43));
|
199
148
|
printf(" assert r_r1.equals(r_r2)\n");
|
200
149
|
puts("end");
|
201
|
-
/*
|
202
|
-
gsl_complex
|
203
|
-
GSL_SET_COMPLEX(&
|
204
|
-
gsl_complex
|
205
|
-
|
206
|
-
puts("def
|
207
|
-
puts("
|
208
|
-
puts("
|
209
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
150
|
+
/* c_test_gsl_complex_arctan */
|
151
|
+
gsl_complex v58;
|
152
|
+
GSL_SET_COMPLEX(&v58, 2.0, 2.0);
|
153
|
+
gsl_complex c_r57;
|
154
|
+
c_r57 = gsl_complex_arctan(v58);
|
155
|
+
puts("def test_gsl_complex_arctan()");
|
156
|
+
puts(" v58 = GSL_Complex.create");
|
157
|
+
puts(" v58.set(2.0,2.0)");
|
158
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctan(v58)");
|
210
159
|
puts(" r_r2 = GSL_Complex.new");
|
211
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
160
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r57),GSL_IMAG(c_r57));
|
212
161
|
printf(" assert r_r1.equals(r_r2)\n");
|
213
162
|
puts("end");
|
214
|
-
/*
|
215
|
-
|
216
|
-
|
217
|
-
|
218
|
-
|
219
|
-
puts("
|
220
|
-
puts("
|
221
|
-
puts("
|
222
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_exp(v39)");
|
163
|
+
/* c_test_gsl_complex_arctanh_real */
|
164
|
+
double v72 = (double)2.0;
|
165
|
+
gsl_complex c_r71;
|
166
|
+
c_r71 = gsl_complex_arctanh_real(v72);
|
167
|
+
puts("def test_gsl_complex_arctanh_real()");
|
168
|
+
puts("");
|
169
|
+
puts(" v72 = 2.0");
|
170
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctanh_real(v72)");
|
223
171
|
puts(" r_r2 = GSL_Complex.new");
|
224
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
172
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r71),GSL_IMAG(c_r71));
|
225
173
|
printf(" assert r_r1.equals(r_r2)\n");
|
226
174
|
puts("end");
|
227
|
-
/*
|
228
|
-
gsl_complex
|
229
|
-
GSL_SET_COMPLEX(&
|
230
|
-
gsl_complex
|
231
|
-
|
232
|
-
puts("def
|
233
|
-
puts("
|
234
|
-
puts("
|
235
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
175
|
+
/* c_test_gsl_complex_inverse */
|
176
|
+
gsl_complex v29;
|
177
|
+
GSL_SET_COMPLEX(&v29, 2.0, 2.0);
|
178
|
+
gsl_complex c_r28;
|
179
|
+
c_r28 = gsl_complex_inverse(v29);
|
180
|
+
puts("def test_gsl_complex_inverse()");
|
181
|
+
puts(" v29 = GSL_Complex.create");
|
182
|
+
puts(" v29.set(2.0,2.0)");
|
183
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_inverse(v29)");
|
236
184
|
puts(" r_r2 = GSL_Complex.new");
|
237
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
185
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r28),GSL_IMAG(c_r28));
|
238
186
|
printf(" assert r_r1.equals(r_r2)\n");
|
239
187
|
puts("end");
|
240
|
-
/*
|
241
|
-
|
242
|
-
|
243
|
-
|
188
|
+
/* c_test_gsl_complex_csc */
|
189
|
+
gsl_complex v47;
|
190
|
+
GSL_SET_COMPLEX(&v47, 2.0, 2.0);
|
191
|
+
gsl_complex c_r46;
|
192
|
+
c_r46 = gsl_complex_csc(v47);
|
193
|
+
puts("def test_gsl_complex_csc()");
|
194
|
+
puts(" v47 = GSL_Complex.create");
|
195
|
+
puts(" v47.set(2.0,2.0)");
|
196
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_csc(v47)");
|
197
|
+
puts(" r_r2 = GSL_Complex.new");
|
198
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r46),GSL_IMAG(c_r46));
|
199
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
200
|
+
puts("end");
|
201
|
+
/* c_test_gsl_complex_arg */
|
202
|
+
gsl_complex v1;
|
203
|
+
GSL_SET_COMPLEX(&v1, 2.0, 2.0);
|
204
|
+
double c_r0;
|
205
|
+
c_r0 = gsl_complex_arg(v1);
|
206
|
+
puts("def test_gsl_complex_arg()");
|
207
|
+
puts(" v1 = GSL_Complex.create");
|
208
|
+
puts(" v1.set(2.0,2.0)");
|
209
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arg(v1)");
|
210
|
+
printf(" assert_in_delta r_r1, %.15g, EPSILON\n", c_r0);
|
211
|
+
puts("end");
|
212
|
+
/* c_test_gsl_complex_cosh */
|
213
|
+
gsl_complex v61;
|
214
|
+
GSL_SET_COMPLEX(&v61, 2.0, 2.0);
|
215
|
+
gsl_complex c_r60;
|
216
|
+
c_r60 = gsl_complex_cosh(v61);
|
217
|
+
puts("def test_gsl_complex_cosh()");
|
218
|
+
puts(" v61 = GSL_Complex.create");
|
219
|
+
puts(" v61.set(2.0,2.0)");
|
220
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_cosh(v61)");
|
221
|
+
puts(" r_r2 = GSL_Complex.new");
|
222
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r60),GSL_IMAG(c_r60));
|
223
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
224
|
+
puts("end");
|
225
|
+
/* c_test_gsl_complex_mul */
|
226
|
+
gsl_complex v8;
|
227
|
+
GSL_SET_COMPLEX(&v8, 2.0, 2.0);
|
228
|
+
gsl_complex v9;
|
229
|
+
GSL_SET_COMPLEX(&v9, 2.0, 2.0);
|
230
|
+
gsl_complex c_r7;
|
231
|
+
c_r7 = gsl_complex_mul(v8,v9);
|
232
|
+
puts("def test_gsl_complex_mul()");
|
233
|
+
puts(" v8 = GSL_Complex.create");
|
234
|
+
puts(" v8.set(2.0,2.0)");
|
235
|
+
puts(" v9 = GSL_Complex.create");
|
236
|
+
puts(" v9.set(2.0,2.0)");
|
237
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul(v8,v9)");
|
238
|
+
puts(" r_r2 = GSL_Complex.new");
|
239
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r7),GSL_IMAG(c_r7));
|
240
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
241
|
+
puts("end");
|
242
|
+
/* c_test_gsl_complex_sqrt_real */
|
243
|
+
double v32 = (double)2.0;
|
244
|
+
gsl_complex c_r31;
|
245
|
+
c_r31 = gsl_complex_sqrt_real(v32);
|
244
246
|
puts("def test_gsl_complex_sqrt_real()");
|
245
247
|
puts("");
|
246
248
|
puts(" v32 = 2.0");
|
@@ -249,29 +251,62 @@ double v32 = (double)2.0;
|
|
249
251
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r31),GSL_IMAG(c_r31));
|
250
252
|
printf(" assert r_r1.equals(r_r2)\n");
|
251
253
|
puts("end");
|
252
|
-
/*
|
253
|
-
|
254
|
-
|
255
|
-
|
256
|
-
|
254
|
+
/* c_test_gsl_complex_arcsin */
|
255
|
+
gsl_complex v50;
|
256
|
+
GSL_SET_COMPLEX(&v50, 2.0, 2.0);
|
257
|
+
gsl_complex c_r49;
|
258
|
+
c_r49 = gsl_complex_arcsin(v50);
|
259
|
+
puts("def test_gsl_complex_arcsin()");
|
260
|
+
puts(" v50 = GSL_Complex.create");
|
261
|
+
puts(" v50.set(2.0,2.0)");
|
262
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsin(v50)");
|
263
|
+
puts(" r_r2 = GSL_Complex.new");
|
264
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r49),GSL_IMAG(c_r49));
|
265
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
266
|
+
puts("end");
|
267
|
+
/* c_test_gsl_complex_tanh */
|
268
|
+
gsl_complex v64;
|
269
|
+
GSL_SET_COMPLEX(&v64, 2.0, 2.0);
|
270
|
+
gsl_complex c_r63;
|
271
|
+
c_r63 = gsl_complex_tanh(v64);
|
272
|
+
puts("def test_gsl_complex_tanh()");
|
273
|
+
puts(" v64 = GSL_Complex.create");
|
274
|
+
puts(" v64.set(2.0,2.0)");
|
275
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_tanh(v64)");
|
276
|
+
puts(" r_r2 = GSL_Complex.new");
|
277
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r63),GSL_IMAG(c_r63));
|
278
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
279
|
+
puts("end");
|
280
|
+
/* c_test_gsl_complex_sub_real */
|
281
|
+
gsl_complex v14;
|
282
|
+
GSL_SET_COMPLEX(&v14, 2.0, 2.0);
|
283
|
+
double v15 = (double)2.0;
|
284
|
+
gsl_complex c_r13;
|
285
|
+
c_r13 = gsl_complex_sub_real(v14,v15);
|
286
|
+
puts("def test_gsl_complex_sub_real()");
|
287
|
+
puts(" v14 = GSL_Complex.create");
|
288
|
+
puts(" v14.set(2.0,2.0)");
|
257
289
|
puts("");
|
258
|
-
puts("
|
259
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
290
|
+
puts(" v15 = 2.0");
|
291
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sub_real(v14,v15)");
|
260
292
|
puts(" r_r2 = GSL_Complex.new");
|
261
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
293
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r13),GSL_IMAG(c_r13));
|
262
294
|
printf(" assert r_r1.equals(r_r2)\n");
|
263
295
|
puts("end");
|
264
|
-
/*
|
265
|
-
gsl_complex
|
266
|
-
GSL_SET_COMPLEX(&
|
267
|
-
|
268
|
-
|
269
|
-
|
270
|
-
puts("
|
271
|
-
puts("
|
272
|
-
puts("
|
296
|
+
/* c_test_gsl_complex_pow_real */
|
297
|
+
gsl_complex v37;
|
298
|
+
GSL_SET_COMPLEX(&v37, 2.0, 2.0);
|
299
|
+
double v38 = (double)2.0;
|
300
|
+
gsl_complex c_r36;
|
301
|
+
c_r36 = gsl_complex_pow_real(v37,v38);
|
302
|
+
puts("def test_gsl_complex_pow_real()");
|
303
|
+
puts(" v37 = GSL_Complex.create");
|
304
|
+
puts(" v37.set(2.0,2.0)");
|
305
|
+
puts("");
|
306
|
+
puts(" v38 = 2.0");
|
307
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_pow_real(v37,v38)");
|
273
308
|
puts(" r_r2 = GSL_Complex.new");
|
274
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
309
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r36),GSL_IMAG(c_r36));
|
275
310
|
printf(" assert r_r1.equals(r_r2)\n");
|
276
311
|
puts("end");
|
277
312
|
/* c_test_gsl_complex_arccos_real */
|
@@ -286,6 +321,19 @@ double v53 = (double)2.0;
|
|
286
321
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r52),GSL_IMAG(c_r52));
|
287
322
|
printf(" assert r_r1.equals(r_r2)\n");
|
288
323
|
puts("end");
|
324
|
+
/* c_test_gsl_complex_arccosh */
|
325
|
+
gsl_complex v67;
|
326
|
+
GSL_SET_COMPLEX(&v67, 2.0, 2.0);
|
327
|
+
gsl_complex c_r66;
|
328
|
+
c_r66 = gsl_complex_arccosh(v67);
|
329
|
+
puts("def test_gsl_complex_arccosh()");
|
330
|
+
puts(" v67 = GSL_Complex.create");
|
331
|
+
puts(" v67.set(2.0,2.0)");
|
332
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccosh(v67)");
|
333
|
+
puts(" r_r2 = GSL_Complex.new");
|
334
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r66),GSL_IMAG(c_r66));
|
335
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
336
|
+
puts("end");
|
289
337
|
/* c_test_gsl_complex_add_imag */
|
290
338
|
gsl_complex v20;
|
291
339
|
GSL_SET_COMPLEX(&v20, 2.0, 2.0);
|
@@ -302,17 +350,30 @@ double v21 = (double)2.0;
|
|
302
350
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r19),GSL_IMAG(c_r19));
|
303
351
|
printf(" assert r_r1.equals(r_r2)\n");
|
304
352
|
puts("end");
|
305
|
-
/*
|
306
|
-
gsl_complex
|
307
|
-
GSL_SET_COMPLEX(&
|
308
|
-
gsl_complex
|
309
|
-
|
310
|
-
puts("def
|
311
|
-
puts("
|
312
|
-
puts("
|
313
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
353
|
+
/* c_test_gsl_complex_log10 */
|
354
|
+
gsl_complex v41;
|
355
|
+
GSL_SET_COMPLEX(&v41, 2.0, 2.0);
|
356
|
+
gsl_complex c_r40;
|
357
|
+
c_r40 = gsl_complex_log10(v41);
|
358
|
+
puts("def test_gsl_complex_log10()");
|
359
|
+
puts(" v41 = GSL_Complex.create");
|
360
|
+
puts(" v41.set(2.0,2.0)");
|
361
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log10(v41)");
|
314
362
|
puts(" r_r2 = GSL_Complex.new");
|
315
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
363
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r40),GSL_IMAG(c_r40));
|
364
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
365
|
+
puts("end");
|
366
|
+
/* c_test_gsl_complex_arccsc */
|
367
|
+
gsl_complex v56;
|
368
|
+
GSL_SET_COMPLEX(&v56, 2.0, 2.0);
|
369
|
+
gsl_complex c_r55;
|
370
|
+
c_r55 = gsl_complex_arccsc(v56);
|
371
|
+
puts("def test_gsl_complex_arccsc()");
|
372
|
+
puts(" v56 = GSL_Complex.create");
|
373
|
+
puts(" v56.set(2.0,2.0)");
|
374
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsc(v56)");
|
375
|
+
puts(" r_r2 = GSL_Complex.new");
|
376
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r55),GSL_IMAG(c_r55));
|
316
377
|
printf(" assert r_r1.equals(r_r2)\n");
|
317
378
|
puts("end");
|
318
379
|
/* c_test_gsl_complex_arccsch */
|
@@ -328,20 +389,20 @@ double v21 = (double)2.0;
|
|
328
389
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r69),GSL_IMAG(c_r69));
|
329
390
|
printf(" assert r_r1.equals(r_r2)\n");
|
330
391
|
puts("end");
|
331
|
-
/*
|
332
|
-
gsl_complex
|
333
|
-
GSL_SET_COMPLEX(&
|
334
|
-
double
|
335
|
-
gsl_complex
|
336
|
-
|
337
|
-
puts("def
|
338
|
-
puts("
|
339
|
-
puts("
|
392
|
+
/* c_test_gsl_complex_div_imag */
|
393
|
+
gsl_complex v26;
|
394
|
+
GSL_SET_COMPLEX(&v26, 2.0, 2.0);
|
395
|
+
double v27 = (double)2.0;
|
396
|
+
gsl_complex c_r25;
|
397
|
+
c_r25 = gsl_complex_div_imag(v26,v27);
|
398
|
+
puts("def test_gsl_complex_div_imag()");
|
399
|
+
puts(" v26 = GSL_Complex.create");
|
400
|
+
puts(" v26.set(2.0,2.0)");
|
340
401
|
puts("");
|
341
|
-
puts("
|
342
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
402
|
+
puts(" v27 = 2.0");
|
403
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div_imag(v26,v27)");
|
343
404
|
puts(" r_r2 = GSL_Complex.new");
|
344
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
405
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r25),GSL_IMAG(c_r25));
|
345
406
|
printf(" assert r_r1.equals(r_r2)\n");
|
346
407
|
puts("end");
|
347
408
|
/* c_test_gsl_complex_cos */
|
@@ -357,17 +418,30 @@ double v13 = (double)2.0;
|
|
357
418
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r44),GSL_IMAG(c_r44));
|
358
419
|
printf(" assert r_r1.equals(r_r2)\n");
|
359
420
|
puts("end");
|
360
|
-
/*
|
361
|
-
gsl_complex
|
362
|
-
GSL_SET_COMPLEX(&
|
363
|
-
gsl_complex
|
364
|
-
|
365
|
-
puts("def
|
366
|
-
puts("
|
367
|
-
puts("
|
368
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
421
|
+
/* c_test_gsl_complex_arccot */
|
422
|
+
gsl_complex v59;
|
423
|
+
GSL_SET_COMPLEX(&v59, 2.0, 2.0);
|
424
|
+
gsl_complex c_r58;
|
425
|
+
c_r58 = gsl_complex_arccot(v59);
|
426
|
+
puts("def test_gsl_complex_arccot()");
|
427
|
+
puts(" v59 = GSL_Complex.create");
|
428
|
+
puts(" v59.set(2.0,2.0)");
|
429
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccot(v59)");
|
369
430
|
puts(" r_r2 = GSL_Complex.new");
|
370
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
431
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r58),GSL_IMAG(c_r58));
|
432
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
433
|
+
puts("end");
|
434
|
+
/* c_test_gsl_complex_arccoth */
|
435
|
+
gsl_complex v73;
|
436
|
+
GSL_SET_COMPLEX(&v73, 2.0, 2.0);
|
437
|
+
gsl_complex c_r72;
|
438
|
+
c_r72 = gsl_complex_arccoth(v73);
|
439
|
+
puts("def test_gsl_complex_arccoth()");
|
440
|
+
puts(" v73 = GSL_Complex.create");
|
441
|
+
puts(" v73.set(2.0,2.0)");
|
442
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccoth(v73)");
|
443
|
+
puts(" r_r2 = GSL_Complex.new");
|
444
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r72),GSL_IMAG(c_r72));
|
371
445
|
printf(" assert r_r1.equals(r_r2)\n");
|
372
446
|
puts("end");
|
373
447
|
/* c_test_gsl_complex_add */
|
@@ -387,25 +461,49 @@ double v13 = (double)2.0;
|
|
387
461
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r3),GSL_IMAG(c_r3));
|
388
462
|
printf(" assert r_r1.equals(r_r2)\n");
|
389
463
|
puts("end");
|
390
|
-
/*
|
391
|
-
gsl_complex
|
392
|
-
GSL_SET_COMPLEX(&
|
393
|
-
gsl_complex
|
394
|
-
|
395
|
-
puts("def
|
396
|
-
puts("
|
397
|
-
puts("
|
398
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
464
|
+
/* c_test_gsl_complex_negative */
|
465
|
+
gsl_complex v30;
|
466
|
+
GSL_SET_COMPLEX(&v30, 2.0, 2.0);
|
467
|
+
gsl_complex c_r29;
|
468
|
+
c_r29 = gsl_complex_negative(v30);
|
469
|
+
puts("def test_gsl_complex_negative()");
|
470
|
+
puts(" v30 = GSL_Complex.create");
|
471
|
+
puts(" v30.set(2.0,2.0)");
|
472
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_negative(v30)");
|
399
473
|
puts(" r_r2 = GSL_Complex.new");
|
400
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
474
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r29),GSL_IMAG(c_r29));
|
401
475
|
printf(" assert r_r1.equals(r_r2)\n");
|
402
476
|
puts("end");
|
403
|
-
/*
|
404
|
-
gsl_complex
|
405
|
-
GSL_SET_COMPLEX(&
|
406
|
-
gsl_complex
|
407
|
-
|
408
|
-
puts("def
|
477
|
+
/* c_test_gsl_complex_tan */
|
478
|
+
gsl_complex v48;
|
479
|
+
GSL_SET_COMPLEX(&v48, 2.0, 2.0);
|
480
|
+
gsl_complex c_r47;
|
481
|
+
c_r47 = gsl_complex_tan(v48);
|
482
|
+
puts("def test_gsl_complex_tan()");
|
483
|
+
puts(" v48 = GSL_Complex.create");
|
484
|
+
puts(" v48.set(2.0,2.0)");
|
485
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_tan(v48)");
|
486
|
+
puts(" r_r2 = GSL_Complex.new");
|
487
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r47),GSL_IMAG(c_r47));
|
488
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
489
|
+
puts("end");
|
490
|
+
/* c_test_gsl_complex_abs */
|
491
|
+
gsl_complex v2;
|
492
|
+
GSL_SET_COMPLEX(&v2, 2.0, 2.0);
|
493
|
+
double c_r1;
|
494
|
+
c_r1 = gsl_complex_abs(v2);
|
495
|
+
puts("def test_gsl_complex_abs()");
|
496
|
+
puts(" v2 = GSL_Complex.create");
|
497
|
+
puts(" v2.set(2.0,2.0)");
|
498
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs(v2)");
|
499
|
+
printf(" assert_in_delta r_r1, %.15g, EPSILON\n", c_r1);
|
500
|
+
puts("end");
|
501
|
+
/* c_test_gsl_complex_sech */
|
502
|
+
gsl_complex v62;
|
503
|
+
GSL_SET_COMPLEX(&v62, 2.0, 2.0);
|
504
|
+
gsl_complex c_r61;
|
505
|
+
c_r61 = gsl_complex_sech(v62);
|
506
|
+
puts("def test_gsl_complex_sech()");
|
409
507
|
puts(" v62 = GSL_Complex.create");
|
410
508
|
puts(" v62.set(2.0,2.0)");
|
411
509
|
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sech(v62)");
|
@@ -413,6 +511,23 @@ double v13 = (double)2.0;
|
|
413
511
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r61),GSL_IMAG(c_r61));
|
414
512
|
printf(" assert r_r1.equals(r_r2)\n");
|
415
513
|
puts("end");
|
514
|
+
/* c_test_gsl_complex_div */
|
515
|
+
gsl_complex v10;
|
516
|
+
GSL_SET_COMPLEX(&v10, 2.0, 2.0);
|
517
|
+
gsl_complex v11;
|
518
|
+
GSL_SET_COMPLEX(&v11, 2.0, 2.0);
|
519
|
+
gsl_complex c_r9;
|
520
|
+
c_r9 = gsl_complex_div(v10,v11);
|
521
|
+
puts("def test_gsl_complex_div()");
|
522
|
+
puts(" v10 = GSL_Complex.create");
|
523
|
+
puts(" v10.set(2.0,2.0)");
|
524
|
+
puts(" v11 = GSL_Complex.create");
|
525
|
+
puts(" v11.set(2.0,2.0)");
|
526
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_div(v10,v11)");
|
527
|
+
puts(" r_r2 = GSL_Complex.new");
|
528
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r9),GSL_IMAG(c_r9));
|
529
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
530
|
+
puts("end");
|
416
531
|
/* c_test_gsl_complex_polar */
|
417
532
|
double v33 = (double)2.0;
|
418
533
|
double v34 = (double)2.0;
|
@@ -428,30 +543,58 @@ double v34 = (double)2.0;
|
|
428
543
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r32),GSL_IMAG(c_r32));
|
429
544
|
printf(" assert r_r1.equals(r_r2)\n");
|
430
545
|
puts("end");
|
431
|
-
/*
|
432
|
-
|
433
|
-
|
434
|
-
|
435
|
-
|
436
|
-
puts("
|
437
|
-
puts("
|
438
|
-
puts("
|
439
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctan(v58)");
|
546
|
+
/* c_test_gsl_complex_arcsin_real */
|
547
|
+
double v51 = (double)2.0;
|
548
|
+
gsl_complex c_r50;
|
549
|
+
c_r50 = gsl_complex_arcsin_real(v51);
|
550
|
+
puts("def test_gsl_complex_arcsin_real()");
|
551
|
+
puts("");
|
552
|
+
puts(" v51 = 2.0");
|
553
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsin_real(v51)");
|
440
554
|
puts(" r_r2 = GSL_Complex.new");
|
441
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
555
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r50),GSL_IMAG(c_r50));
|
442
556
|
printf(" assert r_r1.equals(r_r2)\n");
|
443
557
|
puts("end");
|
444
|
-
/*
|
445
|
-
gsl_complex
|
446
|
-
GSL_SET_COMPLEX(&
|
447
|
-
gsl_complex
|
448
|
-
|
449
|
-
puts("def
|
450
|
-
puts("
|
451
|
-
puts("
|
452
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
558
|
+
/* c_test_gsl_complex_coth */
|
559
|
+
gsl_complex v65;
|
560
|
+
GSL_SET_COMPLEX(&v65, 2.0, 2.0);
|
561
|
+
gsl_complex c_r64;
|
562
|
+
c_r64 = gsl_complex_coth(v65);
|
563
|
+
puts("def test_gsl_complex_coth()");
|
564
|
+
puts(" v65 = GSL_Complex.create");
|
565
|
+
puts(" v65.set(2.0,2.0)");
|
566
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_coth(v65)");
|
453
567
|
puts(" r_r2 = GSL_Complex.new");
|
454
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
568
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r64),GSL_IMAG(c_r64));
|
569
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
570
|
+
puts("end");
|
571
|
+
/* c_test_gsl_complex_mul_real */
|
572
|
+
gsl_complex v16;
|
573
|
+
GSL_SET_COMPLEX(&v16, 2.0, 2.0);
|
574
|
+
double v17 = (double)2.0;
|
575
|
+
gsl_complex c_r15;
|
576
|
+
c_r15 = gsl_complex_mul_real(v16,v17);
|
577
|
+
puts("def test_gsl_complex_mul_real()");
|
578
|
+
puts(" v16 = GSL_Complex.create");
|
579
|
+
puts(" v16.set(2.0,2.0)");
|
580
|
+
puts("");
|
581
|
+
puts(" v17 = 2.0");
|
582
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_real(v16,v17)");
|
583
|
+
puts(" r_r2 = GSL_Complex.new");
|
584
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r15),GSL_IMAG(c_r15));
|
585
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
586
|
+
puts("end");
|
587
|
+
/* c_test_gsl_complex_exp */
|
588
|
+
gsl_complex v39;
|
589
|
+
GSL_SET_COMPLEX(&v39, 2.0, 2.0);
|
590
|
+
gsl_complex c_r38;
|
591
|
+
c_r38 = gsl_complex_exp(v39);
|
592
|
+
puts("def test_gsl_complex_exp()");
|
593
|
+
puts(" v39 = GSL_Complex.create");
|
594
|
+
puts(" v39.set(2.0,2.0)");
|
595
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_exp(v39)");
|
596
|
+
puts(" r_r2 = GSL_Complex.new");
|
597
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r38),GSL_IMAG(c_r38));
|
455
598
|
printf(" assert r_r1.equals(r_r2)\n");
|
456
599
|
puts("end");
|
457
600
|
/* c_test_gsl_complex_arcsec */
|
@@ -467,6 +610,18 @@ double v34 = (double)2.0;
|
|
467
610
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r53),GSL_IMAG(c_r53));
|
468
611
|
printf(" assert r_r1.equals(r_r2)\n");
|
469
612
|
puts("end");
|
613
|
+
/* c_test_gsl_complex_arccosh_real */
|
614
|
+
double v68 = (double)2.0;
|
615
|
+
gsl_complex c_r67;
|
616
|
+
c_r67 = gsl_complex_arccosh_real(v68);
|
617
|
+
puts("def test_gsl_complex_arccosh_real()");
|
618
|
+
puts("");
|
619
|
+
puts(" v68 = 2.0");
|
620
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccosh_real(v68)");
|
621
|
+
puts(" r_r2 = GSL_Complex.new");
|
622
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r67),GSL_IMAG(c_r67));
|
623
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
624
|
+
puts("end");
|
470
625
|
/* c_test_gsl_complex_sub_imag */
|
471
626
|
gsl_complex v22;
|
472
627
|
GSL_SET_COMPLEX(&v22, 2.0, 2.0);
|
@@ -483,29 +638,34 @@ double v23 = (double)2.0;
|
|
483
638
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r21),GSL_IMAG(c_r21));
|
484
639
|
printf(" assert r_r1.equals(r_r2)\n");
|
485
640
|
puts("end");
|
486
|
-
/*
|
487
|
-
gsl_complex
|
488
|
-
GSL_SET_COMPLEX(&
|
489
|
-
gsl_complex
|
490
|
-
|
491
|
-
|
492
|
-
|
493
|
-
puts("
|
494
|
-
puts("
|
641
|
+
/* c_test_gsl_complex_log_b */
|
642
|
+
gsl_complex v42;
|
643
|
+
GSL_SET_COMPLEX(&v42, 2.0, 2.0);
|
644
|
+
gsl_complex v43;
|
645
|
+
GSL_SET_COMPLEX(&v43, 2.0, 2.0);
|
646
|
+
gsl_complex c_r41;
|
647
|
+
c_r41 = gsl_complex_log_b(v42,v43);
|
648
|
+
puts("def test_gsl_complex_log_b()");
|
649
|
+
puts(" v42 = GSL_Complex.create");
|
650
|
+
puts(" v42.set(2.0,2.0)");
|
651
|
+
puts(" v43 = GSL_Complex.create");
|
652
|
+
puts(" v43.set(2.0,2.0)");
|
653
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log_b(v42,v43)");
|
495
654
|
puts(" r_r2 = GSL_Complex.new");
|
496
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
655
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r41),GSL_IMAG(c_r41));
|
497
656
|
printf(" assert r_r1.equals(r_r2)\n");
|
498
657
|
puts("end");
|
499
|
-
/*
|
500
|
-
|
501
|
-
|
502
|
-
|
503
|
-
|
504
|
-
puts("
|
505
|
-
puts("
|
506
|
-
puts("
|
507
|
-
puts("
|
508
|
-
printf("
|
658
|
+
/* c_test_gsl_complex_arccsc_real */
|
659
|
+
double v57 = (double)2.0;
|
660
|
+
gsl_complex c_r56;
|
661
|
+
c_r56 = gsl_complex_arccsc_real(v57);
|
662
|
+
puts("def test_gsl_complex_arccsc_real()");
|
663
|
+
puts("");
|
664
|
+
puts(" v57 = 2.0");
|
665
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccsc_real(v57)");
|
666
|
+
puts(" r_r2 = GSL_Complex.new");
|
667
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r56),GSL_IMAG(c_r56));
|
668
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
509
669
|
puts("end");
|
510
670
|
/* c_test_gsl_complex_arctanh */
|
511
671
|
gsl_complex v71;
|
@@ -520,20 +680,17 @@ double v23 = (double)2.0;
|
|
520
680
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r70),GSL_IMAG(c_r70));
|
521
681
|
printf(" assert r_r1.equals(r_r2)\n");
|
522
682
|
puts("end");
|
523
|
-
/*
|
524
|
-
gsl_complex
|
525
|
-
GSL_SET_COMPLEX(&
|
526
|
-
|
527
|
-
|
528
|
-
|
529
|
-
puts("
|
530
|
-
puts("
|
531
|
-
puts("
|
532
|
-
puts("");
|
533
|
-
puts(" v15 = 2.0");
|
534
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sub_real(v14,v15)");
|
683
|
+
/* c_test_gsl_complex_conjugate */
|
684
|
+
gsl_complex v28;
|
685
|
+
GSL_SET_COMPLEX(&v28, 2.0, 2.0);
|
686
|
+
gsl_complex c_r27;
|
687
|
+
c_r27 = gsl_complex_conjugate(v28);
|
688
|
+
puts("def test_gsl_complex_conjugate()");
|
689
|
+
puts(" v28 = GSL_Complex.create");
|
690
|
+
puts(" v28.set(2.0,2.0)");
|
691
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_conjugate(v28)");
|
535
692
|
puts(" r_r2 = GSL_Complex.new");
|
536
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
693
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r27),GSL_IMAG(c_r27));
|
537
694
|
printf(" assert r_r1.equals(r_r2)\n");
|
538
695
|
puts("end");
|
539
696
|
/* c_test_gsl_complex_sec */
|
@@ -549,17 +706,17 @@ double v15 = (double)2.0;
|
|
549
706
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r45),GSL_IMAG(c_r45));
|
550
707
|
printf(" assert r_r1.equals(r_r2)\n");
|
551
708
|
puts("end");
|
552
|
-
/*
|
553
|
-
gsl_complex
|
554
|
-
GSL_SET_COMPLEX(&
|
555
|
-
gsl_complex
|
556
|
-
|
557
|
-
puts("def
|
558
|
-
puts("
|
559
|
-
puts("
|
560
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
709
|
+
/* c_test_gsl_complex_sinh */
|
710
|
+
gsl_complex v60;
|
711
|
+
GSL_SET_COMPLEX(&v60, 2.0, 2.0);
|
712
|
+
gsl_complex c_r59;
|
713
|
+
c_r59 = gsl_complex_sinh(v60);
|
714
|
+
puts("def test_gsl_complex_sinh()");
|
715
|
+
puts(" v60 = GSL_Complex.create");
|
716
|
+
puts(" v60.set(2.0,2.0)");
|
717
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sinh(v60)");
|
561
718
|
puts(" r_r2 = GSL_Complex.new");
|
562
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
719
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r59),GSL_IMAG(c_r59));
|
563
720
|
printf(" assert r_r1.equals(r_r2)\n");
|
564
721
|
puts("end");
|
565
722
|
/* c_test_gsl_complex_sub */
|
@@ -579,19 +736,43 @@ double v15 = (double)2.0;
|
|
579
736
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r5),GSL_IMAG(c_r5));
|
580
737
|
printf(" assert r_r1.equals(r_r2)\n");
|
581
738
|
puts("end");
|
582
|
-
/*
|
583
|
-
gsl_complex
|
584
|
-
GSL_SET_COMPLEX(&
|
585
|
-
gsl_complex
|
586
|
-
|
587
|
-
puts("def
|
588
|
-
puts("
|
589
|
-
puts("
|
590
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::
|
739
|
+
/* c_test_gsl_complex_sqrt */
|
740
|
+
gsl_complex v31;
|
741
|
+
GSL_SET_COMPLEX(&v31, 2.0, 2.0);
|
742
|
+
gsl_complex c_r30;
|
743
|
+
c_r30 = gsl_complex_sqrt(v31);
|
744
|
+
puts("def test_gsl_complex_sqrt()");
|
745
|
+
puts(" v31 = GSL_Complex.create");
|
746
|
+
puts(" v31.set(2.0,2.0)");
|
747
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_sqrt(v31)");
|
591
748
|
puts(" r_r2 = GSL_Complex.new");
|
592
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(
|
749
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r30),GSL_IMAG(c_r30));
|
750
|
+
printf(" assert r_r1.equals(r_r2)\n");
|
751
|
+
puts("end");
|
752
|
+
/* c_test_gsl_complex_cot */
|
753
|
+
gsl_complex v49;
|
754
|
+
GSL_SET_COMPLEX(&v49, 2.0, 2.0);
|
755
|
+
gsl_complex c_r48;
|
756
|
+
c_r48 = gsl_complex_cot(v49);
|
757
|
+
puts("def test_gsl_complex_cot()");
|
758
|
+
puts(" v49 = GSL_Complex.create");
|
759
|
+
puts(" v49.set(2.0,2.0)");
|
760
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_cot(v49)");
|
761
|
+
puts(" r_r2 = GSL_Complex.new");
|
762
|
+
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r48),GSL_IMAG(c_r48));
|
593
763
|
printf(" assert r_r1.equals(r_r2)\n");
|
594
764
|
puts("end");
|
765
|
+
/* c_test_gsl_complex_abs2 */
|
766
|
+
gsl_complex v3;
|
767
|
+
GSL_SET_COMPLEX(&v3, 2.0, 2.0);
|
768
|
+
double c_r2;
|
769
|
+
c_r2 = gsl_complex_abs2(v3);
|
770
|
+
puts("def test_gsl_complex_abs2()");
|
771
|
+
puts(" v3 = GSL_Complex.create");
|
772
|
+
puts(" v3.set(2.0,2.0)");
|
773
|
+
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs2(v3)");
|
774
|
+
printf(" assert_in_delta r_r1, %.15g, EPSILON\n", c_r2);
|
775
|
+
puts("end");
|
595
776
|
/* c_test_gsl_complex_csch */
|
596
777
|
gsl_complex v63;
|
597
778
|
GSL_SET_COMPLEX(&v63, 2.0, 2.0);
|
@@ -605,187 +786,6 @@ double v15 = (double)2.0;
|
|
605
786
|
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r62),GSL_IMAG(c_r62));
|
606
787
|
printf(" assert r_r1.equals(r_r2)\n");
|
607
788
|
puts("end");
|
608
|
-
/* c_test_gsl_complex_pow */
|
609
|
-
gsl_complex v35;
|
610
|
-
GSL_SET_COMPLEX(&v35, 2.0, 2.0);
|
611
|
-
gsl_complex v36;
|
612
|
-
GSL_SET_COMPLEX(&v36, 2.0, 2.0);
|
613
|
-
gsl_complex c_r34;
|
614
|
-
c_r34 = gsl_complex_pow(v35,v36);
|
615
|
-
puts("def test_gsl_complex_pow()");
|
616
|
-
puts(" v35 = GSL_Complex.create");
|
617
|
-
puts(" v35.set(2.0,2.0)");
|
618
|
-
puts(" v36 = GSL_Complex.create");
|
619
|
-
puts(" v36.set(2.0,2.0)");
|
620
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_pow(v35,v36)");
|
621
|
-
puts(" r_r2 = GSL_Complex.new");
|
622
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r34),GSL_IMAG(c_r34));
|
623
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
624
|
-
puts("end");
|
625
|
-
/* c_test_gsl_complex_arccot */
|
626
|
-
gsl_complex v59;
|
627
|
-
GSL_SET_COMPLEX(&v59, 2.0, 2.0);
|
628
|
-
gsl_complex c_r58;
|
629
|
-
c_r58 = gsl_complex_arccot(v59);
|
630
|
-
puts("def test_gsl_complex_arccot()");
|
631
|
-
puts(" v59 = GSL_Complex.create");
|
632
|
-
puts(" v59.set(2.0,2.0)");
|
633
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccot(v59)");
|
634
|
-
puts(" r_r2 = GSL_Complex.new");
|
635
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r58),GSL_IMAG(c_r58));
|
636
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
637
|
-
puts("end");
|
638
|
-
/* c_test_gsl_complex_negative */
|
639
|
-
gsl_complex v30;
|
640
|
-
GSL_SET_COMPLEX(&v30, 2.0, 2.0);
|
641
|
-
gsl_complex c_r29;
|
642
|
-
c_r29 = gsl_complex_negative(v30);
|
643
|
-
puts("def test_gsl_complex_negative()");
|
644
|
-
puts(" v30 = GSL_Complex.create");
|
645
|
-
puts(" v30.set(2.0,2.0)");
|
646
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_negative(v30)");
|
647
|
-
puts(" r_r2 = GSL_Complex.new");
|
648
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r29),GSL_IMAG(c_r29));
|
649
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
650
|
-
puts("end");
|
651
|
-
/* c_test_gsl_complex_arcsec_real */
|
652
|
-
double v55 = (double)2.0;
|
653
|
-
gsl_complex c_r54;
|
654
|
-
c_r54 = gsl_complex_arcsec_real(v55);
|
655
|
-
puts("def test_gsl_complex_arcsec_real()");
|
656
|
-
puts("");
|
657
|
-
puts(" v55 = 2.0");
|
658
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsec_real(v55)");
|
659
|
-
puts(" r_r2 = GSL_Complex.new");
|
660
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r54),GSL_IMAG(c_r54));
|
661
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
662
|
-
puts("end");
|
663
|
-
/* c_test_gsl_complex_mul_imag */
|
664
|
-
gsl_complex v24;
|
665
|
-
GSL_SET_COMPLEX(&v24, 2.0, 2.0);
|
666
|
-
double v25 = (double)2.0;
|
667
|
-
gsl_complex c_r23;
|
668
|
-
c_r23 = gsl_complex_mul_imag(v24,v25);
|
669
|
-
puts("def test_gsl_complex_mul_imag()");
|
670
|
-
puts(" v24 = GSL_Complex.create");
|
671
|
-
puts(" v24.set(2.0,2.0)");
|
672
|
-
puts("");
|
673
|
-
puts(" v25 = 2.0");
|
674
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_imag(v24,v25)");
|
675
|
-
puts(" r_r2 = GSL_Complex.new");
|
676
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r23),GSL_IMAG(c_r23));
|
677
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
678
|
-
puts("end");
|
679
|
-
/* c_test_gsl_complex_arcsin_real */
|
680
|
-
double v51 = (double)2.0;
|
681
|
-
gsl_complex c_r50;
|
682
|
-
c_r50 = gsl_complex_arcsin_real(v51);
|
683
|
-
puts("def test_gsl_complex_arcsin_real()");
|
684
|
-
puts("");
|
685
|
-
puts(" v51 = 2.0");
|
686
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arcsin_real(v51)");
|
687
|
-
puts(" r_r2 = GSL_Complex.new");
|
688
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r50),GSL_IMAG(c_r50));
|
689
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
690
|
-
puts("end");
|
691
|
-
/* c_test_gsl_complex_abs */
|
692
|
-
gsl_complex v2;
|
693
|
-
GSL_SET_COMPLEX(&v2, 2.0, 2.0);
|
694
|
-
double c_r1;
|
695
|
-
c_r1 = gsl_complex_abs(v2);
|
696
|
-
puts("def test_gsl_complex_abs()");
|
697
|
-
puts(" v2 = GSL_Complex.create");
|
698
|
-
puts(" v2.set(2.0,2.0)");
|
699
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_abs(v2)");
|
700
|
-
printf(" assert_in_delta r_r1, %.15g, EPSILON\n", c_r1);
|
701
|
-
puts("end");
|
702
|
-
/* c_test_gsl_complex_arctanh_real */
|
703
|
-
double v72 = (double)2.0;
|
704
|
-
gsl_complex c_r71;
|
705
|
-
c_r71 = gsl_complex_arctanh_real(v72);
|
706
|
-
puts("def test_gsl_complex_arctanh_real()");
|
707
|
-
puts("");
|
708
|
-
puts(" v72 = 2.0");
|
709
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arctanh_real(v72)");
|
710
|
-
puts(" r_r2 = GSL_Complex.new");
|
711
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r71),GSL_IMAG(c_r71));
|
712
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
713
|
-
puts("end");
|
714
|
-
/* c_test_gsl_complex_mul_real */
|
715
|
-
gsl_complex v16;
|
716
|
-
GSL_SET_COMPLEX(&v16, 2.0, 2.0);
|
717
|
-
double v17 = (double)2.0;
|
718
|
-
gsl_complex c_r15;
|
719
|
-
c_r15 = gsl_complex_mul_real(v16,v17);
|
720
|
-
puts("def test_gsl_complex_mul_real()");
|
721
|
-
puts(" v16 = GSL_Complex.create");
|
722
|
-
puts(" v16.set(2.0,2.0)");
|
723
|
-
puts("");
|
724
|
-
puts(" v17 = 2.0");
|
725
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul_real(v16,v17)");
|
726
|
-
puts(" r_r2 = GSL_Complex.new");
|
727
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r15),GSL_IMAG(c_r15));
|
728
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
729
|
-
puts("end");
|
730
|
-
/* c_test_gsl_complex_csc */
|
731
|
-
gsl_complex v47;
|
732
|
-
GSL_SET_COMPLEX(&v47, 2.0, 2.0);
|
733
|
-
gsl_complex c_r46;
|
734
|
-
c_r46 = gsl_complex_csc(v47);
|
735
|
-
puts("def test_gsl_complex_csc()");
|
736
|
-
puts(" v47 = GSL_Complex.create");
|
737
|
-
puts(" v47.set(2.0,2.0)");
|
738
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_csc(v47)");
|
739
|
-
puts(" r_r2 = GSL_Complex.new");
|
740
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r46),GSL_IMAG(c_r46));
|
741
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
742
|
-
puts("end");
|
743
|
-
/* c_test_gsl_complex_arccosh_real */
|
744
|
-
double v68 = (double)2.0;
|
745
|
-
gsl_complex c_r67;
|
746
|
-
c_r67 = gsl_complex_arccosh_real(v68);
|
747
|
-
puts("def test_gsl_complex_arccosh_real()");
|
748
|
-
puts("");
|
749
|
-
puts(" v68 = 2.0");
|
750
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_arccosh_real(v68)");
|
751
|
-
puts(" r_r2 = GSL_Complex.new");
|
752
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r67),GSL_IMAG(c_r67));
|
753
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
754
|
-
puts("end");
|
755
|
-
/* c_test_gsl_complex_mul */
|
756
|
-
gsl_complex v8;
|
757
|
-
GSL_SET_COMPLEX(&v8, 2.0, 2.0);
|
758
|
-
gsl_complex v9;
|
759
|
-
GSL_SET_COMPLEX(&v9, 2.0, 2.0);
|
760
|
-
gsl_complex c_r7;
|
761
|
-
c_r7 = gsl_complex_mul(v8,v9);
|
762
|
-
puts("def test_gsl_complex_mul()");
|
763
|
-
puts(" v8 = GSL_Complex.create");
|
764
|
-
puts(" v8.set(2.0,2.0)");
|
765
|
-
puts(" v9 = GSL_Complex.create");
|
766
|
-
puts(" v9.set(2.0,2.0)");
|
767
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_mul(v8,v9)");
|
768
|
-
puts(" r_r2 = GSL_Complex.new");
|
769
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r7),GSL_IMAG(c_r7));
|
770
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
771
|
-
puts("end");
|
772
|
-
/* c_test_gsl_complex_log_b */
|
773
|
-
gsl_complex v42;
|
774
|
-
GSL_SET_COMPLEX(&v42, 2.0, 2.0);
|
775
|
-
gsl_complex v43;
|
776
|
-
GSL_SET_COMPLEX(&v43, 2.0, 2.0);
|
777
|
-
gsl_complex c_r41;
|
778
|
-
c_r41 = gsl_complex_log_b(v42,v43);
|
779
|
-
puts("def test_gsl_complex_log_b()");
|
780
|
-
puts(" v42 = GSL_Complex.create");
|
781
|
-
puts(" v42.set(2.0,2.0)");
|
782
|
-
puts(" v43 = GSL_Complex.create");
|
783
|
-
puts(" v43.set(2.0,2.0)");
|
784
|
-
puts(" r_r1 = ::GSL4r::Complex::Methods::gsl_complex_log_b(v42,v43)");
|
785
|
-
puts(" r_r2 = GSL_Complex.new");
|
786
|
-
printf(" r_r2.set(%.15g,%.15g)\n",GSL_REAL(c_r41),GSL_IMAG(c_r41));
|
787
|
-
printf(" assert r_r1.equals(r_r2)\n");
|
788
|
-
puts("end");
|
789
789
|
puts("end");
|
790
790
|
return(0);
|
791
791
|
}
|