OpenMP vs trigonmetric transforms

OpenMP vs trigonmetric transforms

Hi, 

 

!$omp parallel   private(iy,iz,bx,by) shared(ll,nn,mm,b,im1)
!$omp do
  do iz = 1,ll
    do iy = 1,nn
  
        bx=.0
        by=.0
        !fetch data 
        bx(1:mm)=real(b(1:mm,iy,iz,ik))       
        by(1:mm)= aimag(b(1:mm,iy,iz,ik))  
        !compute forward sine transform
        call dst_fwd(im1,bx(0:im1))          
        call dst_fwd(im1,by(0:im1))

       !set data
        b(1:mm,iy,iz,ik) =cmplx( bx(1:mm),by(1:mm))
    enddo
  enddo
!$omp end do
!$omp end parallel

 

The code snippet above crashes on runtime with no error message when I put OpenMP commands around the the two loops in order to compute Sine Transforms (subroutine dst_fwd ) in parallel. 

Does anyone could give me a hint what may cause the program to crash? 

 

Thanks,

Nicolas

 

 

 

2 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

I suspect we would need a full working example.

Leave a Comment

Please sign in to add a comment. Not a member? Join today