Matrix.i

   1 %module Matrix
   2
   3 %apply int *OUTPUT { size_t *imin, size_t *imax, size_t *jmin, size_t *jmax };
   4
   5 %apply double *OUTPUT { double * min_out, double * max_out };
   6 %{
   7     #include "/usr/local/include/gsl/gsl_matrix.h"
   8     #include "/usr/local/include/gsl/gsl_vector_double.h"
   9     #include "/usr/local/include/gsl/gsl_matrix_double.h"
  10     #include "/usr/local/include/gsl/gsl_matrix_int.h"
  11     #include "/usr/local/include/gsl/gsl_matrix_complex_double.h"
  12     #include "/usr/local/include/gsl/gsl_matrix_char.h"
  13 %}
  14
  15 %include "/usr/local/include/gsl/gsl_matrix.h"
  16 %include "/usr/local/include/gsl/gsl_vector_double.h"
  17 %include "/usr/local/include/gsl/gsl_matrix_double.h"
  18 %include "/usr/local/include/gsl/gsl_matrix_int.h"
  19 %include "/usr/local/include/gsl/gsl_matrix_complex_double.h"
  20 %include "/usr/local/include/gsl/gsl_matrix_char.h"
  21
  22
  23 %perlcode %{
  24
  25 no warnings 'redefine';
  26
  27 @EXPORT_OK = qw/fopen fclose
  28         gsl_matrix_alloc gsl_matrix_calloc gsl_matrix_alloc_from_block
  29         gsl_matrix_alloc_from_matrix gsl_vector_alloc_row_from_matrix
  30         gsl_vector_alloc_col_from_matrix gsl_matrix_free gsl_matrix_submatrix
  31         gsl_matrix_row gsl_matrix_column gsl_matrix_diagonal
  32         gsl_matrix_subdiagonal gsl_matrix_superdiagonal gsl_matrix_subrow
  33         gsl_matrix_subcolumn gsl_matrix_view_array
  34         gsl_matrix_view_array_with_tda gsl_matrix_view_vector
  35         gsl_matrix_view_vector_with_tda gsl_matrix_const_submatrix
  36         gsl_matrix_const_row gsl_matrix_const_column gsl_matrix_const_diagonal
  37         gsl_matrix_const_subdiagonal gsl_matrix_const_superdiagonal
  38         gsl_matrix_const_subrow gsl_matrix_const_subcolumn
  39         gsl_matrix_const_view_array gsl_matrix_const_view_array_with_tda
  40         gsl_matrix_const_view_vector gsl_matrix_const_view_vector_with_tda
  41         gsl_matrix_get gsl_matrix_set gsl_matrix_ptr gsl_matrix_const_ptr
  42         gsl_matrix_set_zero gsl_matrix_set_identity gsl_matrix_set_all
  43         gsl_matrix_fread gsl_matrix_fwrite gsl_matrix_fscanf gsl_matrix_fprintf
  44         gsl_matrix_memcpy gsl_matrix_swap gsl_matrix_swap_rows
  45         gsl_matrix_swap_columns gsl_matrix_swap_rowcol gsl_matrix_transpose
  46         gsl_matrix_transpose_memcpy gsl_matrix_max gsl_matrix_min
  47         gsl_matrix_minmax gsl_matrix_max_index gsl_matrix_min_index
  48         gsl_matrix_minmax_index gsl_matrix_isnull gsl_matrix_ispos
  49         gsl_matrix_isneg gsl_matrix_isnonneg gsl_matrix_add gsl_matrix_sub
  50         gsl_matrix_mul_elements gsl_matrix_div_elements gsl_matrix_scale
  51         gsl_matrix_add_constant gsl_matrix_add_diagonal
  52         gsl_matrix_char_alloc gsl_matrix_char_calloc  gsl_matrix_char_alloc_from_block
  53         gsl_matrix_char_alloc_from_matrix gsl_vector_char_alloc_row_from_matrix gsl_vector_char_alloc_col_from_matrix
  54         gsl_matrix_char_free gsl_matrix_char_submatrix
  55         gsl_matrix_char_row  gsl_matrix_char_column
  56         gsl_matrix_char_diagonal gsl_matrix_char_subdiagonal gsl_matrix_char_superdiagonal
  57         gsl_matrix_char_subrow gsl_matrix_char_subcolumn gsl_matrix_char_view_array
  58         gsl_matrix_char_view_array_with_tda gsl_matrix_char_view_vector gsl_matrix_char_view_vector_with_tda
  59         gsl_matrix_char_const_submatrix gsl_matrix_char_const_row gsl_matrix_char_const_column
  60         gsl_matrix_char_const_diagonal gsl_matrix_char_const_subdiagonal gsl_matrix_char_const_superdiagonal
  61         gsl_matrix_char_const_subrow gsl_matrix_char_const_subcolumn gsl_matrix_char_const_view_array
  62         gsl_matrix_char_const_view_array_with_tda gsl_matrix_char_const_view_vector gsl_matrix_char_const_view_vector_with_tda
  63         gsl_matrix_char_get gsl_matrix_char_set gsl_matrix_char_ptr gsl_matrix_char_const_ptr
  64         gsl_matrix_char_set_zero gsl_matrix_char_set_identity
  65         gsl_matrix_char_set_all  gsl_matrix_char_fread
  66         gsl_matrix_char_fwrite gsl_matrix_char_fscanf gsl_matrix_char_fprintf
  67         gsl_matrix_char_memcpy gsl_matrix_char_swap
  68         gsl_matrix_char_swap_rows gsl_matrix_char_swap_columns
  69         gsl_matrix_char_swap_rowcol gsl_matrix_char_transpose gsl_matrix_char_transpose_memcpy
  70         gsl_matrix_char_max gsl_matrix_char_min
  71         gsl_matrix_char_minmax  gsl_matrix_char_max_index
  72         gsl_matrix_char_min_index gsl_matrix_char_minmax_index
  73         gsl_matrix_char_isnull gsl_matrix_char_ispos gsl_matrix_char_isneg
  74         gsl_matrix_char_isnonneg  gsl_matrix_char_add
  75         gsl_matrix_char_sub gsl_matrix_char_mul_elements gsl_matrix_char_div_elements
  76         gsl_matrix_char_scale gsl_matrix_char_add_constant gsl_matrix_char_add_diagonal
  77         gsl_matrix_int_alloc gsl_matrix_int_calloc gsl_matrix_int_alloc_from_block
  78         gsl_matrix_int_alloc_from_matrix gsl_vector_int_alloc_row_from_matrix gsl_vector_int_alloc_col_from_matrix
  79         gsl_matrix_int_free gsl_matrix_int_submatrix gsl_matrix_int_row
  80         gsl_matrix_int_column gsl_matrix_int_diagonal gsl_matrix_int_subdiagonal
  81         gsl_matrix_int_superdiagonal gsl_matrix_int_subrow gsl_matrix_int_subcolumn gsl_matrix_int_view_array
  82         gsl_matrix_int_view_array_with_tda gsl_matrix_int_view_vector gsl_matrix_int_view_vector_with_tda
  83         gsl_matrix_int_const_submatrix gsl_matrix_int_const_row gsl_matrix_int_const_column
  84         gsl_matrix_int_const_diagonal gsl_matrix_int_const_subdiagonal gsl_matrix_int_const_superdiagonal
  85         gsl_matrix_int_const_subrow gsl_matrix_int_const_subcolumn gsl_matrix_int_const_view_array
  86         gsl_matrix_int_const_view_array_with_tda gsl_matrix_int_const_view_vector gsl_matrix_int_const_view_vector_with_tda
  87         gsl_matrix_int_get gsl_matrix_int_set
  88         gsl_matrix_int_ptr gsl_matrix_int_const_ptr
  89         gsl_matrix_int_set_zero gsl_matrix_int_set_identity gsl_matrix_int_set_all
  90         gsl_matrix_int_fread gsl_matrix_int_fwrite
  91         gsl_matrix_int_fscanf gsl_matrix_int_fprintf
  92         gsl_matrix_int_memcpy gsl_matrix_int_swap
  93         gsl_matrix_int_swap_rows gsl_matrix_int_swap_columns gsl_matrix_int_swap_rowcol
  94         gsl_matrix_int_transpose gsl_matrix_int_transpose_memcpy
  95         gsl_matrix_int_max gsl_matrix_int_min gsl_matrix_int_minmax
  96         gsl_matrix_int_max_index gsl_matrix_int_min_index
  97         gsl_matrix_int_minmax_index  gsl_matrix_int_isnull
  98         gsl_matrix_int_ispos gsl_matrix_int_isneg gsl_matrix_int_isnonneg
  99         gsl_matrix_int_add gsl_matrix_int_sub
 100         gsl_matrix_int_mul_elements gsl_matrix_int_div_elements gsl_matrix_int_scale
 101         gsl_matrix_int_add_constant gsl_matrix_int_add_diagonal
 102         gsl_matrix_get_row gsl_matrix_get_col gsl_matrix_set_row gsl_matrix_set_col
 103         gsl_matrix_complex_alloc gsl_matrix_complex_calloc gsl_matrix_complex_alloc_from_block
 104         gsl_matrix_complex_alloc_from_matrix gsl_vector_complex_alloc_row_from_matrix gsl_vector_complex_alloc_col_from_matrix
 105         gsl_matrix_complex_free gsl_matrix_complex_submatrix gsl_matrix_complex_row
 106         gsl_matrix_complex_column gsl_matrix_complex_diagonal gsl_matrix_complex_subdiagonal
 107         gsl_matrix_complex_superdiagonal gsl_matrix_complex_subrow gsl_matrix_complex_subcolumn
 108         gsl_matrix_complex_view_array gsl_matrix_complex_view_array_with_tda gsl_matrix_complex_view_vector
 109         gsl_matrix_complex_view_vector_with_tda gsl_matrix_complex_const_submatrix gsl_matrix_complex_const_row
 110         gsl_matrix_complex_const_column gsl_matrix_complex_const_diagonal gsl_matrix_complex_const_subdiagonal
 111         gsl_matrix_complex_const_superdiagonal gsl_matrix_complex_const_subrow gsl_matrix_complex_const_subcolumn
 112         gsl_matrix_complex_const_view_array gsl_matrix_complex_const_view_array_with_tda gsl_matrix_complex_const_view_vector
 113         gsl_matrix_complex_const_view_vector_with_tda gsl_matrix_complex_get gsl_matrix_complex_set
 114         gsl_matrix_complex_ptr gsl_matrix_complex_const_ptr gsl_matrix_complex_set_zero
 115         gsl_matrix_complex_set_identity gsl_matrix_complex_set_all gsl_matrix_complex_fread
 116         gsl_matrix_complex_fwrite gsl_matrix_complex_fscanf gsl_matrix_complex_fprintf
 117         gsl_matrix_complex_memcpy gsl_matrix_complex_swap gsl_matrix_complex_swap_rows
 118         gsl_matrix_complex_swap_columns gsl_matrix_complex_swap_rowcol gsl_matrix_complex_transpose
 119         gsl_matrix_complex_transpose_memcpy gsl_matrix_complex_isnull gsl_matrix_complex_ispos
 120         gsl_matrix_complex_isneg gsl_matrix_complex_add gsl_matrix_complex_sub
 121         gsl_matrix_complex_mul_elements gsl_matrix_complex_div_elements gsl_matrix_complex_scale
 122         gsl_matrix_complex_add_constant gsl_matrix_complex_add_diagonal gsl_matrix_complex_get_row
 123         gsl_matrix_complex_get_col gsl_matrix_complex_set_row gsl_matrix_complex_set_col /;
 124
 125
 126
 127 %EXPORT_TAGS = ( all => [ @EXPORT_OK ],
 128                  char => [ qw/
 129                     gsl_matrix_char_alloc
 130                     gsl_matrix_char_calloc
 131                     gsl_matrix_char_alloc_from_block
 132                     gsl_matrix_char_alloc_from_matrix
 133                     gsl_vector_char_alloc_row_from_matrix
 134                     gsl_vector_char_alloc_col_from_matrix
 135                     gsl_matrix_char_free
 136                     gsl_matrix_char_submatrix
 137                     gsl_matrix_char_row
 138                     gsl_matrix_char_column
 139                     gsl_matrix_char_diagonal
 140                     gsl_matrix_char_subdiagonal
 141                     gsl_matrix_char_superdiagonal
 142                     gsl_matrix_char_subrow
 143                     gsl_matrix_char_subcolumn
 144                     gsl_matrix_char_view_array
 145                     gsl_matrix_char_view_array_with_tda
 146                     gsl_matrix_char_view_vector
 147                     gsl_matrix_char_view_vector_with_tda
 148                     gsl_matrix_char_const_submatrix
 149                     gsl_matrix_char_const_row
 150                     gsl_matrix_char_const_column
 151                     gsl_matrix_char_const_diagonal
 152                     gsl_matrix_char_const_subdiagonal
 153                     gsl_matrix_char_const_superdiagonal
 154                     gsl_matrix_char_const_subrow
 155                     gsl_matrix_char_const_subcolumn
 156                     gsl_matrix_char_const_view_array
 157                     gsl_matrix_char_const_view_array_with_tda
 158                     gsl_matrix_char_const_view_vector
 159                     gsl_matrix_char_const_view_vector_with_tda
 160                     gsl_matrix_char_get
 161                     gsl_matrix_char_set
 162                     gsl_matrix_char_ptr
 163                     gsl_matrix_char_const_ptr
 164                     gsl_matrix_char_set_zero
 165                     gsl_matrix_char_set_identity
 166                     gsl_matrix_char_set_all
 167                     gsl_matrix_char_fread
 168                     gsl_matrix_char_fwrite
 169                     gsl_matrix_char_fscanf
 170                     gsl_matrix_char_fprintf
 171                     gsl_matrix_char_memcpy
 172                     gsl_matrix_char_swap
 173                     gsl_matrix_char_swap_rows
 174                     gsl_matrix_char_swap_columns
 175                     gsl_matrix_char_swap_rowcol
 176                     gsl_matrix_char_transpose
 177                     gsl_matrix_char_transpose_memcpy
 178                     gsl_matrix_char_max
 179                     gsl_matrix_char_min
 180                     gsl_matrix_char_minmax
 181                     gsl_matrix_char_max_index
 182                     gsl_matrix_char_min_index
 183                     gsl_matrix_char_minmax_index
 184                     gsl_matrix_char_isnull
 185                     gsl_matrix_char_ispos
 186                     gsl_matrix_char_isneg
 187                     gsl_matrix_char_isnonneg
 188                     gsl_matrix_char_add
 189                     gsl_matrix_char_sub
 190                     gsl_matrix_char_mul_elements
 191                     gsl_matrix_char_div_elements
 192                     gsl_matrix_char_scale
 193                     gsl_matrix_char_add_constant
 194                     gsl_matrix_char_add_diagonal
 195                      /],
 196
 197                  double => [ qw/
 198                     gsl_matrix_alloc
 199                     gsl_matrix_calloc
 200                     gsl_matrix_alloc_from_block
 201                     gsl_matrix_alloc_from_matrix
 202                     gsl_vector_alloc_row_from_matrix
 203                     gsl_vector_alloc_col_from_matrix
 204                     gsl_matrix_free
 205                     gsl_matrix_submatrix
 206                     gsl_matrix_row
 207                     gsl_matrix_column
 208                     gsl_matrix_diagonal
 209                     gsl_matrix_subdiagonal
 210                     gsl_matrix_superdiagonal
 211                     gsl_matrix_subrow
 212                     gsl_matrix_subcolumn
 213                     gsl_matrix_view_array
 214                     gsl_matrix_view_array_with_tda
 215                     gsl_matrix_view_vector
 216                     gsl_matrix_view_vector_with_tda
 217                     gsl_matrix_const_submatrix
 218                     gsl_matrix_const_row
 219                     gsl_matrix_const_column
 220                     gsl_matrix_const_diagonal
 221                     gsl_matrix_const_subdiagonal
 222                     gsl_matrix_const_superdiagonal
 223                     gsl_matrix_const_subrow
 224                     gsl_matrix_const_subcolumn
 225                     gsl_matrix_const_view_array
 226                     gsl_matrix_const_view_array_with_tda
 227                     gsl_matrix_const_view_vector
 228                     gsl_matrix_const_view_vector_with_tda
 229                     gsl_matrix_get
 230                     gsl_matrix_set
 231                     gsl_matrix_ptr
 232                     gsl_matrix_const_ptr
 233                     gsl_matrix_set_zero
 234                     gsl_matrix_set_identity
 235                     gsl_matrix_set_all
 236                     gsl_matrix_fread
 237                     gsl_matrix_fwrite
 238                     gsl_matrix_fscanf
 239                     gsl_matrix_fprintf
 240                     gsl_matrix_memcpy
 241                     gsl_matrix_swap
 242                     gsl_matrix_swap_rows
 243                     gsl_matrix_swap_columns
 244                     gsl_matrix_swap_rowcol
 245                     gsl_matrix_transpose
 246                     gsl_matrix_transpose_memcpy
 247                     gsl_matrix_max
 248                     gsl_matrix_minmax
 249                     gsl_matrix_max_index
 250                     gsl_matrix_min_index
 251                     gsl_matrix_minmax_index
 252                     gsl_matrix_isnull
 253                     gsl_matrix_ispos
 254                     gsl_matrix_isneg
 255                     gsl_matrix_isnonneg
 256                     gsl_matrix_add
 257                     gsl_matrix_mul_elements
 258                     gsl_matrix_div_elements
 259                     gsl_matrix_scale
 260                     gsl_matrix_add_constant
 261                     gsl_matrix_add_diagonal
 262                     /],
 263 int => [ qw/
 264                  gsl_matrix_int_alloc
 265                  gsl_matrix_int_alloc_from_matrix
 266                  gsl_matrix_int_free
 267                  gsl_matrix_int_column
 268                  gsl_matrix_int_superdiagonal
 269                  gsl_matrix_int_view_array_with_tda
 270                  gsl_matrix_int_const_submatrix
 271                  gsl_matrix_int_const_diagonal
 272                  gsl_matrix_int_const_subrow
 273                  gsl_matrix_int_const_view_array_with_tda
 274                  gsl_matrix_int_get
 275                  gsl_matrix_int_ptr
 276                  gsl_matrix_int_set_zero
 277                  gsl_matrix_int_fread
 278                  gsl_matrix_int_fscanf
 279                  gsl_matrix_int_memcpy
 280                  gsl_matrix_int_swap_rows
 281                  gsl_matrix_int_transpose
 282                  gsl_matrix_int_max
 283                  gsl_matrix_int_max_index
 284                  gsl_matrix_int_minmax_index
 285                  gsl_matrix_int_ispos
 286                  gsl_matrix_int_add
 287                  gsl_matrix_int_mul_elements
 288                  gsl_matrix_int_add_constant
 289                  /],
 290   complex => [ qw/
 291                  gsl_matrix_complex_alloc
 292                  gsl_matrix_complex_calloc
 293                  gsl_matrix_complex_alloc_from_block
 294                  gsl_matrix_complex_alloc_from_matrix
 295                  gsl_vector_complex_alloc_row_from_matrix
 296                  gsl_vector_complex_alloc_col_from_matrix
 297                  gsl_matrix_complex_free
 298                  gsl_matrix_complex_submatrix
 299                  gsl_matrix_complex_row
 300                  gsl_matrix_complex_column
 301                  gsl_matrix_complex_diagonal
 302                  gsl_matrix_complex_subdiagonal
 303                  gsl_matrix_complex_superdiagonal
 304                  gsl_matrix_complex_subrow
 305                  gsl_matrix_complex_subcolumn
 306                  gsl_matrix_complex_view_array
 307                  gsl_matrix_complex_view_array_with_tda
 308                  gsl_matrix_complex_view_vector
 309                  gsl_matrix_complex_view_vector_with_tda
 310                  gsl_matrix_complex_const_submatrix
 311                  gsl_matrix_complex_const_row
 312                  gsl_matrix_complex_const_column
 313                  gsl_matrix_complex_const_diagonal
 314                  gsl_matrix_complex_const_subdiagonal
 315                  gsl_matrix_complex_const_superdiagonal
 316                  gsl_matrix_complex_const_subrow
 317                  gsl_matrix_complex_const_subcolumn
 318                  gsl_matrix_complex_const_view_array
 319                  gsl_matrix_complex_const_view_array_with_tda
 320                  gsl_matrix_complex_const_view_vector
 321                  gsl_matrix_complex_const_view_vector_with_tda
 322                  gsl_matrix_complex_get
 323                  gsl_matrix_complex_set
 324                  gsl_matrix_complex_ptr
 325                  gsl_matrix_complex_const_ptr
 326                  gsl_matrix_complex_set_zero gsl_matrix_complex_set_identity
 327                  gsl_matrix_complex_set_all
 328                  gsl_matrix_complex_fread
 329                  gsl_matrix_complex_fwrite
 330                  gsl_matrix_complex_fscanf
 331                  gsl_matrix_complex_fprintf
 332                  gsl_matrix_complex_memcpygsl_matrix_complex_swap
 333                  gsl_matrix_complex_swap_rows
 334                  gsl_matrix_complex_swap_columns
 335                  gsl_matrix_complex_swap_rowcol
 336                  gsl_matrix_complex_transpose
 337                  gsl_matrix_complex_transpose_memcpy
 338                  gsl_matrix_complex_isnull
 339                  gsl_matrix_complex_ispos
 340                  gsl_matrix_complex_isneg
 341                  gsl_matrix_complex_add
 342                  gsl_matrix_complex_sub
 343                  gsl_matrix_complex_mul_elements
 344                  gsl_matrix_complex_div_elements
 345                  gsl_matrix_complex_scale
 346                  gsl_matrix_complex_add_constant
 347                  gsl_matrix_complex_add_diagonal
 348                  gsl_matrix_complex_get_row
 349                  gsl_matrix_complex_get_col
 350                  gsl_matrix_complex_set_row
 351                  gsl_matrix_complex_set_col
 352                  /]);
 353
 354 sub new {
 355     my ($class, $rows, $cols) = @_;
 356     my $this = {};
 357     my $matrix;
 358     if ( defined $rows       && defined $cols &&
 359         $rows > 0            && $cols > 0     &&
 360         (int $rows == $rows) && (int $cols == $cols)){
 361
 362         $matrix  = gsl_matrix_alloc($rows,$cols);
 363     } else {
 364         die __PACKAGE__.'::new($x,$y) - $x and $y must be positive integers';
 365     }
 366     $this->{_matrix} = $matrix;
 367     ($this->{_rows}, $this->{_cols}) = ($rows,$cols);
 368     bless $this, $class;
 369 }
 370 sub raw  { (shift)->{_matrix} }
 371 sub rows { (shift)->{_rows}   }
 372 sub cols { (shift)->{_cols}   }
 373
 374 __END__
 375
 376 =head1 NAME
 377
 378 Math::GSL::Matrix - Mathematical functions concerning Matrices
 379
 380 =head1 SYPNOPSIS
 381
 382 use Math::GSL::Matrix qw/:all/;
 383
 384 =head1 DESCRIPTION
 385
 386
 387 Here is a list of all the functions included in this module :
 388
 389  gsl_matrix_alloc($i, $j) - Return a gsl_matrix of $i rows and $j columns
 390  gsl_matrix_calloc($i, $j) - Return a gsl_matrix of $i rows and $j columns and initialize all of the elements of the matrix to zero
 391  gsl_matrix_alloc_from_block
 392  gsl_matrix_alloc_from_matrix
 393  gsl_vector_alloc_row_from_matrix
 394  gsl_vector_alloc_col_from_matrix
 395  gsl_matrix_free
 396  gsl_matrix_submatrix($m, $k1, $k2, $n1, $n2) - Return a matrix view of the matrix $m. The upper-left element of the submatrix is the element ($k1,$k2) of the original matrix. The submatrix has $n1 rows and $n2 columns.
 397  gsl_matrix_row($m , $i) - Return a vector view of the $i-th row of the matrix $m
 398  gsl_matrix_column($m, $j) - Return a vector view of the $j-th column of the matrix $m
 399  gsl_matrix_diagonal($m) - Return a vector view of the diagonal of the vector. The matrix doesn't have to be square.
 400  gsl_matrix_subdiagonal($m, $k) - Return a vector view of the $k-th subdiagonal of the matrix $m. The diagonal of the matrix corresponds to k=0.
 401  gsl_matrix_superdiagonal($m, $k) - Return a vector view of the $k-th superdiagonal of the matrix $m. The matrix doesn't have to be square.
 402  gsl_matrix_subrow($m, $i, $offset, $n) - Return a vector view of the $i-th row of the matrix $m beginning at offset elements and containing n elements.
 403  gsl_matrix_subcolumn($m, $j, $offset, $n) - Return a vector view of the $j-th column of the matrix $m beginning at offset elements and containing n elements.
 404  gsl_matrix_view_array
 405  gsl_matrix_view_array_with_tda
 406  gsl_matrix_view_vector
 407  gsl_matrix_view_vector_with_tda
 408  gsl_matrix_const_submatrix
 409  gsl_matrix_get($m, $i, $j) - Return the (i,j)-th element of the matrix $m
 410  gsl_matrix_set($m, $i, $j, $x) - Set the value of the (i,j)-th element of the matrix $m to $x
 411  gsl_matrix_ptr
 412  gsl_matrix_const_ptr
 413  gsl_matrix_set_zero($m) - Set all the elements of the matrix $m to zero
 414  gsl_matrix_set_identity($m) - Set the elements of the matrix $m to the corresponding elements of the identity matrix
 415  gsl_matrix_set_all($m, $x) - Set all the elements of the matrix $m to the value $x
 416  gsl_matrix_fread($fh, $m) - Read a file which has been written with gsl_matrix_fwrite from the stream $fh opened with the fopen function and stores the data inside the matrix $m
 417  gsl_matrix_fwrite($fh, $m) - Write the elements of the matrix $m in binary format to a stream $fh opened with the fopen function
 418  gsl_matrix_fscanf($fh, $m) - Read a file which has been written with gsl_matrix_fprintf from the stream $fh opened with the fopen function and stores the data inside the matrix $m
 419  gsl_matrix_fprintf($fh, $m, $format) - Write the elements of the matrix $m in the format $format (for example "%f" is the format for double) to a stream $fh opened with the fopen function
 420  gsl_matrix_memcpy($dest, $src) - Copy the elements of the matrix $src to the matrix $dest. The two matrices must have the same size.
 421  gsl_matrix_swap($m1, $m2) - Exchange the elements of the matrices $m1 and $m2 by copying. The two matrices must have the same size.
 422  gsl_matrix_swap_rows($m, $i, $j) - Exchange the $i-th and $j-th row of the matrix $m. The function returns 0 if the operation suceeded, 1 otherwise.
 423  gsl_matrix_swap_columns($m, $i, $j) - Exchange the $i-th and $j-th column of the matrix $m. The function returns 0 if the operation suceeded, 1 otherwise.
 424  gsl_matrix_swap_rowcol($m, $i, $j) - Exchange the $i-th row and the $j-th column of the matrix $m. The matrix must be square. The function returns 0 if the operation suceeded, 1 otherwise.
 425  gsl_matrix_transpose($m) - This function replaces the matrix m by its transpose by copying the elements of the matrix in-place. The matrix must be square for this operation to be possible.
 426  gsl_matrix_transpose_memcpy($dest, $src) - Make the matrix $dest the transpose of the matrix $src. This function works for all matrices provided that the dimensions of the matrix dest match the transposed dimensions of the matrix src.
 427  gsl_matrix_max($m) - Return the maximum value in the matrix $m
 428  gsl_matrix_min($m) - Return the minimum value in the matrix $m
 429  gsl_matrix_minmax($m) - Return two values, the first is the minimum value of the Matrix $m and the second is the maximum of the same the same matrix.
 430  gsl_matrix_max_index($m) - Return two values, the first is the the i indice of the maximum value of the matrix $m and the second is the j indice of the same value.
 431  gsl_matrix_min_index - Return two values, the first is the the i indice of the minimum value of the matrix $m and the second is the j indice of the same value.
 432  gsl_matrix_minmax_index($m) - Return four values, the first is the i indice of the minimum of the matrix $m, the second is the j indice of the same value, the third is the i indice of the maximum of the matrix $m and the fourth is the j indice of the same value
 433  gsl_matrix_isnull($m) - Return 1 if all the elements of the matrix $m are zero, 0 otherwise
 434  gsl_matrix_ispos($m) - Return 1 if all the elements of the matrix $m are strictly positve, 0 otherwise
 435  gsl_matrix_isneg($m) - Return 1 if all the elements of the matrix $m are strictly negative, 0 otherwise
 436  gsl_matrix_isnonneg($m) - Return 1 if all the elements of the matrix $m are non-negatuive, 0 otherwise
 437  gsl_matrix_add($a, $b) - Add the elements of matrix $b to the elements of matrix $a
 438  gsl_matrix_sub($a, $b) - Subtract the elements of matrix $b from the elements of matrix $a
 439  gsl_matrix_mul_elements($a, $b) - Multiplie the elements of matrix $a by the elements of matrix $b
 440  gsl_matrix_div_elements($a, $b) - Divide the elements of matrix $a by the elements of matrix $b
 441  gsl_matrix_scale($a, $x) - Multiplie the elements of matrix $a by the constant factor $x
 442  gsl_matrix_add_constant($a, $x) - Add the constant value $x to the elements of the matrix $a
 443  gsl_matrix_add_diagonal($a, $x) - Add the constant value $x to the elements of the diagonal of the matrix $a
 444  gsl_matrix_get_row($v, $m, $i) - Copy the elements of the $i-th row of the matrix $m into the vector $v. The lenght of the vector must be of the same as the lenght of the row. The function returns 0 if it succeded, 1 otherwise.
 445  gsl_matrix_get_col($v, $m, $i) - Copy the elements of the $j-th column of the matrix $m into the vector $v. The lenght of the vector must be of the same as the lenght of the column. The function returns 0 if it succeded, 1 otherwise.
 446  gsl_matrix_set_row($m, $i, $v) - Copy the elements of vector $v into the $i-th row of the matrix $m The lenght of the vector must be of the same as the lenght of the row. The function returns 0 if it succeded, 1 otherwise.
 447  gsl_matrix_set_col($m, $j, $v) - Copy the elements of vector $v into the $j-th row of the matrix $m The lenght of the vector must be of the same as the lenght of the column. The function returns 0 if it succeded, 1 otherwise.
 448
 449  The following functions are specific to matrices containing complex numbers :
 450
 451  gsl_matrix_complex_alloc
 452  gsl_matrix_complex_calloc
 453  gsl_matrix_complex_alloc_from_block
 454  gsl_matrix_complex_alloc_from_matrix
 455  gsl_vector_complex_alloc_row_from_matrix
 456  gsl_vector_complex_alloc_col_from_matrix
 457  gsl_matrix_complex_free
 458  gsl_matrix_complex_submatrix
 459  gsl_matrix_complex_row
 460  gsl_matrix_complex_column
 461  gsl_matrix_complex_diagonal
 462  gsl_matrix_complex_subdiagonal
 463  gsl_matrix_complex_superdiagonal
 464  gsl_matrix_complex_subrow
 465  gsl_matrix_complex_subcolumn
 466  gsl_matrix_complex_view_array
 467  gsl_matrix_complex_view_array_with_tda
 468  gsl_matrix_complex_view_vector
 469  gsl_matrix_complex_view_vector_with_tda
 470  gsl_matrix_complex_const_submatrix
 471  gsl_matrix_complex_const_row
 472  gsl_matrix_complex_const_column
 473  gsl_matrix_complex_const_diagonal
 474  gsl_matrix_complex_const_subdiagonal
 475  gsl_matrix_complex_const_superdiagonal
 476  gsl_matrix_complex_const_subrow
 477  gsl_matrix_complex_const_subcolumn
 478  gsl_matrix_complex_const_view_array
 479  gsl_matrix_complex_const_view_array_with_tda
 480  gsl_matrix_complex_const_view_vector
 481  gsl_matrix_complex_const_view_vector_with_tda
 482  gsl_matrix_complex_get
 483  gsl_matrix_complex_set
 484  gsl_matrix_complex_ptr
 485  gsl_matrix_complex_const_ptr
 486  gsl_matrix_complex_set_zero
 487  gsl_matrix_complex_set_identity
 488  gsl_matrix_complex_set_all
 489  gsl_matrix_complex_fread
 490  gsl_matrix_complex_fwrite
 491  gsl_matrix_complex_fscanf
 492  gsl_matrix_complex_fprintf
 493  gsl_matrix_complex_memcpy
 494  gsl_matrix_complex_swap
 495  gsl_matrix_complex_swap_rows
 496  gsl_matrix_complex_swap_columns
 497  gsl_matrix_complex_swap_rowcol
 498  gsl_matrix_complex_transpose
 499  gsl_matrix_complex_transpose_memcpy
 500  gsl_matrix_complex_isnull
 501  gsl_matrix_complex_ispos
 502  gsl_matrix_complex_isneg
 503  gsl_matrix_complex_add
 504  gsl_matrix_complex_sub
 505  gsl_matrix_complex_mul_elements
 506  gsl_matrix_complex_div_elements
 507  gsl_matrix_complex_scale
 508  gsl_matrix_complex_add_constant
 509  gsl_matrix_complex_add_diagonal
 510  gsl_matrix_complex_get_row
 511  gsl_matrix_complex_get_col
 512  gsl_matrix_complex_set_row
 513  gsl_matrix_complex_set_col
 514
 515
 516  The following functions are the same as the previous enonced ones except that they work with other data types than double.
 517
 518  gsl_matrix_const_row gsl_matrix_const_column gsl_matrix_const_diagonal
 519  gsl_matrix_const_subdiagonal gsl_matrix_const_superdiagonal
 520  gsl_matrix_const_subrow gsl_matrix_const_subcolumn
 521  gsl_matrix_const_view_array gsl_matrix_const_view_array_with_tda
 522  gsl_matrix_const_view_vector gsl_matrix_const_view_vector_with_tda
 523  gsl_matrix_char_alloc gsl_matrix_char_calloc  gsl_matrix_char_alloc_from_block
 524  gsl_matrix_char_alloc_from_matrix gsl_vector_char_alloc_row_from_matrix gsl_vector_char_alloc_col_from_matrix
 525  gsl_matrix_char_free gsl_matrix_char_submatrix
 526  gsl_matrix_char_row  gsl_matrix_char_column
 527  gsl_matrix_char_diagonal gsl_matrix_char_subdiagonal gsl_matrix_char_superdiagonal
 528  gsl_matrix_char_subrow gsl_matrix_char_subcolumn gsl_matrix_char_view_array
 529  gsl_matrix_char_view_array_with_tda gsl_matrix_char_view_vector gsl_matrix_char_view_vector_with_tda
 530  gsl_matrix_char_const_submatrix gsl_matrix_char_const_row gsl_matrix_char_const_column
 531  gsl_matrix_char_const_diagonal gsl_matrix_char_const_subdiagonal gsl_matrix_char_const_superdiagonal
 532  gsl_matrix_char_const_subrow gsl_matrix_char_const_subcolumn gsl_matrix_char_const_view_array
 533  gsl_matrix_char_const_view_array_with_tda gsl_matrix_char_const_view_vector gsl_matrix_char_const_view_vector_with_tda
 534  gsl_matrix_char_get gsl_matrix_char_set gsl_matrix_char_ptr gsl_matrix_char_const_ptr
 535  gsl_matrix_char_set_zero gsl_matrix_char_set_identity
 536  gsl_matrix_char_set_all  gsl_matrix_char_fread
 537  gsl_matrix_char_fwrite gsl_matrix_char_fscanf gsl_matrix_char_fprintf
 538  gsl_matrix_char_memcpy gsl_matrix_char_swap
 539  gsl_matrix_char_swap_rows gsl_matrix_char_swap_columns
 540  gsl_matrix_char_swap_rowcol gsl_matrix_char_transpose gsl_matrix_char_transpose_memcpy
 541  gsl_matrix_char_max gsl_matrix_char_min
 542  gsl_matrix_char_minmax  gsl_matrix_char_max_index
 543  gsl_matrix_char_min_index gsl_matrix_char_minmax_index
 544  gsl_matrix_char_isnull gsl_matrix_char_ispos gsl_matrix_char_isneg
 545  gsl_matrix_char_isnonneg  gsl_matrix_char_add
 546  gsl_matrix_char_sub gsl_matrix_char_mul_elements gsl_matrix_char_div_elements
 547  gsl_matrix_char_scale gsl_matrix_char_add_constant gsl_matrix_char_add_diagonal
 548  gsl_matrix_int_alloc gsl_matrix_int_calloc gsl_matrix_int_alloc_from_block
 549  gsl_matrix_int_alloc_from_matrix gsl_vector_int_alloc_row_from_matrix gsl_vector_int_alloc_col_from_matrix
 550  gsl_matrix_int_free gsl_matrix_int_submatrix gsl_matrix_int_row
 551  gsl_matrix_int_column gsl_matrix_int_diagonal gsl_matrix_int_subdiagonal
 552  gsl_matrix_int_superdiagonal gsl_matrix_int_subrow gsl_matrix_int_subcolumn gsl_matrix_int_view_array
 553  gsl_matrix_int_view_array_with_tda gsl_matrix_int_view_vector gsl_matrix_int_view_vector_with_tda
 554  gsl_matrix_int_const_submatrix gsl_matrix_int_const_row gsl_matrix_int_const_column
 555  gsl_matrix_int_const_diagonal gsl_matrix_int_const_subdiagonal gsl_matrix_int_const_superdiagonal
 556  gsl_matrix_int_const_subrow gsl_matrix_int_const_subcolumn gsl_matrix_int_const_view_array
 557  gsl_matrix_int_const_view_array_with_tda gsl_matrix_int_const_view_vector gsl_matrix_int_const_view_vector_with_tda
 558  gsl_matrix_int_get gsl_matrix_int_set
 559  gsl_matrix_int_ptr gsl_matrix_int_const_ptr
 560  gsl_matrix_int_set_zero gsl_matrix_int_set_identity gsl_matrix_int_set_all
 561  gsl_matrix_int_fread gsl_matrix_int_fwrite
 562  gsl_matrix_int_fscanf gsl_matrix_int_fprintf
 563  gsl_matrix_int_memcpy gsl_matrix_int_swap
 564  gsl_matrix_int_swap_rows gsl_matrix_int_swap_columns gsl_matrix_int_swap_rowcol
 565  gsl_matrix_int_transpose gsl_matrix_int_transpose_memcpy
 566  gsl_matrix_int_max gsl_matrix_int_min gsl_matrix_int_minmax
 567  gsl_matrix_int_max_index gsl_matrix_int_min_index
 568  gsl_matrix_int_minmax_index  gsl_matrix_int_isnull
 569  gsl_matrix_int_ispos gsl_matrix_int_isneg gsl_matrix_int_isnonneg
 570  gsl_matrix_int_add gsl_matrix_int_sub
 571  gsl_matrix_int_mul_elements gsl_matrix_int_div_elements gsl_matrix_int_scale
 572  gsl_matrix_int_add_constant gsl_matrix_int_add_diagonal
 573
 574  You have to add the functions you want to use inside the qw /put_funtion_here /.
 575  You can also write use Math::GSL::PowInt qw/:name_of_tag/ to use all avaible functions of the module.
 576  Other tags are also avaible, here is a complete list of all tags for this module :
 577
 578  all
 579  int
 580  double
 581  char
 582  complex
 583
 584  For more informations on the functions, we refer you to the GSL offcial documentation: http://www.gnu.org/software/gsl/manual/html_node/
 585  Tip : search on google: site:http://www.gnu.org/software/gsl/manual/html_node/name_of_the_function_you_want
 586
 587 =head1 EXAMPLES
 588
 589  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
 590
 591  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.
 592
 593  use Math::GSL::Matrix qw/:all/;
 594  $m = gsl_matrix_alloc (10,3);
 595  for($i = 0; $i< 10; $i++){
 596     for($j = 0; $j<3; $j++){
 597         gsl_matrix_set($m, $i, $j, 0.23 + 100*$i + $j);
 598     }
 599  }
 600
 601  for($i = 0; $i< 100; $i++){ # OUT OF RANGE ERROR
 602      for($j=0; $j<3; $j++){
 603         print "m($i, $j) = " . gsl_matrix_get ($m, $i, $j) . "\n";
 604     }
 605  }
 606  gsl_matrix_free ($m);
 607
 608
 609  use Math::GSL::Matrix qw/:all/;
 610
 611  $m = gsl_matrix_alloc (100, 100);
 612  $a = gsl_matrix_alloc (100, 100);
 613
 614  for($i = 0; $i < 100; $i++){
 615      for($j = 0; $j < 100; $j++){
 616          gsl_matrix_set ($m, $i, $j, 0.23 + $i + $j);
 617      }
 618  }
 619
 620  The next program shows how to write a matrix to a file.
 621
 622  $out = fopen("test.dat", "wb");
 623  gsl_matrix_fwrite ($out, $m);
 624  fclose ($out);
 625
 626  $in = fopen("test.dat", "rb");
 627  gsl_matrix_fread ($in, $a);
 628  fclose($in);
 629
 630  $k=0;
 631  for($i = 0; $i < 100; $i++){
 632      for($j = 0; $j < 100; $j++){
 633          $mij = gsl_matrix_get ($m, $i, $j);
 634          $aij = gsl_matrix_get ($a, $i, $j);
 635          if ($mij != $aij){ $k++ };
 636      }
 637  }
 638
 639  gsl_matrix_free($m);
 640  gsl_matrix_free($a);
 641
 642  print "differences = $k (should be zero)\n";
 643
 644 =head1 AUTHOR
 645
 646 Jonathan Leto <jonathan@leto.net> and Thierry Moisan <thierry.moisan@gmail.com>
 647
 648 =head1 COPYRIGHT AND LICENSE
 649
 650 Copyright (C) 2008 Jonathan Leto and Thierry Moisan
 651
 652 This program is free software; you can redistribute it and/or modify it
 653 under the same terms as Perl itself.
 654
 655 =cut
 656 %}