GriddedWake.h

#pragma once
 
#include <marenv/grid/UniformGrid3.h>
#include <marenv/marenv.h>
#include <marenv/wake/WakeBlend.h>
#include <marenv/wake/WakeField.h>
#include <marenv/wake/WakeTypes.h>
 
#include <memory>
#include <string>
 
namespace marenv::wake
{
struct GriddedSnapshot
{
    WakeFrame frame; 
    grid::UniformGrid3<double> logRatio; 
    double sourceVelocity[3]; 
    std::string warning; 
 
    bool Deficit(double t, const double pos[3], double& deficit) const;
 
    bool Sample(double t, const double pos[3], WakeSample& out) const;
 
    void ToFrame(double t, const double pos[3], double xi[3]) const;
};
 
class GriddedWake : public WakeField
{
public:
    GriddedWake()           = default;
    ~GriddedWake() override = default;
 
    void Publish(std::shared_ptr<const GriddedSnapshot> snapshot, double blendTime = 0.0);
 
    Status Sample(double t, const double pos[3], WakeSample& out) const override;
 
    WakeSampleParts SampleParts(double t, const double pos[3]) const override;
 
private:
    using Blend = WakeBlend<GriddedSnapshot>;
 
    static Status SampleState(const Blend::State* state, double t, const double pos[3], WakeSample& out);
 
    Blend m_blend; 
};
} // namespace marenv::wake