da_matmultiover.inc

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