decomposition_precision.hpp

Go to the documentation of this file.

#pragma once
 
#include <optional>
 
#include <uni20/core/scalar_concepts.hpp>
#include <uni20/linalg/precision.hpp>
#include <uni20/tensor/precision.hpp>
 
namespace uni20::linalg::detail
{
template <class S>
concept DecompositionScalar = RealOrComplex<S>;
 
template <TensorView A> Precision decomposition_precision(A const& a, std::optional<Precision> requested)
{
  auto p = requested ? *requested : precision_of(a);
  require_precision<tensor_element_t<A>>(p);
  (void)p.bit_count();
  return p;
}
} // namespace uni20::linalg::detail