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

    if kwargs['accelerated']:
        return c_accel_mlem(tomo, center, recon, theta, **kwargs)