@@ -2460,12 +2460,6 @@ def RBsymChk(RBsym,cubic,coefNames,L=18):
24602460 :param list coefNames: sp. harm coefficient names to be checked/converted
24612461 :param int L: maximum spherical harmonic order no. for cubic generation if needed
24622462 '''
2463- # cubicsigns = {'C(3,1)c':[-1,],'C(4,1)c':[1,1,],'C(6,1)c':[1,-1],'C(6,2)c':[1,-1],'C(7,1)c':[1,-1],'C(8,1)c':[1,1,1],
2464- # 'C(9,1)c':[1,-1],'C(9,2)c':[1,-1],'C(10,1)c':[1,-1,-1],'C(10,2)c':[1,1,-1]}
2465- # cubicnames = {'C(3,1)c':['C(3,2)',],'C(4,1)c':['C(4,0)','C(4,4)'],'C(6,1)c':['C(6,0)','C(6,4)'],
2466- # 'C(6,2)c':['C(6,2)','C(6,6)'],'C(7,1)c':['C(7,2)','C(7,6)'],'C(8,1)c':['C(8,0)','C(8,4)','C(8,8)'],
2467- # 'C(9,1)c':['C(9,2)','C((9,6)'],'C(9,2)c':['C(9,4)','C(9,8)'],
2468- # 'C(10,1)c':['C(10,0)','C(10,4)','C(10,8)'],'C(10,2)c':['C(10,2)','C(10,6)','C(10,10)']}
24692463 newNames = []
24702464 newSgns = []
24712465 if cubic : #sytsym is a cubic site
@@ -3021,7 +3015,8 @@ def SphHarmAng(L,M,P,Th,Ph):
30213015
30223016 :returns ylmp value/array: as reals
30233017 '''
3024- ylmp = spsp .sph_harm (M ,L ,rpd * Th ,rpd * Ph ) #wants radians; order then degree
3018+ ylmp = spsp .sph_harm (M ,L ,rpd * Th ,rpd * Ph ) #wants radians; order then degree; not normalized
3019+ #### TODO: this will be deprecated in future scipy; new one sph_harm_y in scipy 1.15.1
30253020
30263021 if M > 0 :
30273022 return (- 1 )** M * P * np .real (ylmp )* SQ2
0 commit comments