So changed the error over to a tsi, eta representation before the chi squared is calculated to prevent the under-detemined error matrices that were causing the inversion problem. Shapes of lone ip fraction much better, but still getting a hard shoulder on the C++ at sqrt(1000) microns!
Investigating....