Kernel.h

#pragma once
 
#include <marenv/wake/WakeField.h>
#include <marenv/wake/WakeTypes.h>
 
namespace marenv::wake
{
double ElementaryDeficit(const ElementarySource& source, const KernelParams& params, const double pos[3]);
 
double OnsetRamp(double x, double xStart, double length);
 
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
{
    static constexpr int kMaxImages = 4; 
 
    int count; 
    ImageMap maps[kMaxImages]; 
};
 
ImageSet MirrorImages(const KernelParams& params);
 
struct ImageQueries
{
    int count; 
    double point[1 + ImageSet::kMaxImages][3]; 
    double sign[1 + ImageSet::kMaxImages]; 
};
 
ImageQueries MakeImageQueries(const ImageSet& images, const double pos[3], const double shift[3]);
 
struct ImageRegions
{
    int count; 
    WakeQueryRegion region[1 + ImageSet::kMaxImages]; 
};
 
ImageRegions MakeImageRegions(const ImageSet& images, const WakeQueryRegion& region);
 
ElementarySource ImageOf(const ImageMap& map, const ElementarySource& source);
 
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 UnrampedAmplitude(double x) const;
 
    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_rampLength; 
    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