Commit 24ea3e51 authored by Philippe Helluy's avatar Philippe Helluy

up

parent 278a84ec
......@@ -299,7 +299,7 @@ __kernel void time_step(__global const real *fn, __global real *fnp1) {
f2w(fnow + 0, wnow + 0);
real fnext[_N];
// first order relaxation
//first order relaxation
w2f(wnow + 0, fnext + 0);
real om = 2;
......@@ -321,8 +321,8 @@ __kernel void time_step(__global const real *fn, __global real *fnp1) {
// save
for (int ik = 0; ik < _N; ik++) {
int imem = i + j * _NX + ik * ngrid;
fnp1[imem] = fnext[ik];
// fnp1[imem] = fnow[ik];
//fnp1[imem] = fnext[ik];
fnp1[imem] = fnow[ik];
}
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment