| Exec Time (%) | Codelet Name | Error (%) |
|---|---|---|
| 99.8 | shell_lam_shell__159 | 0.1 |
Exec Time : 1.657044e+10 cycles
Matching : 99.8%
__invivo__block_solver_mat_times_vec__167
Real time: 8.521994e+09 -- Predicted time: 1.107973e+10
Callcount: 851
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 393 | ### 1 | 404.635172490719 | 1.309768e+07 | 1.022026e+07 | 21.9 |
| 793 | ### 2 | 399.496157122142 | 1.281759e+07 | 9.627160e+06 | 24.8 |
| 564 | ### 3 | 47.7306514716025 | 1.381433e+07 | 1.160517e+07 | 15.9 |
__invivo__jacobian_lam_jacobian__30
Real time: 4.851621e+09 -- Predicted time: 4.818960e+09
Callcount: 75
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 67 | ### 1 | 58.5640267331168 | 1.667783e+07 | 1.707306e+07 | 2.3 |
| 8 | ### 2 | 3.00225848911494 | 3.147402e+08 | 3.153064e+08 | 0.1 |
| 11 | ### 3 | 3.00331428049316 | 3.441076e+08 | 3.443595e+08 | 0.0 |
| 14 | ### 4 | 3.00234592734919 | 2.707702e+08 | 2.710894e+08 | 0.1 |
| 17 | ### 5 | 3.0018995683394 | 2.113549e+08 | 2.121556e+08 | 0.3 |
| 20 | ### 6 | 3.00521274160711 | 9.259428e+07 | 9.296912e+07 | 0.4 |
| 23 | ### 7 | 3.0141263061733 | 4.583992e+07 | 4.619410e+07 | 0.7 |
__invivo__block_solver_bi_cgstab_block__52
Real time: 9.346048e+09 -- Predicted time: 9.781475e+09
Callcount: 25
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 3.09020560490978 | 8.348555e+07 | 9.862624e+07 | 15.3 |
| 5 | ### 2 | 4.288745255644 | 1.926579e+08 | 1.601882e+08 | 16.8 |
| 9 | ### 3 | 2.99999967970211 | 3.244180e+08 | 2.747442e+08 | 15.3 |
| 12 | ### 4 | 4.04350372155428 | 3.677495e+08 | 3.687185e+08 | 0.2 |
| 15 | ### 5 | 4.14890823579954 | 4.689913e+08 | 4.584233e+08 | 2.2 |
| 19 | ### 6 | 2.99843673850849 | 5.739560e+08 | 5.731351e+08 | 0.1 |
| 20 | ### 7 | 1.99802514917976 | 6.656611e+08 | 6.657009e+08 | 0.0 |
| 22 | ### 8 | 1.99711849508498 | 6.209930e+08 | 6.197234e+08 | 0.2 |
__invivo__flux_lam_flux__58
Real time: 2.336126e+08 -- Predicted time: 2.326995e+08
Callcount: 25
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 1 | 8.372797e+06 | 1.017800e+07 | 17.7 |
| 13 | ### 2 | 23.9173139456854 | 9.379260e+06 | 9.352484e+06 | 0.2 |
x
1
subroutine shell(Re,Pr,nx,ny,nz,2
$nuim,nuex2,nuex4,cfl,scheme,conf,ni,maxit) 3
implicit none4
COMMON /scale/ ro1,p1,ro2,p25
integer nx,ny,nz,ni,n,i,j,k,scheme,conf,si,sj,sk,l,m6
integer im2,im1,ip1,ip2,jm2,jm1,jp1,jp2,km2,km1,kp1,kp27
integer maxit8
real*8 nuim,nuex2,nuex4,cfl,epsilon,dx,dy,dz,dt,rad9
real*8 gm,u1,v1,w1,p1,ro1,u2,v2,w2,p2,ro2,Re,Pr,cfll10
$,time,t1,t2,t3,t4,t5,t6,t7,t8,mu11
real*8 dqnorm12
13
real*8 q(5,nx,ny,nz),dq(5,nx,ny,nz),rhs(5,nx,ny,nz),e(5,nx,ny,nz),14
1 f(5,nx,ny,nz),g(5,nx,ny,nz),ev(5,nx,ny,nz),fv(5,nx,ny,nz),15
2 gv(5,nx,ny,nz),diss(5,nx,ny,nz)16
17
real*8 a(5,5,nx,ny,nz),be(5,5,nx,ny,nz),ce(5,5,nx,ny,nz),18
1 av(5,5,nx,ny,nz),bv(5,5,nx,ny,nz),cv(5,5,nx,ny,nz)19
20
real*8 ae(5,5,nx,ny,nz),axp(5,5,nx,ny,nz),ayp(5,5,nx,ny,nz),21
1 azp(5,5,nx,ny,nz),axm(5,5,nx,ny,nz),aym(5,5,nx,ny,nz),22
2 azm(5,5,nx,ny,nz),ident(5,5,nx,ny,nz)23
24
25
real*8 u(nx,ny,nz),v(nx,ny,nz),w(nx,ny,nz),p(nx,ny,nz),26
1 ro(nx,ny,nz),at(nx,ny,nz)27
28
real*8 xyz(3,nx,ny,nz)29
30
31
C This particular problem is periodic only32
33
34
35
36
C Data initialization37
epsilon=1.d-1638
39
si=nx/840
sj=ny/841
sk=nz/842
dx=1.0d0/(nx-1)43
dy=1.0d0/(ny-1)44
dz=1.0d0/(nz-1)45
time=0.0d046
47
48
49
do k = 1,nz50
t3=(k-1)*dz51
do j=1,ny52
t2=(j-1)*dy53
do i=1,nx54
t1=(i-1)*dx55
xyz(1,i,j,k)=t156
xyz(2,i,j,k)=t257
xyz(3,i,j,k)=t358
do l=1,559
rhs(l,i,j,k)=0.0d060
do m=1,561
a(m,l,i,j,k)=0.0d062
axp(m,l,i,j,k)=0.0d063
ayp(m,l,i,j,k)=0.0d064
azp(m,l,i,j,k)=0.0d065
axm(m,l,i,j,k)=0.0d066
aym(m,l,i,j,k)=0.0d067
azm(m,l,i,j,k)=0.0d068
ident(m,l,i,j,k)=0.0d069
enddo70
ident(l,l,i,j,k)=1.0d071
rhs(l,i,j,k)=0.0d072
dq(l,i,j,k)=0.0d073
enddo74
enddo75
enddo76
enddo 77
78
C Flow Initialization79
80
gm=1.4d081
u1=0.0d082
v1=0.0d083
w1=0.0d084
p1=1.0d085
ro1=1.0d086
u2=0.0d087
v2=0.0d088
w2=0.0d089
p2=0.1d090
ro2=0.1d091
92
t1=ro293
t2=ro2*u294
t3=ro2*v295
t4=ro2*w296
t5=p2/(gm-1.0d0)+0.5d0*ro2*(u2**2+v2**2+w2**2)97
98
99
do k = 1,nz100
do j=1,ny101
do i=1,nx102
q(1,i,j,k)=t1103
q(2,i,j,k)=t2104
q(3,i,j,k)=t3105
q(4,i,j,k)=t4106
q(5,i,j,k)=t5107
enddo108
enddo109
enddo110
111
t1=ro1112
t2=ro1*u1113
t3=ro1*v1114
t4=ro1*w1115
t5=p1/(gm-1.0d0)+0.5d0*ro1*(u1**2+v1**2+w1**2)116
117
if (conf.EQ.0) then118
119
do k=nz/2-sk,nz/2+sk120
do j=ny/2-sj,ny/2+sj121
do i=nx/2-si,nx/2+si122
q(1,i,j,k)=t1123
q(2,i,j,k)=t2124
q(3,i,j,k)=t3125
q(4,i,j,k)=t4126
q(5,i,j,k)=t5127
enddo128
enddo129
enddo130
else131
rad=DMIN1((si*dx)**2,(sj*dy)**2,(sk*dz)**2)132
133
do k=1,nz134
t8=((k-nz/2)*dz)**2135
do j=1,ny136
t7=((j-ny/2)*dy)**2137
do i=1,nx138
t6=((i-nx/2)*dx)**2 +t7 + t8139
if (t6.LE.rad) then140
q(1,i,j,k)=t1141
q(2,i,j,k)=t2142
q(3,i,j,k)=t3143
q(4,i,j,k)=t4144
q(5,i,j,k)=t5145
endif146
enddo147
enddo148
enddo149
endif150
151
152
153
C Propagation in time154
155
156
157
time = 0.0d0158
159
do n=1,ni160
C Time step definition161
162
cfll=0.1d0+(n-1.0d0)*cfl/20.0d0163
if (cfll.ge.cfl) cfll=cfl164
t8=0.0d0165
166
do k=1,nz167
do j=1,ny168
do i=1,nx169
t1=q(1,i,j,k)170
t2=q(2,i,j,k)/t1171
t3=q(3,i,j,k)/t1172
t4=q(4,i,j,k)/t1173
t5=(gm-1.0d0)*(q(5,i,j,k)-0.5d0*t1*(t2*t2+t3*t3+t4*t4))174
t6=dSQRT(gm*t5/t1)175
mu=gm*Pr*(gm*t5/t1)**0.75d0*2.0d0/Re/t1176
t7=((dabs(t2)+t6)/dx+mu/dx**2)**2 +177
1 ((dabs(t3)+t6)/dy+mu/dy**2)**2 +178
2 ((dabs(t4)+t6)/dz+mu/dz**2)**2179
t7=DSQRT(t7)180
t8=max(t8,t7)181
enddo182
enddo183
enddo184
dt=cfll / t8185
186
187
C Left hand side188
call jacobian(q,ae,av,0.0d0,1.0d0,0.0d0,0.0d0,Re,Pr,gm,189
1 nx,ny,nz,dx,1)190
call jacobian(q,be,bv,0.0d0,0.0d0,1.0d0,0.0d0,Re,Pr,gm,191
1 nx,ny,nz,dy,2)192
call jacobian(q,ce,cv,0.0d0,0.0d0,0.0d0,1.0d0,Re,Pr,gm,193
1 nx,ny,nz,dz,3)194
195
196
do k=1,nz197
km1=mod(k-2+nz,nz)+1198
kp1=mod(k,nz)+1199
do j=1,ny200
jm1=mod(j-2+ny,ny)+1201
jp1=mod(j,ny)+1202
do i=1,nx203
im1=mod(i-2+nx,nx)+1204
ip1=mod(i,nx)+1205
do l=1,5206
do m=1,5207
t1=ident(m,l,i,j,k)208
a(m,l,i,j,k)=t1 - 0.5d0*dt*(209
1 (av(m,l,i,j,k)-av(m,l,im1,j,k))/dx+210
2 (bv(m,l,i,j,k)-bv(m,l,i,jm1,k))/dy+211
3 (cv(m,l,i,j,k)-cv(m,l,i,j,km1))/dz)/Re212
4 + 2.0d0*nuim*dt*(1.0d0/dx+1.0d0/dy+1.0d0/dz)*t1 213
214
axp(m,l,i,j,k)=0.5d0*dt/dx*215
1 (ae(m,l,ip1,j,k)-av(m,l,i,j,k)/Re) 216
2 - nuim*dt/dx * t1217
axm(m,l,i,j,k)=-0.5d0*dt/dx*218
1 (ae(m,l,im1,j,k)-av(m,l,im1,j,k)/Re)219
2 - nuim*dt/dx * t1 220
221
ayp(m,l,i,j,k)=0.5d0*dt/dy*222
1 (be(m,l,i,jp1,k)-bv(m,l,i,j,k)/Re)223
2 - nuim*dt/dy * t1224
aym(m,l,i,j,k)=-0.5d0*dt/dy*225
1 (be(m,l,i,jm1,k)-bv(m,l,i,jm1,k)/Re)226
2 - nuim*dt/dy * t1227
228
azp(m,l,i,j,k)=0.5d0*dt/dz*229
1 (ce(m,l,i,j,kp1)-cv(m,l,i,j,k)/Re)230
2 - nuim*dt/dz * t1231
azm(m,l,i,j,k)=-0.5d0*dt/dz*232
1 (ce(m,l,i,j,km1)-cv(m,l,i,j,km1)/Re)233
2 - nuim*dt/dz * t1234
235
enddo236
enddo237
enddo238
enddo239
enddo240
241
242
243
C Right hand side244
call flux(q,e,f,g,ev,fv,gv,Re,Pr,gm,nx,ny,nz,dx,dy,dz)245
246
247
do k=1,nz248
km2=mod(k+nz-3,nz)+1249
km1=mod(k+nz-2,nz)+1250
kp1=mod(k,nz)+1251
kp2=mod(k+1,nz)+1252
do j=1,ny253
jm2=mod(j+ny-3,ny)+1254
jm1=mod(j+ny-2,ny)+1255
jp1=mod(j,ny)+1256
jp2=mod(j+1,ny)+1257
do i=1,nx258
im2=mod(i+nx-3,nx)+1259
im1=mod(i+nx-2,nx)+1260
ip1=mod(i,nx)+1261
ip2=mod(i+1,nx)+1262
do l=1,5263
t1= -0.5d0*dt*(264
1 (e(l,ip1,j,k)-e(l,im1,j,k))/dx +265
2 (f(l,i,jp1,k)-f(l,i,jm1,k))/dy +266
3 (g(l,i,j,kp1)-g(l,i,j,km1))/dz) +267
4 dt/Re*((ev(l,i,j,k)-ev(l,im1,j,k))/dx +268
5 (fv(l,i,j,k)-fv(l,i,jm1,k))/dy +269
6 (gv(l,i,j,k)-gv(l,i,j,km1))/dz)270
271
C Artificial Viscosity272
C Explicit dissipation - second order273
t2=dt*nuex2*(274
1 (q(l,ip1,j,k)-2.0d0*q(l,i,j,k)+q(l,im1,j,k))/dx +275
2 (q(l,i,jp1,k)-2.0d0*q(l,i,j,k)+q(l,i,jm1,k))/dy +276
3 (q(l,i,j,kp1)-2.0d0*q(l,i,j,k)+q(l,i,j,km1))/dz)277
278
C Explicit dissipation - fourth order279
t2=t2-dt*nuex4*(280
1 (q(l,ip2,j,k)-4.0d0*q(l,ip1,j,k)+6.0d0*q(l,i,j,k)-281
2 4.0d0*q(l,im1,j,k)+q(l,im2,j,k))/dx +282
3 (q(l,i,jp2,k)-4.0d0*q(l,i,jp1,k)+6.0d0*q(l,i,j,k)-283
4 4.0d0*q(l,i,jm1,k)+q(l,i,jm2,k))/dy +284
5 (q(l,i,j,kp2)-4.0d0*q(l,i,j,kp1)+6.0d0*q(l,i,j,k)-285
6 4.0d0*q(l,i,j,km1)+q(l,i,j,km2))/dx )286
287
rhs(l,i,j,k)=t1+t2288
C initial guess is rhs289
dq(l,i,j,k)=rhs(l,i,j,k)290
enddo291
enddo292
enddo293
enddo294
295
c write(10,*) 'Time step: ',n,' dt: ',dt 296
297
if (scheme.eq.1) call bi_cgstab_block(dq,rhs,a,axp,298
$ ayp,azp,axm,aym,azm,epsilon,5,nx,ny,nz,maxit)299
300
301
do k=1,nz302
do j=1,ny303
do i=1,nx304
do l=1,5305
q(l,i,j,k)=q(l,i,j,k)+dq(l,i,j,k)306
enddo307
enddo308
enddo309
enddo310
311
c write(10,*) n, q(1,2,2,2)312
c write(10,*) q(2,2,2,2),q(3,2,2,2),q(4,2,2,2)313
c write(10,*) q(5,2,2,2)314
315
316
317
time=time+dt318
enddo319
320
dqnorm = 0.0d0321
do k=1,nz322
do j=1,ny323
do i=1,nx324
do l=1,5325
dqnorm = dqnorm + dq(l,i,j,k)*dq(l,i,j,k)326
enddo327
enddo328
enddo329
enddo330
331
write(30,*) 'dqnorm ==', dqnorm 332
333
return334
endOne invocation only
__invivo__shell_lam_shell__159
Real time: 1.654916e+10 -- Predicted time: 1.651822e+10
Callcount: 1
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 1 | 1.651822e+10 | 1.659105e+10 | 0.4 |