def test_calc_lfp_pointsource_anisotropic(self): sigma = [0.6, 0.3, 0.45] cell = TestCell() cell.xmid = cell.ymid = cell.zmid = np.array([1.2]) sigma_r = np.sqrt(sigma[1] * sigma[2] * 1.2**2 + sigma[0] * sigma[2] * 1.2**2 + sigma[0] * sigma[1] * 1.2**2) phi_analytic = 1./(4*np.pi * sigma_r) np.testing.assert_equal(phi_analytic, lfpcalc.calc_lfp_pointsource_anisotropic(cell, x=0, y=0, z=0, sigma=sigma, r_limit=cell.diam/2))
def test_calc_lfp_pointsource_anisotropic(self): sigma = [0.6, 0.3, 0.45] cell = DummyCell() cell.xmid = cell.ymid = cell.zmid = np.array([1.2]) sigma_r = np.sqrt(sigma[1] * sigma[2] * 1.2**2 + sigma[0] * sigma[2] * 1.2**2 + sigma[0] * sigma[1] * 1.2**2) phi_analytic = 1./(4*np.pi * sigma_r) np.testing.assert_equal(phi_analytic, lfpcalc.calc_lfp_pointsource_anisotropic(cell, x=0, y=0, z=0, sigma=sigma, r_limit=cell.diam/2))