Kernel.h

#pragma once
 
#include <marenv/wake/WakeTypes.h>
 
namespace marenv::wake
{
double ElementaryDeficit(const ElementarySource& source, const KernelParams& params, const double pos[3]);
 
void Validate(const KernelParams& params);
 
struct ImagePlanes
{
    int count; 
    double z[2]; 
};
 
ImagePlanes MirrorPlanes(const KernelParams& params);
 
struct ImageMap
{
    double offset; 
    double sign; 
 
    double QueryDepth(double z) const { return offset + sign * z; }
 
    double ImageDepth(double z) const { return sign * (z - offset); }
};
 
struct ImageSet
{
    int count; 
    ImageMap maps[4]; 
};
 
ImageSet MirrorImages(const KernelParams& params);
 
struct KernelCentreline
{
    double amplitude; 
    double sigma; 
};
 
class KernelProfile
{
public:
    KernelProfile(const ElementarySource& source, const KernelParams& params);
 
    double XStart() const { return m_xStart; }
 
    KernelCentreline At(double x) const;
 
    double PeakAmplitude(double xFrom, double xTo) const;
 
private:
    double CentrelineAmplitude(double x, double pinnedAmplitude) const;
 
    double LinearSigma(double x) const;
 
    double Sigma(double x) const;
 
    double m_diameter; 
    double m_expansion; 
    double m_xStart; 
    double m_thrustCoefficient; 
    double m_nearWakeDeficit; 
    double m_epsilon; 
    double m_plateau; 
    double m_plateauEnd; 
    double m_blendEnd; 
    bool m_porous; 
    double m_farWakeStart; 
    double m_farWakeSigma; 
};
 
double ElementaryDeficit(const ElementarySource& source, const KernelProfile& profile, const double pos[3]);
} // namespace marenv::wake