numerical_config.hpp

Go to the documentation of this file.

/*
 * Copyright © 2025-2026 Akreon LLC
 * SPDX-License-Identifier: Apache-2.0
 */
 
#pragma once
#include "akreon/linalg/space/space_concepts.hpp"
#include "akreon/linalg/vector/algebra_util.hpp"
namespace config {
template <linalg::Real T>
struct FactorConfig {
  std::size_t steps{64}; 
  T tol{1e-8};           
  T floor_eps{1e-10};    
};
 
} // namespace config
namespace linalg {
 
enum struct NumericalEvent : uint8_t {
  NONE = 0,
 
  // detection / observation
  RANK_DEFICIENCY,
  INDEFINITE_FACTORIZATION,
  NEAR_SINGULAR_TRIANGULAR_SOLVE,
  NUMERICAL_SAFEGUARD,
  DB_MAX_ITERS_REACHED,
 
  // failure
  FACTORIZATION_FAILED,
  EVD_FAILED
};
 
inline constexpr bool is_fatal_event(NumericalEvent event) {
  return event == NumericalEvent::FACTORIZATION_FAILED ||
         event == NumericalEvent::EVD_FAILED;
}
 
} // namespace linalg