+void update_graph_v2(std::vector< std::vector< node > > &ugwi_graph, std::vector < RVec > ugwi_x, long double ugwi_epsi) {
+ rvec ugwi_temp;
+ int ugwi_for = ugwi_graph.size();
+ for (int i = 0; i < ugwi_for; i++) {
+ for (int j = i; j < ugwi_for; j++) {
+ rvec_sub(ugwi_x[i].as_vec(), ugwi_x[j].as_vec(), ugwi_temp);
+ rvec_dec(ugwi_temp, ugwi_graph[i][j].r.as_vec());
+ if (norm(ugwi_temp) <= ugwi_epsi) {
+ if (i == j) {
+ ugwi_graph[i][j].n++;
+ }
+ else {
+ ugwi_graph[i][j].n++;
+ ugwi_graph[j][i].n++;
+ }
+ }
+ }
+ }
+}
+