Got width sections to approach constant signal as width sections incr…

…ease

bash/calc_base_signal

@@ -0,0 +1,10 @@
+#!/usr/bin/env bash
+
+rm ~/Desktop/intensities; 
+for n in `seq 30`; do 
+    generate_image /home/tclose/fouts/params/fibre/tract/test_configurations/x-large.tct ~/Desktop/intens_calc/$n.mif -clean --exp_num_width $n; 
+    maxb0 ~/Desktop/intens_calc/$n.mif >> ~/Desktop/intensities; 
+    echo " " >> ~/Desktop/intensities; 
+done
+python /home/tclose/fouts/python/plot_base_signal.py
+

python/plot_base_signal.py

@@ -0,0 +1,6 @@
+#!/usr/bin/env python
+import numpy as np
+import matplotlib.pyplot as plt
+d = np.loadtxt('/home/tclose/Desktop/intensities')
+plt.plot(d)
+plt.show()

src/bts/fibre/tractlet.cpp

@@ -353,7 +353,7 @@ namespace BTS {
             double length_fraction = 1.0 / (double)num_length_sections;
             double width_fraction = 1.0 / (double)num_width_sections;
 
-            double intensity_scale = 1.0 / (double)(MR::Math::pow2(num_width_sections));
+            double intensity_scale = 1.0 / (double)(width_fractions.rows());
 
             const MR::Math::Matrix<double>& position_matrix = Strand::position_matrix(
                     num_length_sections, this->degree());