diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-03-15 09:15:51 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-03-15 09:15:51 -0400 |
| commit | 7fca0f9c8fb9e116fda63ff153557ec3176c16c4 (patch) | |
| tree | a803e9820be72de4bf45cbb49df6485bf883d8ef /phasematching | |
| parent | 13103836b6fd82c47f493860105a0693de5d59f8 (diff) | |
| download | wgmr-7fca0f9c8fb9e116fda63ff153557ec3176c16c4.tar.gz wgmr-7fca0f9c8fb9e116fda63ff153557ec3176c16c4.zip | |
more fixes from Dmitry
Diffstat (limited to 'phasematching')
| -rw-r--r-- | phasematching/WGM_lib.py | 2 |
1 files changed, 1 insertions, 1 deletions
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): ''' |