def METHOD_NAME(self):
    # Radial distribution function is obtained for correlations of unity.
    # This is not a particularly good test.
    # Tilt unit cell to check results for non-orthogonal cells.
    for i in range(-1,2):
        a = io.read('aC.cfg')
        a1, a2, a3 = a.cell
        a3[0] += i*a1[0]
        a.set_cell([a1, a2, a3])
        a.set_scaled_positions(a.get_scaled_positions())
        v = np.ones(len(a))
        SCF1, r1 = spatial_correlation_function(a, v, 10.0, 0.1, 8.0, 0.1,
                                                norm=False)
        SCF2, r2 = spatial_correlation_function(a, v, 10.0, 0.1, 2.0, 0.1,
                                                norm=False)
        #print(np.abs(SCF1-SCF2).max())
        #print(r1[np.argmax(np.abs(SCF1-SCF2))])
        #plt.plot(r1, SCF1, label='$8$')
        #plt.plot(r2, SCF2, label='$2$')
        #plt.legend(loc='best')
        #plt.show()
        self.assertTrue(np.abs(SCF1-SCF2).max() < 0.31)