Differentially Private Approximation of the John Ellipsoid

We study the problem of approximating the John ellipsoid (JE) of a given (centrally symmetric) polytope of $n$ constraints in a Euclidean space under differential privacy (DP). We give the first differentially private algorithm for this problem under the standard model, where neighboring datasets may differ arbitrarily in one a single constraint. Our work also extends to the complimentary problem of Minimum Enclosing Ellipsoid of $n$ points in the Euclidean space. Our approach is based on the recent non-private multiplicative-weights algorithm of pmlr-v99-cohen19a. First we introduce a non-private generalization of the Cohen et al algorithm, yielding a $(1+\gamma )$-approximation of the JE problem while violating at most $\kappa n$ constraints in $O(\log (1/\kappa )/\gamma )$ iterations. This variant works by projecting the intermediate weights assigned to the constraints onto the set of $\kappa$-dense distributions, similarly to bun2020efficientnoisetolerantprivatelearning. We then design a $\rho$-zCDP variant of this algorithm by adding Gaussian noise to the weighted covariance matrix aggregated in each step of the algorithm. Under a mild goodness assumption on the data we can assert that the resulting noisy matrix is close to the true matrix, thereby achieving essentially the same guarantee as the non-private algorithm provided sufficiently many input points. Thus our method achieves an efficient DP poly-time algorithm under concrete sample complexity bounds.