lar::util::details::FastMatrixOperations< T, 4 > Struct Template Reference

Routines for 4x4 matrices. More...

#include <FastMatrixMathHelper.h>

Inheritance diagram for lar::util::details::FastMatrixOperations< T, 4 >:

Public Types
using	Base_t = FastMatrixOperationsBase< T, 4 >
using	Data_t = typename Base_t::Data_t
using	Matrix_t = typename Base_t::Matrix_t
Public Types inherited from lar::util::details::FastMatrixOperationsBase< T, 4 >
using	Data_t = T
using	Matrix_t = std::array< Data_t, Dim *Dim >
using	Vector_t = std::array< Data_t, Dim >

Static Public Member Functions
static Data_t	Determinant (Matrix_t const &mat)
	Computes the determinant of a matrix. More...
static Matrix_t	InvertMatrix (Matrix_t const &mat, Data_t det)
	Computes the determinant of a matrix, using the provided determinant. More...
static Matrix_t	InvertSymmetricMatrix (Matrix_t const &mat, Data_t det)
static Matrix_t	InvertMatrix (Matrix_t const &mat)
	Computes the determinant of a matrix. More...
static Matrix_t	InvertSymmetricMatrix (Matrix_t const &mat)
	Computes the determinant of a matrix. More...
Static Public Member Functions inherited from lar::util::details::FastMatrixOperationsBase< T, 4 >
static Vector_t	MatrixVectorProduct (Matrix_t const &mat, Vector_t const &vec)
	Returns the product of a square matrix times a column vector. More...
static constexpr Data_t	sqr (Data_t v)

Static Public Attributes
static constexpr unsigned int	Dim = Base_t::Dim
Static Public Attributes inherited from lar::util::details::FastMatrixOperationsBase< T, 4 >
static constexpr unsigned int	Dim
	matrix dimensions More...

Detailed Description

template<typename T>
struct lar::util::details::FastMatrixOperations< T, 4 >

Routines for 4x4 matrices.

Definition at line 321 of file FastMatrixMathHelper.h.

Member Typedef Documentation

template<typename T >

using lar::util::details::FastMatrixOperations< T, 4 >::Base_t = FastMatrixOperationsBase<T, 4>

Definition at line 324 of file FastMatrixMathHelper.h.

template<typename T >

using lar::util::details::FastMatrixOperations< T, 4 >::Data_t = typename Base_t::Data_t

Definition at line 326 of file FastMatrixMathHelper.h.

template<typename T >

using lar::util::details::FastMatrixOperations< T, 4 >::Matrix_t = typename Base_t::Matrix_t

Definition at line 327 of file FastMatrixMathHelper.h.

Member Function Documentation

template<typename T >

static Data_t lar::util::details::FastMatrixOperations< T, 4 >::Determinant ( Matrix_t const &  mat )

inlinestatic

Computes the determinant of a matrix.

Definition at line 330 of file FastMatrixMathHelper.h.

          {
            return DeterminantHelper<T, Dim, 0, 1, 2, 3, 0, 1, 2, 3>::compute
              (mat.data());
          }

template<typename T >

auto lar::util::details::FastMatrixOperations< T, 4 >::InvertMatrix ( Matrix_t const &  mat, Data_t  det  )

static

Computes the determinant of a matrix, using the provided determinant.

Definition at line 455 of file FastMatrixMathHelper.h.

{
  //
  // Basically using Cramer's rule;
  // each element [r,c] gets assigned the determinant of the submatrix
  //   after removing c from the rows and r from the columns
  //   (effectively assigning the transpose of the minor matrix)
  //   with the usual sign -1^(r+c)
  //
  //
  Data_t const* data = mat.data();
  Matrix_t Inverse;
  Inverse[0 * Dim + 0] =  DeterminantHelper<T, 4, 1, 2, 3, 1, 2, 3>::compute(data)/det;
  Inverse[0 * Dim + 1] = -DeterminantHelper<T, 4, 0, 2, 3, 1, 2, 3>::compute(data)/det;
  Inverse[0 * Dim + 2] =  DeterminantHelper<T, 4, 0, 1, 3, 1, 2, 3>::compute(data)/det;
  Inverse[0 * Dim + 3] = -DeterminantHelper<T, 4, 0, 1, 2, 1, 2, 3>::compute(data)/det;
  Inverse[1 * Dim + 0] = -DeterminantHelper<T, 4, 1, 2, 3, 0, 2, 3>::compute(data)/det;
  Inverse[1 * Dim + 1] =  DeterminantHelper<T, 4, 0, 2, 3, 0, 2, 3>::compute(data)/det;
  Inverse[1 * Dim + 2] = -DeterminantHelper<T, 4, 0, 1, 3, 0, 2, 3>::compute(data)/det;
  Inverse[1 * Dim + 3] =  DeterminantHelper<T, 4, 0, 1, 2, 0, 2, 3>::compute(data)/det;
  Inverse[2 * Dim + 0] =  DeterminantHelper<T, 4, 1, 2, 3, 0, 1, 3>::compute(data)/det;
  Inverse[2 * Dim + 1] = -DeterminantHelper<T, 4, 0, 2, 3, 0, 1, 3>::compute(data)/det;
  Inverse[2 * Dim + 2] =  DeterminantHelper<T, 4, 0, 1, 3, 0, 1, 3>::compute(data)/det;
  Inverse[2 * Dim + 3] = -DeterminantHelper<T, 4, 0, 1, 2, 0, 1, 3>::compute(data)/det;
  Inverse[3 * Dim + 0] = -DeterminantHelper<T, 4, 1, 2, 3, 0, 1, 2>::compute(data)/det;
  Inverse[3 * Dim + 1] =  DeterminantHelper<T, 4, 0, 2, 3, 0, 1, 2>::compute(data)/det;
  Inverse[3 * Dim + 2] = -DeterminantHelper<T, 4, 0, 1, 3, 0, 1, 2>::compute(data)/det;
  Inverse[3 * Dim + 3] =  DeterminantHelper<T, 4, 0, 1, 2, 0, 1, 2>::compute(data)/det;
  return Inverse;
} // FastMatrixOperations<T, 4>::InvertMatrix()

