matvec.f

Go to the documentation of this file.

        subroutine matvec(xm,v1,v2)
c  computes 3-vector v2 as product of 3x3 matrix xm and 3-vector v1
 
c       Program written by:     frederick s. patt
c                               general sciences corporation
c                               october 19, 1992
 
c  calling arguments
 
c  name         Type    i/o     description
c
c  xm(3,3)      r*4      i      input matrix
c  v1(3)        r*4      i      input vector
c  v2(3)        r*4      i      output vector
 
        real*4 v1(3),v2(3),xm(3,3)
        do i=1,3
            v2(i) = 0.0
            do j=1,3
                v2(i) = v2(i) + xm(i,j)*v1(j) 
            end do
        end do
        return
        end