def test_ip_theta_pi_at_sigma_is_p_irr( self ): import math D = 1e-12 sigma = 1e-8 kf = 0 t = 1e-5 r0 = sigma r = r0 + math.sqrt( 6 * D * t ) gf = mod.BasicPairGreensFunction( D, kf, sigma ) ip = gf.ip_theta( numpy.pi, r, r0, t ) * ( 2 * numpy.pi * r * r ) pirr = mod.p_irr( r, t, r0, kf, D, sigma ) pcorr = gf.ip_corr( numpy.pi, r, r0, t ) * ( 2 * numpy.pi * r * r ) pfree = gf.ip_free( numpy.pi, r, r0, t ) * ( 2 * numpy.pi * r * r ) self.assertNotAlmostEqual( pirr, pfree, 7, 'pcorr estimated to be too small.' + \ ' test may not be valid.' ) #print 'PP', pirr, ip, pcorr, pfree self.assertNotEqual( 0.0, ip ) self.assertAlmostEqual( ip/pirr, 1 )
def p_irr( r, t, r0, kf, D, sigma ): return _gfrd.p_irr( r, t, r0, kf, D, sigma )
def p_irr(r, t, r0, kf, D, sigma): return _gfrd.p_irr(r, t, r0, kf, D, sigma)