summaryrefslogtreecommitdiff
path: root/benchmarks/CUDA/DG/buildsource.m
diff options
context:
space:
mode:
authorTor Aamodt <[email protected]>2010-10-01 08:55:28 -0800
committerTor Aamodt <[email protected]>2010-10-01 08:55:28 -0800
commit11b308e7363e937966b035b4891db32b4eece3bf (patch)
tree50ca4c9ad6f163ac4acb2bf505e64dfebed66947 /benchmarks/CUDA/DG/buildsource.m
parentbb820c116764d7a1b8e071137d32b74e7f34dd2f (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.m106
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)