BathymetryField.h

#pragma once
 
#include <marenv/marenv.h>
 
#include <cmath>
 
namespace marenv::seafloor
{
 
class BathymetryField : private marenv::NoCopyNoMove
{
public:
 
    virtual ~BathymetryField() = default;
 
    virtual Status GetDepth(const double pos[2], double& depthOut) = 0;
 
    virtual Status ComputeSurfaceUnitNormal(const double pos[2], double unitNormalOut[3]) = 0;
 
    virtual Status GetSeabedSigma(const double pos[3], double& sigmaOut) = 0;
    virtual Status GetSeabedTau(const double pos[3], double& tauOut) = 0;
    virtual Status GetSeabedRho(const double pos[3], double& rhoOut) = 0;
 
    virtual Status ComputeRayIntersection(const double origin[3], const double direction[3], double& distanceOut, double maxSearchDistance = 100, double startSearchDistance = -1.0, double stepSize = 1.0, double convergenceTolerance = 0.01);
 
protected:
 
    BathymetryField() = default;
 
    static bool IsValidRay(const double origin[3], double directionNorm, double maxSearchDistance);
 
    static double ValidTolerance(double convergenceTolerance);
 
private:
    static double FlatSeabedDistance(double altitude, const double unitDirection[3]);
};
} // namespace marenv::seafloor