#include #include "inversion.h" #include "storage_adaptors.hpp" using namespace boost::numeric::ublas; int main(int , char * []) { double initialValues[3][3] = { { 1, 2, 3 }, { 5, 1, 4 }, { 6, 7, 1 } }; matrix A(3, 3), Z(3, 3); A = make_matrix_from_pointer(initialValues); InvertMatrix(A, Z); std::cout << "A=" << A << std::endl << "Z=" << Z << std::endl; matrix result(3,3); result = prod(A,Z); std::cout << "Result=" << result << std::endl; identity_matrix identity(3); result = result-identity; double sum = 0; for (unsigned long i=0; i