More robust calculation for higher sampling rates

Adding a c++ wrapper for faster determing the number of reflection paths used in shoebox rooms, this may saves the memory usage for large image matrics