52         subroutine fMtxVdimVecMult( a, b, c, na, nb, nc, m, n )  argument
57         integer na,     nb,     nc,     m,      n  local
66         if ( m .eq. 1 ) then
72         else if ( m .eq. 2 ) then
79         else if ( m .eq. 3 ) then
87         else if ( m .eq. 4 ) then
99                 do j = 1, m 
113 	subroutine flesZero ( a, m, n )  argument
116         integer  m, n, i, j  local
117         real*8   a(n,m)
120           do j = 1, m
132 	subroutine flesCp ( a, b, m, n )  argument
135         integer  m, n, i, j  local
136         real*8   a(n,m), b(n,m)
139           do j = 1, m
151  	subroutine flesScale ( a, s, m, n )  argument
154         integer  m, n, i, j     local
155         real*8   a(n,m), s
158           do j = 1, m
170 	subroutine flesScaleCp ( a, b, s, m, n )  argument
173         integer  m, n, i, j  local
174         real*8   a(n,m), b(n,m), s
177           do j = 1, m
189 	subroutine flesAdd ( a, b, m, n )  argument
192         integer  m, n, i, j  local
193         real*8   a(n,m), b(n,m)
196           do j = 1, m
208 	subroutine flesSub ( a, b, m, n )  argument
211         integer  m, n, i, j  local
212         real*8   a(n,m), b(n,m)
215           do j = 1, m
227  	real*8 function flesDot1 ( a, m, n )
230         integer  m, n, i, j  local
231         real*8   a(n,m)
235           do j = 1, m
247 	real*8 function flesDot2 ( a, b, m, n )
250         integer  m, n, i, j  local
251         real*8   a(n,m), b(n,m)
255           do j = 1, m
267 	subroutine flesDaxpy ( x, y, a, m, n )  argument
270         integer  m, n, i, j  local
271         real*8   x(n,m), y(n,m)
275           do j= 1, m
287 	subroutine flesDxpay ( x, y, a, m, n )  argument
290         integer  m, n, i, j  local
291         real*8   x(n,m), y(n,m)
295           do j = 1, m
307 	subroutine flesInv ( x, m, n )  argument
310         integer  m, n, i, j  local
311         real*8   x(n,m)
314           do j = 1, m
328         subroutine fMtxBlkDot2( x, y, c, m, n )  argument
333         integer m,      n  local
334         real*8  x(n,m), y(n),   c(m)
342         m1 = mod(m,8) + 1
472         do j = m1, m, 8
510         subroutine fMtxBlkDaxpy( x, y, c, m, n )  argument
515         integer m,      n  local
516         real*8  x(n,m), y(n),   c(m)
524         m1 = mod(m,8) + 1
621         do j = m1, m, 8
648         subroutine fMtxBlkDyeax( x, y, c, m, n )  argument
653         integer m,      n  local
654         real*8  x(n,m), y(n),   c(m)
662         m1 = mod(m,8) + 1
764         do j = m1, m, 8
791        subroutine fMtxBlkDmaxpy( x, y, c, m, n )  argument
796         integer m,      n  local
797         real*8  x(n,m), y(n),   c(m)
805         m1 = mod(m,8) + 1
901         do j = m1, m, 8
928         subroutine fMtxVdimVecCp( a, b, na, nb, m, n )  argument
933         integer na,     nb,     m,      n  local
940         if ( m .eq. 1 ) then
946         else if ( m .eq. 2 ) then
953         else if ( m .eq. 3 ) then
961         else if ( m .eq. 4 ) then
973                 do j = 1, m 
989         subroutine fMtxVdimVecDot2( a, b, c, na, nb, m, n )  argument
994         integer na,     nb,     m,      n  local
995         real*8  a(n,na),        b(n,nb),        c(m)
1001         if ( m .eq. 1 ) then
1008         else if ( m .eq. 2 ) then
1017         else if ( m .eq. 3 ) then
1028         else if ( m .eq. 4 ) then
1043             do j = 1, m 
1061         subroutine fMtxVdimVecDaxpy( a, b, c, na, nb, m, n )  argument
1066         integer na,     nb,     m,      n  local
1067         real*8  a(n,na),        b(n,nb),        c(m)
1073         if ( m .eq. 1 ) then
1079         else if ( m .eq. 2 ) then
1086         else if ( m .eq. 3 ) then
1094         else if ( m .eq. 4 ) then
1105             do j = 1, m