gsl4r 0.0.1 → 0.0.2
Sign up to get free protection for your applications and to get access to all the features.
- 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
|
}
|