+// Normalize so that the sum becomes one. Note that we do it twice;
+// this sometimes helps a tiny little bit when we have many samples.
+template<class T>
+void normalize_sum(Tap<T>* vals, unsigned num)
+{
+ for (int normalize_pass = 0; normalize_pass < 2; ++normalize_pass) {
+ double sum = 0.0;
+ for (unsigned i = 0; i < num; ++i) {
+ sum += to_fp64(vals[i].weight);
+ }
+ for (unsigned i = 0; i < num; ++i) {
+ vals[i].weight = from_fp64<T>(to_fp64(vals[i].weight) / sum);
+ }
+ }
+}
+