calibration.h

#pragma once

#include <cmath>
#include <iostream>


namespace wayverb {
namespace waveguide {

inline double distance_for_unit_intensity() {
    //  Intensity for distance is given by 1 / (4 * pi * distance * distance)
    using std::sqrt;
    return sqrt(1.0 / (4 * M_PI));
}

inline double distance_for_unit_pressure(double acoustic_impedance) {
    using std::sqrt;
    return sqrt(acoustic_impedance / (4 * M_PI));
}

template <typename T, typename U>
inline auto rectilinear_calibration_factor(T grid_spacing,
                                           U acoustic_impedance) {
    return distance_for_unit_pressure(acoustic_impedance) /
                     (0.3405 * grid_spacing);
}

}  // namespace waveguide
}  // namespace wayverb