SourceSetWake.h

#pragma once
 
#include <marenv/marenv.h>
#include <marenv/wake/Kernel.h>
#include <marenv/wake/LineKernel.h>
#include <marenv/wake/WakeBlend.h>
#include <marenv/wake/WakeField.h>
#include <marenv/wake/WakeTypes.h>
 
#include <memory>
#include <vector>
 
namespace marenv::wake
{
class SourceSetWake : public WakeField
{
public:
    SourceSetWake() = default;
 
    explicit SourceSetWake(const KernelParams& params);
 
    ~SourceSetWake() override = default;
 
    void Publish(SourceSet sources, const double velocity[3], double t0, double blendTime, double seabedDepth);
 
    Status Sample(double t, const double pos[3], WakeSample& out) const override;
 
    bool MayReach(const WakeQueryRegion& region) const override;
 
    WakeSampleParts SampleParts(double t, const double pos[3]) const override;
 
private:
    struct Snapshot
    {
        std::vector<ElementarySource> sources; 
        double velocity[3]; 
        KernelParams params; 
        double t0; 
        ImageSet images; 
        std::vector<KernelProfile> profiles; 
        std::vector<LineKernelProfile> lineProfiles; 
 
        bool Sample(double t, const double pos[3], WakeSample& out) const;
 
        bool MayReach(const WakeQueryRegion& region) const;
    };
 
    static Status SampleState(const WakeBlend<Snapshot>::State* state, double t, const double pos[3], WakeSample& out);
 
    static std::shared_ptr<const Snapshot> MakeSnapshot(SourceSet sources, const double velocity[3], const KernelParams& params, double t0);
 
    KernelParams m_params; 
    WakeBlend<Snapshot> m_blend; 
};
} // namespace marenv::wake