cone_block.hpp

Go to the documentation of this file.

/*
 * Copyright © 2025-2026 Akreon LLC
 * SPDX-License-Identifier: Apache-2.0
 */
 
#pragma once
 
#include "cone_traits.hpp"
 
namespace conic {
 
template <affine::AffineMapLike Map, typename... Cs>
struct ConeBlock {
  using MatT = affine::affine_mat_t<Map>;
  using T = linalg::matrix_scalar_t<MatT>;
 
  static constexpr std::size_t M = linalg::mat_rows_v<MatT>;
  static constexpr std::size_t N = linalg::mat_cols_v<MatT>;
 
  constexpr explicit ConeBlock(const Map& map_in, Cs... cs)
      : map(map_in), cones(std::forward<Cs>(cs)...) {}
 
  constexpr explicit ConeBlock(const Map& map_in, std::tuple<Cs...> cs)
      : map(map_in), cones(std::move(cs)) {}
 
  constexpr bool validate() const {
    // (dakre) lambda functor to perform the soft check
    std::size_t n_dims = 0;
    auto perform_check = [&](const auto& arg) {
      using ArgT = cone_scalar_t<decltype(arg)>;
      static_assert(
          std::same_as<ArgT, T>,
          "Cone scalar type is mismatched with cone block scalar type");
      n_dims += cone_dim_t<decltype(arg)>;
    };
    // (dakre) apply fold expression over comma operator on stored cones
    std::apply([&](const auto&... args) { (perform_check(args), ...); }, cones);
    return n_dims == M;
  }
 
  Map map; 
  std::tuple<Cs...> cones; 
};
 
} // namespace conic