Skip to content
Snippets Groups Projects

Update PrintPDF.py

Merged Julian Garcia Pardinas requested to merge patch-1 into master
1 file
+ 3
5
Compare changes
  • Side-by-side
  • Inline
@@ -375,14 +375,12 @@ def hj1j2j1pj2p(j1,j2,j1p,j2p):
def comp_num_fit(j1,j2,h,j1p,j2p,hp):
factor = one/three/three/PhSp*two**(3)*pi
factor = one#one/three/three/PhSp*two**(3)*pi
if (j1==j1p) and (j2==j2p) and (h==hp): return (factor*re(IT_j1j2hj1pj2php(j1,j2,h,j1p,j2p,hp)*Nj1j2hj1pj2php(j1,j2,h,j1p,j2p,hp)*hj1j2j1pj2p(j1,j2,j1p,j2p))*ghhp(phi,h,hp)*fjjphhp(cos1,j1,j1p,h,hp)*fjjphhp(cos2,j2,j2p,h,hp)).simplify()
else: return (factor*two*re(IT_j1j2hj1pj2php(j1,j2,h,j1p,j2p,hp)*Nj1j2hj1pj2php(j1,j2,h,j1p,j2p,hp)*hj1j2j1pj2p(j1,j2,j1p,j2p))*ghhp(phi,h,hp)*fjjphhp(cos1,j1,j1p,h,hp)*fjjphhp(cos2,j2,j2p,h,hp)).simplify()
term_list = [[1,1,0,1,1,0],[1,1,1,1,1,1],[1,1,2,1,1,2],[1,1,1,1,1,0],[1,1,2,1,1,0],[1,1,2,1,1,1],[1,0,0,1,0,0],[1,0,0,1,1,0],[1,0,0,1,1,1],[1,0,0,1,1,2],[1,0,0,0,0,0],[0,1,0,0,1,0],[0,1,0,1,1,0],[0,1,0,1,1,1],[0,1,0,1,1,2],[0,1,0,0,0,0],[0,1,0,1,0,0],[0,0,0,0,0,0],[0,0,0,1,1,0],[0,0,0,1,1,1],[0,0,0,1,1,2]]
def angular_comp(j1,j2,h,j1p,j2p,hp):
def printi(i):
return comp_num_fit(*term_list[i])
return (Nj1j2hj1pj2php(j1,j2,h,j1p,j2p,hp)*ghhp(phi,h,hp)*fjjphhp(cos1,j1,j1p,h,hp)*fjjphhp(cos2,j2,j2p,h,hp)).simplify()
\ No newline at end of file
Loading