Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Feature/issue 49 openmp direct #51

Merged
merged 3 commits into from
Mar 19, 2024
Merged

Feature/issue 49 openmp direct #51

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
d6aab33
added openmp directives to helmkernels.f
oneilm Mar 18, 2024
bff5d85
added openmp directives on helmholtz direct kernels
oneilm Mar 19, 2024
96d6497
cleaned up testing code
oneilm Mar 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
21 changes: 21 additions & 0 deletions src/Helmholtz/helmkernels.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -90,6 +90,8 @@ subroutine h3ddirectcp(nd,zk,sources,charge,ns,ztarg,nt,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, zdiff, j, dd, d, ztmp, idim)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand All @@ -107,6 +109,7 @@ subroutine h3ddirectcp(nd,zk,sources,charge,ns,ztarg,nt,
1000 continue
enddo
enddo
c$omp end parallel do


return
Expand Down Expand Up @@ -183,6 +186,9 @@ subroutine h3ddirectcg(nd,zk,sources,charge,ns,ztarg,nt,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, j, zdiff, dd, d, cd, cd1, ztmp1)
c$omp$ private(ztmp2, ztmp3, idim)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand All @@ -206,6 +212,7 @@ subroutine h3ddirectcg(nd,zk,sources,charge,ns,ztarg,nt,
1000 continue
enddo
enddo
c$omp end parallel do


return
Expand Down Expand Up @@ -285,6 +292,8 @@ subroutine h3ddirectdp(nd,zk,sources,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, j, zdiff, dd, d, dinv, cd, cd1, idim, dotprod)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand All @@ -309,6 +318,7 @@ subroutine h3ddirectdp(nd,zk,sources,
1000 continue
enddo
enddo
c$omp end parallel do


return
Expand Down Expand Up @@ -394,6 +404,9 @@ subroutine h3ddirectdg(nd,zk,sources,dipvec,ns,ztarg,nt,pot,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, j, zdiff, dd, d, dinv, dinv2, cd, cd2)
c$omp$ private(cd3, idim, dotprod, cd4)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand Down Expand Up @@ -428,6 +441,7 @@ subroutine h3ddirectdg(nd,zk,sources,dipvec,ns,ztarg,nt,pot,
1000 continue
enddo
enddo
c$omp end parallel do


return
Expand Down Expand Up @@ -506,6 +520,8 @@ subroutine h3ddirectcdp(nd,zk,sources,charge,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, j, zdiff, dd, d, dinv, cd, cd1, idim, dotprod)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand All @@ -532,6 +548,7 @@ subroutine h3ddirectcdp(nd,zk,sources,charge,
1000 continue
enddo
enddo
c$omp end parallel do



Expand Down Expand Up @@ -619,6 +636,9 @@ subroutine h3ddirectcdg(nd,zk,sources,charge,

zkeye = zk*eye

c$omp parallel do default(shared)
c$omp$ private(i, j, zdiff, dd, d, dinv, dinv2, cd, cd2, cd3)
c$omp$ private(idim, dotprod, cd4)
do i=1,nt
do j=1,ns
zdiff(1) = ztarg(1,i)-sources(1,j)
Expand Down Expand Up @@ -657,6 +677,7 @@ subroutine h3ddirectcdg(nd,zk,sources,charge,
1000 continue
enddo
enddo
c$omp end parallel do


return
Expand Down
Loading
Loading