As Multiple Inputs Multiple Outputs (MIMO) is becoming one of the enable techniques in modern wireless communication like 5G/6G and beyond, it is important to design efficient controls for the MIMO antenna arrays to realize critical functions such as high-throughput communication and beamforming. Efficient and rigorous Computational Electromagnetics (CEM) algorithms are key for such control design problem, especially for massive MIMO antenna arrays and beyond. Here we present such universal Fast Fourier Transform (FFT) based iterative CEM algorithm for efficient control design of massive MIMO antenna array and beyond, for example Reconfigurable Intelligent Surface and Large Intelligent Surface (RIS/LIS) MIMO. Our FFT-based iterative CEM algorithm is universal and works for both discrete MIMO and quasi-