add examples

Author: perazz

example/linalg/CMakeLists.txt

@@ -12,6 +12,9 @@ ADD_EXAMPLE(is_skew_symmetric)
 ADD_EXAMPLE(is_square)
 ADD_EXAMPLE(is_symmetric)
 ADD_EXAMPLE(is_triangular)
+ADD_EXAMPLE(inverse1)
+ADD_EXAMPLE(inverse2)
+ADD_EXAMPLE(inverse3)
 ADD_EXAMPLE(outer_product)
 ADD_EXAMPLE(trace)
 ADD_EXAMPLE(state1)

example/linalg/example_inverse1.f90

@@ -0,0 +1,22 @@
+! Matrix inversion example: operator interface
+program example_inverse1
+  use stdlib_linalg_constants, only: dp
+  use stdlib_linalg, only: operator(.inv.),eye
+  implicit none
+
+  real(dp) :: A(2,2), Am1(2,2)
+
+  ! Input matrix (NB Fortran is column major! input columns then transpose)
+  A = transpose(reshape( [4, 3, &
+                          3, 2], [2,2] ))
+
+  ! Invert matrix
+  Am1 = .inv.A
+
+  print *, '         |',Am1(1,:),'|' ! | -2  3 |
+  print *, ' inv(A)= |',Am1(2,:),'|' ! |  3 -4 |
+
+  ! Final check 
+  print *, 'CHECK passed? ',matmul(A,Am1)==eye(2)
+
+end program example_inverse1

example/linalg/example_inverse2.f90

@@ -0,0 +1,22 @@
+! Matrix inversion example: function interface
+program example_inverse2
+  use stdlib_linalg_constants, only: dp
+  use stdlib_linalg, only: inv,eye
+  implicit none
+
+  real(dp) :: A(2,2), Am1(2,2)
+
+  ! Input matrix (NB Fortran is column major! input columns then transpose)
+  A = transpose(reshape( [4, 3, &
+                          3, 2], [2,2] ))
+
+  ! Invert matrix
+  Am1 = inv(A)
+
+  print *, '         |',Am1(1,:),'|' ! | -2  3 |
+  print *, ' inv(A)= |',Am1(2,:),'|' ! |  3 -4 |
+
+  ! Final check 
+  print *, 'CHECK passed? ',matmul(A,Am1)==eye(2)
+
+end program example_inverse2

example/linalg/example_inverse3.f90

@@ -0,0 +1,23 @@
+! Matrix inversion example: in-place inversion
+program example_inverse3
+  use stdlib_linalg_constants, only: dp
+  use stdlib_linalg, only: invert,eye
+  implicit none
+
+  real(dp) :: A(2,2), Am1(2,2)
+
+  ! Input matrix (NB Fortran is column major! input columns then transpose)
+  A = transpose(reshape( [4, 3, &
+                          3, 2], [2,2] ))
+  Am1 = A
+
+  ! Invert matrix
+  call invert(Am1)
+
+  print *, '         |',Am1(1,:),'|' ! | -2  3 |
+  print *, ' inv(A)= |',Am1(2,:),'|' ! |  3 -4 |
+
+  ! Final check 
+  print *, 'CHECK passed? ',matmul(A,Am1)==eye(2)
+
+end program example_inverse3