MPI и Фортран

      include 'mpif.h'
      real A(10,10), B(10,10), R(10,10), R1(10,10)
      integer h, rank, ierr, size
      call MPI_INIT(ierr)
      call MPI_COMM_SIZE(MPI_COMM_WORLD, size, ierr)
      call MPI_COMM_RANK(MPI_COMM_WORLD, rank, ierr)
      
       h = 10/size
 
      DO I=1,10
        DO J=1,10
            A(i,j) = 1
            B(i,j) = 2
        enddo
      enddo  
      
      DO I=h*rank + 1, h*rank + h       
        DO J=1,10
            R(i,j) = 0.0
            DO L=1,10
                R(I,J)=R(I,J)+A(I,L)*B(L,J)
            enddo
        enddo
      enddo  
 
      call MPI_GATHER(R(h*rank + 1,h*rank + 1), h, MPI_DOUBLE_PRECISION, R1(1,1), h, MPI_DOUBLE_PRECISION, 0, MPI_COMM_WORLD, ierr)
      
      open(5, file='matrix')
  6  FORMAT(10F10.3)
      do l=1,10
        write(5,6)(R1(l, k), k=1,10)
      enddo
      
      call MPI_FINALIZE(ierr)
      end

DO I=rank + 1, 10, size