Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <array>
       7              : 
       8              : #include "DataStructures/Tensor/Metafunctions.hpp"
       9              : #include "DataStructures/Tensor/TypeAliases.hpp"
      10              : #include "Utilities/Gsl.hpp"
      11              : 
      12              : /*!
      13              :  * \brief Orthonormalize a set of vectors using the modified Gram-Schmidt
      14              :  * process.
      15              :  *
      16              :  * The first vector in `basis` is normalized, then each subsequent vector is
      17              :  * orthogonalized against all previous vectors and then normalized.
      18              :  */
      19              : template <typename DataType, typename Index, size_t NumVectors>
      20            1 : void gram_schmidt_orthonormalize(
      21              :     const std::array<
      22              :         gsl::not_null<Tensor<DataType, Symmetry<1>, index_list<Index>>*>,
      23              :         NumVectors>& basis,
      24              :     const Tensor<DataType, Symmetry<1, 1>,
      25              :                  index_list<change_index_up_lo<Index>,
      26              :                             change_index_up_lo<Index>>>& metric);