52         subroutine fMtxVdimVecMult( a, b, c, na, nb, nc, m, n )  argument
58         real*8  a(n,na),        b(n,nb),        c(n,nc)  local
69                 c(i,1) = a(i,1) * b(i,1)
75                 c(i,1) = a(i,1) * b(i,1)
76                 c(i,2) = a(i,2) * b(i,2)
82                 c(i,1) = a(i,1) * b(i,1)
83                 c(i,2) = a(i,2) * b(i,2)
84                 c(i,3) = a(i,3) * b(i,3)
90                 c(i,1) = a(i,1) * b(i,1)
91                 c(i,2) = a(i,2) * b(i,2)
92                 c(i,3) = a(i,3) * b(i,3)
93                 c(i,4) = a(i,4) * b(i,4)
100                     c(i,j) = a(i,j) * b(i,j)
132 	subroutine flesCp ( a, b, m, n )  argument
136         real*8   a(n,m), b(n,m)  local
140             b(i,j) = a(i,j)
170 	subroutine flesScaleCp ( a, b, s, m, n )  argument
174         real*8   a(n,m), b(n,m), s  local
178             b(i,j) = a(i,j) * s
189 	subroutine flesAdd ( a, b, m, n )  argument
193         real*8   a(n,m), b(n,m)  local
197             b(i,j) = b(i,j) + a(i,j)
208 	subroutine flesSub ( a, b, m, n )  argument
212         real*8   a(n,m), b(n,m)  local
216             b(i,j) = b(i,j) - a(i,j)
247 	real*8 function flesDot2 ( a, b, m, n )
251         real*8   a(n,m), b(n,m)  local
256             flesDot2 = flesDot2 + a(i,j) * b(i,j)
928         subroutine fMtxVdimVecCp( a, b, na, nb, m, n )  argument
934         real*8  a(n,na),        b(n,nb)  local
943                 b(i,1) = a(i,1)
949                 b(i,1) = a(i,1)
950                 b(i,2) = a(i,2)
956                 b(i,1) = a(i,1)
957                 b(i,2) = a(i,2)
958                 b(i,3) = a(i,3)
964                 b(i,1) = a(i,1)
965                 b(i,2) = a(i,2)
966                 b(i,3) = a(i,3)
967                 b(i,4) = a(i,4)
974                     b(i,j) = a(i,j)
989         subroutine fMtxVdimVecDot2( a, b, c, na, nb, m, n )  argument
995         real*8  a(n,na),        b(n,nb),        c(m)  local
1005                 c(1) = c(1) + a(i,1) * b(i,1)
1013                 c(1) = c(1) + a(i,1) * b(i,1)
1014                 c(2) = c(2) + a(i,2) * b(i,2)
1023                 c(1) = c(1) + a(i,1) * b(i,1)
1024                 c(2) = c(2) + a(i,2) * b(i,2)
1025                 c(3) = c(3) + a(i,3) * b(i,3)
1035                 c(1) = c(1) + a(i,1) * b(i,1)
1036                 c(2) = c(2) + a(i,2) * b(i,2)
1037                 c(3) = c(3) + a(i,3) * b(i,3)
1038                 c(4) = c(4) + a(i,4) * b(i,4)
1046                     c(j) = c(j) + a(i,j) * b(i,j)
1061         subroutine fMtxVdimVecDaxpy( a, b, c, na, nb, m, n )  argument
1067         real*8  a(n,na),        b(n,nb),        c(m)  local
1076                 b(i,1) = b(i,1) + c(1) * a(i,1)
1082                 b(i,1) = b(i,1) + c(1) * a(i,1)
1083                 b(i,2) = b(i,2) + c(2) * a(i,2)
1089                 b(i,1) = b(i,1) + c(1) * a(i,1)
1090                 b(i,2) = b(i,2) + c(2) * a(i,2)
1091                 b(i,3) = b(i,3) + c(3) * a(i,3)
1097                 b(i,1) = b(i,1) + c(1) * a(i,1)
1098                 b(i,2) = b(i,2) + c(2) * a(i,2)
1099                 b(i,3) = b(i,3) + c(3) * a(i,3)
1100                 b(i,4) = b(i,4) + c(4) * a(i,4)
1107                     b(i,j) = b(i,j) + c(j) * a(i,j)