da_matmulti.inc

1 subroutine da_matmulti(mata,matb,matc,ni,nj,nab)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer, intent(in) :: ni, nj, nab
10    real, intent(in) :: mata(ni,nab), matb(nab, nj)
11    real, intent(out) :: matc(ni,nj)
12 
13    integer             :: i, j, k            ! Loop counters
14 
15    matc = 0.
16 
17    do i=1,ni
18       do j=1,nj 
19          do k=1,nab
20             matc(i,j) = matc(i,j) + mata(i,k)*matb(k,j)
21          end do
22       end do
23    end do
24 
25 end subroutine da_matmulti
26 
27