pod/Histogram2D.pod

   1 %perlcode %{
   2 @EXPORT_OK = qw/
   3                gsl_histogram2d_alloc
   4                gsl_histogram2d_calloc
   5                gsl_histogram2d_calloc_uniform
   6                gsl_histogram2d_free
   7                gsl_histogram2d_increment
   8                gsl_histogram2d_accumulate
   9                gsl_histogram2d_find
  10                gsl_histogram2d_get
  11                gsl_histogram2d_get_xrange
  12                gsl_histogram2d_get_yrange
  13                gsl_histogram2d_xmax
  14                gsl_histogram2d_xmin
  15                gsl_histogram2d_nx
  16                gsl_histogram2d_ymax
  17                gsl_histogram2d_ymin
  18                gsl_histogram2d_ny
  19                gsl_histogram2d_reset
  20                gsl_histogram2d_calloc_range
  21                gsl_histogram2d_set_ranges_uniform
  22                gsl_histogram2d_set_ranges
  23                gsl_histogram2d_memcpy
  24                gsl_histogram2d_clone
  25                gsl_histogram2d_max_val
  26                gsl_histogram2d_max_bin
  27                gsl_histogram2d_min_val
  28                gsl_histogram2d_min_bin
  29                gsl_histogram2d_xmean
  30                gsl_histogram2d_ymean
  31                gsl_histogram2d_xsigma
  32                gsl_histogram2d_ysigma
  33                gsl_histogram2d_cov
  34                gsl_histogram2d_sum
  35                gsl_histogram2d_equal_bins_p
  36                gsl_histogram2d_add
  37                gsl_histogram2d_sub
  38                gsl_histogram2d_mul
  39                gsl_histogram2d_div
  40                gsl_histogram2d_scale
  41                gsl_histogram2d_shift
  42                gsl_histogram2d_fwrite
  43                gsl_histogram2d_fread
  44                gsl_histogram2d_fprintf
  45                gsl_histogram2d_fscanf
  46                gsl_histogram2d_pdf_alloc
  47                gsl_histogram2d_pdf_init
  48                gsl_histogram2d_pdf_free
  49                gsl_histogram2d_pdf_sample
  50              /;
  51 %EXPORT_TAGS = ( all => [ @EXPORT_OK ] );
  52
  53 =head1 NAME
  54
  55 Math::GSL::Histogram2D - Create and manipulate histograms of data in 2 dimensions
  56
  57 =head1 SYNOPSIS
  58
  59     use Math::GSL::Histogram2D qw/:all/;
  60
  61     my $H = gsl_histogram2d_alloc(100, 100);
  62     gsl_histogram2d_set_ranges_uniform($H,0,101, 0, 101);
  63     gsl_histogram2d_increment($H, -50, -12);  # ignored
  64     gsl_histogram2d_increment($H, 70 );
  65     gsl_histogram2d_increment($H, 85.2 );
  66
  67     my $G = gsl_histogram2d_clone($H);
  68     my $value = gsl_histogram2d_get($G, 70, 33);
  69     my ($max,$min) = (gsl_histogram2d_min_val($H), gsl_histogram2d_max_val($H) );
  70     my $sum = gsl_histogram2d_sum($H);
  71
  72 =cut
  73
  74 =head1 DESCRIPTION
  75
  76 Here is a list of all the functions included in this module :
  77
  78 =over
  79
  80 =item C<gsl_histogram2d_alloc($nx, $ny)> - This function allocates memory for a two-dimensional histogram with $nx bins in the x direction and $ny bins in the y direction. The function returns a pointer to a newly created gsl_histogram2d struct. If insufficient memory is available a null pointer is returned and the error handler is invoked with an error code of $GSL_ENOMEM. The bins and ranges must be initialized with one of the functions below before the histogram is ready for use.
  81
  82 =item C<gsl_histogram2d_calloc >
  83
  84 =item C<gsl_histogram2d_calloc_uniform >
  85
  86 =item C<gsl_histogram2d_free($h)> - This function frees the 2D histogram h and all of the memory associated with it.
  87
  88 =item C<gsl_histogram2d_increment($h, $xmin, $xmax, $ymin, $ymax)> - This function sets the ranges of the existing histogram $h to cover the ranges $xmin to $xmax and $ymin to $ymax uniformly. The values of the histogram bins are reset to zero.
  89
  90 =item C<gsl_histogram2d_accumulate($h, $x, $y, $weight)> - This function is similar to gsl_histogram2d_increment but increases the value of the appropriate bin in the histogram $h by the floating-point number $weight.
  91
  92 =item C<gsl_histogram2d_find($h, $x, $y)> - This function finds indices i and j to the to the bin which covers the coordinates ($x,$y). The bin is located using a binary search. The search includes an optimization for histograms with uniform ranges, and will return the correct bin immediately in this case. If ($x,$y) is found then the function return GSL_SUCCESS, i and j in this order. If ($x,$y) lies outside the valid range of the histogram then the function returns $GSL_EDOM and the error handler is invoked.
  93
  94 =item C<gsl_histogram2d_get($h, $i, $j)> - This function returns the contents of the ($i,$j)-th bin of the histogram $h. If ($i,$j) lies outside the valid range of indices for the histogram then the error handler is called with an error code of $GSL_EDOM and the function returns 0.
  95
  96 =item C<gsl_histogram2d_get_xrange($h, $i)> - This functions finds the upper and lower range limits of the $i-th in the x direction of the histogram $h. The range limits are stored in returned after 0 or 1 in this order : xlower and xupper. The lower limits are inclusive (i.e. events with these coordinates are included in the bin) and the upper limits are exclusive (i.e. events with the value of the upper limit are not included and fall in the neighboring higher bin, if it exists). The functions returns 0 has first value to indicate success. If $i lies outside the valid range of indices for the histogram then the error handler is called with an error code of GSL_EDOM.
  97
  98 =item C<gsl_histogram2d_get_yrange($h, $j)> - This functions finds the upper and lower range limits of the $j-th in the y direction of the histogram $h. The range limits are stored in returned after 0 or 1 in this order : ylower and yupper. The lower limits are inclusive (i.e. events with these coordinates are included in the bin) and the upper limits are exclusive (i.e. events with the value of the upper limit are not included and fall in the neighboring higher bin, if it exists). The functions returns 0 has first value to indicate success. If $j lies outside the valid range of indices for the histogram then the error handler is called with an error code of GSL_EDOM.
  99
 100 =item C<gsl_histogram2d_xmax($h)> - This functions returns the maximum upper range limit for the x direction of the histogram $h.
 101
 102 =item C<gsl_histogram2d_xmin($h)> - This functions returns the minimum lower range limit for the x direction of the histogram $h.
 103
 104 =item C<gsl_histogram2d_nx($h)> - This functions the number of bins for the x direction.
 105
 106 =item C<gsl_histogram2d_ymax($h)> - This functions returns the maximum upper range limit for the y direction of the histogram $h.
 107
 108 =item C<gsl_histogram2d_ymin($h)> - This functions returns the minimum lower range limit for the y direction of the histogram $h.
 109
 110 =item C<gsl_histogram2d_ny($h)> - This functions the number of bins for the y direction.
 111
 112 =item C<gsl_histogram2d_reset($h)> - This function resets all the bins of the histogram $h to zero.
 113
 114 =item C<gsl_histogram2d_calloc_range >
 115
 116 =item C<gsl_histogram2d_set_ranges($h, $xrange, $xsize, $yrange, $ysize)> - This function sets the ranges of the existing histogram $h using the arrays $xrange and $yrange of size $xsize and $ysize respectively. The values of the histogram bins are reset to zero.
 117
 118 =item C<gsl_histogram2d_set_ranges_uniform($h, $xmin, $xmax, $ymin, $ymax)> - This function sets the ranges of the existing histogram $h to cover the ranges $xmin to $xmax and $ymin to $ymax uniformly. The values of the histogram bins are reset to zero.
 119
 120 =item C<gsl_histogram2d_memcpy($dest, $src)> - This function copies the histogram $src into the pre-existing histogram $dest, making $dest into an exact copy of $src. The two histograms must be of the same size.
 121
 122 =item C<gsl_histogram2d_clone($src)> - This function returns a pointer to a newly created histogram which is an exact copy of the histogram $src.
 123
 124 =item C<gsl_histogram2d_max_val($h)> - This function returns the maximum value contained in the histogram bins.
 125
 126 =item C<gsl_histogram2d_max_bin($h)> - This function finds the indices of the bin containing the maximum value in the histogram $h and returns the result in this order : 0 if the operation succeded, 1 otherwise, i and j. In the case where several bins contain the same maximum value the first bin found is returned.
 127
 128 =item C<gsl_histogram2d_min_val($h)> - This function returns the minimum value contained in the histogram bins.
 129
 130 =item C<gsl_histogram2d_min_bin($h)> - This function finds the indices of the bin containing the minimum value in the histogram $h and returns the result in this order : 0 if the operation succeded, 1 otherwise, i and j. In the case where several bins contain the same minimum value the first bin found is returned.
 131
 132 =item C<gsl_histogram2d_xmean($h)> - This function returns the mean of the histogrammed x variable, where the histogram is regarded as a probability distribution. Negative bin values are ignored for the purposes of this calculation.
 133
 134 =item C<gsl_histogram2d_ymean($h)> - This function returns the mean of the histogrammed y variable, where the histogram is regarded as a probability distribution. Negative bin values are ignored for the purposes of this calculation.
 135
 136 =item C<gsl_histogram2d_xsigma($h)> - This function returns the standard deviation of the histogrammed x variable, where the histogram is regarded as a probability distribution. Negative bin values are ignored for the purposes of this calculation.
 137
 138 =item C<gsl_histogram2d_ysigma($h)> - This function returns the standard deviation of the histogrammed y variable, where the histogram is regarded as a probability distribution. Negative bin values are ignored for the purposes of this calculation.
 139
 140 =item C<gsl_histogram2d_cov($h)> - This function returns the covariance of the histogrammed x and y variables, where the histogram is regarded as a probability distribution. Negative bin values are ignored for the purposes of this calculation.
 141
 142 =item C<gsl_histogram2d_sum($h)> - This function returns the sum of all bin values. Negative bin values are included in the sum.
 143
 144 =item C<gsl_histogram2d_equal_bins_p($h1, $h2)> - This function returns 1 if all the individual bin ranges of the two histograms are identical, and 0 otherwise.
 145
 146 =item C<gsl_histogram2d_add($h1, $h2)> - This function adds the contents of the bins in histogram $h2 to the corresponding bins of histogram $h1, i.e. h'_1(i,j) = h_1(i,j) + h_2(i,j). The two histograms must have identical bin ranges.
 147
 148 =item C<gsl_histogram2d_sub($h1, $h2)> - This function subtracts the contents of the bins in histogram $h2 from the corresponding bins of histogram $h1, i.e. h'_1(i,j) = h_1(i,j) - h_2(i,j). The two histograms must have identical bin ranges.
 149
 150 =item C<gsl_histogram2d_mul($h1, $h2)> - This function multiplies the contents of the bins of histogram $h1 by the contents of the corresponding bins in histogram $h2, i.e. h'_1(i,j) = h_1(i,j) * h_2(i,j). The two histograms must have identical bin ranges.
 151
 152 =item C<gsl_histogram2d_div($h1, $h2)> - This function divides the contents of the bins of histogram $h1 by the contents of the corresponding bins in histogram $h2, i.e. h'_1(i,j) = h_1(i,j) / h_2(i,j). The two histograms must have identical bin ranges.
 153
 154 =item C<gsl_histogram2d_scale($h, $scale)> - This function multiplies the contents of the bins of histogram $h by the constant scale, i.e. h'_1(i,j) = h_1(i,j) $scale.
 155
 156 =item C<gsl_histogram2d_shift($h, $offset)> - This function shifts the contents of the bins of histogram $h by the constant offset, i.e. h'_1(i,j) = h_1(i,j) + $offset.
 157
 158 =item C<gsl_histogram2d_fwrite($stream $h)> - This function writes the ranges and bins of the histogram $h to the stream $stream (opened with the gsl_fopen function from the Math::GSL module) in binary format. The return value is 0 for success and $GSL_EFAILED if there was a problem writing to the file. Since the data is written in the native binary format it may not be portable between different architectures.
 159
 160 =item C<gsl_histogram2d_fread($stream $h)> - This function reads into the histogram $h from the stream $stream (opened with the gsl_fopen function from the Math::GSL module) in binary format. The histogram $h must be preallocated with the correct size since the function uses the number of x and y bins in $h to determine how many bytes to read. The return value is 0 for success and $GSL_EFAILED if there was a problem reading from the file. The data is assumed to have been written in the native binary format on the same architecture.
 161
 162 =item C<gsl_histogram2d_fprintf($stream, $h, $range_format, $bin_format)> - This function writes the ranges and bins of the histogram $h line-by-line to the stream $stream (opened with the gsl_fopen function from the Math::GSL module) using the format specifiers $range_format and $bin_format. These should be one of the %g, %e or %f formats for floating point numbers. The function returns 0 for success and $GSL_EFAILED if there was a problem writing to the file. The histogram output is formatted in five columns, and the columns are separated by spaces, like this,
 163
 164 =over
 165
 166 =item xrange[0] xrange[1] yrange[0] yrange[1] bin(0,0)
 167
 168 =item xrange[0] xrange[1] yrange[1] yrange[2] bin(0,1)
 169
 170 =item xrange[0] xrange[1] yrange[2] yrange[3] bin(0,2)
 171
 172 =item ....
 173
 174 =item xrange[0] xrange[1] yrange[ny-1] yrange[ny] bin(0,ny-1)
 175
 176 =item
 177
 178 =item xrange[1] xrange[2] yrange[0] yrange[1] bin(1,0)
 179
 180 =item xrange[1] xrange[2] yrange[1] yrange[2] bin(1,1)
 181
 182 =item xrange[1] xrange[2] yrange[1] yrange[2] bin(1,2)
 183
 184 =item ....
 185
 186 =item xrange[1] xrange[2] yrange[ny-1] yrange[ny] bin(1,ny-1)
 187
 188 =item
 189
 190 =item ....
 191
 192 =item
 193
 194 =item xrange[nx-1] xrange[nx] yrange[0] yrange[1] bin(nx-1,0)
 195
 196 =item xrange[nx-1] xrange[nx] yrange[1] yrange[2] bin(nx-1,1)
 197
 198 =item xrange[nx-1] xrange[nx] yrange[1] yrange[2] bin(nx-1,2)
 199
 200 =item ....
 201
 202 =item xrange[nx-1] xrange[nx] yrange[ny-1] yrange[ny] bin(nx-1,ny-1)
 203
 204 =back
 205
 206 Each line contains the lower and upper limits of the bin and the contents of the bin. Since the upper limits of the each bin are the lower limits of the neighboring bins there is duplication of these values but this allows the histogram to be manipulated with line-oriented tools.
 207
 208 =item C<gsl_histogram2d_fscanf($stream, $h)> - This function reads formatted data from the stream $stream (opened with the gsl_fopen function from the Math::GSL module) into the histogram $h. The data is assumed to be in the five-column format used by gsl_histogram2d_fprintf. The histogram $h must be preallocated with the correct lengths since the function uses the sizes of $h to determine how many numbers to read. The function returns 0 for success and $GSL_EFAILED if there was a problem reading from the file.
 209
 210 =item C<gsl_histogram2d_pdf_alloc($nx, $ny)> - This function allocates memory for a two-dimensional probability distribution of size $nx-by-$ny and returns a pointer to a newly initialized gsl_histogram2d_pdf struct. If insufficient memory is available a null pointer is returned and the error handler is invoked with an error code of $GSL_ENOMEM.
 211
 212 =item C<gsl_histogram2d_pdf_init($p, $h)> - This function initializes the two-dimensional probability distribution calculated $p from the histogram $h. If any of the bins of $h are negative then the error handler is invoked with an error code of GSL_EDOM because a probability distribution cannot contain negative values.
 213
 214 =item C<gsl_histogram2d_pdf_free($p)> - This function frees the two-dimensional probability distribution function $p and all of the memory associated with it.
 215
 216 =item C<gsl_histogram2d_pdf_sample($p, $r1, $r2)> - This function uses two uniform random numbers between zero and one, $r1 and $r2, to compute a single random sample from the two-dimensional probability distribution p. The function returns 0 if the operation succeeded, 1 otherwise followed by the x and y values of the sample.
 217
 218 =back
 219
 220 =head1 AUTHORS
 221
 222 Jonathan Leto <jonathan@leto.net> and Thierry Moisan <thierry.moisan@gmail.com>
 223
 224 =head1 COPYRIGHT AND LICENSE
 225
 226 Copyright (C) 2008-2009 Jonathan Leto and Thierry Moisan
 227
 228 This program is free software; you can redistribute it and/or modify it
 229 under the same terms as Perl itself.
 230
 231 =cut
 232
 233 %}
 234