40 |
|
* |
41 |
|
* Created by J. Daniel Gezelter on 09/26/06. |
42 |
|
* @author J. Daniel Gezelter |
43 |
< |
* @version $Id: BondOrderParameter.cpp,v 1.23 2009-11-25 20:01:59 gezelter Exp $ |
43 |
> |
* @version $Id$ |
44 |
|
* |
45 |
|
*/ |
46 |
|
|
107 |
|
Wigner3jm(&lPass, &lPass, &lPass, |
108 |
|
&m1Pass, &m2m, &m2M, |
109 |
|
THRCOF, &mSize, &error); |
110 |
< |
|
110 |
> |
|
111 |
|
m2Min[lm] = (int)floor(m2m); |
112 |
|
m2Max[lm] = (int)floor(m2M); |
113 |
|
|
250 |
|
for(int m = -l; m <= l; m++){ |
251 |
|
sphericalHarmonic.setM(m); |
252 |
|
q[std::make_pair(l,m)] += sphericalHarmonic.getValueAt(costheta, phi); |
253 |
+ |
|
254 |
|
} |
255 |
|
} |
256 |
|
nBonds++; |
264 |
|
q2[l] = 0.0; |
265 |
|
for (int m = -l; m <= l; m++){ |
266 |
|
q[std::make_pair(l,m)] /= (RealType)nBonds; |
267 |
+ |
|
268 |
|
q2[l] += norm(q[std::make_pair(l,m)]); |
269 |
|
} |
270 |
|
q_l[l] = sqrt(q2[l] * 4.0 * NumericConstant::PI / (RealType)(2*l + 1)); |