File: example_qr.f90

package info (click to toggle)
fortran-stdlib 0.8.1-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 34,008 kB
  • sloc: f90: 24,178; ansic: 1,244; cpp: 623; python: 119; makefile: 13
file content (15 lines) | stat: -rw-r--r-- 321 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
program example_qr  
  use stdlib_linalg, only: qr
  implicit none
  real :: A(104, 32), Q(104,32), R(32,32)
  
  ! Create a random matrix
  call random_number(A)

  ! Compute its QR factorization (reduced)
  call qr(A,Q,R)

  ! Test factorization: Q*R = A 
  print *, maxval(abs(matmul(Q,R)-A)) 

end program example_qr