www.pudn.com > Wang_Lecture_2012.rar > axb.f, change:2001-11-22,size:634b


c*******************************************************************************
c*******************************************************************************
        subroutine axb(a,b,n,l,m,c)
        implicit none
c
c	First implemented in Potsdam, Feb, 1999
c	Last modified: Potsdam, Nov, 2001, by R. Wang
c
c
c       calculate c=a*b
c
        integer n,l,m
        double precision a(n,l),b(l,m),c(n,m)
c
        integer i,j,k
c
        do j=1,m
          do i=1,n
            c(i,j)=0.d0
            do k=1,l
              c(i,j)=c(i,j)+a(i,k)*b(k,j)
            enddo
          enddo
        enddo
        return
        end