separator =
' ';
22 for(
size_trow = 0; row < v.size(); row++) {
23 for(
size_tcol = 0; col < v[row].size(); col++)
24out << std::left << std::setw(width) << std::setfill(separator)
39 const intrange2 = range >> 1;
45 for(
inti = 0; i <
mat_size; i++) {
49A[i][j] =
static_cast<double>(std::rand() % range - range2);
52std::clock_t start_t = std::clock();
54std::clock_t end_t = std::clock();
55std::cout <<
"Time taken: " 56<<
static_cast<double>(end_t - start_t) / CLOCKS_PER_SEC <<
"\n";
58std::cout <<
"A = \n"<< A <<
"\n";
59std::cout <<
"L = \n"<< L <<
"\n";
60std::cout <<
"U = \n"<< U <<
"\n";
67std::cout <<
"Determinant test 1...";
70std::cout <<
"passed\n";
72std::cout <<
"Determinant test 2...";
75std::cout <<
"passed\n";
77std::cout <<
"Determinant test 3...";
78 matrix<float>A3({{1.2, 2.3, 3.4}, {4.5, 5.6, 6.7}, {7.8, 8.9, 9.0}});
80std::cout <<
"passed\n";
85std::srand(std::time(NULL));
std::ostream & operator<<(std::ostream &out, matrix< T > const &v)
Functions associated with LU Decomposition of a square matrix.
double determinant_lu(const matrix< T > &A)
int lu_decomposition(const matrix< T > &A, matrix< double > *L, matrix< double > *U)
std::vector< std::valarray< T > > matrix
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4