template<typename T >

static Matrix_t lar::util::details::FastMatrixOperations< T, 4 >::InvertMatrix ( Matrix_t const &  mat )

inlinestatic

Computes the determinant of a matrix.

Definition at line 344 of file FastMatrixMathHelper.h.

          { return InvertMatrix(mat, Determinant(mat)); }

template<typename T >

auto lar::util::details::FastMatrixOperations< T, 4 >::InvertSymmetricMatrix ( Matrix_t const &  mat, Data_t  det  )

static

Computes the determinant of a symmatric matrix, using the provided determinant

Definition at line 489 of file FastMatrixMathHelper.h.

{
  //
  // Same algorithm as InvertMatrix(), but use the fact that the result is
  // also symmetric
  //
  Data_t const* data = mat.data();
  Matrix_t Inverse;
  Inverse[0 * Dim + 0] =  DeterminantHelper<T, 4, 1, 2, 3, 1, 2, 3>::compute(data)/det;
  Inverse[0 * Dim + 1] =
  Inverse[1 * Dim + 0] = -DeterminantHelper<T, 4, 1, 2, 3, 0, 2, 3>::compute(data)/det;
  Inverse[0 * Dim + 2] =
  Inverse[2 * Dim + 0] =  DeterminantHelper<T, 4, 1, 2, 3, 0, 1, 3>::compute(data)/det;
  Inverse[0 * Dim + 3] =
  Inverse[3 * Dim + 0] = -DeterminantHelper<T, 4, 1, 2, 3, 0, 1, 2>::compute(data)/det;
  Inverse[1 * Dim + 1] =  DeterminantHelper<T, 4, 0, 2, 3, 0, 2, 3>::compute(data)/det;
  Inverse[1 * Dim + 2] =
  Inverse[2 * Dim + 1] = -DeterminantHelper<T, 4, 0, 2, 3, 3, 0, 1>::compute(data)/det;
  Inverse[1 * Dim + 3] =
  Inverse[3 * Dim + 1] =  DeterminantHelper<T, 4, 0, 2, 3, 0, 1, 2>::compute(data)/det;
  Inverse[2 * Dim + 2] =  DeterminantHelper<T, 4, 0, 1, 3, 0, 1, 3>::compute(data)/det;
  Inverse[2 * Dim + 3] =
  Inverse[3 * Dim + 2] = -DeterminantHelper<T, 4, 0, 1, 3, 0, 1, 2>::compute(data)/det;
  Inverse[3 * Dim + 3] =  DeterminantHelper<T, 4, 0, 1, 2, 0, 1, 2>::compute(data)/det;
  return Inverse;
} // FastMatrixOperations<T, 4>::InvertSymmetricMatrix()

template<typename T >

static Matrix_t lar::util::details::FastMatrixOperations< T, 4 >::InvertSymmetricMatrix ( Matrix_t const &  mat )

inlinestatic

Computes the determinant of a matrix.

Definition at line 348 of file FastMatrixMathHelper.h.

          { return InvertSymmetricMatrix(mat, Determinant(mat)); }

Member Data Documentation

template<typename T >

constexpr unsigned int lar::util::details::FastMatrixOperations< T, 4 >::Dim = Base_t::Dim

static

Definition at line 325 of file FastMatrixMathHelper.h.

---

The documentation for this struct was generated from the following file:

• lardata/lardata/Utilities/FastMatrixMathHelper.h