Sparse matrix approximation method for an active optical control system

We develop a sparse matrix approximation method to decompose a wave front into a basis set of actuator influence functions for an active optical system consisting of a deformable mirror and a segmented primary mirror. The wave front used is constructed by Zernike polynomials to simulate the output of a phase-retrieval algorithm. Results of a Monte Carlo simulation of the optical control loop are compared with the standard, nonsparse approach in terms of accuracy and precision, as well as computational speed and memory. The sparse matrix approximation method can yield more than a 50-fold increase in the speed and a 20-fold reduction in matrix size and a commensurate decrease in required memory, with less than 10% degradation in solution accuracy. Our method is also shown to be better than when elements are selected for the sparse matrix on a magnitude basis alone. We show that the method developed is a viable alternative to use of the full control matrix in a phase-retrieval-based active optical control system.