CAMLprim value ml_gsl_histogram_equal_bins_p(value vh1, value vh2) { gsl_histogram h1, h2; histo_of_val(&h1, vh1); histo_of_val(&h2, vh2); return Val_bool(gsl_histogram_equal_bins_p(&h1, &h2)); }
/* * gsl_histogram_add: * add two histograms */ int gsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2) { size_t i; if (!gsl_histogram_equal_bins_p (h1, h2)) { GSL_ERROR ("histograms have different binning", GSL_EINVAL); } for (i = 0; i < h1->n; i++) { h1->bin[i] += h2->bin[i]; } return GSL_SUCCESS; }