def test_Fz(self):
        r = np.random.rand(self.M.nFz)
        P = self.M.getInterpolationMat(self.M.gridFz, 'Fz')
        assert np.abs(P[:,(self.M.nFx+self.M.nFy):]*r - r).max() < TOL
dataset/ETHPy150Open simpeg/simpeg/tests/mesh/test_TreeMesh.py/Test3DInterpolation.test_Fz/original