symmetric_matrix_space.hpp

Go to the documentation of this file.

/*
 * Copyright © 2025-2026 Akreon LLC
 * SPDX-License-Identifier: Apache-2.0
 */
 
#pragma once
 
#include "coordinate_map.hpp"
#include "akreon/linalg/operator/symmetric_matrix.hpp"
 
namespace linalg {
 
template <Real T, std::size_t Dim>
struct SymmetricMatrixSpace {};
 
template <Real T, std::size_t Dim>
struct space_traits<SymmetricMatrixSpace<T, Dim>> {
  using scalar_type = T;
  using element_type = SymmetricMatrix<Matrix<T, Dim, Dim>>;
  using coordinate_type = Vector<T, Dim*(Dim + 1) / 2>;
  static constexpr std::size_t dim = Dim * (Dim + 1) / 2;
  static constexpr std::size_t DIM = dim;
};
 
template <Real T, std::size_t Dim>
struct coordinate_map<SymmetricMatrixSpace<T, Dim>> {
  using Space = SymmetricMatrixSpace<T, Dim>;
 
  static constexpr auto to_coordinates(const space_element_t<Space>& element) {
    return element.svec;
  }
 
  static constexpr auto from_coordinates(
      const space_coordinate_t<Space>& coordinate) {
    return space_element_t<Space>{make_smat<Dim>(coordinate)};
  }
};
 
} // namespace linalg