transform.hpp

// Copyright Take Vos 2021.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt)
 
#pragma once
 
#include "matrix2.hpp"
#include "matrix3.hpp"
#include "translate2.hpp"
#include "translate3.hpp"
#include "rotate2.hpp"
#include "rotate3.hpp"
#include "scale2.hpp"
#include "scale3.hpp"
#include "perspective.hpp"
#include "vector2.hpp"
#include "vector3.hpp"
#include "point2.hpp"
#include "point3.hpp"
#include "extent2.hpp"
#include "extent3.hpp"
#include "aarectangle.hpp"
#include "rectangle.hpp"
#include "quad.hpp"
#include "circle.hpp"
#include "corner_radii.hpp"
#include "line_segment.hpp"
#include "../macros.hpp"
#include <type_traits>
#include <functional>
#include <exception>
#include <compare>
 
hi_export_module(hikogui.geometry : transform);
 
hi_export namespace hi { inline namespace v1 {
 
template<typename T>
struct transform2 : public std::false_type {};
template<typename T>
struct transform3 : public std::false_type {};
 
// clang-format off
template<> struct transform2<matrix2> : public std::true_type {};
template<> struct transform3<matrix3> : public std::true_type {};
template<> struct transform2<translate2> : public std::true_type {};
template<> struct transform3<translate3> : public std::true_type {};
template<> struct transform2<rotate2> : public std::true_type {};
template<> struct transform3<rotate3> : public std::true_type {};
template<> struct transform2<scale2> : public std::true_type {};
template<> struct transform3<scale3> : public std::true_type {};
template<> struct transform3<perspective> : public std::true_type {};
// clang-format on
 
template<typename T>
constexpr bool transform2_v = transform2<T>::value;
template<typename T>
constexpr bool transform3_v = transform3<T>::value;
 
template<typename T>
concept transformer2 = transform2_v<T>;
template<typename T>
concept transformer3 = transform3_v<T>;
template<typename T>
concept transformer = transformer2<T> or transformer3<T>;

[[nodiscard]] constexpr matrix2 operator*(matrix2 const& lhs, matrix2 const& rhs) noexcept
{
    return matrix2{lhs * get<0>(rhs), lhs * get<1>(rhs), lhs * get<2>(rhs), lhs * get<3>(rhs)};
}

[[nodiscard]] constexpr matrix3 operator*(matrix3 const& lhs, matrix3 const& rhs) noexcept
{
    return matrix3{lhs * get<0>(rhs), lhs * get<1>(rhs), lhs * get<2>(rhs), lhs * get<3>(rhs)};
}
 
[[nodiscard]] constexpr matrix2 operator*(translate2 const& lhs, matrix2 const& rhs) noexcept
{
    return matrix2{get<0>(rhs), get<1>(rhs), get<2>(rhs), get<3>(rhs) + f32x4{lhs}};
}
 
[[nodiscard]] constexpr matrix3 operator*(translate3 const& lhs, matrix3 const& rhs) noexcept
{
    return matrix3{get<0>(rhs), get<1>(rhs), get<2>(rhs), get<3>(rhs) + f32x4{lhs}};
}
 
[[nodiscard]] constexpr translate2 operator*(translate2 const& lhs, translate2 const& rhs) noexcept
{
    return translate2{f32x4{lhs} + f32x4{rhs}};
}
 
[[nodiscard]] constexpr translate3 operator*(translate3 const& lhs, translate3 const& rhs) noexcept
{
    return translate3{f32x4{lhs} + f32x4{rhs}};
}
 
[[nodiscard]] constexpr matrix2 operator*(translate2 const& lhs, scale2 const& rhs) noexcept
{
    hi_axiom(lhs.holds_invariant() && rhs.holds_invariant());
    return matrix2{f32x4{rhs}.x000(), f32x4{rhs}._0y00(), f32x4{rhs}._00z0(), f32x4{lhs}.xyz1()};
}
 
[[nodiscard]] constexpr matrix3 operator*(translate3 const& lhs, scale3 const& rhs) noexcept
{
    return matrix3{f32x4{rhs}.x000(), f32x4{rhs}._0y00(), f32x4{rhs}._00z0(), f32x4{lhs}.xyz1()};
}
 
[[nodiscard]] constexpr matrix2 operator*(translate2 const& lhs, rotate2 const& rhs) noexcept
{
    return lhs * matrix2(rhs);
}
 
[[nodiscard]] constexpr matrix3 operator*(translate3 const& lhs, rotate3 const& rhs) noexcept
{
    return lhs * matrix3(rhs);
}
 
[[nodiscard]] constexpr matrix2 operator*(scale2 const& lhs, translate2 const& rhs) noexcept
{
    hi_axiom(lhs.holds_invariant() && rhs.holds_invariant());
    return matrix2{f32x4{lhs}.x000(), f32x4{lhs}._0y00(), f32x4{lhs}._00z0(), f32x4{lhs} * f32x4{rhs}.xyz1()};
}
 
[[nodiscard]] constexpr matrix3 operator*(scale3 const& lhs, translate3 const& rhs) noexcept
{
    hi_axiom(lhs.holds_invariant() && rhs.holds_invariant());
    return matrix3{f32x4{lhs}.x000(), f32x4{lhs}._0y00(), f32x4{lhs}._00z0(), f32x4{lhs} * f32x4{rhs}.xyz1()};
}
 
[[nodiscard]] constexpr scale2 operator*(scale2 const& lhs, scale2 const& rhs) noexcept
{
    return scale2{f32x4{lhs} * f32x4{rhs}};
}
 
[[nodiscard]] constexpr scale3 operator*(scale3 const& lhs, scale3 const& rhs) noexcept
{
    return scale3{f32x4{lhs} * f32x4{rhs}};
}
 
/******************* Shape operations ******************/

[[nodiscard]] inline float operator*(matrix2 const& lhs, float const& rhs) noexcept
{
    // As if _col0 * rhs.xxxx() in operator*(f32x4 rhs)
    auto const abs_scale = hypot<0b0011>(get<0>(lhs) * f32x4::broadcast(rhs)).x();
 
    // We want to keep the sign of the original scaler, even if the matrix has rotation.
    return std::copysign(abs_scale, rhs);
}

[[nodiscard]] inline float operator*(matrix3 const &lhs, float const& rhs) noexcept
{
    // As if _col0 * rhs.xxxx() in operator*(f32x4 rhs)
    auto const abs_scale = hypot<0b0111>(get<0>(lhs) * f32x4::broadcast(rhs)).x();
 
    // We want to keep the sign of the original scaler, even if the matrix has rotation.
    return std::copysign(abs_scale, rhs);
}
 
[[nodiscard]] constexpr float operator*(translate2 const&, float const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr float operator*(translate3 const&, float const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr float operator*(scale2 const& lhs, float const& rhs) noexcept
{
    return lhs.x() * rhs;
}
 
[[nodiscard]] constexpr float operator*(scale3 const& lhs, float const& rhs) noexcept
{
    return lhs.x() * rhs;
}
 
[[nodiscard]] constexpr float operator*(rotate2 const&, float const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr float operator*(rotate3 const&, float const& rhs) noexcept
{
    return rhs;
}

[[nodiscard]] constexpr vector2 operator*(matrix2 const& lhs, vector2 const& rhs) noexcept
{
    return vector2{get<0>(lhs) * f32x4{rhs}.xxxx() + get<1>(lhs) * f32x4{rhs}.yyyy()};
}
 
[[nodiscard]] constexpr vector2 operator*(translate2 const&, vector2 const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr vector2 operator*(translate3 const&, vector2 const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr vector2 operator*(scale2 const& lhs, vector2 const& rhs) noexcept
{
    return vector2{f32x4{lhs} * f32x4{rhs}};
}
 
[[nodiscard]] constexpr vector2 operator*(scale3 const& lhs, vector2 const& rhs) noexcept
{
    return vector2{f32x4{lhs} * f32x4{rhs}};
}

[[nodiscard]] constexpr vector3 operator*(matrix3 const& lhs, vector3 const& rhs) noexcept
{
    return vector3{get<0>(lhs) * f32x4{rhs}.xxxx() + get<1>(lhs) * f32x4{rhs}.yyyy() + get<2>(lhs) * f32x4{rhs}.zzzz()};
}
 
[[nodiscard]] constexpr vector3 operator*(translate3 const&, vector3 const& rhs) noexcept
{
    return rhs;
}
 
[[nodiscard]] constexpr vector3 operator*(scale3 const& lhs, vector3 const& rhs) noexcept
{
    return vector3{f32x4{lhs} * f32x4{rhs}};
}

[[nodiscard]] constexpr point2 operator*(matrix2 const& lhs, point2 const& rhs) noexcept
{
    return point2{get<0>(lhs) * f32x4{rhs}.xxxx() + get<1>(lhs) * f32x4{rhs}.yyyy() + get<3>(lhs)};
}
 
[[nodiscard]] constexpr point2 operator*(translate2 const& lhs, point2 const& rhs) noexcept
{
    return point2{f32x4{lhs} + f32x4{rhs}};
}
 
[[nodiscard]] constexpr point2 operator*(scale2 const& lhs, point2 const& rhs) noexcept
{
    return point2{f32x4{lhs} * f32x4{rhs}};
}
 
[[nodiscard]] constexpr point2 operator*(scale3 const& lhs, point2 const& rhs) noexcept
{
    return point2{f32x4{lhs} * f32x4{rhs}};
}

[[nodiscard]] constexpr point3 operator*(matrix3 const& lhs, point3 const& rhs) noexcept
{
    return point3{
        get<0>(lhs) * f32x4{rhs}.xxxx() + get<1>(lhs) * f32x4{rhs}.yyyy() + get<2>(lhs) * f32x4{rhs}.zzzz() +
        get<3>(lhs) * f32x4{rhs}.wwww()};
}
 
[[nodiscard]] constexpr point3 operator*(translate3 const& lhs, point3 const& rhs) noexcept
{
    return point3{f32x4{lhs} + f32x4{rhs}};
}
 
[[nodiscard]] constexpr point3 operator*(scale3 const& lhs, point3 const& rhs) noexcept
{
    return point3{f32x4{lhs} * f32x4{rhs}};
}
 
[[nodiscard]] constexpr extent2 operator*(transformer2 auto const& lhs, extent2 const& rhs) noexcept
{
    return extent2{lhs *vector2{rhs}};
}
 
[[nodiscard]] constexpr extent3 operator*(transformer3 auto const& lhs, std::convertible_to<extent3> auto const& rhs) noexcept
{
    return extent3{lhs *vector3{rhs}};
}
 
[[nodiscard]] constexpr aarectangle operator*(translate2 const& lhs, aarectangle const& rhs) noexcept
{
    return aarectangle{lhs * get<0>(rhs), lhs * get<3>(rhs)};
}

[[nodiscard]] constexpr aarectangle operator*(scale2 const& lhs, aarectangle const& rhs) noexcept
{
    return aarectangle{lhs * get<0>(rhs), lhs * get<3>(rhs)};
}

[[nodiscard]] constexpr rectangle operator*(transformer auto const& lhs, aarectangle const& rhs) noexcept
{
    auto const rhs_ = rectangle{rhs};
    return rectangle{lhs * rhs_.origin, lhs * rhs_.right, lhs * rhs_.up};
}

[[deprecated("Use full_mul() or fast_mul() instead.")]] [[nodiscard]] constexpr rectangle
operator*(matrix2 const& lhs, aarectangle const& rhs) noexcept
{
    auto const rhs_ = rectangle{rhs};
    return rectangle{lhs * rhs_.origin, lhs * rhs_.right, lhs * rhs_.up};
}

[[nodiscard]] constexpr rectangle full_mul(matrix2 const& lhs, aarectangle const& rhs) noexcept
{
    auto const rhs_ = rectangle{rhs};
    return rectangle{lhs * rhs_.origin, lhs * rhs_.right, lhs * rhs_.up};
}

[[nodiscard]] constexpr aarectangle fast_mul(matrix2 const& lhs, aarectangle const& rhs) noexcept
{
    return aarectangle{lhs * get<0>(rhs), lhs * get<3>(rhs)};
}

[[nodiscard]] constexpr rectangle operator*(transformer auto const& lhs, rectangle const& rhs) noexcept
{
    return rectangle{lhs * rhs.origin, lhs * rhs.right, lhs * rhs.up};
}

[[nodiscard]] constexpr quad operator*(transformer auto const& lhs, quad const& rhs) noexcept
{
    return quad{lhs * rhs.p0, lhs * rhs.p1, lhs * rhs.p2, lhs * rhs.p3};
}

[[nodiscard]] constexpr circle operator*(transformer auto const& lhs, circle const& rhs) noexcept
{
    return circle{lhs * midpoint(rhs), lhs * rhs.radius()};
}

[[nodiscard]] constexpr line_segment operator*(transformer auto const& lhs, line_segment const& rhs) noexcept
{
    return line_segment{lhs * rhs.origin(), lhs * rhs.direction()};
}

[[nodiscard]] constexpr corner_radii operator*(matrix3 const& lhs, corner_radii const& rhs) noexcept
{
    return {lhs * get<0>(rhs), lhs * get<1>(rhs), lhs * get<2>(rhs), lhs * get<3>(rhs)};
}

template<typename Lhs, transformer Rhs>
constexpr Lhs& operator*=(Lhs& lhs, Rhs const& rhs) noexcept
    requires requires(Lhs& a, Rhs const& b) { a = b * a; }
{
    return lhs = rhs * lhs;
}
 
[[nodiscard]] constexpr aarectangle fit(aarectangle const& bounds, aarectangle const& rectangle) noexcept
{
    auto const resized_rectangle = aarectangle{
        rectangle.left(),
        rectangle.bottom(),
        std::min(rectangle.width(), bounds.width()),
        std::min(rectangle.height(), bounds.height())};
 
    auto const translate_from_p0 = max(vector2{}, get<0>(bounds) - get<0>(resized_rectangle));
    auto const translate_from_p3 = min(vector2{}, get<3>(bounds) - get<3>(resized_rectangle));
    return translate2{translate_from_p0 + translate_from_p3} * resized_rectangle;
}

[[nodiscard]] constexpr quad scale_from_center(quad const& lhs, scale2 const& rhs) noexcept
{
    auto const top_extra = (lhs.top() * rhs.x() - lhs.top()) * 0.5f;
    auto const bottom_extra = (lhs.bottom() * rhs.x() - lhs.bottom()) * 0.5f;
    auto const left_extra = (lhs.left() * rhs.y() - lhs.left()) * 0.5f;
    auto const right_extra = (lhs.right() * rhs.y() - lhs.right()) * 0.5f;
 
    return {
        lhs.p0 - bottom_extra - left_extra,
        lhs.p1 + bottom_extra - right_extra,
        lhs.p2 - top_extra + left_extra,
        lhs.p3 + top_extra + right_extra};
}
 
}} // namespace hi::v1