doc/html/test_sytrf_singular.html

 // ========================================================================= //
 // This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
 // and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
 // ========================================================================= //

 // Containers.
 #include <myramath/dense/Matrix.h>
 #include <myramath/dense/transpose.h>

 // Algorithms.
 #include <myramath/dense/sytrf.h>

 // Reporting.
 #include <tests/myratest.h>

 using namespace myra;

 ADD_TEST("sytrf_singular","[dense][lapack]")
   {
   // Make random A ..
   auto A = Matrix<double>::random(5,5);
   A = A+transpose(A);
   // .. make singular, third row/column is all zero.
   A.column(3).zero();
   A.row(3).zero();
   REQUIRE_EXCEPTION( sytrf_inplace('L','R',A); );
   }

myra::Matrix::random
static Matrix< Number > random(int I, int J)
Generates a random Matrix of specified size.
Definition: Matrix.cpp:353

myra
Definition: syntax.dox:1

sytrf.h
LDL&#39; decompositions for real symmetric Matrix A (indefinite is fine).

transpose.h
Returns a transposed copy of a Matrix. The inplace version only works on a square operand...

Matrix.h
General purpose dense matrix container, O(i*j) storage.