GET_MINVALUES.f90

Go to the documentation of this file.

!    Copyright (C) 2012 The SPEED FOUNDATION
!    Author: Ilario Mazzieri
!
!    This file is part of SPEED.
!
!    SPEED is free software; you can redistribute it and/or modify it
!    under the terms of the GNU Affero General Public License as
!    published by the Free Software Foundation, either version 3 of the
!    License, or (at your option) any later version.
!
!    SPEED is distributed in the hope that it will be useful, but
!    WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
!    Affero General Public License for more details.
!
!    You should have received a copy of the GNU Affero General Public License
!    along with SPEED.  If not, see <http://www.gnu.org/licenses/>.
 
 
 
     subroutine  get_minvalues(i_glo, v_glo, n_glo, i_loc, n_loc, np)
     
  
     
     implicit none
     
     integer*4 :: n_glo, n_loc, np
     integer*4 :: i,j
 
     integer*4, dimension(1) :: pos
     integer*4, dimension(np) :: ind     
     integer*4, dimension(np) :: pos_glo
     integer*4, dimension(n_glo) :: i_glo
     integer*4, dimension(n_loc) :: i_loc
          
     real*8, dimension(n_glo) :: v_glo
     real*8, dimension(n_loc) :: v_loc
     real*8, dimension(np) :: val
     
    ! write(*,*) v_glo
    ! read(*,*)
    ! write(*,*) i_glo
    ! read(*,*)
     
 
     do i = 1, n_loc
        do j = 1, np
          
          pos_glo(j) = n_loc*(j-1) + i
          ind(j)     = i_glo(n_loc*(j-1) + i)
          val(j)     = v_glo(n_loc*(j-1) + i)
 
    !      write(*,*) j, val(j), ind(j)
    !      read(*,*)
        enddo  
        
       pos = minloc(val)
    !   write(*,*) i, val 
    !   write(*,*) pos(1),pos_glo(pos(1))
    !   write(*,*) '====================='
       !i_loc(i) = ind(pos(1))
       i_loc(i) = pos_glo(pos(1))
            
     enddo         
 
     end subroutine  get_minvalues