refinement_info.hpp

Go to the documentation of this file.

#pragma once
 
#include <uni20/linalg/solve_info.hpp>
 
namespace uni20::linalg
{
template <RealOrComplex Factor, RealOrComplex Working = Factor, RealOrComplex Residual = Working>
  requires(Complex<Factor> == Complex<Working> && Complex<Working> == Complex<Residual>)
struct LuRefinementPrecision
{
    Precision factor;
    Precision working;
    Precision residual;
 
    static constexpr LuRefinementPrecision native()
      requires(!has_runtime_precision_v<Factor> && !has_runtime_precision_v<Working> &&
               !has_runtime_precision_v<Residual>)
    {
      return {native_precision<Factor>(), native_precision<Working>(), native_precision<Residual>()};
    }
};
 
template <Real R> struct SolveError
{
    R residual_infinity_norm;
    R normwise_backward_error;
    R componentwise_backward_error;
};
 
enum class RefinementStatus
{
  converged,
  stagnated,
  iteration_limit,
  factorization_failed,
  correction_failed,
  nonfinite_input,
  nonfinite_result
};
 
struct RefinementInfo
{
    RefinementStatus status = RefinementStatus::converged;
    std::optional<SolveInfo> solve_info = {};
    bool succeeded() const noexcept { return this->status == RefinementStatus::converged; }
};
 
template <Real R> struct RefinementOptions
{
    std::optional<R> backward_error_tolerance = {};
    std::size_t max_iterations = 10;
};
 
template <Real R> struct RefinementColumn
{
    SolveError<R> error;
    std::size_t corrections = 0;
    bool converged = false;
    bool stagnated = false;
};
} // namespace uni20::linalg