compute_AR_CI_optimized function

Joint Anderson-Rubin CI computation