3 #include
"/usr/local/include/gsl/gsl_matrix.h"
4 #include
"/usr/local/include/gsl/gsl_matrix_double.h"
5 #include
"/usr/local/include/gsl/gsl_matrix_int.h"
6 #include
"/usr/local/include/gsl/gsl_matrix_complex_double.h"
7 #include
"/usr/local/include/gsl/gsl_matrix_char.h"
10 %include
"/usr/local/include/gsl/gsl_matrix.h"
11 %include
"/usr/local/include/gsl/gsl_matrix_double.h"
12 %include
"/usr/local/include/gsl/gsl_matrix_int.h"
13 %include
"/usr/local/include/gsl/gsl_matrix_complex_double.h"
14 %include
"/usr/local/include/gsl/gsl_matrix_char.h"
16 %typemap
(in
) gsl_vector const
* {
17 printf
("gsl_vector_view const *\n");
21 FILE *fopen
(char
*, char
*);
26 no warnings 'redefine'
;
28 @EXPORT_OK
= qw
/fopen fclose
29 gsl_matrix_alloc gsl_matrix_calloc gsl_matrix_alloc_from_block
30 gsl_matrix_alloc_from_matrix gsl_vector_alloc_row_from_matrix
31 gsl_vector_alloc_col_from_matrix gsl_matrix_free gsl_matrix_submatrix
32 gsl_matrix_row gsl_matrix_column gsl_matrix_diagonal
33 gsl_matrix_subdiagonal gsl_matrix_superdiagonal gsl_matrix_subrow
34 gsl_matrix_subcolumn gsl_matrix_view_array
35 gsl_matrix_view_array_with_tda gsl_matrix_view_vector
36 gsl_matrix_view_vector_with_tda gsl_matrix_const_submatrix
37 gsl_matrix_const_row gsl_matrix_const_column gsl_matrix_const_diagonal
38 gsl_matrix_const_subdiagonal gsl_matrix_const_superdiagonal
39 gsl_matrix_const_subrow gsl_matrix_const_subcolumn
40 gsl_matrix_const_view_array gsl_matrix_const_view_array_with_tda
41 gsl_matrix_const_view_vector gsl_matrix_const_view_vector_with_tda
42 gsl_matrix_get gsl_matrix_set gsl_matrix_ptr gsl_matrix_const_ptr
43 gsl_matrix_set_zero gsl_matrix_set_identity gsl_matrix_set_all
44 gsl_matrix_fread gsl_matrix_fwrite gsl_matrix_fscanf gsl_matrix_fprintf
45 gsl_matrix_memcpy gsl_matrix_swap gsl_matrix_swap_rows
46 gsl_matrix_swap_columns gsl_matrix_swap_rowcol gsl_matrix_transpose
47 gsl_matrix_transpose_memcpy gsl_matrix_max gsl_matrix_min
48 gsl_matrix_minmax gsl_matrix_max_index gsl_matrix_min_index
49 gsl_matrix_minmax_index gsl_matrix_isnull gsl_matrix_ispos
50 gsl_matrix_isneg gsl_matrix_isnonneg gsl_matrix_add gsl_matrix_sub
51 gsl_matrix_mul_elements gsl_matrix_div_elements gsl_matrix_scale
52 gsl_matrix_add_constant gsl_matrix_add_diagonal
53 gsl_matrix_char_alloc gsl_matrix_char_calloc gsl_matrix_char_alloc_from_block
54 gsl_matrix_char_alloc_from_matrix gsl_vector_char_alloc_row_from_matrix gsl_vector_char_alloc_col_from_matrix
55 gsl_matrix_char_free gsl_matrix_char_submatrix
56 gsl_matrix_char_row gsl_matrix_char_column
57 gsl_matrix_char_diagonal gsl_matrix_char_subdiagonal gsl_matrix_char_superdiagonal
58 gsl_matrix_char_subrow gsl_matrix_char_subcolumn gsl_matrix_char_view_array
59 gsl_matrix_char_view_array_with_tda gsl_matrix_char_view_vector gsl_matrix_char_view_vector_with_tda
60 gsl_matrix_char_const_submatrix gsl_matrix_char_const_row gsl_matrix_char_const_column
61 gsl_matrix_char_const_diagonal gsl_matrix_char_const_subdiagonal gsl_matrix_char_const_superdiagonal
62 gsl_matrix_char_const_subrow gsl_matrix_char_const_subcolumn gsl_matrix_char_const_view_array
63 gsl_matrix_char_const_view_array_with_tda gsl_matrix_char_const_view_vector gsl_matrix_char_const_view_vector_with_tda
64 gsl_matrix_char_get gsl_matrix_char_set gsl_matrix_char_ptr gsl_matrix_char_const_ptr
65 gsl_matrix_char_set_zero gsl_matrix_char_set_identity
66 gsl_matrix_char_set_all gsl_matrix_char_fread
67 gsl_matrix_char_fwrite gsl_matrix_char_fscanf gsl_matrix_char_fprintf
68 gsl_matrix_char_memcpy gsl_matrix_char_swap
69 gsl_matrix_char_swap_rows gsl_matrix_char_swap_columns
70 gsl_matrix_char_swap_rowcol gsl_matrix_char_transpose gsl_matrix_char_transpose_memcpy
71 gsl_matrix_char_max gsl_matrix_char_min
72 gsl_matrix_char_minmax gsl_matrix_char_max_index
73 gsl_matrix_char_min_index gsl_matrix_char_minmax_index
74 gsl_matrix_char_isnull gsl_matrix_char_ispos gsl_matrix_char_isneg
75 gsl_matrix_char_isnonneg gsl_matrix_char_add
76 gsl_matrix_char_sub gsl_matrix_char_mul_elements gsl_matrix_char_div_elements
77 gsl_matrix_char_scale gsl_matrix_char_add_constant gsl_matrix_char_add_diagonal
78 gsl_matrix_int_alloc gsl_matrix_int_calloc gsl_matrix_int_alloc_from_block
79 gsl_matrix_int_alloc_from_matrix gsl_vector_int_alloc_row_from_matrix gsl_vector_int_alloc_col_from_matrix
80 gsl_matrix_int_free gsl_matrix_int_submatrix gsl_matrix_int_row
81 gsl_matrix_int_column gsl_matrix_int_diagonal gsl_matrix_int_subdiagonal
82 gsl_matrix_int_superdiagonal gsl_matrix_int_subrow gsl_matrix_int_subcolumn gsl_matrix_int_view_array
83 gsl_matrix_int_view_array_with_tda gsl_matrix_int_view_vector gsl_matrix_int_view_vector_with_tda
84 gsl_matrix_int_const_submatrix gsl_matrix_int_const_row gsl_matrix_int_const_column
85 gsl_matrix_int_const_diagonal gsl_matrix_int_const_subdiagonal gsl_matrix_int_const_superdiagonal
86 gsl_matrix_int_const_subrow gsl_matrix_int_const_subcolumn gsl_matrix_int_const_view_array
87 gsl_matrix_int_const_view_array_with_tda gsl_matrix_int_const_view_vector gsl_matrix_int_const_view_vector_with_tda
88 gsl_matrix_int_get gsl_matrix_int_set
89 gsl_matrix_int_ptr gsl_matrix_int_const_ptr
90 gsl_matrix_int_set_zero gsl_matrix_int_set_identity gsl_matrix_int_set_all
91 gsl_matrix_int_fread gsl_matrix_int_fwrite
92 gsl_matrix_int_fscanf gsl_matrix_int_fprintf
93 gsl_matrix_int_memcpy gsl_matrix_int_swap
94 gsl_matrix_int_swap_rows gsl_matrix_int_swap_columns gsl_matrix_int_swap_rowcol
95 gsl_matrix_int_transpose gsl_matrix_int_transpose_memcpy
96 gsl_matrix_int_max gsl_matrix_int_min gsl_matrix_int_minmax
97 gsl_matrix_int_max_index gsl_matrix_int_min_index
98 gsl_matrix_int_minmax_index gsl_matrix_int_isnull
99 gsl_matrix_int_ispos gsl_matrix_int_isneg gsl_matrix_int_isnonneg
100 gsl_matrix_int_add gsl_matrix_int_sub
101 gsl_matrix_int_mul_elements gsl_matrix_int_div_elements gsl_matrix_int_scale
102 gsl_matrix_int_add_constant gsl_matrix_int_add_diagonal
103 gsl_matrix_get_row gsl_matrix_get_col gsl_matrix_set_row gsl_matrix_set_col
105 %EXPORT_TAGS
= ( all
=> [ @EXPORT_OK
], char
=> [gsl_matrix_char_alloc
, gsl_matrix_char_calloc
, gsl_matrix_char_alloc_from_block
, gsl_matrix_char_alloc_from_matrix
, gsl_vector_char_alloc_row_from_matrix
, gsl_vector_char_alloc_col_from_matrix
, gsl_matrix_char_free
, gsl_matrix_char_submatrix
, gsl_matrix_char_row
, gsl_matrix_char_column
, gsl_matrix_char_diagonal
, gsl_matrix_char_subdiagonal
, gsl_matrix_char_superdiagonal
, gsl_matrix_char_subrow
, gsl_matrix_char_subcolumn
, gsl_matrix_char_view_array
, gsl_matrix_char_view_array_with_tda
, gsl_matrix_char_view_vector
, gsl_matrix_char_view_vector_with_tda
, gsl_matrix_char_const_submatrix
, gsl_matrix_char_const_row
, gsl_matrix_char_const_column
, gsl_matrix_char_const_diagonal
, gsl_matrix_char_const_subdiagonal
, gsl_matrix_char_const_superdiagonal
, gsl_matrix_char_const_subrow
, gsl_matrix_char_const_subcolumn
, gsl_matrix_char_const_view_array
, gsl_matrix_char_const_view_array_with_tda
, gsl_matrix_char_const_view_vector
, gsl_matrix_char_const_view_vector_with_tda
, gsl_matrix_char_get
, gsl_matrix_char_set
, gsl_matrix_char_ptr
, gsl_matrix_char_const_ptr
, gsl_matrix_char_set_zero
, gsl_matrix_char_set_identity
, gsl_matrix_char_set_all
, gsl_matrix_char_fread
, gsl_matrix_char_fwrite
, gsl_matrix_char_fscanf
, gsl_matrix_char_fprintf
, gsl_matrix_char_memcpy
, gsl_matrix_char_swap
, gsl_matrix_char_swap_rows
, gsl_matrix_char_swap_columns
, gsl_matrix_char_swap_rowcol
, gsl_matrix_char_transpose
, gsl_matrix_char_transpose_memcpy
, gsl_matrix_char_max
, gsl_matrix_char_min
, gsl_matrix_char_minmax
, gsl_matrix_char_max_index
, gsl_matrix_char_min_index
, gsl_matrix_char_minmax_index
, gsl_matrix_char_isnull
, gsl_matrix_char_ispos
, gsl_matrix_char_isneg
, gsl_matrix_char_isnonneg
, gsl_matrix_char_add
, gsl_matrix_char_sub
, gsl_matrix_char_mul_elements
, gsl_matrix_char_div_elements
, gsl_matrix_char_scale
, gsl_matrix_char_add_constant
, gsl_matrix_char_add_diagonal
], double
=> [ gsl_matrix_alloc
, gsl_matrix_calloc
, gsl_matrix_alloc_from_block
, gsl_matrix_alloc_from_matrix
, gsl_vector_alloc_row_from_matrix
, gsl_vector_alloc_col_from_matrix
, gsl_matrix_free
, gsl_matrix_submatrix
, gsl_matrix_row
, gsl_matrix_column
, gsl_matrix_diagonal
, gsl_matrix_subdiagonal
, gsl_matrix_superdiagonal
, gsl_matrix_subrow
, gsl_matrix_subcolumn
, gsl_matrix_view_array
, gsl_matrix_view_array_with_tda
, gsl_matrix_view_vector
, gsl_matrix_view_vector_with_tda
, gsl_matrix_const_submatrix
, gsl_matrix_const_row
, gsl_matrix_const_column
, gsl_matrix_const_diagonal
, gsl_matrix_const_subdiagonal
, gsl_matrix_const_superdiagonal
, gsl_matrix_const_subrow
, gsl_matrix_const_subcolumn
, gsl_matrix_const_view_array
, gsl_matrix_const_view_array_with_tda
, gsl_matrix_const_view_vector
, gsl_matrix_const_view_vector_with_tda
, gsl_matrix_get
, gsl_matrix_set
, gsl_matrix_ptr
, gsl_matrix_const_ptr
, gsl_matrix_set_zero
, gsl_matrix_set_identity
, gsl_matrix_set_all
, gsl_matrix_fread
, gsl_matrix_fwrite
, gsl_matrix_fscanf
, gsl_matrix_fprintf
, gsl_matrix_memcpy
, gsl_matrix_swap
, gsl_matrix_swap_rows
, gsl_matrix_swap_columns
, gsl_matrix_swap_rowcol
, gsl_matrix_transpose
, gsl_matrix_transpose_memcpy
, gsl_matrix_max
, gsl_matrix_minmax
, gsl_matrix_max_index
, gsl_matrix_min_index
, gsl_matrix_minmax_index
, gsl_matrix_isnull
, gsl_matrix_ispos
, gsl_matrix_isneg
, gsl_matrix_isnonneg
, gsl_matrix_add
, gsl_matrix_mul_elements
, gsl_matrix_div_elements
, gsl_matrix_scale
, gsl_matrix_add_constant
, gsl_matrix_add_diagonal
], int
=> [gsl_matrix_int_alloc
, gsl_matrix_int_alloc_from_matrix
, gsl_matrix_int_free
, gsl_matrix_int_column
, gsl_matrix_int_superdiagonal
, gsl_matrix_int_view_array_with_tda
, gsl_matrix_int_const_submatrix
, gsl_matrix_int_const_diagonal
, gsl_matrix_int_const_subrow
, gsl_matrix_int_const_view_array_with_tda
, gsl_matrix_int_get
, gsl_matrix_int_ptr
, gsl_matrix_int_set_zero
, gsl_matrix_int_fread
, gsl_matrix_int_fscanf
, gsl_matrix_int_memcpy
, gsl_matrix_int_swap_rows
, gsl_matrix_int_transpose
, gsl_matrix_int_max
, gsl_matrix_int_max_index
, gsl_matrix_int_minmax_index
, gsl_matrix_int_ispos
, gsl_matrix_int_add
, gsl_matrix_int_mul_elements
, gsl_matrix_int_add_constant
]);
112 Mathematical functions concerning Matrices
116 use Math
::GSL
::Matrix qw
/:all
/;
120 Here is a list of all the functions included in this module
:
122 gsl_matrix_alloc gsl_matrix_calloc gsl_matrix_alloc_from_block
123 gsl_matrix_alloc_from_matrix gsl_vector_alloc_row_from_matrix
124 gsl_vector_alloc_col_from_matrix gsl_matrix_free gsl_matrix_submatrix
125 gsl_matrix_row gsl_matrix_column gsl_matrix_diagonal
126 gsl_matrix_subdiagonal gsl_matrix_superdiagonal gsl_matrix_subrow
127 gsl_matrix_subcolumn gsl_matrix_view_array
128 gsl_matrix_view_array_with_tda gsl_matrix_view_vector
129 gsl_matrix_view_vector_with_tda gsl_matrix_const_submatrix
130 gsl_matrix_const_row gsl_matrix_const_column gsl_matrix_const_diagonal
131 gsl_matrix_const_subdiagonal gsl_matrix_const_superdiagonal
132 gsl_matrix_const_subrow gsl_matrix_const_subcolumn
133 gsl_matrix_const_view_array gsl_matrix_const_view_array_with_tda
134 gsl_matrix_const_view_vector gsl_matrix_const_view_vector_with_tda
135 gsl_matrix_get gsl_matrix_set gsl_matrix_ptr gsl_matrix_const_ptr
136 gsl_matrix_set_zero gsl_matrix_set_identity gsl_matrix_set_all
137 gsl_matrix_fread gsl_matrix_fwrite gsl_matrix_fscanf gsl_matrix_fprintf
138 gsl_matrix_memcpy gsl_matrix_swap gsl_matrix_swap_rows
139 gsl_matrix_swap_columns gsl_matrix_swap_rowcol gsl_matrix_transpose
140 gsl_matrix_transpose_memcpy gsl_matrix_max gsl_matrix_min
141 gsl_matrix_minmax gsl_matrix_max_index gsl_matrix_min_index
142 gsl_matrix_minmax_index gsl_matrix_isnull gsl_matrix_ispos
143 gsl_matrix_isneg gsl_matrix_isnonneg gsl_matrix_add gsl_matrix_sub
144 gsl_matrix_mul_elements gsl_matrix_div_elements gsl_matrix_scale
145 gsl_matrix_add_constant gsl_matrix_add_diagonal
146 gsl_matrix_char_alloc gsl_matrix_char_calloc gsl_matrix_char_alloc_from_block
147 gsl_matrix_char_alloc_from_matrix gsl_vector_char_alloc_row_from_matrix gsl_vector_char_alloc_col_from_matrix
148 gsl_matrix_char_free gsl_matrix_char_submatrix
149 gsl_matrix_char_row gsl_matrix_char_column
150 gsl_matrix_char_diagonal gsl_matrix_char_subdiagonal gsl_matrix_char_superdiagonal
151 gsl_matrix_char_subrow gsl_matrix_char_subcolumn gsl_matrix_char_view_array
152 gsl_matrix_char_view_array_with_tda gsl_matrix_char_view_vector gsl_matrix_char_view_vector_with_tda
153 gsl_matrix_char_const_submatrix gsl_matrix_char_const_row gsl_matrix_char_const_column
154 gsl_matrix_char_const_diagonal gsl_matrix_char_const_subdiagonal gsl_matrix_char_const_superdiagonal
155 gsl_matrix_char_const_subrow gsl_matrix_char_const_subcolumn gsl_matrix_char_const_view_array
156 gsl_matrix_char_const_view_array_with_tda gsl_matrix_char_const_view_vector gsl_matrix_char_const_view_vector_with_tda
157 gsl_matrix_char_get gsl_matrix_char_set gsl_matrix_char_ptr gsl_matrix_char_const_ptr
158 gsl_matrix_char_set_zero gsl_matrix_char_set_identity
159 gsl_matrix_char_set_all gsl_matrix_char_fread
160 gsl_matrix_char_fwrite gsl_matrix_char_fscanf gsl_matrix_char_fprintf
161 gsl_matrix_char_memcpy gsl_matrix_char_swap
162 gsl_matrix_char_swap_rows gsl_matrix_char_swap_columns
163 gsl_matrix_char_swap_rowcol gsl_matrix_char_transpose gsl_matrix_char_transpose_memcpy
164 gsl_matrix_char_max gsl_matrix_char_min
165 gsl_matrix_char_minmax gsl_matrix_char_max_index
166 gsl_matrix_char_min_index gsl_matrix_char_minmax_index
167 gsl_matrix_char_isnull gsl_matrix_char_ispos gsl_matrix_char_isneg
168 gsl_matrix_char_isnonneg gsl_matrix_char_add
169 gsl_matrix_char_sub gsl_matrix_char_mul_elements gsl_matrix_char_div_elements
170 gsl_matrix_char_scale gsl_matrix_char_add_constant gsl_matrix_char_add_diagonal
171 gsl_matrix_int_alloc gsl_matrix_int_calloc gsl_matrix_int_alloc_from_block
172 gsl_matrix_int_alloc_from_matrix gsl_vector_int_alloc_row_from_matrix gsl_vector_int_alloc_col_from_matrix
173 gsl_matrix_int_free gsl_matrix_int_submatrix gsl_matrix_int_row
174 gsl_matrix_int_column gsl_matrix_int_diagonal gsl_matrix_int_subdiagonal
175 gsl_matrix_int_superdiagonal gsl_matrix_int_subrow gsl_matrix_int_subcolumn gsl_matrix_int_view_array
176 gsl_matrix_int_view_array_with_tda gsl_matrix_int_view_vector gsl_matrix_int_view_vector_with_tda
177 gsl_matrix_int_const_submatrix gsl_matrix_int_const_row gsl_matrix_int_const_column
178 gsl_matrix_int_const_diagonal gsl_matrix_int_const_subdiagonal gsl_matrix_int_const_superdiagonal
179 gsl_matrix_int_const_subrow gsl_matrix_int_const_subcolumn gsl_matrix_int_const_view_array
180 gsl_matrix_int_const_view_array_with_tda gsl_matrix_int_const_view_vector gsl_matrix_int_const_view_vector_with_tda
181 gsl_matrix_int_get gsl_matrix_int_set
182 gsl_matrix_int_ptr gsl_matrix_int_const_ptr
183 gsl_matrix_int_set_zero gsl_matrix_int_set_identity gsl_matrix_int_set_all
184 gsl_matrix_int_fread gsl_matrix_int_fwrite
185 gsl_matrix_int_fscanf gsl_matrix_int_fprintf
186 gsl_matrix_int_memcpy gsl_matrix_int_swap
187 gsl_matrix_int_swap_rows gsl_matrix_int_swap_columns gsl_matrix_int_swap_rowcol
188 gsl_matrix_int_transpose gsl_matrix_int_transpose_memcpy
189 gsl_matrix_int_max gsl_matrix_int_min gsl_matrix_int_minmax
190 gsl_matrix_int_max_index gsl_matrix_int_min_index
191 gsl_matrix_int_minmax_index gsl_matrix_int_isnull
192 gsl_matrix_int_ispos gsl_matrix_int_isneg gsl_matrix_int_isnonneg
193 gsl_matrix_int_add gsl_matrix_int_sub
194 gsl_matrix_int_mul_elements gsl_matrix_int_div_elements gsl_matrix_int_scale
195 gsl_matrix_int_add_constant gsl_matrix_int_add_diagonal
196 gsl_matrix_get_row gsl_matrix_get_col gsl_matrix_set_row gsl_matrix_set_col
198 You have to add the functions you want to use inside the qw
/put_funtion_here
/. You can also write use Math
::GSL
::PowInt qw
/:name_of_tag
/ to use all avaible functions of the module. Other tags are also avaible
, here is a complete list of all tags for this module
:
205 For more informations on the functions
, we refer you to the GSL offcial documentation
: http
://www.gnu.org
/software
/gsl
/manual
/html_node
/
206 Tip
: search on google
: site
:http
://www.gnu.org
/software
/gsl
/manual
/html_node
/name_of_the_function_you_want
210 Most of the examples from this section are perl versions of the examples there
: http
://www.gnu.org
/software
/gsl
/manual
/html_node
/Example-programs-for-matrices.html
213 The program below shows how to allocate
, initialize and read from a matrix using the functions gsl_matrix_alloc
, gsl_matrix_set and gsl_matrix_get.
215 use Math
::GSL
::Matrix qw
/:all
/;
216 $m
= gsl_matrix_alloc
(10,3);
217 for
($i
= 0; $i
< 10; $i
++){
218 for
($j
= 0; $j
<3; $j
++){
219 gsl_matrix_set
($m
, $i
, $j
, 0.23 + 100*$i
+ $j
);
222 for
($i
= 0; $i
< 100; $i
++){ # OUT
OF RANGE
ERROR
223 for
($j
=0; $j
<3; $j
++){
224 print
"m($i, $j) = " . gsl_matrix_get
($m
, $i
, $j
) .
"\n";
227 gsl_matrix_free
($m
);
230 use Math
::GSL
::Matrix qw
/:all
/;
232 $m
= gsl_matrix_alloc
(100, 100);
233 $a
= gsl_matrix_alloc
(100, 100);
235 for
($i
= 0; $i
< 100; $i
++){
236 for
($j
= 0; $j
< 100; $j
++){
237 gsl_matrix_set
($m
, $i
, $j
, 0.23 + $i
+ $j
);
241 The next program shows how to write a matrix to a file.
243 $out
= fopen
("test.dat", "wb");
244 gsl_matrix_fwrite
($out
, $m
);
247 $in
= fopen
("test.dat", "rb");
248 gsl_matrix_fread
($in
, $a
);
252 for
($i
= 0; $i
< 100; $i
++){
253 for
($j
= 0; $j
< 100; $j
++){
254 $mij
= gsl_matrix_get
($m
, $i
, $j
);
255 $aij
= gsl_matrix_get
($a
, $i
, $j
);
256 if
($mij
!= $aij
){ $k
++ };
263 print
"differences = $k (should be zero)\n";
267 Jonathan Leto
<jonathan@leto.net
> and Thierry Moisan
<thierry.moisan@gmail.com
>
269 =head1 COPYRIGHT
AND LICENSE
271 Copyright
(C
) 2008 Jonathan Leto and Thierry Moisan
273 This program is free software
; you can redistribute it and
/or modify it
274 under the same terms as Perl itself.