- reset_x(index.size(), NULL, index.size(), NULL, etalon, w_rls[j]);
- reset_x(index.size(), NULL, index.size(), NULL, traj, w_rls[j]);
- do_fit(index.size(), w_rls[j], etalon, traj);
- update_graph(graph[j], traj, epsi);
- sfree(etalon);
- sfree(traj);
+ domains_local[a].push_back(i);
+ flags[i] = false;
+ }
+ }
+
+ snew(w_rls, domains_local.size());
+ for (int i = 0; i < domains_local.size(); i++) {
+ snew(w_rls[i], index.size());
+ for (int j = 0; j < index.size(); j++) {
+ w_rls[i][j] = 0;
+ }
+ for (int j = 0; j < domains_local[i].size(); j++) {
+ w_rls[i][domains_local[i][j]] = 1;
+ }
+ }
+
+ frankenstein_trajectory.resize(trajectory.size());
+ frankenstein_trajectory.back().resize(index.size());
+
+ for (int i = 0; i < domains_local.size(); i++) {
+ rvec *basic, *traj;
+ snew(basic, index.size());
+ for (int k = 0; k < index.size(); k++) {
+ copy_rvec(trajectory[basic_frame][k].as_vec(), basic[k]);
+ }
+ snew(traj, index.size());
+ for (int k = 0; k < index.size(); k++) {
+ copy_rvec(trajectory.back()[k].as_vec(), traj[k]);