void Histogram1D::multiply_with_ratio_exponented(const Histogram1D & nominator, const Histogram1D & denominator, double exponent){ check_compatibility(nominator); check_compatibility(denominator); const double * n_data = nominator.get_data(); const double * d_data = denominator.get_data(); double * data = get_data(); const size_t n = size(); for(size_t i=0; i<n; i++){ if(d_data[i]>0.0) data[i] *= pow(n_data[i] / d_data[i], exponent); } }
void Histogram1D::operator*=(const Histogram1D & h) { check_compatibility(h); const double * hdata = h.get_data(); double * data = get_data(); const size_t n = size(); for (size_t i = 0; i < n; ++i) { data[i] *= hdata[i]; } }