Vector.i

   1 %module Vector
   2
   3 %include "typemaps.i"
   4
   5 %apply int *INOUT { size_t *imin, size_t *imax };
   6
   7 %apply double *INOUT { double * min_out, double * max_out };
   8
   9 %typemap(in) double *v {
  10     AV *tempav;
  11     I32 len;
  12     int i;
  13     SV **tv;
  14     if (!SvROK($input))
  15         croak("Math::GSL : $input is not a reference!");
  16     if (SvTYPE(SvRV($input)) != SVt_PVAV)
  17         croak("Math::GSL : $input is not an array ref!");
  18
  19     tempav = (AV*)SvRV($input);
  20     len = av_len(tempav);
  21     $1 = (double *) malloc((len+1)*sizeof(double));
  22     for (i = 0; i <= len; i++) {
  23         tv = av_fetch(tempav, i, 0);
  24         $1[i] = (double) SvNV(*tv);
  25     }
  26     $1[i] = GSL_NAN;
  27 }
  28
  29
  30 FILE * fopen(char *, char *);
  31 int fclose(FILE *);
  32
  33 %{
  34     #include "/usr/include/stdio.h"
  35     #include "/usr/local/include/gsl/gsl_nan.h"
  36     #include "/usr/local/include/gsl/gsl_vector.h"
  37     #include "/usr/local/include/gsl/gsl_vector_char.h"
  38     #include "/usr/local/include/gsl/gsl_vector_complex.h"
  39     #include "/usr/local/include/gsl/gsl_vector_complex_double.h"
  40     #include "/usr/local/include/gsl/gsl_vector_double.h"
  41     #include "/usr/local/include/gsl/gsl_vector_float.h"
  42     #include "/usr/local/include/gsl/gsl_vector_int.h"
  43 %}
  44
  45 %include "/usr/local/include/gsl/gsl_nan.h"
  46 %include "/usr/local/include/gsl/gsl_vector.h"
  47 %include "/usr/local/include/gsl/gsl_vector_char.h"
  48 %include "/usr/local/include/gsl/gsl_vector_complex.h"
  49 %include "/usr/local/include/gsl/gsl_vector_complex_double.h"
  50 %include "/usr/local/include/gsl/gsl_vector_double.h"
  51 %include "/usr/local/include/gsl/gsl_vector_int.h"
  52
  53
  54
  55 %perlcode %{
  56
  57
  58 @EXPORT_OK  = qw/fopen fclose
  59                  gsl_vector_alloc gsl_vector_calloc gsl_vector_alloc_from_b gsl_vector_alloc_from_v
  60                  gsl_vector_free gsl_vector_view_array gsl_vector_view_array_w
  61                  gsl_vector_const_view_a gsl_vector_subvector gsl_vector_subvector_wi gsl_vector_subvector_with_stride
  62                  gsl_vector_const_subvec gsl_vector_const_subvec gsl_vector_get gsl_vector_set
  63                  gsl_vector_ptr gsl_vector_const_ptr gsl_vector_set_zero gsl_vector_set_all
  64                  gsl_vector_set_basis gsl_vector_fread gsl_vector_fwrite gsl_vector_fscanf
  65                  gsl_vector_fprintf gsl_vector_memcpy gsl_vector_reverse gsl_vector_swap
  66                  gsl_vector_swap_elements gsl_vector_max gsl_vector_min gsl_vector_minmax
  67                  gsl_vector_max_index gsl_vector_min_index gsl_vector_minmax_index
  68                  gsl_vector_add gsl_vector_sub gsl_vector_mul gsl_vector_div
  69                  gsl_vector_scale gsl_vector_add_constant gsl_vector_isnull
  70                  gsl_vector_ispos gsl_vector_isneg gsl_vector_isnonneg
  71                  /;
  72 %EXPORT_TAGS = ( all => [ @EXPORT_OK ] );
  73
  74 sub new {
  75     my ($class, $values) = @_;
  76     my $length  = $#$values;
  77     my $vector;
  78     die __PACKAGE__.'::new($x) - $x must be an int or nonempty array reference'
  79         if( !(defined $values) || ($length == -1));
  80
  81     if ( ref $values eq 'ARRAY' ){
  82         $vector  = gsl_vector_alloc($length+1);
  83         map { gsl_vector_set($vector, $_, $values->[$_] ) }  (0 .. $length);
  84     } elsif ( int $values == $values && $values > 0) {
  85         $vector  = gsl_vector_alloc($length);
  86     } else {
  87         die __PACKAGE__.'::new($x) - $x must be an int or array reference';
  88     }
  89     my $self = {};
  90     $self->{_vector} = $vector;
  91     bless $self, $class;
  92 }
  93
  94 sub get {
  95     my ($self, $indices) = @_;
  96     return  map {  gsl_vector_get($self->{_vector}, $_ ) } @$indices ;
  97 }
  98
  99 sub set {
 100     my ($self, $indices, $values) = @_;
 101     die (__PACKAGE__.'::set($indices, $values) - $indices and $values must be array references of the same length')
 102         unless ( ref $indices eq 'ARRAY' && ref $values eq 'ARRAY' &&  $#$indices == $#$values );
 103     eval {
 104         map {  gsl_vector_set($self->{_vector}, $indices->[$_], $values->[$_] ) } (0..$#$indices);
 105     };
 106     return;
 107 }
 108
 109 __END__
 110
 111 =head1 NAME
 112
 113 Math::GSL::Vector - Functions concerning vectors
 114
 115 =head1 SYPNOPSIS
 116
 117     use Math::GSL::Vector qw/:all/;
 118
 119 =head1 DESCRIPTION
 120
 121 Here is a list of all the functions included in this module :
 122
 123     gsl_vector_alloc($x) - create a vector of size $x
 124
 125     gsl_vector_calloc($x) - create a vector of size $x and initializes all the elements of the vector to zero
 126
 127     gsl_vector_alloc_from_b
 128
 129     gsl_vector_alloc_from_v
 130
 131     gsl_vector_free($v) - free a previously allocated vector $v
 132
 133     gsl_vector_view_array
 134
 135     gsl_vector_view_array_w
 136
 137     gsl_vector_const_view_a
 138
 139     gsl_vector_subvector - return a vector_view type which contains a subvector of $v, with a size of $size, starting from the $offset position
 140
 141     gsl_vector_subvector_wi
 142
 143     gsl_vector_subvector_with_stride($v, $offset, $stride, $size) - return a vector_view type which contains a subvector of $v, with a size of $size, starting from the $offset position and with a $stride step between each element of $v
 144
 145     gsl_vector_const_subvec
 146
 147     gsl_vector_get($v, $i) - return the $i-th element of a vector $v
 148
 149     gsl_vector_set($v, $i, $x) - return the vector $v with his $i-th element set to $x
 150
 151     gsl_vector_ptr
 152
 153     gsl_vector_const_ptr
 154
 155     gsl_vector_set_zero($v) - set all the elements of $v to 0
 156
 157     gsl_vector_set_all($v, $x) - set all the elements of $v to $x
 158
 159     gsl_vector_set_basis($v, $i) - set all the elements of $v to 0 except for the $i-th element which is set to 1 and return 0 if the operation succeded, 1 otherwise.
 160
 161     gsl_vector_fread($file, $v)
 162
 163     gsl_vector_fwrite
 164
 165     gsl_vector_fscanf
 166
 167     gsl_vector_fprintf
 168
 169     gsl_vector_memcpy
 170
 171     gsl_vector_reverse($v) - reverse the order of the elements of the vector $v and return 0 if the opertaion succeded, 1 otherwise
 172
 173     gsl_vector_swap($v, $v2) - swap the values of the vectors $v and $v2 and return 0 if the opertaion succeded, 1 otherwise
 174
 175     gsl_vector_swap_elements($v, $i, $j) - permute the elements at position $i and $j in the vector $v and return 0 if the operation succeded, 1 otherwise.
 176
 177     gsl_vector_max($v) - return the maximum value in the vector $v
 178
 179     gsl_vector_min($v) - return the minimum value in the vector $v
 180
 181     gsl_vector_minmax
 182
 183     gsl_vector_max_index($v) - return the position of the maximum value in the vector $v
 184
 185     gsl_vector_min_index - return the position of the minimum value in the vector $v
 186
 187     gsl_vector_minmax_index
 188
 189     gsl_vector_add($v, $v2) - add the elements of $v2 to the elements of $v, the two vectors must have the same lenght and return 0 if the operation succeded, 1 otherwise.
 190
 191     gsl_vector_sub($v, $v2) - substract the elements of $v2 from the elements of $v, the two vectors must have the same lenght and return 0 if the operation succeded, 1 otherwise.
 192
 193     gsl_vector_mul($v, $v2) - multiply the elements of $v by the elements of $v2, the two vectors must have the same lenght and return 0 if the operation succeded, 1 otherwise.
 194
 195     gsl_vector_div($v, $v2) - divides the elements of $v by the elements of $v2, the two vectors must have the same lenght and return 0 if the operation succeded, 1 otherwise.
 196
 197     gsl_vector_scale($v, $x) - multiplty the elements of the vector $v by a constant $x and return 0 if the operation succeded, 1 otherwise.
 198
 199     gsl_vector_add_constant($v, $x) - add a constant $x to the elements of the vector $v and return 0 if the operation succeded, 1 otherwise.
 200
 201     gsl_vector_isnull($v) - verify if all the elements of the vector $v are null, return 0 if it's the case, 1 otherwise.
 202
 203     gsl_vector_ispos($v) - verify if all the elements of the vector $v are positive, return 0 if it's the case, 1 otherwise.
 204
 205     gsl_vector_isneg($v) - verify if all the elements of the vector $v are negative, return 0 if it's the case, 1 otherwise.
 206
 207     gsl_vector_isnonneg($v) - verify if all the elements the vector $v are not negative, return 0 if it's the case, 1 otherwise.
 208
 209
 210 You have to add the functions you want to use inside the qw /put_funtion_here / with spaces between each function. You can also write use Math::GSL::Complex qw/:all/ to use all avaible functions of the module.
 211
 212
 213 Precision on the vector_view type : every modification you'll make on a vector_view will also modify the original vector.
 214 For example, the following code will zero the even elements of the vector v of length n, while leaving the odd elements untouched :
 215
 216 $v_even= gsl_vector_subvector_with_stride ($v, 0, 2, $size/2);
 217 gsl_vector_set_zero ($v_even->{vector});
 218
 219
 220 For more informations on the functions, we refer you to the GSL offcial documentation: http://www.gnu.org/software/gsl/manual/html_node/
 221 Tip : search on google: site:http://www.gnu.org/software/gsl/manual/html_node/ name_of_the_function_you_want
 222
 223 =head1 EXAMPLES
 224
 225
 226 =head1 AUTHOR
 227
 228 Jonathan Leto <jonathan@leto.net> and Thierry Moisan <thierry.moisan@gmail.com>
 229
 230 =head1 COPYRIGHT AND LICENSE
 231
 232 Copyright (C) 2008 Jonathan Leto and Thierry Moisan
 233
 234 This program is free software; you can redistribute it and/or modify it
 235 under the same terms as Perl itself.
 236
 237 =cut
 238 %}