GVCSP в Фортране 90

Question

Я хочу точно определить собственное значение и собственный вектор матрицы, который является симметричным c положительно определенным. Однако я не могу использовать функцию GVCSP в Фортране в Linux. Есть ли в Фортране еще одна функция, эквивалентная функции GVCSP.

Вот мой код:

program main

  integer , parameter :: n=4
  real*8 A(n,n), eigval(n), eigvec(n,n),eigvalmatrix(n,n)

  integer i,j

  do i=1,n
      eigval(i)=0
      do j=1,n
          eigvec(i,j)=0
          eigvalmatrix(i,j)=0
      end do
  end do

A(1,1)=-1
A(1,2)=0
A(1,3)=4
A(1,4)=2

A(2,1)=0
A(2,2)=2
A(2,3)=3
A(2,4)=-2

A(3,1)=4
A(3,2)=3
A(3,3)=0
A(3,4)=-4

A(4,1)=2
A(4,2)=-2
A(4,3)=-4
A(4,4)=1

call GVCSP(A, eigvalmatrix, eigval,eigvec)
  end program