| 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)); |