algorithm_selector.hpp

// This file is part of the ACTS project.
//
// Copyright (C) 2016 CERN for the benefit of the ACTS project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
#pragma once
 
// Project include(s).
#include "detray/algebra/generic/algorithms/matrix/determinant/hard_coded.hpp"
#include "detray/algebra/generic/algorithms/matrix/determinant/partial_pivot_lud.hpp"
#include "detray/algebra/generic/algorithms/matrix/inverse/hard_coded.hpp"
#include "detray/algebra/generic/algorithms/matrix/inverse/partial_pivot_lud.hpp"
 
namespace detray::algebra::generic {

 
// Default algorithm
template <std::size_t N, typename... Args>
struct determinant_selector {
  using type = matrix::determinant::partial_pivot_lud<Args...>;
};

template <typename... Args>
struct determinant_selector<2, Args...> {
  using type = matrix::determinant::hard_coded<Args...>;
};

template <concepts::square_matrix M>
using determinant_t =
    typename determinant_selector<detray::traits::max_rank<M>, M>::type;

template <std::size_t N, typename... Args>
struct inversion_selector {
  using type = matrix::inverse::partial_pivot_lud<Args...>;
};

template <typename... Args>
struct inversion_selector<2, Args...> {
  using type = matrix::inverse::hard_coded<Args...>;
};

template <concepts::square_matrix M>
using inversion_t =
    typename inversion_selector<detray::traits::max_rank<M>, M>::type;
 
}  // namespace detray::algebra::generic