Commit b90c52e8 authored by Josef Hrabal's avatar Josef Hrabal
Browse files

Added OpenCoarrays: docs.it4i/salomon/software/numerical-languages/

Added 'Fortran' section under salomon/software/numerical-languages
Added 'OpenCoarrays' page under /salomon/software/numerical-language/fortran
parent ab77b1a2
Pipeline #2120 passed with stages
in 1 minute and 10 seconds
# OpenCoarrays
## Introduction
Coarray Fortran (CAF) is an extension of Fortran language and offers a simple interface for parallel processing and memory sharing.
The advantage is that only small changes are required to convert existing Fortran code to support a robust and potentially efficient parallelism.
A CAF program is interpreted as if it was replicated a number of times and all copies were executed asynchronously.
The number of copies is decided at execution time. Each copy (called *image*) has its own private variables.
The variable syntax of Fortran language is extended with indexes in square brackets (called *co-dimension*), which represents a reference to data distributed across images.
By default, the CAF is using Message Passing Interface (MPI) for lower-level communication, so there are some similarities with MPI.
Read more on <>
### Indexing of coarray images
Indexing of individual images can be shown on the simple *Hello World* program:
program hello_world
implicit none
print *, 'Hello world from image ', this_image() , 'of', num_images()
end program hello_world
* num_images() - returns the number of all images
* this_image() - returns the image index - numbered from 1 to num_images()
### Co-dimension variables declaration
Coarray variables can be declared with the **codimension[*]** attribute or by adding trailing index **[*]** after the variable name.
Notice, the ***** character always has to be in the square brackets.
integer, codimension[*] :: scalar
integer :: scalar[*]
real, dimension(64), codimension[*] :: vector
real :: vector(64)[*]
### Images synchronization
Because each image is running on its own, the image synchronization is needed to ensure, that all altered data are distributed to all images.
Synchronization can be done across all images or only between selected images. Be aware, that selective synchronization can lead to the race condition problems like deadlock.
Example program:
program synchronization_test
implicit none
integer :: i ! Local variable
integer :: numbers[*] ! Scalar coarray
! Genereate random number on image 1
if (this_image() == 1) then
numbers = floor(rand(1) * 1000)
! Distribute information to other images
do i = 2, num_images()
numbers[i] = numbers
end do
end if
sync all ! Barrier to synchronize all images
print *, 'The random number is', numbers
end program synchronization_test
* sync all - Synchronize all images between each other
* sync images(*) - Synchronize this image to all other
* sync images(*index*) - Synchronize this image to image with *index*
!!! note
**number** is the local variable while **number[*index*]** accesses the variable in a specific image.
**number[this_image()]** is the same as **number**.
### Compiling CAF program
The preferred method for compiling a CAF program is by invoking the *caf* compiler wrapper.
The above mentioned *Hello World* program can be compiled as follows:
$ caf hello_world.f90 -o hello_world.x
!!! warning
The input file extension has to be **.f90** or **.F90** to be interpreted as *Fortran 90*.
If the input file extension is **.f** or **.F** the source code will be interpreted as *Fortran 77*.
Another method for compiling is by invoking the *mpif90* compiler wrapper directly:
$ mpif90 hello_world.f90 -o hello_world.x -L $CAF_LIBRARY_PATH -fcoarray=lib -lcaf_mpi
### Running CAF program
The preferred method for running a CAF program is by invoking the *cafrun* wrapper:
$ cafrun -np 4 ./hello_world.x
Hello world from image 1 of 4
Hello world from image 2 of 4
Hello world from image 3 of 4
Hello world from image 4 of 4
where **-np 4** is number of images to run.
Another method is using the *mpiexec* directly:
$ mpiexec -np 4 ./synchronization_test.x
The random number is 242
The random number is 242
The random number is 242
The random number is 242
The parameters of **cafrun** and **mpiexec** are the same.
\ No newline at end of file
......@@ -116,6 +116,8 @@ pages:
- Matlab: salomon/software/numerical-languages/
- Octave: salomon/software/numerical-languages/
- R: salomon/software/numerical-languages/
- 'Fortran':
- OpenCoarrays: salomon/software/numerical-languages/
- Operating System: salomon/software/
- Anselm Software:
- Available Modules:
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment