test_tomopy.py

import tomopy as tp
import numpy as np
import time
#timing functions
def tic():
    global startTime_for_tictoc
    startTime_for_tictoc = time.time()

def toc():
    if 'startTime_for_tictoc' in globals():
       return time.time() - startTime_for_tictoc

# projections sizes [Ntheta,Nz,N]
Nz = 2048
N = 2448
Ntheta = 1500
# init some data
proj = np.ones([Ntheta,Nz,N],dtype='float32')
# init angles
theta = np.linspace(0,np.pi,Ntheta,endpoint=False).astype('float32')
# rotation center
rot_center = N/2

tic()
rec = tp.recon(proj, theta, center = N/2, sinogram_order = False, algorithm = 'gridrec', filter_name = 'hann')
print('Time for reconstruction:', toc())