// ---------------------------------------------------------------------
//
// Copyright (C) 2022 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------


for (deal_II_dimension : DIMENSIONS;
     deal_II_scalar_vectorized : REAL_SCALARS_VECTORIZED)
  {
    template void select_vmult<1>(
      deal_II_scalar_vectorized * dst,
      const deal_II_scalar_vectorized *         src,
      AlignedVector<deal_II_scalar_vectorized> &tmp,
      const unsigned int                        n_rows,
      const std::array<const deal_II_scalar_vectorized *, deal_II_dimension>
        &mass_matrix,
      const std::array<const deal_II_scalar_vectorized *, deal_II_dimension>
        &derivative_matrix);

    template void select_apply_inverse<1>(
      deal_II_scalar_vectorized * dst,
      const deal_II_scalar_vectorized *src,
      const unsigned int               n_rows,
      const std::array<const deal_II_scalar_vectorized *, deal_II_dimension>
        &eigenvectors,
      const std::array<const deal_II_scalar_vectorized *, deal_II_dimension>
        &                              eigenvalues,
      const deal_II_scalar_vectorized *inverted_eigenvalues);
  }

for (deal_II_dimension : DIMENSIONS; deal_II_scalar : REAL_SCALARS)
  {
    template void select_vmult<1>(
      deal_II_scalar * dst,
      const deal_II_scalar *                                       src,
      AlignedVector<deal_II_scalar> &                              tmp,
      const unsigned int                                           n_rows,
      const std::array<const deal_II_scalar *, deal_II_dimension> &mass_matrix,
      const std::array<const deal_II_scalar *, deal_II_dimension>
        &derivative_matrix);

    template void select_apply_inverse<1>(
      deal_II_scalar * dst,
      const deal_II_scalar *                                       src,
      const unsigned int                                           n_rows,
      const std::array<const deal_II_scalar *, deal_II_dimension> &eigenvectors,
      const std::array<const deal_II_scalar *, deal_II_dimension> &eigenvalues,
      const deal_II_scalar *inverted_eigenvalues);
  }
