diff --git a/src/control-variates/regions-integrator-parallel-control-variates.h b/src/control-variates/regions-integrator-parallel-control-variates.h
index bbedfb1..a5997f8 100644
--- a/src/control-variates/regions-integrator-parallel-control-variates.h
+++ b/src/control-variates/regions-integrator-parallel-control-variates.h
@@ -40,7 +40,7 @@ class RegionsIntegratorParallelControlVariates {
                 Range<Float,DIM> binrange = range;
                 for (std::size_t i=0;i<DIMBINS;++i)
                     binrange = binrange.subrange_dimension(i,range.min(i)+pos[i]*drange[i],range.min(i)+(pos[i]+1)*drange[i]);
-                std::vector<Range<Float,DIM>> region_bin_ranges(perbin[pos].size(),range_primary<DIM>());
+                std::vector<Range<Float,DIM>> region_bin_ranges(perbin[pos].size(),range_primary<DIM,Float>());
                 //This is for the control variate and for the region_bin_ranges
                 for (std::size_t i = 0; i<region_bin_ranges.size(); ++i) {
                     region_bin_ranges[i] = binrange.intersection(perbin[pos][i]->range());