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)
328         subroutine fMtxBlkDot2( x, y, c, m, n )  argument
334         real*8  x(n,m), y(n),   c(m)  local
354         c(1) = tmp1
364         c(1) = tmp1
365         c(2) = tmp2
377         c(1) = tmp1
378         c(2) = tmp2
379         c(3) = tmp3
393         c(1) = tmp1
394         c(2) = tmp2
395         c(3) = tmp3
396         c(4) = tmp4
412         c(1) = tmp1
413         c(2) = tmp2
414         c(3) = tmp3
415         c(4) = tmp4
416         c(5) = tmp5
434         c(1) = tmp1
435         c(2) = tmp2
436         c(3) = tmp3
437         c(4) = tmp4
438         c(5) = tmp5
439         c(6) = tmp6
459         c(1) = tmp1
460         c(2) = tmp2
461         c(3) = tmp3
462         c(4) = tmp4
463         c(5) = tmp5
464         c(6) = tmp6
465         c(7) = tmp7
491             c(j+0) = tmp1
492             c(j+1) = tmp2
493             c(j+2) = tmp3
494             c(j+3) = tmp4
495             c(j+4) = tmp5
496             c(j+5) = tmp6
497             c(j+6) = tmp7
498             c(j+7) = tmp8
510         subroutine fMtxBlkDaxpy( x, y, c, m, n )  argument
516         real*8  x(n,m), y(n),   c(m)  local
531         tmp1 = c(1)
539         tmp1 = c(1)
540         tmp2 = c(2)
548         tmp1 = c(1)
549         tmp2 = c(2)
550         tmp3 = c(3)
560         tmp1 = c(1)
561         tmp2 = c(2)
562         tmp3 = c(3)
563         tmp4 = c(4)
572         tmp1 = c(1)
573         tmp2 = c(2)
574         tmp3 = c(3)
575         tmp4 = c(4)
576         tmp5 = c(5)
586         tmp1 = c(1)
587         tmp2 = c(2)
588         tmp3 = c(3)
589         tmp4 = c(4)
590         tmp5 = c(5)
591         tmp6 = c(6)
601         tmp1 = c(1)
602         tmp2 = c(2)
603         tmp3 = c(3)
604         tmp4 = c(4)
605         tmp5 = c(5)
606         tmp6 = c(6)
607         tmp7 = c(7)
622             tmp1 = c(j+0)
623             tmp2 = c(j+1)
624             tmp3 = c(j+2)
625             tmp4 = c(j+3)
626             tmp5 = c(j+4)
627             tmp6 = c(j+5)
628             tmp7 = c(j+6)
629             tmp8 = c(j+7)
648         subroutine fMtxBlkDyeax( x, y, c, m, n )  argument
654         real*8  x(n,m), y(n),   c(m)  local
669         tmp1 = c(1)
677         tmp1 = c(1)
678         tmp2 = c(2)
686         tmp1 = c(1)
687         tmp2 = c(2)
688         tmp3 = c(3)
697         tmp1 = c(1)
698         tmp2 = c(2)
699         tmp3 = c(3)
700         tmp4 = c(4)
709         tmp1 = c(1)
710         tmp2 = c(2)
711         tmp3 = c(3)
712         tmp4 = c(4)
713         tmp5 = c(5)
723         tmp1 = c(1)
724         tmp2 = c(2)
725         tmp3 = c(3)
726         tmp4 = c(4)
727         tmp5 = c(5)
728         tmp6 = c(6)
738         tmp1 = c(1)
739         tmp2 = c(2)
740         tmp3 = c(3)
741         tmp4 = c(4)
742         tmp5 = c(5)
743         tmp6 = c(6)
744         tmp7 = c(7)
765             tmp1 = c(j+0)
766             tmp2 = c(j+1)
767             tmp3 = c(j+2)
768             tmp4 = c(j+3)
769             tmp5 = c(j+4)
770             tmp6 = c(j+5)
771             tmp7 = c(j+6)
772             tmp8 = c(j+7)
791        subroutine fMtxBlkDmaxpy( x, y, c, m, n )  argument
797         real*8  x(n,m), y(n),   c(m)  local
812         tmp1 = c(1)
820         tmp1 = c(1)
821         tmp2 = c(2)
829         tmp1 = c(1)
830         tmp2 = c(2)
831         tmp3 = c(3)
840         tmp1 = c(1)
841         tmp2 = c(2)
842         tmp3 = c(3)
843         tmp4 = c(4)
852         tmp1 = c(1)
853         tmp2 = c(2)
854         tmp3 = c(3)
855         tmp4 = c(4)
856         tmp5 = c(5)
866         tmp1 = c(1)
867         tmp2 = c(2)
868         tmp3 = c(3)
869         tmp4 = c(4)
870         tmp5 = c(5)
871         tmp6 = c(6)
881         tmp1 = c(1)
882         tmp2 = c(2)
883         tmp3 = c(3)
884         tmp4 = c(4)
885         tmp5 = c(5)
886         tmp6 = c(6)
887         tmp7 = c(7)
902             tmp1 = c(j+0)
903             tmp2 = c(j+1)
904             tmp3 = c(j+2)
905             tmp4 = c(j+3)
906             tmp5 = c(j+4)
907             tmp6 = c(j+5)
908             tmp7 = c(j+6)
909             tmp8 = c(j+7)
989         subroutine fMtxVdimVecDot2( a, b, c, na, nb, m, n )  argument
995         real*8  a(n,na),        b(n,nb),        c(m)  local
1003             c(1) = 0
1005                 c(1) = c(1) + a(i,1) * b(i,1)
1010             c(1) = 0
1011             c(2) = 0
1013                 c(1) = c(1) + a(i,1) * b(i,1)
1014                 c(2) = c(2) + a(i,2) * b(i,2)
1019             c(1) = 0
1020             c(2) = 0
1021             c(3) = 0
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)
1030             c(1) = 0
1031             c(2) = 0
1032             c(3) = 0
1033             c(4) = 0
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)
1044                 c(j) = 0
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)