diff options
| author | Tor Aamodt <[email protected]> | 2010-10-01 08:55:28 -0800 |
|---|---|---|
| committer | Tor Aamodt <[email protected]> | 2010-10-01 08:55:28 -0800 |
| commit | 11b308e7363e937966b035b4891db32b4eece3bf (patch) | |
| tree | 50ca4c9ad6f163ac4acb2bf505e64dfebed66947 /benchmarks/CUDA/DG/buildsource.m | |
| parent | bb820c116764d7a1b8e071137d32b74e7f34dd2f (diff) | |
integrating recent changes from fermi-test into fermi
(i'll use "fermi" for more disruptive changes to the pipeline model such
as updating the MSHRs and getting rid of the warp tracker, ripping out DWF, etc...)
[git-p4: depot-paths = "//depot/gpgpu_sim_research/fermi/distribution/": change = 7805]
Diffstat (limited to 'benchmarks/CUDA/DG/buildsource.m')
| -rw-r--r-- | benchmarks/CUDA/DG/buildsource.m | 106 |
1 files changed, 0 insertions, 106 deletions
diff --git a/benchmarks/CUDA/DG/buildsource.m b/benchmarks/CUDA/DG/buildsource.m deleted file mode 100644 index c93d363..0000000 --- a/benchmarks/CUDA/DG/buildsource.m +++ /dev/null @@ -1,106 +0,0 @@ - -function buildsource(in_N) - -Globals2D; - - -NODETOL = 1e-8; - -N = in_N; - -Np = (N+1)*(N+2)/2; Nfp = N+1; Nfaces = 3; - -[x,y] = Nodes2D(N); -[r,s] = xytors(x,y); - -% find all the nodes that lie on each edge -fmask1 = find( abs(s+1) < NODETOL)'; -fmask2 = find( abs(r+s) < NODETOL)'; -fmask3 = find( abs(r+1) < NODETOL)'; -Fmask = [fmask1;fmask2;fmask3]'; - -% vandermonde matrix -V = Vandermonde2D(N, r, s); - -% derivative matrices -[Dr, Ds] = Dmatrices2D(N, r, s, V); - -% lift matrix -LIFT = Lift2D(); - -% correct to 0-index -Fmask = Fmask-1; - - -fid = fopen(sprintf('dataN%02d.h', N), 'w'); - -fprintf(fid, sprintf('#ifndef DATAN%02d \n\n', N)); -fprintf(fid, sprintf('#define DATAN%02d 1\n\n', N)); -fprintf(fid, 'double p_r[%d] = { ', Np); -fprintf(fid, '%17.15g ', r(1)); -for n=2:Np - fprintf(fid, ', %17.15g ', r(n)); -end -fprintf(fid, '};\n') - -fprintf(fid, 'double p_s[%d] = {', Np); -fprintf(fid, '%17.15g ', s(1)); -for n=2:Np - fprintf(fid, ', %17.15g ', s(n)); -end -fprintf(fid, '};\n') - -fprintf(fid, 'double p_Dr[%d][%d] = {', Np, Np); -for m=1:Np -fprintf(fid, '{%17.15g ', Dr(m,1)); -for n=2:Np - fprintf(fid, ', %17.15g ', Dr(m,n)); -end -if(m<Np) -fprintf(fid, '},\n') -else -fprintf(fid, '}};\n') -end -end - -fprintf(fid, 'double p_Ds[%d][%d] = {', Np, Np); -for m=1:Np -fprintf(fid, '{%17.15g ', Ds(m,1)); -for n=2:Np - fprintf(fid, ', %17.15g ', Ds(m,n)); -end -if(m<Np) -fprintf(fid, '},\n') -else -fprintf(fid, '}};\n') -end -end - -fprintf(fid, 'double p_LIFT[%d][%d] = {', Np, Nfp*Nfaces); -for m=1:Np -fprintf(fid, '{%17.15g ', LIFT(m,1)); -for n=2:Nfp*Nfaces - fprintf(fid, ', %17.15g ', LIFT(m,n)); -end -if(m<Np) -fprintf(fid, '},\n') -else -fprintf(fid, '}};\n') -end -end - -fprintf(fid, 'int p_Fmask[%d][%d] = {', Nfaces, Nfp); -for m=1:Nfaces -fprintf(fid, '{%d ', Fmask(1,m)); -for n=2:Nfp - fprintf(fid, ', %d ', Fmask(n,m)); -end -if(m<Nfaces) -fprintf(fid, '},\n') -else -fprintf(fid, '}};\n') -end -end - -fprintf(fid, sprintf('\n#endif\n')); -fclose(fid) |
