da_innerprod.inc

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