Parallel Computing with OpenMP

This example demonstrates using OpenMP for parallel computing.

fortranCopy codeprogram parallel_sum
    use omp_lib
    implicit none
    integer, parameter :: n = 1000
    integer :: i
    real :: sum, local_sum
    real, dimension(n) :: array

    ! Initialize the array
    array = 1.0

    ! Initialize sum
    sum = 0.0

    ! Parallel region
    !$omp parallel private(local_sum, i) reduction(+:sum)
    local_sum = 0.0
    !$omp do
    do i = 1, n
        local_sum = local_sum + array(i)
    end do
    !$omp end do
    sum = sum + local_sum
    !$omp end parallel

    print *, 'Sum of array: ', sum
end program parallel_sum