--- a/LibDBlasL3F77.mod
+++ b/LibDBlasL3F77.mod
@@ -7,6 +7,7 @@
   (* Letzte Veraenderung                                                    *)
   (*                                                                        *)
   (* 29.10.17, MRi: Erstellen der ersten Version nur mit dgemm              *)
+  (* 23.06.18, MRi: Einfuegen von zgemm                                     *)
   (*------------------------------------------------------------------------*)
   (* Offene Punkte                                                          *)
   (*                                                                        *)
@@ -19,8 +20,10 @@
 
   (* $Id: LibDBlasL3F77.mod,v 1.1 2017/11/01 10:42:09 mriedl Exp mriedl $ *)
 
-FROM LibDBlasL1F77 IMPORT CHAR1,INTEGER4,REAL4,DOUBLEPRECISION;
-FROM LibDBlasLxF77 IMPORT sgemm_,dgemm_;
+FROM LibDBlasL1F77 IMPORT CHAR1,INTEGER4,REAL4,DOUBLEPRECISION,DOUBLECOMPLEX;
+FROM LibDBlasLxF77 IMPORT sgemm_,dgemm_,dgemmomp,dgemmomp2,zgemm_;
+
+IMPORT TIO;
 
 PROCEDURE sgemm(    TA    : CHAR1;
                     TB    : CHAR1;
@@ -68,4 +71,53 @@
       dgemm_(TB,TA,N,M,K,Alpha,B[0][0],ldb,A[0][0],lda,Beta,C[0][0],ldc);
 END dgemm;
 
+PROCEDURE dgemmOMP(    TA    : CHAR1;
+                       TB    : CHAR1;
+                       M,N,K : INTEGER4;
+                       Alpha : DOUBLEPRECISION;
+                   VAR A     : ARRAY OF ARRAY OF DOUBLEPRECISION;
+                       lda   : INTEGER4;
+                   VAR B     : ARRAY OF ARRAY OF DOUBLEPRECISION;
+                       ldb   : INTEGER4;
+                       Beta  : DOUBLEPRECISION;
+                   VAR C     : ARRAY OF ARRAY OF DOUBLEPRECISION;
+                       ldc   : INTEGER4);
+
+          (*---------------------------------------------------------------*)
+          (* Aufruf der Fortran Version von BLAS3 subroutine dgemm         *)
+          (*---------------------------------------------------------------*)
+
+BEGIN
+      lda:=HIGH(A[0]) + 1;
+      ldb:=HIGH(B[0]) + 1;
+      ldc:=HIGH(C[0]) + 1;
+(*
+      dgemmomp2(TB,TA,N,M,K,Alpha,B[0][0],ldb,A[0][0],lda,Beta,C[0][0],ldc);
+ *)
+      dgemmomp (TB,TA,N,M,K,Alpha,B[0][0],ldb,A[0][0],lda,Beta,C[0][0],ldc);
+END dgemmOMP;
+
+PROCEDURE zgemm(    TA    : CHAR1;
+                    TB    : CHAR1;
+                    M,N,K : INTEGER4;
+                    Alpha : DOUBLECOMPLEX;
+                VAR A     : ARRAY OF ARRAY OF DOUBLECOMPLEX;
+                    lda   : INTEGER4;
+                VAR B     : ARRAY OF ARRAY OF DOUBLECOMPLEX;
+                    ldb   : INTEGER4;
+                    Beta  : DOUBLECOMPLEX;
+                VAR C     : ARRAY OF ARRAY OF DOUBLECOMPLEX;
+                    ldc   : INTEGER4);
+
+          (*---------------------------------------------------------------*)
+          (* Aufruf der Fortran Version von BLAS3 subroutine zgemm         *)
+          (*---------------------------------------------------------------*)
+
+BEGIN
+      lda:=HIGH(A[0]) + 1;
+      ldb:=HIGH(B[0]) + 1;
+      ldc:=HIGH(C[0]) + 1;
+      zgemm_(TB,TA,N,M,K,Alpha,B[0][0],ldb,A[0][0],lda,Beta,C[0][0],ldc);
+END zgemm;
+
 END LibDBlasL3F77.