[139f46]: LibDBlasL3F77.mod History

Parent: [3b207b] (diff)

Child: [28b809] (diff)

Download this file

LibDBlasL3F77.mod    123 lines (105 with data), 5.3 kB 

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
IMPLEMENTATION MODULE LibDBlasL3F77;



  (*------------------------------------------------------------------------*)

  (* Schnittstell zu Level 3 Blas FORTRAN77 Subroutinen.                    *)

  (* Interface to level 3 Blas FORTRAN77 subroutines.                       *)

  (*------------------------------------------------------------------------*)

  (* Letzte Veraenderung                                                    *)

  (*                                                                        *)

  (* 29.10.17, MRi: Erstellen der ersten Version nur mit dgemm              *)

  (* 23.06.18, MRi: Einfuegen von zgemm                                     *)

  (*------------------------------------------------------------------------*)

  (* Offene Punkte                                                          *)

  (*                                                                        *)

  (* - Testen                                                               *)

  (* - Lot of other routines ;-)                                            *)

  (*------------------------------------------------------------------------*)

  (* Implementation : Michael Riedl                                         *)

  (* Licence        : GNU Lesser General Public License (LGPL)              *)

  (*------------------------------------------------------------------------*)



  (* $Id: LibDBlasL3F77.mod,v 1.3 2018/09/12 13:20:49 mriedl Exp mriedl $ *)



FROM LibDBlasL1F77 IMPORT CHAR1,INTEGER4,REAL4,DOUBLEPRECISION,DOUBLECOMPLEX;

FROM LibDBlasLxF77 IMPORT sgemm_,dgemm_,dgemmomp,dgemmomp2,zgemm_;





PROCEDURE sgemm(    TA    : CHAR1;

                    TB    : CHAR1;

                    M,N,K : INTEGER4;

                    Alpha : REAL4;

                VAR A     : ARRAY OF ARRAY OF REAL4;

                    lda   : INTEGER4;

                VAR B     : ARRAY OF ARRAY OF REAL4;

                    ldb   : INTEGER4;

                    Beta  : REAL4;

                VAR C     : ARRAY OF ARRAY OF REAL4;

                    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;

      sgemm_(TB,TA,N,M,K,Alpha,B[0][0],ldb,A[0][0],lda,Beta,C[0][0],ldc);

END sgemm;



PROCEDURE dgemm(    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;

      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.