105 |
|
|
106 |
|
Vector3d pos1 = sd1->getPos(); |
107 |
|
Vector3d pos2 = sd2->getPos(); |
108 |
< |
Vector3d r12 = pos1 - pos2; |
108 |
> |
Vector3d r12 = pos2 - pos1; |
109 |
|
currentSnapshot_->wrapVector(r12); |
110 |
|
|
111 |
|
double distance = r12.length(); |
152 |
|
double GofRTheta::evaluateAngle(StuntDouble* sd1, StuntDouble* sd2) { |
153 |
|
Vector3d pos1 = sd1->getPos(); |
154 |
|
Vector3d pos2 = sd2->getPos(); |
155 |
< |
Vector3d r12 = pos1 - pos2; |
155 |
> |
Vector3d r12 = pos2 - pos1; |
156 |
|
currentSnapshot_->wrapVector(r12); |
157 |
|
r12.normalize(); |
158 |
|
Vector3d dipole = sd1->getElectroFrame().getColumn(2); |
162 |
|
|
163 |
|
double GofROmega::evaluateAngle(StuntDouble* sd1, StuntDouble* sd2) { |
164 |
|
Vector3d v1 = sd1->getElectroFrame().getColumn(2); |
165 |
< |
Vector3d v2 = sd1->getElectroFrame().getColumn(2); |
165 |
> |
Vector3d v2 = sd2->getElectroFrame().getColumn(2); |
166 |
|
v1.normalize(); |
167 |
|
v2.normalize(); |
168 |
|
return dot(v1, v2); |