We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 21fa91f commit 7d8ca7dCopy full SHA for 7d8ca7d
example/linalg/example_schur_eigvals.f90
@@ -4,13 +4,12 @@ program example_schur_eigenvalues
4
use stdlib_linalg, only: schur
5
use stdlib_linalg_constants, only: dp
6
implicit none
7
- real(dp), allocatable :: A(:,:), T(:,:), Z(:,:)
8
- complex(dp), allocatable :: eigenvalues(:)
9
- integer :: n
10
-
+
+ integer, parameter :: n = 5
+ real(dp), dimension(n,n) :: A, T, Z
+ complex(dp), dimension(n) :: eigenvalues
11
12
! Create a random real-valued square matrix
- n = 5
13
- allocate(A(n,n), T(n,n), Z(n,n), eigenvalues(n))
14
call random_number(A)
15
16
! Compute the Schur decomposition and eigenvalues
@@ -29,6 +28,5 @@ program example_schur_eigenvalues
29
28
! Test factorization: Z*T*Z^T = A
30
print *, "Max error in reconstruction:", maxval(abs(matmul(Z, matmul(T, transpose(Z))) - A))
31
32
- deallocate(A, T, Z, eigenvalues)
33
end program example_schur_eigenvalues
34
0 commit comments