! 
! Copyright (c) 2017, NVIDIA CORPORATION.  All rights reserved.
!
! Licensed under the Apache License, Version 2.0 (the "License");
! you may not use this file except in compliance with the License.
! You may obtain a copy of the License at
!
!     http://www.apache.org/licenses/LICENSE-2.0
!
! Unless required by applicable law or agreed to in writing, software
! distributed under the License is distributed on an "AS IS" BASIS,
! WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
! See the License for the specific language governing permissions and
! limitations under the License.
! 


#include "mmul_dir.h"

subroutine F90_matmul_cplx16_contmxm(dest, src1, src2, k_extent,m_extent,n_extent)

  COMPLEX*16, dimension(k_extent,m_extent) :: src1
  COMPLEX*16, dimension(m_extent,n_extent) :: src2
  COMPLEX*16, dimension(k_extent,n_extent) :: dest
  DESC_INT k_extent
  DESC_INT n_extent
  DESC_INT m_extent

  DESC_INT k
  DESC_INT n
  DESC_INT m

  do n=1,n_extent
    do k=1,k_extent
      dest(k,n) = 0
    end do
  end do
  do n=1,n_extent
    do m=1,m_extent
      do k=1,k_extent
        dest(k,n) = dest(k,n) + src1(k,m) * src2(m,n)
      end do
   end do
  end do

end subroutine