Beispiel-Programm in Fortran 90

Sourcen:

• lutest.f90 , kinds.f90

benutzte LAPACK-Routinen:

• dgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
  Lösen eines linearen Gleichungssystems $a x = b$ mit mehreren rechten Seiten

• dlamch(name)
  div. rechnerabhängige Parameter,
  name = 'epsilon' $\rightarrow$ Maschinengenauigkeit

• dlange(norm, m, n, a, lda, work)
  div. Normen der Matrix
  norm = 'i' $\rightarrow \Vert a\Vert _\infty$    (max. Zeilensumme)

• dgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
  BLAS Matrix-Multiplikation: $c = \alpha a b + \beta c$

Geschwindigkeits-Messungen:

• Vergleich mit Matlab:

  
  

  

  
  

• parallele Performance: