QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Minor fixes 

Former-commit-id: b421a45d8dbc3990ebebecbcea55be588c08c4af
This commit is contained in:
Marek Nečada 2016-07-06 12:00:39 +03:00
parent 16c98b8c8f
commit 3018a3999f
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Scattering and Shit.lyx
View File

@ -532,7 +532,7 @@ Definition [T](2.40);
\begin_inset Formula \begin_inset Formula
\begin{eqnarray*} \begin{eqnarray*}
\widetilde{\vect N}_{mn}^{(j)} & = & \frac{n(n+1)}{kr}\sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}P_{n}^{m}\left(\cos\theta\right)e^{im\phi}z_{n}^{j}\left(kr\right)\hat{\vect r}\\ \widetilde{\vect N}_{mn}^{(j)} & = & \frac{n(n+1)}{kr}\sqrt{\frac{2n+1}{4\pi}\frac{\left(n-m\right)!}{\left(n+m\right)!}}P_{n}^{m}\left(\cos\theta\right)e^{im\phi}z_{n}^{j}\left(kr\right)\hat{\vect r}\\
& & +\left[\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}+i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right)\\ & & +\left[\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}+i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}\frac{1}{kr}\frac{\ud\left(kr\, z_{n}^{j}\left(kr\right)\right)}{\ud(kr)}\\
\widetilde{\vect M}_{mn}^{(j)} & = & \left[i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}-\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right) \widetilde{\vect M}_{mn}^{(j)} & = & \left[i\tilde{\pi}_{mn}\left(\cos\theta\right)\hat{\vect{\theta}}-\tilde{\tau}_{mn}\left(\cos\theta\right)\hat{\vect{\phi}}\right]e^{im\phi}z_{n}^{j}\left(kr\right)
\end{eqnarray*} \end{eqnarray*}

4
qpms/qpms_p.py
View File

@ -423,7 +423,7 @@ def Ã(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J):
+math.lgamma(n+ν+1) - math.lgamma(2*(n+ν)+1)) +math.lgamma(n+ν+1) - math.lgamma(2*(n+ν)+1))
presum = math.exp(exponent) presum = math.exp(exponent)
presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n) / (4*n) presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n) / (4*n)
qmax = floor(q_max(-m,n,μ,ν)) #nemá tu být +m? qmax = math.floor(q_max(-m,n,μ,ν)) #nemá tu být +m?
q = np.arange(qmax+1, dtype=int) q = np.arange(qmax+1, dtype=int)
# N.B. -m !!!!!! # N.B. -m !!!!!!
a1q = a_q(-m,n,μ,ν) # there is redundant calc. of qmax a1q = a_q(-m,n,μ,ν) # there is redundant calc. of qmax
@ -455,7 +455,7 @@ def B̃(m,n,μ,ν,kdlj,θlj,φlj,r_ge_d,J):
presum = math.exp(exponent) presum = math.exp(exponent)
presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n+1) / ( presum = presum * np.exp(1j*(μ-m)*φlj) * (-1)**m * 1j**(ν+n+1) / (
(4*n)*(n+1)*(n+m+1)) (4*n)*(n+1)*(n+m+1))
Qmax = floor(q_max(-m,n+1,μ,ν)) Qmax = math.floor(q_max(-m,n+1,μ,ν))
q = np.arange(Qmax+1, dtype=int) q = np.arange(Qmax+1, dtype=int)
if (μ == ν): # it would disappear in the sum because of the factor (ν-μ) anyway if (μ == ν): # it would disappear in the sum because of the factor (ν-μ) anyway
ã2q = 0 ã2q = 0