So far I have been calculating Spi from the density, but there really is no need to do this.  I have the local mu(r) so I can just go straight ahead for Spi.  The density with NLCE is useful because I can use it to get the global MU for the given N and then just stick with that global Mu