fixed error in recognition of overlaps

example/ALL_test.cpp

@@ -2262,7 +2262,7 @@ int main(int argc, char** argv)
                         for (int n = 0; n < global_dim[curr_dim]; ++n)
                         {
                             /* 
-                            if (local_rank == 1)
+                            if (local_rank == 13)
                                 std::cout << "DEBUG: " << n << " " <<
                                     old_border.at(0) << " " <<
                                     old_border.at(1) << " / " <<
@@ -2277,14 +2277,18 @@ int main(int argc, char** argv)
                                     (borders.at(2*n) <= old_border.at(0) &&
                                      borders.at(2*n+1) > old_border.at(0)) ||
                                     (borders.at(2*n) <= old_border.at(1) &&
-                                     borders.at(2*n+1) >= old_border.at(1))
+                                     borders.at(2*n+1) >= old_border.at(1)) ||
+                                    (borders.at(2*n) > old_border.at(0) &&
+                                     borders.at(2*n+1) < old_border.at(1))
                                )
                                 send_neig.push_back(n);
                             if (
                                     (old_borders.at(2*n) <= local_borders.at(0) &&
                                      old_borders.at(2*n+1) > local_borders.at(0)) ||
                                     (old_borders.at(2*n) <= local_borders.at(1) &&
-                                     old_borders.at(2*n+1) >= local_borders.at(1))
+                                     old_borders.at(2*n+1) >= local_borders.at(1)) ||
+                                    (old_borders.at(2*n) > local_borders.at(0) &&
+                                     old_borders.at(2*n+1) < local_borders.at(1))
                                )
                                 recv_neig.push_back(n);
                         }
@@ -2359,6 +2363,20 @@ int main(int argc, char** argv)
 
                         int recvs = 0;
 
+                        /*
+                        std::cout << local_rank << " "  << send_neig.at(0) << " "
+                                                        << ((send_neig.size() >= 2)?send_neig.at(1):-1) << " "
+                                                        << ((send_neig.size() >= 3)?send_neig.at(2):-1) << " "
+                                                        << send_neig.size() << " | "
+                                                        << recv_neig.at(0) << " "
+                                                        << ((recv_neig.size() >= 2)?recv_neig.at(1):-1) << " "
+                                                        << ((recv_neig.size() >= 3)?recv_neig.at(2):-1) << " "
+                                                        << recv_neig.size() << " | "
+                                                        << old_border.at(0) << " " << old_border.at(1) << " | "
+                                                        << local_borders.at(0) << " " << local_borders.at(1)
+                                                        << std::endl;
+                        */
+
                         std::vector<std::vector<double>> recv_vec(recv_neig.size());
                         while (recvs < recv_neig.size())
                         {