tomography working with hack to call operator directly
added iterative method for speed