|
@@ -140,7 +140,7 @@ STARPUFFT(twist1_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
STARPUFFT(complex) * restrict in = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[0]);
|
|
|
STARPUFFT(complex) * restrict twisted1 = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[1]);
|
|
|
|
|
|
- //printf("twist1 %d %d %g\n", i, j, (double) cabs(plan->in[i+j]));
|
|
|
+ /* printf("twist1 %d %d %g\n", i, j, (double) cabs(plan->in[i+j])); */
|
|
|
|
|
|
for (k = 0; k < n2; k++)
|
|
|
for (l = 0; l < m2; l++)
|
|
@@ -169,7 +169,7 @@ STARPUFFT(fft1_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
_fftw_complex * restrict worker_in1 = (STARPUFFT(complex) *)plan->plans[workerid].in1;
|
|
|
_fftw_complex * restrict worker_out1 = (STARPUFFT(complex) *)plan->plans[workerid].out1;
|
|
|
|
|
|
- //printf("fft1 %d %d %g\n", i, j, (double) cabs(twisted1[0]));
|
|
|
+ /* printf("fft1 %d %d %g\n", i, j, (double) cabs(twisted1[0])); */
|
|
|
|
|
|
memcpy(worker_in1, twisted1, plan->totsize2 * sizeof(*worker_in1));
|
|
|
_FFTW(execute)(plan->plans[workerid].plan1_cpu);
|
|
@@ -198,7 +198,7 @@ STARPUFFT(twist2_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
|
|
|
STARPUFFT(complex) * restrict twisted2 = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[0]);
|
|
|
|
|
|
- //printf("twist2 %d %d %g\n", kk, ll, (double) cabs(plan->fft1[kk+ll]));
|
|
|
+ /* printf("twist2 %d %d %g\n", kk, ll, (double) cabs(plan->fft1[kk+ll])); */
|
|
|
|
|
|
for (kkk = 0; kkk < n3; kkk++) {
|
|
|
int k = kk * n3 + kkk;
|
|
@@ -218,8 +218,8 @@ STARPUFFT(fft2_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
{
|
|
|
struct STARPUFFT(args) *args = _args;
|
|
|
STARPUFFT(plan) plan = args->plan;
|
|
|
- //int kk = args->kk;
|
|
|
- //int ll = args->ll;
|
|
|
+ /* int kk = args->kk; */
|
|
|
+ /* int ll = args->ll; */
|
|
|
int workerid = starpu_worker_get_id();
|
|
|
|
|
|
task_per_worker[workerid]++;
|
|
@@ -227,7 +227,7 @@ STARPUFFT(fft2_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
const STARPUFFT(complex) *twisted2 = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[0]);
|
|
|
STARPUFFT(complex) *fft2 = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[1]);
|
|
|
|
|
|
- //printf("fft2 %d %d %g\n", kk, ll, (double) cabs(twisted2[plan->totsize4-1]));
|
|
|
+ /* printf("fft2 %d %d %g\n", kk, ll, (double) cabs(twisted2[plan->totsize4-1])); */
|
|
|
|
|
|
_fftw_complex * restrict worker_in2 = (STARPUFFT(complex) *)plan->plans[workerid].in2;
|
|
|
_fftw_complex * restrict worker_out2 = (STARPUFFT(complex) *)plan->plans[workerid].out2;
|
|
@@ -259,7 +259,7 @@ STARPUFFT(twist3_2d_kernel_cpu)(void *descr[], void *_args)
|
|
|
|
|
|
const STARPUFFT(complex) * restrict fft2 = (STARPUFFT(complex) *)STARPU_VECTOR_GET_PTR(descr[0]);
|
|
|
|
|
|
- //printf("twist3 %d %d %g\n", kk, ll, (double) cabs(fft2[0]));
|
|
|
+ /* printf("twist3 %d %d %g\n", kk, ll, (double) cabs(fft2[0])); */
|
|
|
|
|
|
for (kkk = 0; kkk < n3; kkk++) {
|
|
|
int k = kk * n3 + kkk;
|
|
@@ -530,7 +530,7 @@ STARPUFFT(plan_dft_2d)(int n, int m, int sign, unsigned flags)
|
|
|
/* Create twist1 task */
|
|
|
plan->twist1_tasks[z] = task = starpu_task_create();
|
|
|
task->cl = &STARPUFFT(twist1_2d_codelet);
|
|
|
- //task->buffers[0].handle = to be filled at execution
|
|
|
+ /* task->buffers[0].handle = to be filled at execution */
|
|
|
task->buffers[0].mode = STARPU_R;
|
|
|
task->buffers[1].handle = plan->twisted1_handle[z];
|
|
|
task->buffers[1].mode = STARPU_W;
|