diff --git a/src/scheme/ScalarNodalScheme.cpp b/src/scheme/ScalarNodalScheme.cpp index 22ecc8509048a45c158247fc492a77ebb0de0f1d..4f9719f639e6181fafecb8eb02f1ae3919572a63 100644 --- a/src/scheme/ScalarNodalScheme.cpp +++ b/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;