Référence du module mod_mat

Fonctions/Subroutines
subroutine	matmul (n, a, b, c)
	matrix-matrix product Plus de détails...

Documentation de la fonction/subroutine

matmul()

subroutine mod_mat::matmul	(	integer, intent(in)	n,
		real(kind=kind(1.d0)), dimension(n,n), intent(in)	a,
		real(kind=kind(1.d0)), dimension(n,n), intent(in)	b,
		real(kind=kind(1.d0)), dimension(n,n), intent(out)	c
	)

matrix-matrix product

Paramètres

[in]	n	dimension
[in]	a	matrix
[in]	b	matrix

Renvoie

c matrix with c = a*b

Définition à la ligne 11 du fichier mod_mat.f90.

  implicit none  ! on force l'initialisation de toutes les variables
  integer, intent(in)       :: n
  real(kind=kind(1.d0)), intent(in)    :: a(n,n),b(n,n)
  real(kind=kind(1.d0)), intent(out)   :: c(n,n)
  integer    :: i,j,k,ii,jj,kk
  integer    :: is

  c = 0.
  is = 20
  do j = 1,n,is
     do i = 1,n,is
        do k = 1,n,is
           do jj = j,min(n,j+is-1)
              do ii = i,min(n,i+is-1)
                 do kk = k,min(n,k+is-1)
                    c(ii,jj) = c(ii,jj) + a(ii,kk)*b(kk,jj)
                 end do
              end do
           end do
        end do
     end do
  end do