Fixes the new ray tracing cpp function after merge.

Author: fakufaku

pyroomacoustics/directivities/base.py

@@ -110,7 +110,7 @@ def sample_rays(self, n_rays, rng=None):
 
         Returns
         -------
-        ray_directions: numpy.ndarray, shape (n_dim, n_rays)
+        ray_directions: numpy.ndarray, shape (n_rays, n_dim)
             An array containing the unit vectors in its columns
         energies: numpy.ndarray, shape (n_rays, n_bands)
             An energy carried per ray so that the expectation over all the rays

pyroomacoustics/libroom_src/room.cpp

@@ -953,7 +953,7 @@ void Room<D>::ray_tracing(
     const Eigen::Matrix<float, Eigen::Dynamic, D> &unit_vectors,
     const Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> &energies,
     const Vectorf<D> &source_pos) {
-  float energy_0 = 2.f / unit_vectors.rows();
+  float energy_0 = energy_0_numerator / unit_vectors.rows();
 
   if (unit_vectors.rows() != energies.rows()) {
     throw std::runtime_error("Error: The same number of rays and energies is expected.");