mirror of https://github.com/mosra/magnum.git
You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
626 lines
29 KiB
626 lines
29 KiB
#ifndef Magnum_Math_RectangularMatrix_h |
|
#define Magnum_Math_RectangularMatrix_h |
|
/* |
|
Copyright © 2010, 2011, 2012 Vladimír Vondruš <mosra@centrum.cz> |
|
|
|
This file is part of Magnum. |
|
|
|
Magnum is free software: you can redistribute it and/or modify |
|
it under the terms of the GNU Lesser General Public License version 3 |
|
only, as published by the Free Software Foundation. |
|
|
|
Magnum 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 Lesser General Public License version 3 for more details. |
|
*/ |
|
|
|
/** @file |
|
* @brief Class Magnum::Math::RectangularMatrix |
|
*/ |
|
|
|
#include <cmath> |
|
#include <limits> |
|
#include <Utility/Debug.h> |
|
#include <Utility/ConfigurationValue.h> |
|
|
|
#include "MathTypeTraits.h" |
|
|
|
#include "magnumVisibility.h" |
|
|
|
namespace Magnum { namespace Math { |
|
|
|
/** @todo Properly test all constexpr */ |
|
|
|
template<std::size_t, std::size_t, class> class RectangularMatrix; |
|
|
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
namespace Implementation { |
|
template<std::size_t ...> struct Sequence {}; |
|
|
|
/* E.g. GenerateSequence<3>::Type is Sequence<0, 1, 2> */ |
|
template<std::size_t N, std::size_t ...sequence> struct GenerateSequence: |
|
GenerateSequence<N-1, N-1, sequence...> {}; |
|
|
|
template<std::size_t ...sequence> struct GenerateSequence<0, sequence...> { |
|
typedef Sequence<sequence...> Type; |
|
}; |
|
|
|
/* Implementation for RectangularMatrix<cols, rows, T>::from(const RectangularMatrix<cols, rows, U>&) */ |
|
template<class T, class U, std::size_t c, std::size_t ...sequence> inline constexpr Math::RectangularMatrix<c, sizeof...(sequence)/c, T> rectangularMatrixFrom(Sequence<sequence...>, const Math::RectangularMatrix<c, sizeof...(sequence)/c, U>& matrix) { |
|
return {T(matrix.data()[sequence])...}; |
|
} |
|
} |
|
#endif |
|
|
|
template<std::size_t size, class T> class Vector; |
|
|
|
/** |
|
@brief Rectangular matrix |
|
@tparam cols Column count |
|
@tparam rows Row count |
|
@tparam T Data type |
|
|
|
See @ref matrix-vector for brief introduction. See also Matrix (square) and |
|
Vector. |
|
|
|
The data are stored in column-major order, to reflect that, all indices in |
|
math formulas are in reverse order (i.e. @f$ \boldsymbol A_{ji} @f$ instead |
|
of @f$ \boldsymbol A_{ij} @f$). |
|
*/ |
|
template<std::size_t cols, std::size_t rows, class T> class RectangularMatrix { |
|
static_assert(cols != 0 && rows != 0, "Matrix cannot have zero elements"); |
|
|
|
public: |
|
typedef T Type; /**< @brief Data type */ |
|
const static std::size_t Cols = cols; /**< @brief %Matrix column count */ |
|
const static std::size_t Rows = rows; /**< @brief %Matrix row count */ |
|
|
|
/** |
|
* @brief %Matrix from array |
|
* @return Reference to the data as if it was Matrix, thus doesn't |
|
* perform any copying. |
|
* |
|
* @attention Use with caution, the function doesn't check whether the |
|
* array is long enough. |
|
*/ |
|
inline constexpr static RectangularMatrix<cols, rows, T>& from(T* data) { |
|
return *reinterpret_cast<RectangularMatrix<cols, rows, T>*>(data); |
|
} |
|
/** @overload */ |
|
inline constexpr static const RectangularMatrix<cols, rows, T>& from(const T* data) { |
|
return *reinterpret_cast<const RectangularMatrix<cols, rows, T>*>(data); |
|
} |
|
|
|
/** |
|
* @brief %Matrix from column vectors |
|
* @param first First column vector |
|
* @param next Next column vectors |
|
* |
|
* @todo Creating matrix from arbitrary combination of matrices with n rows |
|
*/ |
|
template<class ...U> inline constexpr static RectangularMatrix<cols, rows, T> from(const Vector<rows, T>& first, const U&... next) { |
|
static_assert(sizeof...(next)+1 == cols, "Improper number of arguments passed to Matrix from Vector constructor"); |
|
return from(typename Implementation::GenerateSequence<rows>::Type(), first, next...); |
|
} |
|
|
|
/** |
|
* @brief %Matrix from another of different type |
|
* |
|
* Performs only default casting on the values, no rounding or |
|
* anything else. Example usage: |
|
* @code |
|
* RectangularMatrix<4, 1, float> floatingPoint(1.3f, 2.7f, -15.0f, 7.0f); |
|
* RectangularMatrix<4, 1, std::int8_t> integral(RectangularMatrix<4, 1, std::int8_t>::from(floatingPoint)); |
|
* // integral == {1, 2, -15, 7} |
|
* @endcode |
|
*/ |
|
template<class U> inline constexpr static RectangularMatrix<cols, rows, T> from(const RectangularMatrix<cols, rows, U>& other) { |
|
return Implementation::rectangularMatrixFrom<T, U>(typename Implementation::GenerateSequence<cols*rows>::Type(), other); |
|
} |
|
|
|
/** @brief Zero-filled matrix constructor */ |
|
inline constexpr /*implicit*/ RectangularMatrix(): _data() {} |
|
|
|
/** |
|
* @brief Initializer-list constructor |
|
* @param first First value |
|
* @param next Next values |
|
* |
|
* Note that the values are in column-major order. |
|
* @todoc Remove workaround when Doxygen supports uniform initialization |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<class ...U> inline constexpr /*implicit*/ RectangularMatrix(T first, U... next): _data{first, next...} { |
|
static_assert(sizeof...(next)+1 == cols*rows, "Improper number of arguments passed to RectangularMatrix constructor"); |
|
} |
|
#else |
|
template<class ...U> inline constexpr RectangularMatrix(T first, U... next); |
|
#endif |
|
|
|
/** @brief Copy constructor */ |
|
inline constexpr RectangularMatrix(const RectangularMatrix<cols, rows, T>&) = default; |
|
|
|
/** @brief Assignment operator */ |
|
inline RectangularMatrix<cols, rows, T>& operator=(const RectangularMatrix<cols, rows, T>&) = default; |
|
|
|
/** |
|
* @brief Raw data |
|
* @return One-dimensional array of `size*size` length in column-major |
|
* order. |
|
*/ |
|
inline T* data() { return _data; } |
|
inline constexpr const T* data() const { return _data; } /**< @overload */ |
|
|
|
/** |
|
* @brief %Matrix column |
|
* |
|
* For accessing individual elements prefer to use operator(), as it |
|
* is guaranteed to not involve unnecessary conversions. |
|
*/ |
|
inline Vector<rows, T>& operator[](std::size_t col) { |
|
return Vector<rows, T>::from(_data+col*rows); |
|
} |
|
/** @overload */ |
|
inline constexpr const Vector<rows, T>& operator[](std::size_t col) const { |
|
return Vector<rows, T>::from(_data+col*rows); |
|
} |
|
|
|
/** |
|
* @brief Element on given position |
|
* |
|
* Prefer this instead of using `[][]`. |
|
* @see operator[] |
|
*/ |
|
inline T& operator()(std::size_t col, std::size_t row) { |
|
return _data[col*rows+row]; |
|
} |
|
/** @overload */ |
|
inline constexpr const T& operator()(std::size_t col, std::size_t row) const { |
|
return _data[col*rows+row]; |
|
} |
|
|
|
/** |
|
* @brief Equality operator |
|
* |
|
* @see Vector::operator<(), Vector::operator<=(), Vector::operator>=(), |
|
* Vector::operator>() |
|
*/ |
|
inline bool operator==(const RectangularMatrix<cols, rows, T>& other) const { |
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
if(!MathTypeTraits<T>::equals(_data[i], other._data[i])) return false; |
|
|
|
return true; |
|
} |
|
|
|
/** |
|
* @brief Non-equality operator |
|
* |
|
* @see Vector::operator<(), Vector::operator<=(), Vector::operator>=(), |
|
* Vector::operator>() |
|
*/ |
|
inline constexpr bool operator!=(const RectangularMatrix<cols, rows, T>& other) const { |
|
return !operator==(other); |
|
} |
|
|
|
/** |
|
* @brief Add and assign matrix |
|
* |
|
* The computation is done in-place. @f[ |
|
* \boldsymbol A_{ji} = \boldsymbol A_{ji} + \boldsymbol B_{ji} |
|
* @f] |
|
*/ |
|
RectangularMatrix<cols, rows, T>& operator+=(const RectangularMatrix<cols, rows, T>& other) { |
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
_data[i] += other._data[i]; |
|
|
|
return *this; |
|
} |
|
|
|
/** |
|
* @brief Add matrix |
|
* |
|
* @see operator+=() |
|
*/ |
|
inline RectangularMatrix<cols, rows, T> operator+(const RectangularMatrix<cols, rows, T>& other) const { |
|
return RectangularMatrix<cols, rows, T>(*this)+=other; |
|
} |
|
|
|
/** |
|
* @brief Negated matrix |
|
* |
|
* The computation is done in-place. @f[ |
|
* \boldsymbol A_{ji} = -\boldsymbol A_{ji} |
|
* @f] |
|
*/ |
|
RectangularMatrix<cols, rows, T> operator-() const { |
|
RectangularMatrix<cols, rows, T> out; |
|
|
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
out._data[i] = -_data[i]; |
|
|
|
return out; |
|
} |
|
|
|
/** |
|
* @brief Subtract and assign matrix |
|
* |
|
* The computation is done in-place. @f[ |
|
* \boldsymbol A_{ji} = \boldsymbol A_{ji} - \boldsymbol B_{ji} |
|
* @f] |
|
*/ |
|
RectangularMatrix<cols, rows, T>& operator-=(const RectangularMatrix<cols, rows, T>& other) { |
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
_data[i] -= other._data[i]; |
|
|
|
return *this; |
|
} |
|
|
|
/** |
|
* @brief Subtract matrix |
|
* |
|
* @see operator-=() |
|
*/ |
|
inline RectangularMatrix<cols, rows, T> operator-(const RectangularMatrix<cols, rows, T>& other) const { |
|
return RectangularMatrix<cols, rows, T>(*this)-=other; |
|
} |
|
|
|
/** |
|
* @brief Multiply matrix with number |
|
* |
|
* @see operator*=(U), operator*(U, const RectangularMatrix<cols, rows, T>&) |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>>::type operator*(U number) const { |
|
#else |
|
template<class U> inline RectangularMatrix<cols, rows, T> operator*(U number) const { |
|
#endif |
|
return RectangularMatrix<cols, rows, T>(*this)*=number; |
|
} |
|
|
|
/** |
|
* @brief Multiply matrix with number and assign |
|
* |
|
* The computation is done in-place. @f[ |
|
* \boldsymbol A_{ji} = a \boldsymbol A_{ji} |
|
* @f] |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>&>::type operator*=(U number) { |
|
#else |
|
template<class U> RectangularMatrix<cols, rows, T>& operator*=(U number) { |
|
#endif |
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
_data[i] *= number; |
|
|
|
return *this; |
|
} |
|
|
|
/** |
|
* @brief Divide matrix with number |
|
* |
|
* @see operator/=(), operator/(U, const RectangularMatrix<cols, rows, T>&) |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>>::type operator/(U number) const { |
|
#else |
|
template<class U> inline RectangularMatrix<cols, rows, T> operator/(U number) const { |
|
#endif |
|
return RectangularMatrix<cols, rows, T>(*this)/=number; |
|
} |
|
|
|
/** |
|
* @brief Divide matrix with number and assign |
|
* |
|
* The computation is done in-place. @f[ |
|
* \boldsymbol A_{ji} = \frac{\boldsymbol A_{ji}} a |
|
* @f] |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>&>::type operator/=(U number) { |
|
#else |
|
template<class U> RectangularMatrix<cols, rows, T>& operator/=(U number) { |
|
#endif |
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
_data[i] /= number; |
|
|
|
return *this; |
|
} |
|
|
|
/** |
|
* @brief Multiply matrix |
|
* |
|
* @f[ |
|
* (\boldsymbol {AB})_{ji} = \sum_{k=0}^{m-1} \boldsymbol A_{ki} \boldsymbol B_{jk} |
|
* @f] |
|
*/ |
|
template<std::size_t size> RectangularMatrix<size, rows, T> operator*(const RectangularMatrix<size, cols, T>& other) const { |
|
RectangularMatrix<size, rows, T> out; |
|
|
|
for(std::size_t col = 0; col != size; ++col) |
|
for(std::size_t row = 0; row != rows; ++row) |
|
for(std::size_t pos = 0; pos != cols; ++pos) |
|
out(col, row) += (*this)(pos, row)*other(col, pos); |
|
|
|
return out; |
|
} |
|
|
|
/** |
|
* @brief Multiply vector |
|
* |
|
* Internally the same as multiplying with one-column matrix, but |
|
* returns vector. @f[ |
|
* (\boldsymbol {Aa})_i = \sum_{k=0}^{m-1} \boldsymbol A_{ki} \boldsymbol a_k |
|
* @f] |
|
*/ |
|
Vector<rows, T> operator*(const Vector<rows, T>& other) const { |
|
return operator*(static_cast<RectangularMatrix<1, rows, T>>(other)); |
|
} |
|
|
|
/** @brief Transposed matrix */ |
|
RectangularMatrix<rows, cols, T> transposed() const { |
|
RectangularMatrix<rows, cols, T> out; |
|
|
|
for(std::size_t col = 0; col != cols; ++col) |
|
for(std::size_t row = 0; row != rows; ++row) |
|
out(row, col) = (*this)(col, row); |
|
|
|
return out; |
|
} |
|
|
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
protected: |
|
T _data[rows*cols]; |
|
#endif |
|
|
|
private: |
|
template<std::size_t ...sequence, class ...U> inline constexpr static RectangularMatrix<cols, rows, T> from(Implementation::Sequence<sequence...> s, const Vector<rows, T>& first, U... next) { |
|
return from(s, next..., first[sequence]...); |
|
} |
|
template<std::size_t ...sequence, class ...U> inline constexpr static RectangularMatrix<cols, rows, T> from(Implementation::Sequence<sequence...>, T first, U... next) { |
|
return RectangularMatrix<cols, rows, T>(first, next...); |
|
} |
|
}; |
|
|
|
/** @relates RectangularMatrix |
|
@brief Multiply number with matrix |
|
|
|
Same as RectangularMatrix::operator*(U) const. |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<std::size_t cols, std::size_t rows, class T, class U> inline typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>>::type operator*(U number, const RectangularMatrix<cols, rows, T>& matrix) { |
|
#else |
|
template<std::size_t cols, std::size_t rows, class T, class U> inline RectangularMatrix<cols, rows, T> operator*(U number, const RectangularMatrix<cols, rows, T>& matrix) { |
|
#endif |
|
return matrix*number; |
|
} |
|
|
|
/** @relates RectangularMatrix |
|
@brief Divide matrix with number and invert |
|
|
|
@f[ |
|
\boldsymbol B_{ji} = \frac a {\boldsymbol A_{ji}} |
|
@f] |
|
@see RectangularMatrix::operator/() |
|
*/ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
template<std::size_t cols, std::size_t rows, class T, class U> typename std::enable_if<std::is_arithmetic<U>::value, RectangularMatrix<cols, rows, T>>::type operator/(U number, const RectangularMatrix<cols, rows, T>& matrix) { |
|
#else |
|
template<std::size_t cols, std::size_t rows, class T, class U> RectangularMatrix<cols, rows, T> operator/(U number, const RectangularMatrix<cols, rows, T>& matrix) { |
|
#endif |
|
RectangularMatrix<cols, rows, T> out; |
|
|
|
for(std::size_t i = 0; i != cols*rows; ++i) |
|
out.data()[i] = number/matrix.data()[i]; |
|
|
|
return out; |
|
} |
|
|
|
/** @debugoperator{Magnum::Math::RectangularMatrix} */ |
|
template<std::size_t cols, std::size_t rows, class T> Corrade::Utility::Debug operator<<(Corrade::Utility::Debug debug, const Magnum::Math::RectangularMatrix<cols, rows, T>& value) { |
|
debug << "Matrix("; |
|
debug.setFlag(Corrade::Utility::Debug::SpaceAfterEachValue, false); |
|
for(std::size_t row = 0; row != rows; ++row) { |
|
if(row != 0) debug << ",\n "; |
|
for(std::size_t col = 0; col != cols; ++col) { |
|
if(col != 0) debug << ", "; |
|
debug << typename MathTypeTraits<T>::NumericType(value[col][row]); |
|
} |
|
} |
|
debug << ')'; |
|
debug.setFlag(Corrade::Utility::Debug::SpaceAfterEachValue, true); |
|
return debug; |
|
} |
|
|
|
/* Explicit instantiation for types used in OpenGL */ |
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
/* Vectors */ |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 2, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 3, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 4, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 2, int>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 3, int>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 4, int>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 2, unsigned int>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 3, unsigned int>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 4, unsigned int>&); |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 2, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 3, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<1, 4, double>&); |
|
#endif |
|
|
|
/* Square matrices */ |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 2, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 3, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 4, float>&); |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 2, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 3, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 4, double>&); |
|
#endif |
|
|
|
/* Rectangular matrices */ |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 3, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 2, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 4, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 2, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 4, float>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 3, float>&); |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 3, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 2, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<2, 4, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 2, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<3, 4, double>&); |
|
extern template Corrade::Utility::Debug MAGNUM_EXPORT operator<<(Corrade::Utility::Debug, const RectangularMatrix<4, 3, double>&); |
|
#endif |
|
#endif |
|
|
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
#define MAGNUM_RECTANGULARMATRIX_SUBCLASS_IMPLEMENTATION(cols, rows, ...) \ |
|
inline constexpr static __VA_ARGS__& from(T* data) { \ |
|
return *reinterpret_cast<__VA_ARGS__*>(data); \ |
|
} \ |
|
inline constexpr static const __VA_ARGS__& from(const T* data) { \ |
|
return *reinterpret_cast<const __VA_ARGS__*>(data); \ |
|
} \ |
|
template<class ...U> inline constexpr static __VA_ARGS__ from(const Math::Vector<rows, T>& first, const U&... next) { \ |
|
return Math::RectangularMatrix<cols, rows, T>::from(first, next...); \ |
|
} \ |
|
template<class U> inline constexpr static RectangularMatrix<cols, rows, T> from(const Math::RectangularMatrix<cols, rows, U>& other) { \ |
|
return Math::RectangularMatrix<cols, rows, T>::from(other); \ |
|
} \ |
|
\ |
|
inline __VA_ARGS__& operator=(const Math::RectangularMatrix<cols, rows, T>& other) { \ |
|
Math::RectangularMatrix<cols, rows, T>::operator=(other); \ |
|
return *this; \ |
|
} |
|
|
|
#define MAGNUM_RECTANGULARMATRIX_SUBCLASS_OPERATOR_IMPLEMENTATION(cols, rows, ...) \ |
|
inline __VA_ARGS__ operator-() const { \ |
|
return Math::RectangularMatrix<cols, rows, T>::operator-(); \ |
|
} \ |
|
inline __VA_ARGS__ operator+(const Math::RectangularMatrix<cols, rows, T>& other) const { \ |
|
return Math::RectangularMatrix<cols, rows, T>::operator+(other); \ |
|
} \ |
|
inline __VA_ARGS__& operator+=(const Math::RectangularMatrix<cols, rows, T>& other) { \ |
|
Math::RectangularMatrix<cols, rows, T>::operator+=(other); \ |
|
return *this; \ |
|
} \ |
|
inline __VA_ARGS__ operator-(const Math::RectangularMatrix<cols, rows, T>& other) const { \ |
|
return Math::RectangularMatrix<cols, rows, T>::operator-(other); \ |
|
} \ |
|
inline __VA_ARGS__& operator-=(const Math::RectangularMatrix<cols, rows, T>& other) { \ |
|
Math::RectangularMatrix<cols, rows, T>::operator-=(other); \ |
|
return *this; \ |
|
} \ |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, __VA_ARGS__>::type operator*(U number) const { \ |
|
return Math::RectangularMatrix<cols, rows, T>::operator*(number); \ |
|
} \ |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, __VA_ARGS__&>::type operator*=(U number) { \ |
|
Math::RectangularMatrix<cols, rows, T>::operator*=(number); \ |
|
return *this; \ |
|
} \ |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, __VA_ARGS__>::type operator/(U number) const { \ |
|
return Math::RectangularMatrix<cols, rows, T>::operator/(number); \ |
|
} \ |
|
template<class U> inline typename std::enable_if<std::is_arithmetic<U>::value, __VA_ARGS__&>::type operator/=(U number) { \ |
|
Math::RectangularMatrix<cols, rows, T>::operator/=(number); \ |
|
return *this; \ |
|
} |
|
#endif |
|
|
|
}} |
|
|
|
namespace Corrade { namespace Utility { |
|
|
|
/** @configurationvalue{Magnum::Math::RectangularMatrix} */ |
|
template<std::size_t cols, std::size_t rows, class T> struct ConfigurationValue<Magnum::Math::RectangularMatrix<cols, rows, T>> { |
|
ConfigurationValue() = delete; |
|
|
|
/** @brief Writes elements separated with spaces */ |
|
static std::string toString(const Magnum::Math::RectangularMatrix<cols, rows, T>& value, ConfigurationValueFlags flags) { |
|
std::string output; |
|
|
|
for(std::size_t row = 0; row != rows; ++row) { |
|
for(std::size_t col = 0; col != cols; ++col) { |
|
if(!output.empty()) output += ' '; |
|
output += ConfigurationValue<T>::toString(value(col, row), flags); |
|
} |
|
} |
|
|
|
return output; |
|
} |
|
|
|
/** @brief Reads elements separated with whitespace */ |
|
static Magnum::Math::RectangularMatrix<cols, rows, T> fromString(const std::string& stringValue, ConfigurationValueFlags flags) { |
|
Magnum::Math::RectangularMatrix<cols, rows, T> result; |
|
|
|
std::size_t oldpos = 0, pos = std::string::npos, i = 0; |
|
do { |
|
pos = stringValue.find(' ', oldpos); |
|
std::string part = stringValue.substr(oldpos, pos-oldpos); |
|
|
|
if(!part.empty()) { |
|
result(i%cols, i/cols) = ConfigurationValue<T>::fromString(part, flags); |
|
++i; |
|
} |
|
|
|
oldpos = pos+1; |
|
} while(pos != std::string::npos); |
|
|
|
return result; |
|
} |
|
}; |
|
|
|
#ifndef DOXYGEN_GENERATING_OUTPUT |
|
/* Vectors */ |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 2, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 3, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 4, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 2, int>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 3, int>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 4, int>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 2, unsigned int>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 3, unsigned int>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 4, unsigned int>>; |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 2, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 3, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<1, 4, double>>; |
|
#endif |
|
|
|
/* Square matrices */ |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 2, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 3, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 4, float>>; |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 2, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 3, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 4, double>>; |
|
#endif |
|
|
|
/* Rectangular matrices */ |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 3, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 2, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 4, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 2, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 4, float>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 3, float>>; |
|
#ifndef MAGNUM_TARGET_GLES |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 3, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 2, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<2, 4, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 2, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<3, 4, double>>; |
|
extern template struct MAGNUM_EXPORT ConfigurationValue<Magnum::Math::RectangularMatrix<4, 3, double>>; |
|
#endif |
|
#endif |
|
|
|
}} |
|
|
|
/* Include also Vector, so the definition is complete */ |
|
#include "Vector.h" |
|
|
|
#endif
|
|
|