From 7fca0f9c8fb9e116fda63ff153557ec3176c16c4 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 15 Mar 2013 09:15:51 -0400 Subject: more fixes from Dmitry --- phasematching/WGM_lib.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/phasematching/WGM_lib.py b/phasematching/WGM_lib.py index ee03b3d..3beefed 100644 --- a/phasematching/WGM_lib.py +++ b/phasematching/WGM_lib.py @@ -286,7 +286,7 @@ def AngOvlp(L1,p1,L2,p2,L3,p3,limit=3): #def HG(p,L,x): return (2**L*math.factorial(p))**-0.5*hermite(p)(x*math.sqrt(L))*(L/math.pi)**0.25*math.exp(-L*x**2/2.0) def HG(p,L,x): return hermite(p)(x*math.sqrt(L))*math.exp(-L*x**2/2.0) # removed a large factor that would drop out in normalization anyway def HGnorm(p,L,x): return HG(p,L,x)*(4*math.pi*quad(lambda x: HG(p,L,x)**2, 0, limit*((p+1)*math.pi/L/2)**0.5)[0])**(-0.5) - return 2*quad(lambda x: HGnorm(p1,L1,x)*HGnorm(p2,L2,x)*HGnorm(p3,L3,x), 0, limit*((max(p1,p2,p3)+1)*math.pi/min(L1,L2,L3)/2)**0.5)[0] + return 4*math.pi*quad(lambda x: HGnorm(p1,L1,x)*HGnorm(p2,L2,x)*HGnorm(p3,L3,x), 0, limit*((max(p1,p2,p3)+1)*math.pi/min(L1,L2,L3)/2)**0.5)[0] def T_phasematch(L1,p1,q1,L2,p2,q2,L3,p3,q3,R,r,T0=70,MgO=0,MgO_th=5.0,E = 0,Tstep=0.1, n = nLNO1): ''' -- cgit v1.2.3