/**
	This file is a part of the rexy/r0nk/atlas project
	Copyright (C) 2020 rexy712

	This program is free software: you can redistribute it and/or modify
	it under the terms of the GNU General Public License as published by
	the Free Software Foundation, either version 3 of the License, or
	(at your option) any later version.

	This program is distributed in the hope that it will be useful,
	but WITHOUT ANY WARRANTY; without even the implied warranty of
	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.	See the
	GNU General Public License for more details.

	You should have received a copy of the GNU General Public License
	along with this program.	If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef REXY_MAT_HPP
#define REXY_MAT_HPP

#include <cstdlib> //size_t
#include <utility> //integer_sequence
#include <type_traits> //decay_t, is_same, integral_constant
#include "detail/math.hpp"
#include "detail/matrix.hpp"

namespace math{

	template<typename T, size_t C, size_t R>
	class matrix : public detail::matrix_base<T,C,R>
	{
	private:
		using base = detail::matrix_base<T,C,R>;
	public:
		using value_type = typename base::value_type;
		using size_type = typename base::size_type;
		using pointer = typename base::pointer;
		using const_pointer = typename base::const_pointer;
		using reference = typename base::reference;
		using const_reference = typename base::const_reference;
	public:
		using detail::matrix_base<T,C,R>::matrix_base;
		using detail::matrix_base<T,C,R>::operator=;
	};

	template<typename T>
	class matrix<T,3,3> : public detail::matrix_base<T,3,3>
	{
	private:
		using base = detail::matrix_base<T,3,3>;
	public:
		using value_type = typename base::value_type;
		using size_type = typename base::size_type;
		using pointer = typename base::pointer;
		using const_pointer = typename base::const_pointer;
		using reference = typename base::reference;
		using const_reference = typename base::const_reference;
	public:
		using detail::matrix_base<T,3,3>::matrix_base;
		using detail::matrix_base<T,3,3>::operator=;

		static matrix rotation(value_type angle);
		static constexpr matrix rotation(value_type sin, value_type cos);
		static constexpr matrix rotation(value_type angle_x, value_type angle_y, value_type angle_z);
	};

	namespace detail{

		template<typename T>
		struct is_matrix_helper {
			template<typename U, size_t W, size_t H>
			static std::true_type test(matrix<U,W,H>*);
			static std::false_type test(void*);

			static constexpr bool value = std::is_same<std::true_type,decltype(test(static_cast<std::decay_t<T>*>(nullptr)))>::value;
		};

	}

	template<typename... Ms>
	struct is_matrix {
		static constexpr bool value = (detail::is_matrix_helper<Ms>::value && ...);
	};

	namespace detail{

		template<typename M1, typename M2>
		struct are_same_size_matrix {
			using l = std::decay_t<M1>;
			using r = std::decay_t<M2>;
			static constexpr bool value = is_matrix<M1,M2>::value && l::Columns == r::Columns && l::Rows == r::Rows;
		};
		template<typename... Ms>
		using enable_if_matrix = std::enable_if_t<is_matrix<Ms...>::value,int>;

		template<typename M1, typename M2>
		using enable_if_eq_matrix = std::enable_if_t<are_same_size_matrix<M1,M2>::value,int>;


		template<typename T, typename U, size_t W, size_t H>
		constexpr bool operator==(const matrix<T,W,H>& left, const matrix<U,W,H> right);
		template<typename T, typename U, size_t W, size_t H>
		constexpr bool operator!=(const matrix<T,W,H>& left, const matrix<U,W,H> right);

		template<typename T, typename U, size_t R1, size_t C1, size_t C2>
		constexpr auto operator*(const matrix<T,C1,R1>& left, const matrix<U,C2,C1>& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr auto operator*(const matrix<T,C,R>& left, U&& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr auto operator/(const matrix<T,C,R>& left, U&& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr auto operator+(const matrix<T,C,R>& left, const matrix<U,C,R>& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr auto operator-(const matrix<T,C,R>& left, const matrix<U,C,R>& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr auto operator-(const matrix<T,C,R>& left);

		template<typename T, typename U, size_t R1, size_t C1, size_t C2>
		constexpr decltype(auto) operator*=(matrix<T,C1,R1>& left, const matrix<U,C2,C1>& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr decltype(auto) operator*=(matrix<T,C,R>& left, U&& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr decltype(auto) operator/=(matrix<T,C,R>& left, U&& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr decltype(auto) operator+=(matrix<T,C,R>& left, const matrix<U,C,R>& right);
		template<typename T, typename U, size_t C, size_t R>
		constexpr decltype(auto) operator-=(matrix<T,C,R>& left, const matrix<U,C,R>& right);

	}

}

#include "mat.tpp"

#endif