}
for (int j = 0; j < frankenstein_trajectory.size() - i; j++) {
for (int f1 = 0; f1 < index.size(); f1++) {
- for (int f2 = f1; f2 < index.size(); f2++) {
+ for (int f2 = 0; f2 < index.size(); f2++) {
rvec temp1, temp2;
rvec_sub(frankenstein_trajectory[basic_frame][f1], frankenstein_trajectory[j][f1], temp1);
rvec_dec(temp1, medx[f1]);
}
}
for (int j = 0; j < index.size(); j++) {
- for (int f = j; f < index.size(); f++) {
+ for (int f = 0; f < index.size(); f++) {
crltns[i][j][f] = a[j][f] / (std::sqrt(b[j][f] * c[j][f]));
}
}