Forgot to mention that the phase equalisation and matching will be constrained mostly below 1 kHz because it seems to have more perceptual ramifications there.The plan is to create a FIR inverse filter to compensate for both the magnitude and phase response followed by matching of the left and right channels such that both the magnitude and phase looks similar. The spatial measurements (per channel) are constrained to a smal listening zone to derive an average response representative of the zone to be optimised. Because the responses change significantly with a small change in the spatial locations, the equalisation of an average response is required to minimise overfitting and overcompensating by the inverse filter. Granted, the dispersion of the impulse response is an issue in these types of environments and it helps if option exist to minimise it. I think it will be a great idea to also have an averaging function based on an RMS averaged magnitude and a vector averaged phase to see which works better in this case. I will appreciate it if you can consider it.