*=*=*=*= aMXVA_adj.html =*=*=*=*
SUBROUTINE aMXVA_adj (SA,IAC,IAR,SB,IB,SC,IC,NRA,NCA) IMPLICIT NONE INTEGER IAC,IAR,IB,IC,NRA,NCA REAL SA(1), SB(1), SC(1) INTEGER I,J DO 220 J =1,NCA DO 210 I =1, NRA SB( 1 + (J-1)*IB )=SB( 1 + (J-1)*IB )+ . SA( 1 + (I-1)*IAC + (J-1)*IAR )*SC( 1 + (I-1)*IC ) 210 CONTINUE 220 CONTINUE DO 100 I = 1, NRA SC( 1 + (I-1)*IC ) = 0. c ( C(i) := 0. ) 100 CONTINUE c multiply matrix from sa and vector from sb c DO 220 J = 1, NCA c DO 210 I = 1, NRA c SC( 1 + (I-1)*IC ) c & = SC( 1 + (I-1)*IC ) c & + SA( 1 + (I-1)*IAC + (J-1)*IAR ) c & * SB( 1 + (J-1)*IB ) c ( C(i) := C(i) + A(i,j)*B(j) ) c210 CONTINUE c220 CONTINUE RETURN END