THeavisideStep2D.h

Go to the documentation of this file.

#pragma once
 
#include "Math/Function/TMathFunction2D.h"
#include "Math/TVector2.h"
 
#include <cmath>
#include <algorithm>
 
namespace ph::math
{
 
template<typename Value>
class THeavisideStep2D : public TMathFunction2D<Value>
{
public:
    static auto makeHorizontal() -> THeavisideStep2D;
    static auto makeVertical() -> THeavisideStep2D;
 
    explicit THeavisideStep2D(TVector2<Value> normal);
 
    Value evaluate(Value x, Value y) const override;
 
private:
    TVector2<Value> m_unitNormal;
};
 
template<typename Value>
inline auto THeavisideStep2D<Value>::makeHorizontal()
-> THeavisideStep2D
{
    return THeavisideStep2D({1, 0});
}
 
template<typename Value>
inline auto THeavisideStep2D<Value>::makeVertical()
-> THeavisideStep2D
{
    return THeavisideStep2D({0, 1});
}
 
template<typename Value>
inline THeavisideStep2D<Value>::THeavisideStep2D(TVector2<Value> normal)
    : m_unitNormal(normal.normalize())
{}
 
template<typename Value>
inline Value THeavisideStep2D<Value>::evaluate(const Value x, const Value y) const
{
    const Value projectedSignedLengthOnNormal = m_unitNormal.dot({x, y});
    return projectedSignedLengthOnNormal > 0 
        ? 1 
        : (projectedSignedLengthOnNormal < 0 ? 0 : static_cast<Value>(0.5));
}
 
}// end namespace ph::math