|
Packit |
67cb25 |
/* gsl_histogram2d_maxval.c
|
|
Packit |
67cb25 |
* Copyright (C) 2000 Simone Piccardi
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This library is free software; you can redistribute it and/or
|
|
Packit |
67cb25 |
* modify it under the terms of the GNU General Public License as
|
|
Packit |
67cb25 |
* published by the Free Software Foundation; either version 3 of the
|
|
Packit |
67cb25 |
* License, or (at your option) any later version.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is distributed in the hope that it will be useful,
|
|
Packit |
67cb25 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Packit |
67cb25 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Packit |
67cb25 |
* General Public License for more details.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* You should have received a copy of the GNU General Public License along
|
|
Packit |
67cb25 |
* with this library; if not, write to the Free Software Foundation, Inc.,
|
|
Packit |
67cb25 |
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
/***************************************************************
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* File gsl_histogram2d_maxval.c:
|
|
Packit |
67cb25 |
* Routine to find maximum and minumum content of a 2D hisogram.
|
|
Packit |
67cb25 |
* Need GSL library and header.
|
|
Packit |
67cb25 |
* Contains the routines:
|
|
Packit |
67cb25 |
* gsl_histogram2d_max_val find max content values
|
|
Packit |
67cb25 |
* gsl_histogram2d_min_val find min content values
|
|
Packit |
67cb25 |
* gsl_histogram2d_bin_max find coordinates of max contents bin
|
|
Packit |
67cb25 |
* gsl_histogram2d_bin_min find coordinates of min contents bin
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Author: S. Piccardi
|
|
Packit |
67cb25 |
* Jan. 2000
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
***************************************************************/
|
|
Packit |
67cb25 |
#include <config.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_errno.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_histogram2d.h>
|
|
Packit |
67cb25 |
/*
|
|
Packit |
67cb25 |
* Return the maximum contents value of a 2D histogram
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
double
|
|
Packit |
67cb25 |
gsl_histogram2d_max_val (const gsl_histogram2d * h)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
const size_t nx = h->nx;
|
|
Packit |
67cb25 |
const size_t ny = h->ny;
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
double max = h->bin[0 * ny + 0];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < nx * ny; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (h->bin[i] > max)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
max = h->bin[i];
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
return max;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/*
|
|
Packit |
67cb25 |
* Find the bin index for maximum value of a 2D histogram
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
gsl_histogram2d_max_bin (const gsl_histogram2d * h, size_t * imax_out, size_t * jmax_out)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
const size_t nx = h->nx;
|
|
Packit |
67cb25 |
const size_t ny = h->ny;
|
|
Packit |
67cb25 |
size_t imax = 0, jmax = 0, i, j;
|
|
Packit |
67cb25 |
double max = h->bin[0 * ny + 0];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < nx; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < ny; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x = h->bin[i * ny + j];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (x > max)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
max = x;
|
|
Packit |
67cb25 |
imax = i;
|
|
Packit |
67cb25 |
jmax = j;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
*imax_out = imax;
|
|
Packit |
67cb25 |
*jmax_out = jmax;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/*
|
|
Packit |
67cb25 |
* Return the minimum contents value of a 2D histogram
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
double
|
|
Packit |
67cb25 |
gsl_histogram2d_min_val (const gsl_histogram2d * h)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
const size_t nx = h->nx;
|
|
Packit |
67cb25 |
const size_t ny = h->ny;
|
|
Packit |
67cb25 |
size_t i;
|
|
Packit |
67cb25 |
double min = h->bin[0 * ny + 0];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < nx * ny; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
if (h->bin[i] < min)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
min = h->bin[i];
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return min;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
/*
|
|
Packit |
67cb25 |
* Find the bin index for minimum value of a 2D histogram
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
void
|
|
Packit |
67cb25 |
gsl_histogram2d_min_bin (const gsl_histogram2d * h, size_t * imin_out, size_t * jmin_out)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
const size_t nx = h->nx;
|
|
Packit |
67cb25 |
const size_t ny = h->ny;
|
|
Packit |
67cb25 |
size_t imin = 0, jmin = 0, i, j;
|
|
Packit |
67cb25 |
double min = h->bin[0 * ny + 0];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
for (i = 0; i < nx; i++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
for (j = 0; j < ny; j++)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
double x = h->bin[i * ny + j];
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
if (x < min)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
min = x;
|
|
Packit |
67cb25 |
imin = i;
|
|
Packit |
67cb25 |
jmin = j;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
*imin_out = imin;
|
|
Packit |
67cb25 |
*jmin_out = jmin;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|