void QwtHistogram::loadDataFromMatrix() {
if (!d_matrix)
return;
int size = d_matrix->numRows() * d_matrix->numCols();
const double *data = d_matrix->matrixModel()->dataVector();
int n;
gsl_histogram *h;
if (d_autoBin) {
double min, max;
d_matrix->range(&min, &max);
d_begin = floor(min);
d_end = ceil(max);
d_bin_size = 1.0;
n = static_cast<int>(floor((d_end - d_begin) / d_bin_size));
if (!n)
return;
h = gsl_histogram_alloc(n);
if (!h)
return;
gsl_histogram_set_ranges_uniform(h, floor(min), ceil(max));
} else {
n = static_cast<int>((d_end - d_begin) / d_bin_size + 1);
if (!n)
return;
h = gsl_histogram_alloc(n);
if (!h)
return;
double *range = new double[n + 2];
for (int i = 0; i <= n + 1; i++)
range[i] = d_begin + i * d_bin_size;
gsl_histogram_set_ranges(h, range, n + 1);
delete[] range;
}
for (int i = 0; i < size; i++)
gsl_histogram_increment(h, data[i]);
QVarLengthArray<double> X(n), Y(n); // stores ranges (x) and bins (y)
for (int i = 0; i < n; i++) {
Y[i] = gsl_histogram_get(h, i);
double lower, upper;
gsl_histogram_get_range(h, i, &lower, &upper);
X[i] = lower;
}
setData(X.data(), Y.data(), n);
d_mean = gsl_histogram_mean(h);
d_standard_deviation = gsl_histogram_sigma(h);
d_min = gsl_histogram_min_val(h);
d_max = gsl_histogram_max_val(h);
gsl_histogram_free(h);
}
int
main (int argc, char **argv)
{
double a = 0.0, b = 1.0;
size_t n = 10;
if (argc != 3 && argc !=4)
{
printf ("Usage: gsl-histogram xmin xmax [n]\n");
printf (
"Computes a histogram of the data on stdin using n bins from xmin to xmax.\n"
"If n is unspecified then bins of integer width are used.\n");
exit (0);
}
a = atof (argv[1]);
b = atof (argv[2]);
if (argc == 4)
{
n = atoi (argv[3]);
}
else
{
n = (int)(b - a) ;
}
{
double x;
gsl_histogram *h = gsl_histogram_alloc (n);
gsl_histogram_set_ranges_uniform (h, a, b);
while (fscanf(stdin, "%lg", &x) == 1)
{
gsl_histogram_increment(h, x);
}
#ifdef DISPLAY_STATS
{
double mean = gsl_histogram_mean (h);
double sigma = gsl_histogram_sigma (h);
fprintf (stdout, "# mean = %g\n", mean);
fprintf (stdout, "# sigma = %g\n", sigma);
}
#endif
gsl_histogram_fprintf (stdout, h, "%g", "%g");
gsl_histogram_free (h);
}
return 0;
}
void QwtHistogram::initData(double *Y, int size)
{
if(size<2 || (size == 2 && Y[0] == Y[1]))
{//non valid histogram data
double x[2], y[2];
for (int i = 0; i<2; i++ ){
y[i] = 0;
x[i] = 0;
}
setData(x, y, 2);
return;
}
int n = 10;//default value
QVarLengthArray<double> x(n), y(n);//double x[n], y[n]; //store ranges (x) and bins (y)
gsl_histogram * h = gsl_histogram_alloc (n);
if (!h)
return;
gsl_vector *v;
v = gsl_vector_alloc (size);
for (int i = 0; i<size; i++ )
gsl_vector_set (v, i, Y[i]);
double min, max;
gsl_vector_minmax (v, &min, &max);
gsl_vector_free (v);
d_begin = floor(min);
d_end = ceil(max);
gsl_histogram_set_ranges_uniform (h, floor(min), ceil(max));
for (int i = 0; i<size; i++ )
gsl_histogram_increment (h, Y[i]);
for (int i = 0; i<n; i++ ){
y[i] = gsl_histogram_get (h, i);
double lower, upper;
gsl_histogram_get_range (h, i, &lower, &upper);
x[i] = lower;
}
setData(x.data(), y.data(), n);//setData(x, y, n);
d_bin_size = (d_end - d_begin)/(double)n;
d_autoBin = true;
d_mean = gsl_histogram_mean(h);
d_standard_deviation = gsl_histogram_sigma(h);
d_min = gsl_histogram_min_val(h);
d_max = gsl_histogram_max_val(h);
gsl_histogram_free (h);
}
void QwtHistogram::loadData()
{
if (d_matrix){
loadDataFromMatrix();
return;
}
int r = abs(d_end_row - d_start_row) + 1;
QVarLengthArray<double> Y(r);
int ycol = d_table->colIndex(title().text());
int size = 0;
for (int i = 0; i<r; i++ ){
QString yval = d_table->text(i, ycol);
if (!yval.isEmpty()){
bool valid_data = true;
Y[size] = ((Graph *)plot())->locale().toDouble(yval, &valid_data);
if (valid_data)
size++;
}
}
if(size < 2 || (size==2 && Y[0] == Y[1])){//non valid histogram
double X[2];
Y.resize(2);
for (int i = 0; i<2; i++ ){
Y[i] = 0;
X[i] = 0;
}
setData(X, Y.data(), 2);
return;
}
int n;
gsl_histogram *h;
if (d_autoBin){
n = 10;
h = gsl_histogram_alloc (n);
if (!h)
return;
gsl_vector *v = gsl_vector_alloc (size);
for (int i = 0; i<size; i++ )
gsl_vector_set (v, i, Y[i]);
double min, max;
gsl_vector_minmax (v, &min, &max);
gsl_vector_free (v);
d_begin = floor(min);
d_end = ceil(max);
d_bin_size = (d_end - d_begin)/(double)n;
gsl_histogram_set_ranges_uniform (h, floor(min), ceil(max));
} else {
n = int((d_end - d_begin)/d_bin_size + 1);
h = gsl_histogram_alloc (n);
if (!h)
return;
double *range = new double[n+2];
for (int i = 0; i<= n+1; i++ )
range[i] = d_begin + i*d_bin_size;
gsl_histogram_set_ranges (h, range, n+1);
delete[] range;
}
for (int i = 0; i<size; i++ )
gsl_histogram_increment (h, Y[i]);
#ifdef Q_CC_MSVC
QVarLengthArray<double> X(n); //stores ranges (x) and bins (y)
#else
double X[n]; //stores ranges (x) and bins (y)
#endif
Y.resize(n);
for (int i = 0; i<n; i++ ){
Y[i] = gsl_histogram_get (h, i);
double lower, upper;
gsl_histogram_get_range (h, i, &lower, &upper);
X[i] = lower;
}
#ifdef Q_CC_MSVC
setData(X.data(), Y.data(), n);
#else
setData(X, Y.data(), n);
#endif
d_mean = gsl_histogram_mean(h);
d_standard_deviation = gsl_histogram_sigma(h);
d_min = gsl_histogram_min_val(h);
d_max = gsl_histogram_max_val(h);
gsl_histogram_free (h);
}
CAMLprim value ml_gsl_histogram_sigma(value vh)
{
gsl_histogram h;
histo_of_val(&h, vh);
return copy_double(gsl_histogram_sigma(&h));
}