rb-gsl 1.16.0.4 → 1.16.0.5
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/.gitignore +12 -0
- data/.travis.yml +24 -0
- data/AUTHORS +6 -0
- data/COPYING +341 -339
- data/ChangeLog +608 -587
- data/Gemfile +4 -0
- data/README.md +77 -0
- data/Rakefile +14 -33
- data/THANKS +1 -1
- data/examples/blas/blas.rb +1 -1
- data/examples/bspline.rb +3 -3
- data/examples/complex/functions.rb +4 -4
- data/examples/complex/mul.rb +1 -1
- data/examples/const/physconst.rb +4 -4
- data/examples/const/travel.rb +4 -4
- data/examples/eigen/nonsymmv.rb +1 -1
- data/examples/eigen/qhoscillator.rb +3 -3
- data/examples/fft/radix2.rb +1 -1
- data/examples/fft/real-halfcomplex.rb +3 -3
- data/examples/fft/real-halfcomplex2.rb +3 -3
- data/examples/fit/expfit.rb +1 -1
- data/examples/fit/multifit.rb +1 -1
- data/examples/fit/ndlinear.rb +44 -44
- data/examples/fit/nonlinearfit.rb +1 -1
- data/examples/fit/wlinear.rb +3 -3
- data/examples/function/function.rb +1 -1
- data/examples/function/min.rb +1 -1
- data/examples/function/synchrotron.rb +2 -2
- data/examples/gallery/koch.rb +1 -1
- data/examples/histogram/cauchy.rb +2 -2
- data/examples/histogram/exponential.rb +1 -1
- data/examples/histogram/histo3d.rb +1 -1
- data/examples/histogram/histogram-pdf.rb +2 -2
- data/examples/histogram/xexp.rb +1 -1
- data/examples/integration/ahmed.rb +2 -2
- data/examples/integration/cosmology.rb +7 -7
- data/examples/integration/friedmann.rb +4 -4
- data/examples/integration/qng.rb +1 -1
- data/examples/interp/demo.rb +2 -2
- data/examples/linalg/LQ_solve.rb +1 -1
- data/examples/linalg/LU.rb +1 -1
- data/examples/linalg/LU2.rb +1 -1
- data/examples/linalg/LU_narray.rb +1 -1
- data/examples/linalg/PTLQ.rb +1 -1
- data/examples/linalg/QRPT.rb +1 -1
- data/examples/linalg/chol.rb +1 -1
- data/examples/linalg/chol_narray.rb +1 -1
- data/examples/linalg/complex.rb +1 -1
- data/examples/math/elementary.rb +1 -1
- data/examples/math/functions.rb +1 -1
- data/examples/math/inf_nan.rb +1 -1
- data/examples/math/minmax.rb +1 -1
- data/examples/math/power.rb +1 -1
- data/examples/math/test.rb +1 -1
- data/examples/min.rb +1 -1
- data/examples/multimin/bundle.rb +1 -1
- data/examples/multimin/cqp.rb +17 -17
- data/examples/multiroot/fsolver3.rb +1 -1
- data/examples/odeiv/binarysystem.rb +12 -12
- data/examples/odeiv/demo.rb +3 -3
- data/examples/odeiv/frei1.rb +7 -7
- data/examples/odeiv/frei2.rb +4 -4
- data/examples/odeiv/oscillator.rb +1 -1
- data/examples/odeiv/sedov.rb +3 -3
- data/examples/odeiv/whitedwarf.rb +11 -11
- data/examples/permutation/ex1.rb +2 -2
- data/examples/permutation/permutation.rb +1 -1
- data/examples/poly/demo.rb +1 -1
- data/examples/random/diffusion.rb +1 -1
- data/examples/random/generator.rb +2 -2
- data/examples/random/randomwalk.rb +3 -3
- data/examples/random/rng.rb +1 -1
- data/examples/roots/bisection.rb +1 -1
- data/examples/roots/brent.rb +1 -1
- data/examples/roots/demo.rb +1 -1
- data/examples/roots/newton.rb +2 -2
- data/examples/roots/recombination.gp +0 -1
- data/examples/sf/hyperg.rb +1 -1
- data/examples/sf/sphbessel.rb +1 -1
- data/examples/sort/sort.rb +1 -1
- data/examples/tamu_anova.rb +4 -4
- data/examples/vector/add.rb +1 -1
- data/examples/vector/decimate.rb +1 -1
- data/examples/vector/gnuplot.rb +8 -8
- data/examples/vector/vector.rb +2 -2
- data/examples/wavelet/wavelet1.rb +1 -1
- data/ext/{gsl → gsl_native}/alf.c +9 -9
- data/ext/{gsl → gsl_native}/array.c +64 -151
- data/ext/{gsl → gsl_native}/array_complex.c +61 -63
- data/ext/{gsl → gsl_native}/blas.c +0 -0
- data/ext/{gsl → gsl_native}/blas1.c +32 -32
- data/ext/{gsl → gsl_native}/blas2.c +54 -58
- data/ext/{gsl → gsl_native}/blas3.c +54 -54
- data/ext/{gsl → gsl_native}/block.c +3 -6
- data/ext/{gsl/block_source.c → gsl_native/block_source.h} +109 -111
- data/ext/gsl_native/bspline.c +122 -0
- data/ext/{gsl → gsl_native}/bundle.c +0 -0
- data/ext/{gsl → gsl_native}/cdf.c +76 -88
- data/ext/{gsl → gsl_native}/cheb.c +74 -85
- data/ext/{gsl → gsl_native}/combination.c +9 -16
- data/ext/{gsl → gsl_native}/common.c +7 -37
- data/ext/{gsl → gsl_native}/complex.c +113 -115
- data/ext/gsl_native/const.c +331 -0
- data/ext/{gsl → gsl_native}/const_additional.c +12 -33
- data/ext/gsl_native/cqp.c +283 -0
- data/ext/{gsl → gsl_native}/deriv.c +22 -30
- data/ext/{gsl → gsl_native}/dht.c +20 -28
- data/ext/{gsl → gsl_native}/diff.c +23 -25
- data/ext/{gsl → gsl_native}/dirac.c +44 -45
- data/ext/{gsl → gsl_native}/eigen.c +1031 -1082
- data/ext/{gsl → gsl_native}/error.c +13 -13
- data/ext/{gsl → gsl_native}/extconf.rb +5 -13
- data/ext/{gsl → gsl_native}/fft.c +192 -198
- data/ext/{gsl → gsl_native}/fit.c +16 -16
- data/ext/gsl_native/fresnel.c +312 -0
- data/ext/{gsl → gsl_native}/function.c +29 -35
- data/ext/{gsl → gsl_native}/geometry.c +15 -15
- data/ext/{gsl → gsl_native}/graph.c +37 -87
- data/ext/{gsl → gsl_native}/gsl.c +10 -23
- data/ext/{gsl → gsl_native}/gsl_narray.c +15 -24
- data/ext/{gsl → gsl_native}/histogram.c +123 -150
- data/ext/{gsl → gsl_native}/histogram2d.c +71 -96
- data/ext/{gsl → gsl_native}/histogram3d.c +72 -72
- data/ext/{gsl → gsl_native}/histogram3d_source.c +195 -196
- data/ext/{gsl → gsl_native}/histogram_find.c +31 -33
- data/ext/{gsl → gsl_native}/histogram_oper.c +43 -52
- data/ext/{gsl → gsl_native}/ieee.c +7 -15
- data/ext/{gsl → gsl_native}/include/rb_gsl.h +0 -7
- data/ext/{gsl → gsl_native}/include/rb_gsl_array.h +15 -25
- data/ext/{gsl → gsl_native}/include/rb_gsl_cheb.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_common.h +44 -52
- data/ext/{gsl → gsl_native}/include/rb_gsl_complex.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_const.h +0 -5
- data/ext/{gsl → gsl_native}/include/rb_gsl_dirac.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_eigen.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_fft.h +0 -3
- data/ext/{gsl → gsl_native}/include/rb_gsl_fit.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_function.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_graph.h +2 -2
- data/ext/{gsl → gsl_native}/include/rb_gsl_histogram.h +8 -8
- data/ext/{gsl → gsl_native}/include/rb_gsl_histogram3d.h +50 -50
- data/ext/{gsl → gsl_native}/include/rb_gsl_integration.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_interp.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_linalg.h +0 -1
- data/ext/{gsl → gsl_native}/include/rb_gsl_math.h +0 -4
- data/ext/{gsl → gsl_native}/include/rb_gsl_odeiv.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_poly.h +3 -3
- data/ext/{gsl → gsl_native}/include/rb_gsl_rational.h +1 -1
- data/ext/{gsl → gsl_native}/include/rb_gsl_rng.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_root.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_sf.h +39 -43
- data/ext/{gsl → gsl_native}/include/rb_gsl_statistics.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_tensor.h +0 -0
- data/ext/{gsl → gsl_native}/include/rb_gsl_with_narray.h +0 -0
- data/ext/{gsl → gsl_native}/include/templates_off.h +0 -0
- data/ext/{gsl → gsl_native}/include/templates_on.h +1 -1
- data/ext/{gsl → gsl_native}/integration.c +160 -183
- data/ext/{gsl → gsl_native}/interp.c +23 -31
- data/ext/gsl_native/jacobi.c +733 -0
- data/ext/{gsl → gsl_native}/linalg.c +448 -576
- data/ext/{gsl → gsl_native}/linalg_complex.c +84 -98
- data/ext/{gsl → gsl_native}/math.c +45 -63
- data/ext/{gsl → gsl_native}/matrix.c +4 -6
- data/ext/{gsl → gsl_native}/matrix_complex.c +111 -114
- data/ext/{gsl → gsl_native}/matrix_double.c +74 -77
- data/ext/{gsl → gsl_native}/matrix_int.c +49 -50
- data/ext/{gsl/matrix_source.c → gsl_native/matrix_source.h} +280 -302
- data/ext/{gsl → gsl_native}/min.c +43 -74
- data/ext/{gsl → gsl_native}/monte.c +45 -59
- data/ext/{gsl → gsl_native}/multifit.c +138 -147
- data/ext/{gsl → gsl_native}/multimin.c +55 -83
- data/ext/{gsl → gsl_native}/multimin_fsdf.c +15 -15
- data/ext/{gsl → gsl_native}/multiroots.c +61 -64
- data/ext/{gsl → gsl_native}/multiset.c +3 -7
- data/ext/{gsl → gsl_native}/ndlinear.c +320 -321
- data/ext/{gsl → gsl_native}/nmf.c +11 -11
- data/ext/{gsl → gsl_native}/nmf_wrap.c +1 -1
- data/ext/{gsl → gsl_native}/ntuple.c +13 -13
- data/ext/{gsl → gsl_native}/odeiv.c +95 -106
- data/ext/gsl_native/ool.c +879 -0
- data/ext/{gsl → gsl_native}/permutation.c +36 -34
- data/ext/{gsl → gsl_native}/poly.c +2 -5
- data/ext/{gsl → gsl_native}/poly2.c +6 -6
- data/ext/{gsl/poly_source.c → gsl_native/poly_source.h} +232 -275
- data/ext/{gsl → gsl_native}/qrng.c +6 -17
- data/ext/{gsl → gsl_native}/randist.c +218 -243
- data/ext/{gsl → gsl_native}/rational.c +9 -9
- data/ext/{gsl → gsl_native}/rng.c +29 -45
- data/ext/{gsl → gsl_native}/root.c +43 -43
- data/ext/{gsl → gsl_native}/sf.c +163 -197
- data/ext/{gsl → gsl_native}/sf_airy.c +0 -0
- data/ext/{gsl → gsl_native}/sf_bessel.c +6 -6
- data/ext/{gsl → gsl_native}/sf_clausen.c +0 -0
- data/ext/{gsl → gsl_native}/sf_coulomb.c +39 -39
- data/ext/{gsl → gsl_native}/sf_coupling.c +29 -29
- data/ext/{gsl → gsl_native}/sf_dawson.c +0 -0
- data/ext/{gsl → gsl_native}/sf_debye.c +0 -9
- data/ext/{gsl → gsl_native}/sf_dilog.c +0 -0
- data/ext/{gsl → gsl_native}/sf_elementary.c +2 -2
- data/ext/{gsl → gsl_native}/sf_ellint.c +42 -42
- data/ext/{gsl → gsl_native}/sf_elljac.c +2 -2
- data/ext/{gsl → gsl_native}/sf_erfc.c +0 -4
- data/ext/{gsl → gsl_native}/sf_exp.c +2 -2
- data/ext/{gsl → gsl_native}/sf_expint.c +1 -11
- data/ext/{gsl → gsl_native}/sf_fermi_dirac.c +0 -0
- data/ext/{gsl → gsl_native}/sf_gamma.c +1 -5
- data/ext/{gsl → gsl_native}/sf_gegenbauer.c +0 -0
- data/ext/{gsl → gsl_native}/sf_hyperg.c +0 -0
- data/ext/{gsl → gsl_native}/sf_laguerre.c +3 -3
- data/ext/{gsl → gsl_native}/sf_lambert.c +0 -0
- data/ext/{gsl → gsl_native}/sf_legendre.c +0 -0
- data/ext/{gsl → gsl_native}/sf_log.c +3 -3
- data/ext/gsl_native/sf_mathieu.c +235 -0
- data/ext/{gsl → gsl_native}/sf_power.c +0 -0
- data/ext/{gsl → gsl_native}/sf_psi.c +2 -11
- data/ext/{gsl → gsl_native}/sf_synchrotron.c +0 -0
- data/ext/{gsl → gsl_native}/sf_transport.c +0 -0
- data/ext/{gsl → gsl_native}/sf_trigonometric.c +3 -3
- data/ext/{gsl → gsl_native}/sf_zeta.c +0 -4
- data/ext/{gsl → gsl_native}/signal.c +60 -64
- data/ext/{gsl → gsl_native}/siman.c +36 -40
- data/ext/{gsl → gsl_native}/sort.c +3 -4
- data/ext/{gsl → gsl_native}/spline.c +25 -43
- data/ext/{gsl → gsl_native}/stats.c +102 -115
- data/ext/{gsl → gsl_native}/sum.c +32 -32
- data/ext/{gsl → gsl_native}/tamu_anova.c +0 -0
- data/ext/{gsl → gsl_native}/tensor.c +2 -4
- data/ext/{gsl/tensor_source.c → gsl_native/tensor_source.h} +144 -144
- data/ext/{gsl → gsl_native}/vector.c +2 -4
- data/ext/{gsl → gsl_native}/vector_complex.c +165 -171
- data/ext/{gsl → gsl_native}/vector_double.c +167 -174
- data/ext/{gsl → gsl_native}/vector_int.c +23 -25
- data/ext/{gsl/vector_source.c → gsl_native/vector_source.h} +400 -414
- data/ext/{gsl → gsl_native}/wavelet.c +215 -236
- data/lib/gsl.rb +1 -7
- data/lib/gsl/gnuplot.rb +2 -2
- data/lib/gsl/version.rb +1 -3
- data/rb-gsl.gemspec +29 -0
- data/rdoc/blas.rdoc +3 -3
- data/rdoc/bspline.rdoc +7 -7
- data/rdoc/cheb.rdoc +8 -8
- data/rdoc/cholesky_complex.rdoc +1 -1
- data/rdoc/combi.rdoc +9 -9
- data/rdoc/complex.rdoc +11 -11
- data/rdoc/const.rdoc +20 -20
- data/rdoc/dht.rdoc +6 -6
- data/rdoc/diff.rdoc +6 -6
- data/rdoc/ehandling.rdoc +3 -3
- data/rdoc/eigen.rdoc +11 -11
- data/rdoc/fft.rdoc +26 -26
- data/rdoc/fit.rdoc +18 -18
- data/rdoc/graph.rdoc +2 -2
- data/rdoc/hist.rdoc +16 -16
- data/rdoc/hist2d.rdoc +4 -4
- data/rdoc/hist3d.rdoc +3 -3
- data/rdoc/integration.rdoc +16 -16
- data/rdoc/interp.rdoc +11 -11
- data/rdoc/intro.rdoc +3 -3
- data/rdoc/linalg.rdoc +20 -20
- data/rdoc/linalg_complex.rdoc +1 -1
- data/rdoc/math.rdoc +13 -13
- data/rdoc/matrix.rdoc +22 -22
- data/rdoc/min.rdoc +9 -9
- data/rdoc/monte.rdoc +3 -3
- data/rdoc/multimin.rdoc +12 -12
- data/rdoc/multiroot.rdoc +12 -12
- data/rdoc/narray.rdoc +9 -9
- data/rdoc/ndlinear.rdoc +4 -4
- data/rdoc/nonlinearfit.rdoc +17 -17
- data/rdoc/ntuple.rdoc +5 -5
- data/rdoc/odeiv.rdoc +12 -12
- data/rdoc/perm.rdoc +11 -11
- data/rdoc/poly.rdoc +17 -17
- data/rdoc/qrng.rdoc +9 -9
- data/rdoc/randist.rdoc +10 -10
- data/rdoc/ref.rdoc +49 -49
- data/rdoc/rng.rdoc +9 -9
- data/rdoc/roots.rdoc +12 -12
- data/rdoc/sf.rdoc +35 -35
- data/rdoc/siman.rdoc +3 -3
- data/rdoc/sort.rdoc +6 -6
- data/rdoc/start.rdoc +2 -2
- data/rdoc/stats.rdoc +13 -13
- data/rdoc/sum.rdoc +4 -4
- data/rdoc/tensor.rdoc +3 -3
- data/rdoc/tut.rdoc +1 -1
- data/rdoc/use.rdoc +3 -3
- data/rdoc/vector.rdoc +28 -28
- data/rdoc/vector_complex.rdoc +6 -6
- data/rdoc/wavelet.rdoc +8 -8
- data/test/gsl/blas_test.rb +1 -1
- data/test/gsl/index_test.rb +61 -0
- data/test/gsl/matrix_test.rb +21 -0
- data/uncrustify.cfg +1693 -0
- metadata +212 -243
- data/README +0 -28
- data/ext/gsl/bspline.c +0 -130
- data/ext/gsl/const.c +0 -672
- data/ext/gsl/cqp.c +0 -283
- data/ext/gsl/fcmp.c +0 -64
- data/ext/gsl/fresnel.c +0 -312
- data/ext/gsl/jacobi.c +0 -739
- data/ext/gsl/ool.c +0 -879
- data/ext/gsl/oper_complex_source.c +0 -251
- data/ext/gsl/sf_mathieu.c +0 -238
- data/rdoc/index.rdoc +0 -50
@@ -1,251 +0,0 @@
|
|
1
|
-
/* matrix/oper_complex_source.c
|
2
|
-
*
|
3
|
-
* Copyright (C) 1996, 1997, 1998, 1999, 2000 Brian Gough
|
4
|
-
*
|
5
|
-
* This program is free software; you can redistribute it and/or modify
|
6
|
-
* it under the terms of the GNU General Public License as published by
|
7
|
-
* the Free Software Foundation; either version 2 of the License, or (at
|
8
|
-
* your option) any later version.
|
9
|
-
*
|
10
|
-
* This program is distributed in the hope that it will be useful, but
|
11
|
-
* WITHOUT ANY WARRANTY; without even the implied warranty of
|
12
|
-
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
13
|
-
* General Public License for more details.
|
14
|
-
*
|
15
|
-
* You should have received a copy of the GNU General Public License
|
16
|
-
* along with this program; if not, write to the Free Software
|
17
|
-
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
|
18
|
-
*/
|
19
|
-
|
20
|
-
#ifndef GSL_1_2_LATER
|
21
|
-
#include <config.h>
|
22
|
-
#include <stdlib.h>
|
23
|
-
#include <math.h>
|
24
|
-
#include <gsl/gsl_math.h>
|
25
|
-
#include <gsl/gsl_matrix.h>
|
26
|
-
|
27
|
-
#define BASE_GSL_COMPLEX
|
28
|
-
#include "include/templates_on.h"
|
29
|
-
|
30
|
-
int
|
31
|
-
FUNCTION (gsl_matrix, add) (GSL_TYPE (gsl_matrix) * a,
|
32
|
-
const GSL_TYPE (gsl_matrix) * b)
|
33
|
-
{
|
34
|
-
const size_t M = a->size1;
|
35
|
-
const size_t N = a->size2;
|
36
|
-
|
37
|
-
if (b->size1 != M || b->size2 != N)
|
38
|
-
{
|
39
|
-
GSL_ERROR ("matrices must have same dimensions", GSL_EBADLEN);
|
40
|
-
}
|
41
|
-
else
|
42
|
-
{
|
43
|
-
const size_t tda_a = a->tda;
|
44
|
-
const size_t tda_b = b->tda;
|
45
|
-
|
46
|
-
size_t i, j;
|
47
|
-
|
48
|
-
for (i = 0; i < M; i++)
|
49
|
-
{
|
50
|
-
for (j = 0; j < N; j++)
|
51
|
-
{
|
52
|
-
const size_t aij = 2 * (i * tda_a + j);
|
53
|
-
const size_t bij = 2 * (i * tda_b + j);
|
54
|
-
|
55
|
-
a->data[aij] += b->data[bij];
|
56
|
-
a->data[aij + 1] += b->data[bij + 1];
|
57
|
-
}
|
58
|
-
}
|
59
|
-
|
60
|
-
return GSL_SUCCESS;
|
61
|
-
}
|
62
|
-
}
|
63
|
-
|
64
|
-
int
|
65
|
-
FUNCTION (gsl_matrix, sub) (GSL_TYPE (gsl_matrix) * a,
|
66
|
-
const GSL_TYPE (gsl_matrix) * b)
|
67
|
-
{
|
68
|
-
const size_t M = a->size1;
|
69
|
-
const size_t N = a->size2;
|
70
|
-
|
71
|
-
if (b->size1 != M || b->size2 != N)
|
72
|
-
{
|
73
|
-
GSL_ERROR ("matrices must have same dimensions", GSL_EBADLEN);
|
74
|
-
}
|
75
|
-
else
|
76
|
-
{
|
77
|
-
const size_t tda_a = a->tda;
|
78
|
-
const size_t tda_b = b->tda;
|
79
|
-
|
80
|
-
size_t i, j;
|
81
|
-
|
82
|
-
for (i = 0; i < M; i++)
|
83
|
-
{
|
84
|
-
for (j = 0; j < N; j++)
|
85
|
-
{
|
86
|
-
const size_t aij = 2 * (i * tda_a + j);
|
87
|
-
const size_t bij = 2 * (i * tda_b + j);
|
88
|
-
|
89
|
-
a->data[aij] -= b->data[bij];
|
90
|
-
a->data[aij + 1] -= b->data[bij + 1];
|
91
|
-
}
|
92
|
-
}
|
93
|
-
|
94
|
-
return GSL_SUCCESS;
|
95
|
-
}
|
96
|
-
}
|
97
|
-
|
98
|
-
int
|
99
|
-
FUNCTION (gsl_matrix, mul_elements) (GSL_TYPE (gsl_matrix) * a,
|
100
|
-
const GSL_TYPE (gsl_matrix) * b)
|
101
|
-
{
|
102
|
-
const size_t M = a->size1;
|
103
|
-
const size_t N = a->size2;
|
104
|
-
|
105
|
-
if (b->size1 != M || b->size2 != N)
|
106
|
-
{
|
107
|
-
GSL_ERROR ("matrices must have same dimensions", GSL_EBADLEN);
|
108
|
-
}
|
109
|
-
else
|
110
|
-
{
|
111
|
-
const size_t tda_a = a->tda;
|
112
|
-
const size_t tda_b = b->tda;
|
113
|
-
|
114
|
-
size_t i, j;
|
115
|
-
|
116
|
-
for (i = 0; i < M; i++)
|
117
|
-
{
|
118
|
-
for (j = 0; j < N; j++)
|
119
|
-
{
|
120
|
-
const size_t aij = 2 * (i * tda_a + j);
|
121
|
-
const size_t bij = 2 * (i * tda_b + j);
|
122
|
-
|
123
|
-
ATOMIC ar = a->data[aij];
|
124
|
-
ATOMIC ai = a->data[aij + 1];
|
125
|
-
|
126
|
-
ATOMIC br = b->data[bij];
|
127
|
-
ATOMIC bi = b->data[bij + 1];
|
128
|
-
|
129
|
-
a->data[aij] = ar * br - ai * bi;
|
130
|
-
a->data[aij + 1] = ar * bi + ai * br;
|
131
|
-
}
|
132
|
-
}
|
133
|
-
|
134
|
-
return GSL_SUCCESS;
|
135
|
-
}
|
136
|
-
}
|
137
|
-
|
138
|
-
int
|
139
|
-
FUNCTION (gsl_matrix, div_elements) (GSL_TYPE (gsl_matrix) * a,
|
140
|
-
const GSL_TYPE (gsl_matrix) * b)
|
141
|
-
{
|
142
|
-
const size_t M = a->size1;
|
143
|
-
const size_t N = a->size2;
|
144
|
-
|
145
|
-
if (b->size1 != M || b->size2 != N)
|
146
|
-
{
|
147
|
-
GSL_ERROR ("matrices must have same dimensions", GSL_EBADLEN);
|
148
|
-
}
|
149
|
-
else
|
150
|
-
{
|
151
|
-
const size_t tda_a = a->tda;
|
152
|
-
const size_t tda_b = b->tda;
|
153
|
-
|
154
|
-
size_t i, j;
|
155
|
-
|
156
|
-
for (i = 0; i < M; i++)
|
157
|
-
{
|
158
|
-
for (j = 0; j < N; j++)
|
159
|
-
{
|
160
|
-
const size_t aij = 2 * (i * tda_a + j);
|
161
|
-
const size_t bij = 2 * (i * tda_b + j);
|
162
|
-
|
163
|
-
ATOMIC ar = a->data[aij];
|
164
|
-
ATOMIC ai = a->data[aij + 1];
|
165
|
-
|
166
|
-
ATOMIC br = b->data[bij];
|
167
|
-
ATOMIC bi = b->data[bij + 1];
|
168
|
-
|
169
|
-
ATOMIC s = 1.0 / hypot(br, bi);
|
170
|
-
|
171
|
-
ATOMIC sbr = s * br;
|
172
|
-
ATOMIC sbi = s * bi;
|
173
|
-
|
174
|
-
a->data[aij] = (ar * sbr + ai * sbi) * s;
|
175
|
-
a->data[aij + 1] = (ai * sbr - ar * sbi) * s;
|
176
|
-
}
|
177
|
-
}
|
178
|
-
|
179
|
-
return GSL_SUCCESS;
|
180
|
-
}
|
181
|
-
}
|
182
|
-
|
183
|
-
int FUNCTION (gsl_matrix, scale) (GSL_TYPE (gsl_matrix) * a, const BASE x)
|
184
|
-
{
|
185
|
-
const size_t M = a->size1;
|
186
|
-
const size_t N = a->size2;
|
187
|
-
const size_t tda = a->tda;
|
188
|
-
|
189
|
-
size_t i, j;
|
190
|
-
|
191
|
-
ATOMIC xr = GSL_REAL(x);
|
192
|
-
ATOMIC xi = GSL_IMAG(x);
|
193
|
-
|
194
|
-
for (i = 0; i < M; i++)
|
195
|
-
{
|
196
|
-
for (j = 0; j < N; j++)
|
197
|
-
{
|
198
|
-
const size_t aij = 2 * (i * tda + j);
|
199
|
-
|
200
|
-
ATOMIC ar = a->data[aij];
|
201
|
-
ATOMIC ai = a->data[aij + 1];
|
202
|
-
|
203
|
-
a->data[aij] = ar * xr - ai * xi;
|
204
|
-
a->data[aij + 1] = ar * xi + ai * xr;
|
205
|
-
}
|
206
|
-
}
|
207
|
-
|
208
|
-
return GSL_SUCCESS;
|
209
|
-
}
|
210
|
-
|
211
|
-
int FUNCTION (gsl_matrix, add_constant) (GSL_TYPE (gsl_matrix) * a, const BASE x)
|
212
|
-
{
|
213
|
-
const size_t M = a->size1;
|
214
|
-
const size_t N = a->size2;
|
215
|
-
const size_t tda = a->tda;
|
216
|
-
|
217
|
-
size_t i, j;
|
218
|
-
|
219
|
-
for (i = 0; i < M; i++)
|
220
|
-
{
|
221
|
-
for (j = 0; j < N; j++)
|
222
|
-
{
|
223
|
-
a->data[2 * (i * tda + j)] += GSL_REAL (x);
|
224
|
-
a->data[2 * (i * tda + j) + 1] += GSL_IMAG (x);
|
225
|
-
}
|
226
|
-
}
|
227
|
-
|
228
|
-
return GSL_SUCCESS;
|
229
|
-
}
|
230
|
-
|
231
|
-
|
232
|
-
int FUNCTION (gsl_matrix, add_diagonal) (GSL_TYPE (gsl_matrix) * a, const BASE x)
|
233
|
-
{
|
234
|
-
const size_t M = a->size1;
|
235
|
-
const size_t N = a->size2;
|
236
|
-
const size_t tda = a->tda;
|
237
|
-
const size_t loop_lim = (M < N ? M : N);
|
238
|
-
size_t i;
|
239
|
-
for (i = 0; i < loop_lim; i++)
|
240
|
-
{
|
241
|
-
a->data[2 * (i * tda + i)] += GSL_REAL (x);
|
242
|
-
a->data[2 * (i * tda + i) + 1] += GSL_IMAG (x);
|
243
|
-
}
|
244
|
-
|
245
|
-
return GSL_SUCCESS;
|
246
|
-
}
|
247
|
-
|
248
|
-
|
249
|
-
#include "include/templates_off.h"
|
250
|
-
#undef BASE_GSL_COMPLEX
|
251
|
-
#endif
|
data/ext/gsl/sf_mathieu.c
DELETED
@@ -1,238 +0,0 @@
|
|
1
|
-
#include "include/rb_gsl.h"
|
2
|
-
#ifdef GSL_1_9_LATER
|
3
|
-
|
4
|
-
static VALUE cWorkspace;
|
5
|
-
|
6
|
-
static VALUE rb_gsl_sf_mathieu_alloc(VALUE klass, VALUE n, VALUE q)
|
7
|
-
{
|
8
|
-
gsl_sf_mathieu_workspace *w;
|
9
|
-
w = gsl_sf_mathieu_alloc((size_t) FIX2INT(n), NUM2DBL(q));
|
10
|
-
return Data_Wrap_Struct(klass, 0, gsl_sf_mathieu_free, w);
|
11
|
-
}
|
12
|
-
|
13
|
-
static VALUE sf_mathieu_eval(VALUE order, VALUE qq,
|
14
|
-
int (*f)(int, double, gsl_sf_result*))
|
15
|
-
{
|
16
|
-
gsl_sf_result r;
|
17
|
-
(*f)(FIX2INT(order), NUM2DBL(qq), &r);
|
18
|
-
return rb_float_new(r.val);
|
19
|
-
}
|
20
|
-
|
21
|
-
static VALUE sf_mathieu_eval2(VALUE n1, VALUE n2, VALUE q, VALUE x,
|
22
|
-
int (*f)(int, int, double, double, gsl_sf_result*))
|
23
|
-
{
|
24
|
-
gsl_sf_result r;
|
25
|
-
(*f)(FIX2INT(n1),FIX2INT(n2), NUM2DBL(q), NUM2DBL(x), &r);
|
26
|
-
return rb_float_new(r.val);
|
27
|
-
}
|
28
|
-
|
29
|
-
static VALUE sf_mathieu_array_eval(int argc, VALUE *argv,
|
30
|
-
int (*f)(int, int, double, gsl_sf_mathieu_workspace*, double[]))
|
31
|
-
{
|
32
|
-
gsl_sf_mathieu_workspace *w;
|
33
|
-
gsl_vector *v;
|
34
|
-
int n1, n2;
|
35
|
-
double q;
|
36
|
-
switch (argc) {
|
37
|
-
case 4:
|
38
|
-
if (!rb_obj_is_kind_of(argv[3], cWorkspace)) {
|
39
|
-
rb_raise(rb_eTypeError, "Wrong argument type 3 (%s detected, %s expected)",
|
40
|
-
rb_class2name(CLASS_OF(argv[3])), rb_class2name(cWorkspace));
|
41
|
-
}
|
42
|
-
n1 = FIX2INT(argv[0]);
|
43
|
-
n2 = FIX2INT(argv[1]);
|
44
|
-
q = NUM2DBL(argv[2]);
|
45
|
-
Data_Get_Struct(argv[3], gsl_sf_mathieu_workspace, w);
|
46
|
-
break;
|
47
|
-
default:
|
48
|
-
rb_raise(rb_eArgError, "Wrong number of arguments. (%d for 4)", argc);
|
49
|
-
}
|
50
|
-
v = gsl_vector_alloc(n2 - n1 + 1);
|
51
|
-
(*f)(n1, n2, q, w, v->data);
|
52
|
-
return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
|
53
|
-
}
|
54
|
-
|
55
|
-
static VALUE sf_mathieu_array_eval2(int argc, VALUE *argv,
|
56
|
-
int (*f)(int, int, double, double, gsl_sf_mathieu_workspace*, double[]))
|
57
|
-
{
|
58
|
-
gsl_sf_mathieu_workspace *w;
|
59
|
-
gsl_vector *v;
|
60
|
-
int n1, n2;
|
61
|
-
double q, x;
|
62
|
-
switch (argc) {
|
63
|
-
case 5:
|
64
|
-
if (!rb_obj_is_kind_of(argv[4], cWorkspace)) {
|
65
|
-
rb_raise(rb_eTypeError, "Wrong argument type 4 (%s detected, %s expected)",
|
66
|
-
rb_class2name(CLASS_OF(argv[4])), rb_class2name(cWorkspace));
|
67
|
-
}
|
68
|
-
n1 = FIX2INT(argv[0]);
|
69
|
-
n2 = FIX2INT(argv[1]);
|
70
|
-
q = NUM2DBL(argv[2]);
|
71
|
-
x = NUM2DBL(argv[3]);
|
72
|
-
Data_Get_Struct(argv[4], gsl_sf_mathieu_workspace, w);
|
73
|
-
break;
|
74
|
-
default:
|
75
|
-
rb_raise(rb_eArgError, "Wrong number of arguments. (%d for 5)", argc);
|
76
|
-
}
|
77
|
-
v = gsl_vector_alloc(n2 - n1 + 1);
|
78
|
-
(*f)(n1, n2, q, x, w, v->data);
|
79
|
-
return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
|
80
|
-
}
|
81
|
-
static VALUE sf_mathieu_array_eval3(int argc, VALUE *argv,
|
82
|
-
int (*f)(int, int, int, double, double, gsl_sf_mathieu_workspace*, double[]))
|
83
|
-
{
|
84
|
-
gsl_sf_mathieu_workspace *w;
|
85
|
-
gsl_vector *v;
|
86
|
-
int n1, n2, n3;
|
87
|
-
double q, x;
|
88
|
-
switch (argc) {
|
89
|
-
case 6:
|
90
|
-
if (!rb_obj_is_kind_of(argv[5], cWorkspace)) {
|
91
|
-
rb_raise(rb_eTypeError, "Wrong argument type 5 (%s detected, %s expected)",
|
92
|
-
rb_class2name(CLASS_OF(argv[5])), rb_class2name(cWorkspace));
|
93
|
-
}
|
94
|
-
n1 = FIX2INT(argv[0]);
|
95
|
-
n2 = FIX2INT(argv[1]);
|
96
|
-
n3 = FIX2INT(argv[2]);
|
97
|
-
q = NUM2DBL(argv[3]);
|
98
|
-
x = NUM2DBL(argv[4]);
|
99
|
-
Data_Get_Struct(argv[5], gsl_sf_mathieu_workspace, w);
|
100
|
-
break;
|
101
|
-
default:
|
102
|
-
rb_raise(rb_eArgError, "Wrong number of arguments. (%d for 6)", argc);
|
103
|
-
}
|
104
|
-
v = gsl_vector_alloc(n3 - n2 + 1);
|
105
|
-
(*f)(n1, n2, n3, q, x, w, v->data);
|
106
|
-
return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
|
107
|
-
}
|
108
|
-
static VALUE sf_mathieu_eval_int_double2(VALUE order, VALUE qq, VALUE zz,
|
109
|
-
int (*f)(int, double, double, gsl_sf_result*))
|
110
|
-
{
|
111
|
-
gsl_sf_result r;
|
112
|
-
(*f)(FIX2INT(order), NUM2DBL(qq), NUM2DBL(zz), &r);
|
113
|
-
return rb_float_new(r.val);
|
114
|
-
}
|
115
|
-
static VALUE sf_mathieu_eval_e_int_double2(VALUE order, VALUE qq, VALUE zz,
|
116
|
-
int (*f)(int, double, double, gsl_sf_result*))
|
117
|
-
{
|
118
|
-
gsl_sf_result *r;
|
119
|
-
VALUE val;
|
120
|
-
val = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
|
121
|
-
(*f)(FIX2INT(order), NUM2DBL(qq), NUM2DBL(zz), r);
|
122
|
-
return val;
|
123
|
-
}
|
124
|
-
|
125
|
-
static VALUE sf_mathieu_eval_e_int2_double2(VALUE n1, VALUE n2, VALUE qq, VALUE zz,
|
126
|
-
int (*f)(int, int, double, double, gsl_sf_result*))
|
127
|
-
{
|
128
|
-
gsl_sf_result *r;
|
129
|
-
VALUE val;
|
130
|
-
val = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
|
131
|
-
(*f)(FIX2INT(n1), FIX2INT(n2), NUM2DBL(qq), NUM2DBL(zz), r);
|
132
|
-
return val;
|
133
|
-
}
|
134
|
-
/**********/
|
135
|
-
static VALUE rb_gsl_sf_mathieu_a_e(VALUE module, VALUE order, VALUE qq)
|
136
|
-
{
|
137
|
-
return rb_gsl_sf_eval_e_int_double(gsl_sf_mathieu_a, order, qq);
|
138
|
-
}
|
139
|
-
static VALUE rb_gsl_sf_mathieu_a(VALUE module, VALUE order, VALUE qq)
|
140
|
-
{
|
141
|
-
return sf_mathieu_eval(order, qq, gsl_sf_mathieu_a);
|
142
|
-
}
|
143
|
-
static VALUE rb_gsl_sf_mathieu_a_array(VALUE module, int argc, VALUE *argv)
|
144
|
-
{
|
145
|
-
return sf_mathieu_array_eval(argc, argv, gsl_sf_mathieu_a_array);
|
146
|
-
}
|
147
|
-
static VALUE rb_gsl_sf_mathieu_b_e(VALUE module, VALUE order, VALUE qq)
|
148
|
-
{
|
149
|
-
return rb_gsl_sf_eval_e_int_double(gsl_sf_mathieu_b, order, qq);
|
150
|
-
}
|
151
|
-
static VALUE rb_gsl_sf_mathieu_b(VALUE module, VALUE order, VALUE qq)
|
152
|
-
{
|
153
|
-
return sf_mathieu_eval(order, qq, gsl_sf_mathieu_b);
|
154
|
-
}
|
155
|
-
static VALUE rb_gsl_sf_mathieu_b_array(VALUE module, int argc, VALUE *argv)
|
156
|
-
{
|
157
|
-
return sf_mathieu_array_eval(argc, argv, gsl_sf_mathieu_b_array);
|
158
|
-
}
|
159
|
-
static VALUE rb_gsl_sf_mathieu_ce_e(VALUE module, VALUE order, VALUE qq, VALUE zz)
|
160
|
-
{
|
161
|
-
return sf_mathieu_eval_e_int_double2(order, qq, zz, gsl_sf_mathieu_ce);
|
162
|
-
}
|
163
|
-
static VALUE rb_gsl_sf_mathieu_ce(VALUE module, VALUE order, VALUE qq, VALUE zz)
|
164
|
-
{
|
165
|
-
return sf_mathieu_eval_int_double2(order, qq, zz, gsl_sf_mathieu_ce);
|
166
|
-
}
|
167
|
-
static VALUE rb_gsl_sf_mathieu_ce_array(VALUE module, int argc, VALUE *argv)
|
168
|
-
{
|
169
|
-
return sf_mathieu_array_eval2(argc, argv, gsl_sf_mathieu_ce_array);
|
170
|
-
}
|
171
|
-
static VALUE rb_gsl_sf_mathieu_se_e(VALUE module, VALUE order, VALUE qq, VALUE zz)
|
172
|
-
{
|
173
|
-
return sf_mathieu_eval_e_int_double2(order, qq, zz, gsl_sf_mathieu_se);
|
174
|
-
}
|
175
|
-
static VALUE rb_gsl_sf_mathieu_se(VALUE module, VALUE order, VALUE qq, VALUE zz)
|
176
|
-
{
|
177
|
-
return sf_mathieu_eval_int_double2(order, qq, zz, gsl_sf_mathieu_se);
|
178
|
-
}
|
179
|
-
static VALUE rb_gsl_sf_mathieu_se_array(VALUE module, int argc, VALUE *argv)
|
180
|
-
{
|
181
|
-
return sf_mathieu_array_eval2(argc, argv, gsl_sf_mathieu_se_array);
|
182
|
-
}
|
183
|
-
|
184
|
-
/*****/
|
185
|
-
static VALUE rb_gsl_sf_mathieu_Mc_e(VALUE module, VALUE n1, VALUE n2, VALUE q, VALUE x)
|
186
|
-
{
|
187
|
-
return sf_mathieu_eval_e_int2_double2(n1, n2, q, x, gsl_sf_mathieu_Mc);
|
188
|
-
}
|
189
|
-
static VALUE rb_gsl_sf_mathieu_Mc(VALUE module, VALUE n1, VALUE n2, VALUE q, VALUE x)
|
190
|
-
{
|
191
|
-
return sf_mathieu_eval2(n1, n2, q, x, gsl_sf_mathieu_Mc);
|
192
|
-
}
|
193
|
-
static VALUE rb_gsl_sf_mathieu_Mc_array(VALUE module, int argc, VALUE *argv)
|
194
|
-
{
|
195
|
-
return sf_mathieu_array_eval3(argc, argv, gsl_sf_mathieu_Mc_array);
|
196
|
-
}
|
197
|
-
static VALUE rb_gsl_sf_mathieu_Ms_e(VALUE module, VALUE n1, VALUE n2, VALUE q, VALUE x)
|
198
|
-
{
|
199
|
-
return sf_mathieu_eval_e_int2_double2(n1, n2, q, x, gsl_sf_mathieu_Ms);
|
200
|
-
}
|
201
|
-
static VALUE rb_gsl_sf_mathieu_Ms(VALUE module, VALUE n1, VALUE n2, VALUE q, VALUE x)
|
202
|
-
{
|
203
|
-
return sf_mathieu_eval2(n1, n2, q, x, gsl_sf_mathieu_Ms);
|
204
|
-
}
|
205
|
-
static VALUE rb_gsl_sf_mathieu_Ms_array(VALUE module, int argc, VALUE *argv)
|
206
|
-
{
|
207
|
-
return sf_mathieu_array_eval3(argc, argv, gsl_sf_mathieu_Ms_array);
|
208
|
-
}
|
209
|
-
/*****/
|
210
|
-
void Init_sf_mathieu(VALUE module)
|
211
|
-
{
|
212
|
-
VALUE mMathieu;
|
213
|
-
|
214
|
-
mMathieu = rb_define_module_under(module, "Mathieu");
|
215
|
-
cWorkspace = rb_define_class_under(mMathieu, "Workspace", cGSL_Object);
|
216
|
-
rb_define_singleton_method(cWorkspace, "alloc", rb_gsl_sf_mathieu_alloc, 2);
|
217
|
-
|
218
|
-
rb_define_module_function(module, "mathieu_a", rb_gsl_sf_mathieu_a, 2);
|
219
|
-
rb_define_module_function(module, "mathieu_a_e", rb_gsl_sf_mathieu_a_e, 2);
|
220
|
-
rb_define_module_function(module, "mathieu_a_array", rb_gsl_sf_mathieu_a_array, -1);
|
221
|
-
rb_define_module_function(module, "mathieu_b", rb_gsl_sf_mathieu_b, 2);
|
222
|
-
rb_define_module_function(module, "mathieu_b_e", rb_gsl_sf_mathieu_b_e, 2);
|
223
|
-
rb_define_module_function(module, "mathieu_b_array", rb_gsl_sf_mathieu_b_array, -1);
|
224
|
-
rb_define_module_function(module, "mathieu_ce", rb_gsl_sf_mathieu_ce, 3);
|
225
|
-
rb_define_module_function(module, "mathieu_ce_e", rb_gsl_sf_mathieu_ce_e, 3);
|
226
|
-
rb_define_module_function(module, "mathieu_ce_array", rb_gsl_sf_mathieu_ce_array, -1);
|
227
|
-
rb_define_module_function(module, "mathieu_se", rb_gsl_sf_mathieu_se, 3);
|
228
|
-
rb_define_module_function(module, "mathieu_se_e", rb_gsl_sf_mathieu_se_e, 3);
|
229
|
-
rb_define_module_function(module, "mathieu_se_array", rb_gsl_sf_mathieu_se_array, -1);
|
230
|
-
rb_define_module_function(module, "mathieu_Mc", rb_gsl_sf_mathieu_Mc, 4);
|
231
|
-
rb_define_module_function(module, "mathieu_Mc_e", rb_gsl_sf_mathieu_Mc_e, 4);
|
232
|
-
rb_define_module_function(module, "mathieu_Mc_array", rb_gsl_sf_mathieu_Mc_array, -1);
|
233
|
-
rb_define_module_function(module, "mathieu_Ms", rb_gsl_sf_mathieu_Ms, 4);
|
234
|
-
rb_define_module_function(module, "mathieu_Ms_e", rb_gsl_sf_mathieu_Ms_e, 4);
|
235
|
-
rb_define_module_function(module, "mathieu_Ms_array", rb_gsl_sf_mathieu_Ms_array, -1);
|
236
|
-
}
|
237
|
-
|
238
|
-
#endif
|