pressure_intensity.h

#pragma once

#include <cmath>

namespace wayverb {
namespace core {

template <typename T, typename U>
inline auto pressure_to_intensity(T pressure, U Z) {
    using std::copysign;
    return copysign(pressure * pressure / Z, pressure);
}

template <typename T, typename U>
inline auto intensity_to_pressure(T intensity, U Z) {
    using std::copysign;
    using std::sqrt;
    using std::abs;
    return copysign(sqrt(abs(intensity * Z)), intensity);
}

double intensity_for_distance(double distance);
double pressure_for_distance(double distance, double acoustic_impedance);

}  // namespace core
}  // namespace wayverb