PGI User Forum
 SearchSearch   MemberlistMemberlist     RegisterRegister   ProfileProfile    Log inLog in 

CUDA-x86.

WRF on Windows
Goto page Previous  1, 2
 
Post new topic   Reply to topic    PGI User Forum Forum Index -> Accelerator Programming
View previous topic :: View next topic  
Author Message
TheMatt



Joined: 06 Jul 2009
Posts: 317
Location: Greenbelt, MD

PostPosted: Fri Oct 23, 2009 10:02 am    Post subject: Reply with quote

mkcolg wrote:
Hi Matt,

I asked Craig about the BLAS inlining issue you mentioned but he says he didn't perform any. Can you give more details?

Thanks,
Mat

Well, it was BLAS-type, quoth the article:
Quote:
...and two BLAS-type subroutine calls that we manually inlined:
Code:
    do k = kts, kte
       call vsrec( tvec1(its), den(its,k), ite-its+1 )
       do i = its, ite
          tvec1(i) = tvec1(i) * den0
       enddo
       call vssqrt( denfac(its,k), tvec1(its), ite-its+1 )
    enddo

Routine vsrec computes the reciprocal of each vector element, and vssqrt computes the square root, so this is equivalent to:
Code:
    do k = kts, kte
       do i = its, ite
          denfac(i,k) = den0 / den(i,k)
       enddo
    enddo


It just looked to my eye that the sqrt call never made it into the code snippet in the article:
Code:
denfac(i,k) = sqrt(den0 / den(i,k))
Back to top
View user's profile
mkcolg



Joined: 30 Jun 2004
Posts: 6129
Location: The Portland Group Inc.

PostPosted: Fri Oct 23, 2009 10:21 am    Post subject: Reply with quote

Good eye! It's just typo. The actual code does have the sqrt. We'll get the article updated.

Thanks,
Mat
Back to top
View user's profile
Display posts from previous:   
Post new topic   Reply to topic    PGI User Forum Forum Index -> Accelerator Programming All times are GMT - 7 Hours
Goto page Previous  1, 2
Page 2 of 2

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © phpBB Group