minor changes
fix bug allocating too small arrays
Finish forward operator implementation

- Back to double precision NFFT: single did not work correctly
- Implement contrast inversion
- Fix CTF application