27         do i = 1, nshl
28           i0 = (nsd) * (i - 1)
29           Diagl(:,i,1) = xKebe(1:npro,i0+1,i0+1)
30           Diagl(:,i,2) = xKebe(1:npro,i0+2,i0+2)
31           Diagl(:,i,3) = xKebe(1:npro,i0+3,i0+3)
37         do i = 1, nshl
38           Diagl(:,i,4) = xGoC(1:npro,nn2+i,i)
60         integer i,      j  local
68             do i = 1, n
69                 c(i,1) = a(i,1) * b(i,1)
74             do i = 1, n
75                 c(i,1) = a(i,1) * b(i,1)
76                 c(i,2) = a(i,2) * b(i,2)
81             do i = 1, n
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)
89             do i = 1, n
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)
98             do i = 1, n 
100                     c(i,j) = a(i,j) * b(i,j)
116         integer  m, n, i, j  local
119         do i = 1, n
121             a(i,j) = 0.0e-0
135         integer  m, n, i, j  local
138         do i = 1, n
140             b(i,j) = a(i,j)
154         integer  m, n, i, j     local
157         do i = 1, n
159             a(i,j) = a(i,j) * s
173         integer  m, n, i, j  local
176         do i = 1, n
178             b(i,j) = a(i,j) * s
192         integer  m, n, i, j  local
195         do i = 1, n
197             b(i,j) = b(i,j) + a(i,j)
211         integer  m, n, i, j  local
214         do i = 1, n
216             b(i,j) = b(i,j) - a(i,j)
230         integer  m, n, i, j  local
234         do i = 1, n
236             flesDot1 = flesDot1 + a(i,j) * a(i,j)
250         integer  m, n, i, j  local
254         do i = 1, n
256             flesDot2 = flesDot2 + a(i,j) * b(i,j)
270         integer  m, n, i, j  local
274         do i = 1, n
276             y(i,j) = y(i,j) + a * x(i,j)
290         integer  m, n, i, j  local
294         do i = 1, n
296             y(i,j) = a * y(i,j) + x(i,j)
310         integer  m, n, i, j  local
313         do i = 1, n
315             if ( x(i,j) .ne. 0 ) x(i,j) = 1. / x(i,j)
338         integer i,      j,      m1  local
351         do i = 1, n
352             tmp1 = tmp1 + x(i,1) * y(i)
360         do i = 1, n
361             tmp1 = tmp1 + x(i,1) * y(i)
362             tmp2 = tmp2 + x(i,2) * y(i)
372         do i = 1, n
373             tmp1 = tmp1 + x(i,1) * y(i)
374             tmp2 = tmp2 + x(i,2) * y(i)
375             tmp3 = tmp3 + x(i,3) * y(i)
387         do i = 1, n
388             tmp1 = tmp1 + x(i,1) * y(i)
389             tmp2 = tmp2 + x(i,2) * y(i)
390             tmp3 = tmp3 + x(i,3) * y(i)
391             tmp4 = tmp4 + x(i,4) * y(i)
405         do i = 1, n
406             tmp1 = tmp1 + x(i,1) * y(i)
407             tmp2 = tmp2 + x(i,2) * y(i)
408             tmp3 = tmp3 + x(i,3) * y(i)
409             tmp4 = tmp4 + x(i,4) * y(i)
410             tmp5 = tmp5 + x(i,5) * y(i)
426         do i = 1, n
427             tmp1 = tmp1 + x(i,1) * y(i)
428             tmp2 = tmp2 + x(i,2) * y(i)
429             tmp3 = tmp3 + x(i,3) * y(i)
430             tmp4 = tmp4 + x(i,4) * y(i)
431             tmp5 = tmp5 + x(i,5) * y(i)
432             tmp6 = tmp6 + x(i,6) * y(i)
450         do i = 1, n
451             tmp1 = tmp1 + x(i,1) * y(i)
452             tmp2 = tmp2 + x(i,2) * y(i)
453             tmp3 = tmp3 + x(i,3) * y(i)
454             tmp4 = tmp4 + x(i,4) * y(i)
455             tmp5 = tmp5 + x(i,5) * y(i)
456             tmp6 = tmp6 + x(i,6) * y(i)
457             tmp7 = tmp7 + x(i,7) * y(i)
481             do i = 1, n
482                 tmp1 = tmp1 + x(i,j+0) * y(i)
483                 tmp2 = tmp2 + x(i,j+1) * y(i)
484                 tmp3 = tmp3 + x(i,j+2) * y(i)
485                 tmp4 = tmp4 + x(i,j+3) * y(i)
486                 tmp5 = tmp5 + x(i,j+4) * y(i)
487                 tmp6 = tmp6 + x(i,j+5) * y(i)
488                 tmp7 = tmp7 + x(i,j+6) * y(i)
489                 tmp8 = tmp8 + x(i,j+7) * y(i)
520         integer i,      j,      m1  local
532         do i = 1, n
533             y(i) = y(i)
534      1           + tmp1 * x(i,1)
541         do i = 1, n
542             y(i) = y(i)
543      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
552         do i = 1, n
553             y(i) = y(i)
554      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
555      2           + tmp3 * x(i,3)
564         do i = 1, n
565             y(i) = y(i)
566      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
567      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
577         do i = 1, n
578             y(i) = y(i)
579      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
580      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
581      3           + tmp5 * x(i,5)
592         do i = 1, n
593             y(i) = y(i)
594      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
595      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
596      3           + tmp5 * x(i,5) + tmp6 * x(i,6)
608         do i = 1, n
609             y(i) = y(i)
610      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
611      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
612      3           + tmp5 * x(i,5) + tmp6 * x(i,6)
613      4           + tmp7 * x(i,7)
630             do i = 1, n
631                 y(i) = y(i)
632      1               + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
633      2               + tmp3 * x(i,j+2) + tmp4 * x(i,j+3)
634      3               + tmp5 * x(i,j+4) + tmp6 * x(i,j+5)
635      4               + tmp7 * x(i,j+6) + tmp8 * x(i,j+7)
658         integer i,      j,      m1  local
670         do i = 1, n
671             y(i) =
672      1           + tmp1 * x(i,1)
679         do i = 1, n
680             y(i) =
681      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
689         do i = 1, n
690             y(i) =
691      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
692      2           + tmp3 * x(i,3)
701         do i = 1, n
702             y(i) =
703      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
704      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
714         do i = 1, n
715             y(i) =
716      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
717      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
718      3           + tmp5 * x(i,5)
729         do i = 1, n
730             y(i) =
731      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
732      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
733      3           + tmp5 * x(i,5) + tmp6 * x(i,6)
745         do i = 1, n
746             y(i) =
747      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
748      2           + tmp3 * x(i,3) + tmp4 * x(i,4)
749      3           + tmp5 * x(i,5) + tmp6 * x(i,6)
750      4           + tmp7 * x(i,7)
755         do i = 1, n
756             y(i) = 0
773             do i = 1, n
774                 y(i) = y(i)
775      1               + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
776      2               + tmp3 * x(i,j+2) + tmp4 * x(i,j+3)
777      3               + tmp5 * x(i,j+4) + tmp6 * x(i,j+5)
778      4               + tmp7 * x(i,j+6) + tmp8 * x(i,j+7)
801         integer i,      j,      m1  local
813         do i = 1, n
814             y(i) = y(i)
815      1           - tmp1 * x(i,1)
822         do i = 1, n
823             y(i) = y(i)
824      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
832         do i = 1, n
833             y(i) = y(i)
834      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
835      2           - tmp3 * x(i,3)
844         do i = 1, n
845             y(i) = y(i)
846      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
847      2           - tmp3 * x(i,3) - tmp4 * x(i,4)
857         do i = 1, n
858             y(i) = y(i)
859      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
860      2           - tmp3 * x(i,3) - tmp4 * x(i,4)
861      3           - tmp5 * x(i,5)
872         do i = 1, n
873             y(i) = y(i)
874      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
875      2           - tmp3 * x(i,3) - tmp4 * x(i,4)
876      3           - tmp5 * x(i,5) - tmp6 * x(i,6)
888         do i = 1, n
889             y(i) = y(i)
890      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
891      2           - tmp3 * x(i,3) - tmp4 * x(i,4)
892      3           - tmp5 * x(i,5) - tmp6 * x(i,6)
893      4           - tmp7 * x(i,7)
910             do i = 1, n
911                 y(i) = y(i)
912      1               - tmp1 * x(i,j+0) - tmp2 * x(i,j+1)
913      2               - tmp3 * x(i,j+2) - tmp4 * x(i,j+3)
914      3               - tmp5 * x(i,j+4) - tmp6 * x(i,j+5)
915      4               - tmp7 * x(i,j+6) - tmp8 * x(i,j+7)
936         integer i,      j  local
942             do i = 1, n
943                 b(i,1) = a(i,1)
948             do i = 1, n
949                 b(i,1) = a(i,1)
950                 b(i,2) = a(i,2)
955             do i = 1, n
956                 b(i,1) = a(i,1)
957                 b(i,2) = a(i,2)
958                 b(i,3) = a(i,3)
963             do i = 1, n
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)
972             do i = 1, n
974                     b(i,j) = a(i,j)
997         integer i,      j  local
1004             do i = 1, n
1005                 c(1) = c(1) + a(i,1) * b(i,1)
1012             do i = 1, n
1013                 c(1) = c(1) + a(i,1) * b(i,1)
1014                 c(2) = c(2) + a(i,2) * b(i,2)
1022             do i = 1, n
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)
1034             do i = 1, n
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)
1045                 do i = 1, n 
1046                     c(j) = c(j) + a(i,j) * b(i,j)
1069         integer i,      j  local
1075             do i = 1, n
1076                 b(i,1) = b(i,1) + c(1) * a(i,1)
1081             do i = 1, n
1082                 b(i,1) = b(i,1) + c(1) * a(i,1)
1083                 b(i,2) = b(i,2) + c(2) * a(i,2)
1088             do i = 1, n
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)
1096             do i = 1, n
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)
1106                 do i = 1, n 
1107                     b(i,j) = b(i,j) + c(j) * a(i,j)
1139     	do i = 1, nshl
1140            i0 = (nsd) * (i - 1)  
1143              Qtemp(:,i,1) = Qtemp(:,i,1) 
1146              Qtemp(:,i,2) = Qtemp(:,i,2)
1149              Qtemp(:,i,3) = Qtemp(:,i,3)
1188         do i = 1, nshl
1189            i0 = (nsd) * (i - 1)
1193             Qtemp(:,i,1) = Qtemp(:,i,1)
1198             Qtemp(:,i,2) = Qtemp(:,i,2)
1202             Qtemp(:,i,3) = Qtemp(:,i,3)
1212         do i = 1, nshl 
1213            i0 = (nsd) * (i - 1) 
1216             Qtemp(:,i,1) = Qtemp(:,i,1)
1218             Qtemp(:,i,2) = Qtemp(:,i,2)
1220             Qtemp(:,i,3) = Qtemp(1:,i,3)
1257         do i = 1, nshl
1261              Qtemp(:,i,nQs) = Qtemp(:,i,nQs)
1262      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1263      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1264      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1300         do i = 1, nshl
1304              Qtemp(:,i,nQs) = Qtemp(:,i,nQs)
1305      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1306      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1307      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1316         do i = 1, nshl
1317            i0 = nnm2 + i
1320              Qtemp(:,i,nQs) = Qtemp(:,i,nQs)
1359         do i = 1, nshl
1360            i0 = (nsd) * (i - 1)
1364             Qtemp(:,i,1) = Qtemp(:,i,1)
1369             Qtemp(:,i,2) = Qtemp(:,i,2)
1373             Qtemp(:,i,3) = Qtemp(:,i,3)
1383        do i = 1, nshl
1384            i0 = (nsd) * (i - 1)
1387             Qtemp(:,i,1) = Qtemp(:,i,1)
1389             Qtemp(:,i,2) = Qtemp(:,i,2)
1391             Qtemp(:,i,3) = Qtemp(:,i,3)
1399        do i = 1, nshl
1403              Qtemp(:,i,nQs) = Qtemp(:,i,nQs)
1404      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1405      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1406      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1415         do i = 1, nshl
1416            i0 = nnm2 + i
1419              Qtemp(:,i,nQs) = Qtemp(:,i,nQs)
1447         do i = 1, nshl
1448            Diagl(:,i,1) = xTe(1:npro,i,i)
1479         do i = 1, nshl
1482             Qtemp(:,i,nQs) = Qtemp(:,i,nQs) 
1483      &                     + xTe(1:npro,i,j) * Ptemp(:,j,nPs)