Backup before changing the BC

src/scheme/ScalarNodalScheme.cpp

@@ -518,6 +518,7 @@ class ScalarNodalSchemeHandler::ScalarNodalScheme : public ScalarNodalSchemeHand
                   const NodeId node_id = face_nodes[i_node];
                   if (not node_is_neumann[node_id]) {
                     compute_node_boundary_values[node_id] += value_list[i_face] * mes_l[face_id];
+                    // compute_node_boundary_values[node_id] = xr[node_id][0];
                     sum_mes_l[node_id] += mes_l[face_id];
                   } else {
                     compute_node_boundary_values[node_id] = value_list[i_face];
@@ -994,7 +995,7 @@ class ScalarNodalSchemeHandler::ScalarNodalScheme : public ScalarNodalSchemeHand
             ur[node_id] += Pj[cell_id] * Cjr(cell_id, i_node);
           }
           ur[node_id] = -inverse(node_betar[node_id]) * ur[node_id];
-          // std::cout << node_id << "; ur = " << ur[node_id] << "\n";
+          std::cout << node_id << "; ur = " << ur[node_id] << "\n";
         } else if (is_boundary_node[node_id] and node_is_dirichlet[node_id]) {
           for (size_t i_cell = 0; i_cell < node_to_cell.size(); ++i_cell) {
             const CellId cell_id = node_to_cell[i_cell];
@@ -1006,8 +1007,8 @@ class ScalarNodalSchemeHandler::ScalarNodalScheme : public ScalarNodalSchemeHand
           } else {
             ur[node_id] = inverse(corner_betar[node_id]) * ur[node_id];
           }
+          std::cout << "bord : " << node_id << "; ur = " << ur[node_id] << "\n";
         }
-        // std::cout << node_id << "; ur = " << ur[node_id] << "\n";
       }
 
       LinearSolver solver;