53 template<
typename dtype>
64 auto lMatrix = inMatrix.template astype<double>();
70 lMatrix(row, col) = 0.;
76 const double& a_kk = lMatrix(k, k);
84 lMatrix(i, k) /= lMatrix(k, k);
88 lMatrix(i,
j) -= lMatrix(i, k) * lMatrix(
j, k);
#define THROW_RUNTIME_ERROR(msg)
Definition: Error.hpp:40
#define STATIC_ASSERT_ARITHMETIC(dtype)
Definition: StaticAsserts.hpp:39
const Shape & shape() const noexcept
Definition: NdArrayCore.hpp:4511
uint32 rows
Definition: Core/Shape.hpp:44
bool issquare() const noexcept
Definition: Core/Shape.hpp:125
uint32 cols
Definition: Core/Shape.hpp:45
constexpr auto j
Definition: Core/Constants.hpp:42
Definition: cholesky.hpp:41
NdArray< double > cholesky(const NdArray< dtype > &inMatrix)
Definition: cholesky.hpp:54
auto sqrt(dtype inValue) noexcept
Definition: sqrt.hpp:48
Shape shape(const NdArray< dtype > &inArray) noexcept
Definition: Functions/Shape.hpp:42
std::uint32_t uint32
Definition: Types.hpp:40