I have a subroutine as below: subroutine lm(n, x, y, res) integer(8) :: n real(8) :: x(n), y(n), res real(8) :: sum_x, sum_y, sum_xx, sum_xy integer(8) :: i sum_x = 0d0 sum_y = 0d0 sum_xx = 0d0 sum_xy = 0d0 !$omp parallel do reduction(+:sum_x) do i = 1, n sum_x = sum_…

[image] fortran4r: … the subroutine will be much less verbose. res = (n*sum(x*y) - sum(x)*sum(y))/(n*sum(x**2) - sum(x)**2) Edit: If you still want to do it the verbose way, you might consider merging the loops. !$omp parallel do reduction(+:sum_x,sum_y,sum_xx,sum_xy) do i = 1, n su…

Parallel sum function

Help

fortran4r October 15, 2021, 1:42pm 4

Can co_sum be used in shared libs?

Topic		Replies	Views
Coarray nested loop Help	3	542	November 9, 2021
OpenMP question: private vs shared work arrays for reduction Help	8	2428	February 23, 2021
Learning coarrays, collective subroutines and other parallel features of Modern Fortran Help	48	2432	May 11, 2021
Reduction specifier for do concurrent Help	0	663	May 4, 2022
OpenMP reduction with max/min function Help	1	2963	January 2, 2023

Parallel sum function

Related topics