C++ Armadillo :: Inverse
---------------------------------------------------------------------------------------------------------------
Description::
This is the basic C++ program of C++ armadillo library which shows how to find inverse of the given matrix. It can be carried out by the help of inbuilt function "inv" in armadillo library.
- If A is not square sized, a std::logic_error exception is thrown
- If A appears to be singular:
- B = inv(A) resets B and throws a std::runtime_error exception
- inv(B,A) resets B and returns a bool set to false (exception is not thrown)
---------------------------------------------------------------------------------------------------------------
Program::
#include <iostream>
#include <armadillo>using namespace std;
using namespace arma;
int main()
{ //initialize the random generator
//Create a 5 x 5 random matrix and printing it
mat A = randu<mat>(5,5); //random matrix of size 5x5 (generared by this syntax)
mat B = inv(A); //inverse of matrix
cout << "Matrix A::\n"<<endl;
cout<< A << endl;
cout << "Matrix B(Inverse)::\n"<<endl;
cout << B << endl;
return 0;
}
---------------------------------------------------------------------------------------------------------------
Commands to run program on terminal::
$ g++ filename.cpp -o objectname -O2 -larmadillo
$ ./objectname
---------------------------------------------------------------------------------------------------------------
Output::
Matrix A::
0.7868 0.4049 0.2742 0.8571 0.2393
0.2505 0.2513 0.5610 0.4998 0.3201
0.7107 0.0227 0.1400 0.4194 0.9105
0.9467 0.5206 0.5439 0.7443 0.1648
0.0193 0.3447 0.5219 0.2492 0.2455
Matrix B(Inverse)::
-1.5142 -0.2802 0.4407 2.0824 -1.1915
1.8683 -5.3369 0.0782 -0.3022 5.0493
-2.5073 3.0391 -0.2686 1.5127 -1.5379
2.3545 2.1098 -0.6624 -1.9925 -1.2510
0.4364 -1.0874 1.0987 -0.9328 1.6167
Comments