+ domains[i].resize(0);
+ check_domains(delta, frames, graph[i]);
+ std::cout << "finding domains' sizes\n";
+ find_domain_sizes(graph[i], domsizes);
+
+ /*for (int j1 = 0; j1 < domsizes.size(); j1++) {
+ for (int j2 = 0; j2 < domsizes[j1].size(); j2++) {
+ if (domsizes[j1][j2] != 0) {
+ std::cout << domsizes[j1][j2] << " ";
+ }
+ }
+ }*/
+
+ while (check_domsizes(domsizes, domain_min_size)) {
+ domains[i].push_back(a);
+ if (DomainSearchingAlgorythm == 0) {
+ get_maxsized_domain(domains[i].back(), graph[i], domsizes);
+ } else if (DomainSearchingAlgorythm == 1) {
+ get_min_domain(domains[i].back(), graph[i], domsizes, domain_min_size);
+ }
+ std::cout << "new domain: " << domains[i].back().size() << " atoms\n";
+ delete_domain_from_graph(graph[i], domains[i].back());
+ domsizes.resize(0);
+ find_domain_sizes(graph[i], domsizes);
+ }