def c_accel_mlem(tomo, center, recon, theta, **kwargs):

    if LIB_TOMOPY_ACCEL is None:
        _missing_library("MLEM ACCEL")