| 1 |
!
-
call blacs_get (0, 0, val)
-
ictxt = val(1)
|
| 2 |
!
|
| 3 |
! Calculate the USERMAP for the BLACS grid, and set up the grid. The
|
| 4 |
! processor ID assignment is done by the same algorithm used by PGHPF,
|
| 5 |
! and may be different for other HPF implementions.
|
| 6 |
!
-
nprow = size (usermap, 1)
-
npcol = size (usermap, 2)
-
do j = 1, npcol
-
do i = 1, nprow
-
usermap(i,j) = (j - 1) * nprow + (i - 1)
-
end do
-
end do
-
call blacs_gridmap (ictxt, usermap, nprow, nprow, npcol)
|
| 7 |
!
-
ip = 0
-
call pdcft3 (x, y, n1, n2, n3, isign, scale, ictxt, ip)
|
| 8 |
!
-
call blacs_gridexit (ictxt)
|
| 9 |
!
|
