154
edits
Changes
→Ubuntu 12.04 LTS and CUDA 5 Toolkit Installation Guide
==Ubuntu 12.04 LTS and CUDA 5 Toolkit Installation Guide==
[http://zenit.senecac.on.ca/wiki/index.php/GPU610/DPS915_Ubuntu_and_CUDA_Installation See the guide here; work in progress]
= Fortran to C =
Sample code from the TOMO project
== Original Fortran Subroutine ==
<pre>
SUBROUTINE longtrack_self(direction,nrep,yp,xp,turnnow)
!-------------------------------------------------------------------------
! h: principal harmonic number
! eta0: phase slip factor
! E0: energy of synchronous particle m
! beta0: relativistic beta of synchronous particle
! phi0: synchronous phase
! q: charge state of particles
! dphi: phase difference between considered particle and synchronous one
! denergy: energy difference between considered particle and synchronous one
! nrep: pass cavity nrep times before returning data
! direction: to inverse the time advance (rotation in the bucket), 1 or -1
! xp and yp: time and energy in pixels
! dtbin and dEbin: GLOBAL time and energy pixel size in s and MeV
! omegarev0: revolution frequency
! VRF1,VRF2,VRF1dot,VRF2dot: GLOBAL RF voltages and derivatives of volts
! turnnow: present turn
!---------------------------------------------------------------------------
IMPLICIT NONE
REAL(SP), DIMENSION(:), INTENT(INOUT) :: xp,yp
REAL(SP), DIMENSION(SIZE(xp)) :: dphi,denergy,selfvolt
!HPF$ distribute dphi(block)
!HPF$ align with dphi :: denergy,selfvolt,xp
INTEGER :: mm
INTEGER :: i,p,nrep,direction,turnnow
dphi=(xp+xorigin)*h*omegarev0(turnnow)*dtbin-phi0(turnnow)
denergy=(yp-yat0)*dEbin
IF (direction.GT.0) THEN
p=turnnow/dturns+1
DO i=1,nrep
forall(mm=1:size(xp)) dphi(mm)=dphi(mm)-c1(turnnow)*denergy(mm)
turnnow=turnnow+1
forall(mm=1:size(xp)) xp(mm)=dphi(mm)+phi0(turnnow)-&
xorigin*h*omegarev0(turnnow)*dtbin
forall(mm=1:size(xp)) xp(mm)=(xp(mm)-&
phiwrap*FLOOR(xp(mm)/phiwrap))/(h*omegarev0(turnnow)*dtbin)
forall(mm=1:size(xp)) selfvolt(mm)=vself(p,FLOOR(xp(mm))+1)
forall(mm=1:size(xp)) denergy(mm)=denergy(mm)+q*((&
(VRF1+VRF1dot*tatturn(turnnow))*SIN(dphi(mm)+phi0(turnnow))+&
(VRF2+VRF2dot*tatturn(turnnow))*&
SIN(hratio*(dphi(mm)+phi0(turnnow)-phi12)))+selfvolt(mm))-c2(turnnow)
END DO
ELSE
p=turnnow/dturns
DO i=1,nrep
forall(mm=1:size(xp)) selfvolt(mm)=vself(p,FLOOR(xp(mm))+1)
forall(mm=1:size(xp)) denergy(mm)=denergy(mm)-q*((&
(VRF1+VRF1dot*tatturn(turnnow))*SIN(dphi(mm)+phi0(turnnow))+&
(VRF2+VRF2dot*tatturn(turnnow))*&
SIN(hratio*(dphi(mm)+phi0(turnnow)-phi12)))+selfvolt(mm))+c2(turnnow)
turnnow=turnnow-1
forall(mm=1:size(xp)) dphi(mm)=dphi(mm)+c1(turnnow)*denergy(mm)
forall(mm=1:size(xp)) xp(mm)=dphi(mm)+phi0(turnnow)-&
xorigin*h*omegarev0(turnnow)*dtbin
forall(mm=1:size(xp)) xp(mm)=(xp(mm)-&
phiwrap*FLOOR(xp(mm)/phiwrap))/(h*omegarev0(turnnow)*dtbin)
END DO
END IF
yp=denergy/dEbin+yat0
END SUBROUTINE longtrack_self
</pre>
== Modified Fortran Subroutine ==
<pre>
SUBROUTINE longtrack_self(direction,nrep,yp,xp,turnnow)
!-------------------------------------------------------------------------
! h: principal harmonic number
! eta0: phase slip factor
! E0: energy of synchronous particle
! beta0: relativistic beta of synchronous particle
! phi0: synchronous phase
! q: charge state of particles
! dphi: phase difference between considered particle and synchronous one
! denergy: energy difference between considered particle and synchronous one
! nrep: pass cavity nrep times before returning data
! direction: to inverse the time advance (rotation in the bucket), 1 or -1
! xp and yp: time and energy in pixels
! dtbin and dEbin: GLOBAL time and energy pixel size in s and MeV
! omegarev0: revolution frequency
! VRF1,VRF2,VRF1dot,VRF2dot: GLOBAL RF voltages and derivatives of volts
! turnnow: present turn
!---------------------------------------------------------------------------
IMPLICIT NONE
REAL(SP), DIMENSION(:), INTENT(INOUT) :: xp,yp
REAL(SP), DIMENSION(SIZE(xp)) :: dphi,denergy,selfvolt
!HPF$ distribute dphi(block)
!HPF$ align with dphi :: denergy,selfvolt,xp
INTEGER :: mm
INTEGER :: i,p,nrep,direction,turnnow
CALL gputrack_self(direction,nrep,yp,xp,turnnow, &
SIZE(xp),dphi,denergy, &
c1, &
c2, &
dEbin, &
dtbin, &
h, &
hratio, &
omegarev0, &
phi0, &
phi12, &
q, &
tatturn, &
VRF1, &
VRF1dot, &
VRF2, &
VRF2dot, &
xorigin, &
yat0, &
p, &
dturns, &
phiwrap, &
selfvolt, &
profilecount-1, &
wraplength, &
vself )
END SUBROUTINE longtrack_self
</pre>
== New C Function ==
<pre>
#include <stdio.h>
#include <math.h>
void gputrack_self_ ( \
int *direction, \
int *nrep, \
float *yp, \
float *xp, \
int *turnnow, \
int *sizeofarrays, \
float *dphi, \
float *denergy, \
float *c1, \
float *c2, \
float *dEbin, \
float *dtbin, \
float *h, \
float *hratio, \
float *omegarev0, \
float *phi0, \
float *phi12, \
float *q, \
float *tatturn, \
float *VRF1, \
float *VRF1dot, \
float *VRF2, \
float *VRF2dot, \
float *xorigin, \
float *yat0, \
int *p, \
int *dturns, \
float *phiwrap, \
float *selfvolt, \
int *vselfDimRow, \
int *vselfDimCol, \
float *vself \
)
{
/* Local Variables */
int l,i,mm,t;
l = *sizeofarrays;
t = *turnnow;
// longtrack_self specific local variables
int cp;
cp = *p;
/* dphi=(xp+xorigin)*h*omegarev0(turnnow)*dtbin-phi0(turnnow) */
for(mm = 0; mm < l; mm++) {
dphi[mm] = (xp[mm] + *xorigin) * *h * omegarev0[t] * *dtbin - phi0[t];
}
/* denergy=(yp-yat0)*dEbin */
for(mm = 0; mm < l; mm++) {
denergy[mm] = (yp[mm] - *yat0) * *dEbin;
}
/* IF (direction.GT.0) THEN */
if (*direction > 0) {
/* p=turnnow/dturns+1 */
cp = t / *dturns + 1;
/* DO i=1,nrep */
for(i = 1; i <= *nrep; i++) {
/* forall(mm=1:size(xp)) dphi(mm)=dphi(mm)-c1(turnnow)*denergy(mm) */
for(mm=0;mm<l;mm++) {
dphi[mm] = dphi[mm] - c1[t] *denergy[mm];
}
/* turnnow=turnnow+1 */
t=t+1;
/* forall(mm=1:size(xp)) xp(mm)=dphi(mm)+phi0(turnnow)-&
xorigin*h*omegarev0(turnnow)*dtbin */
for(mm=0;mm<l;mm++) {
xp[mm] = dphi[mm] + phi0[t] - \
*xorigin * *h * omegarev0[t] * *dtbin;
}
/* forall(mm=1:size(xp)) xp(mm)=(xp(mm)-&
phiwrap*FLOOR(xp(mm)/phiwrap))/(h*omegarev0(turnnow)*dtbin) */
for(mm = 0; mm < l; mm++) {
xp[mm] = (xp[mm] - \
*phiwrap * floor(xp[mm] / *phiwrap)) / (*h * omegarev0[t] * *dtbin);
}
/* forall(mm=1:size(xp)) selfvolt(mm)=vself(p,FLOOR(xp(mm))+1) */
for(mm = 0; mm < l; mm++) {
int itemp = floor(xp[mm]);
selfvolt[mm] = vself[(*vselfDimRow * (itemp)) + (cp-1)];
}
/* forall(mm=1:size(xp)) denergy(mm)=denergy(mm)+q*((&
(VRF1+VRF1dot*tatturn(turnnow))*SIN(dphi(mm)+phi0(turnnow))+&
(VRF2+VRF2dot*tatturn(turnnow))*&
SIN(hratio*(dphi(mm)+phi0(turnnow)-phi12)))+selfvolt(mm))-c2(turnnow) */
for(mm = 0; mm < l; mm++) {
denergy[mm] = denergy[mm] + *q *(( \
(*VRF1 + *VRF1dot * tatturn[t]) * sin(dphi[mm] + phi0[t]) + \
(*VRF2 + *VRF2dot * tatturn[t]) * \
sin(*hratio * (dphi[mm] + phi0[t] - *phi12))) + selfvolt[mm]) -c2[t];
}
/* END DO */
}
}
else {
// p=turnnow/dturns
cp = t / *dturns;
// DO i=1,nrep
for (i=1;i<=*nrep;i++) {
// forall(mm=1:size(xp)) selfvolt(mm)=vself(p,FLOOR(xp(mm))+1)
for(mm = 0; mm < l; mm++) {
int itemp = (int)floor(xp[mm]);
selfvolt[mm] = vself[(*vselfDimRow*(itemp)) + (cp-1)];
}
/* forall(mm=1:size(xp)) denergy(mm)=denergy(mm)-q*((&
(VRF1+VRF1dot*tatturn(turnnow))*SIN(dphi(mm)+phi0(turnnow))+&
(VRF2+VRF2dot*tatturn(turnnow))*&
SIN(hratio*(dphi(mm)+phi0(turnnow)-phi12)))+selfvolt(mm))+c2(turnnow) */
for(mm = 0; mm < l; mm++) {
denergy[mm]=denergy[mm] - *q *(( \
(*VRF1 + *VRF1dot * tatturn[t]) *sin(dphi[mm] + phi0[t]) + \
(*VRF2 + *VRF2dot * tatturn[t]) * \
sin(*hratio * (dphi[mm] + phi0[t] - *phi12))) + selfvolt[mm]) + c2[t];
}
// turnnow=turnnow-1
t--;
/* forall(mm=1:size(xp)) dphi(mm)=dphi(mm)-c1(turnnow)*denergy(mm) */
for(mm = 0; mm < l; mm++) {
dphi[mm]=dphi[mm] + c1[t] * denergy[mm];
}
/* forall(mm=1:size(xp)) xp(mm)=dphi(mm)+phi0(turnnow)-&
xorigin*h*omegarev0(turnnow)*dtbin */
for(mm = 0; mm < l; mm++) {
xp[mm] = dphi[mm] + phi0[t] - \
*xorigin * *h * omegarev0[t] * *dtbin;
}
/* forall(mm=1:size(xp)) xp(mm)=(xp(mm)-&
phiwrap*FLOOR(xp(mm)/phiwrap))/(h*omegarev0(turnnow)*dtbin) */
for(mm = 0; mm < l; mm++) {
xp[mm] = (xp[mm] - \
*phiwrap * floor(xp[mm] / *phiwrap)) / (*h * omegarev0[t] * *dtbin);
}
}
}
// yp=denergy/dEbin+yat0
for(mm=0; mm<l; mm++) {
yp[mm] = denergy[mm] / *dEbin + *yat0;
}
*turnnow = t;
return;
}
</pre>
