fft.forward(*cseq);
  shared_vector tst_real_to_complex(char dir, const fftbx::triple& N)  {    fftbx::triple M = fftbx::Ncomplex_from_Nreal(N);    shared_vector cseq = init_cseq(M);    fftbx::real_to_complex_3d<std::vector<double> > fft(N);    if (dir == 'f') {      fft.forward(*cseq);    }    else {      fft.backward(*cseq);    }    return cseq;  }