make_constant_diagonal.hpp

Go to the documentation of this file.

/* This file is part of OpenKalman, a header-only C++ library for
 * Kalman filters and other recursive filters.
 *
 * Copyright (c) 2025 Christopher Lee Ogden <ogden@gatech.edu>
 *
 * 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/.
 */

#ifndef OPENKALMAN_MAKE_CONSTANT_DIAGONAL_HPP
#define OPENKALMAN_MAKE_CONSTANT_DIAGONAL_HPP

#include "patterns/patterns.hpp"
#include "linear-algebra/concepts/constant_diagonal_object.hpp"
#include "make_constant.hpp"
#include "to_diagonal.hpp"

namespace OpenKalman
{
#ifdef __cpp_concepts
  template<values::value C, patterns::pattern_collection P>
  constexpr constant_diagonal_object auto
#else
  template<typename C, typename P, std::enable_if_t<
    values::value<C> and patterns::pattern_collection<P>, int> = 0>
  constexpr auto
#endif
  make_constant_diagonal(C c, P&& p)
  {
    auto nested_p = patterns::views::diagonal_of(patterns::to_extents(p));
    return to_diagonal(make_constant(std::move(c), std::move(nested_p)), std::forward<P>(p));
  }


#ifdef __cpp_concepts
  template<values::value C, patterns::pattern...Ps>
  constexpr constant_diagonal_object auto
#else
  template<typename C, typename...Ps, std::enable_if_t<values::value<C> and (... and patterns::pattern<Ps>), int> = 0>
  constexpr auto
#endif
  make_constant_diagonal(C c, Ps&&...ps)
  {
    return make_constant_diagonal(std::move(c), std::tuple{std::forward<Ps>(ps)...});
  }


#ifdef __cpp_concepts
  template<patterns::pattern_collection P, values::value C> requires
    std::default_initializable<P> and
    values::fixed<collections::size_of<P>>
  constexpr constant_diagonal_object auto
#else
  template<typename C, typename P, std::enable_if_t<
    patterns::pattern_collection<P> and
    values::value<C> and
    values::fixed<collections::size_of<P>>, int> = 0>
  constexpr auto
#endif
  make_constant_diagonal(C c)
  {
    return make_constant_diagonal(std::move(c), P{});
  }


}

#endif