fix theta_jr

src/scheme/ScalarNodalScheme.cpp

@@ -513,17 +513,16 @@ class ScalarNodalSchemeHandler::ScalarNodalScheme : public ScalarNodalSchemeHand
 
       const NodeValuePerCell<const double> theta = [&] {
         NodeValuePerCell<double> compute_theta{mesh->connectivity()};
-        // parallel_for(
-        //  mesh->numberOfCells(), PUGS_LAMBDA(CellId cell_id) {
-        for (CellId cell_id = 0; cell_id < mesh->numberOfCells(); ++cell_id) {
+        parallel_for(
+          mesh->numberOfCells(), PUGS_LAMBDA(CellId cell_id) {
             const auto& cell_nodes = cell_to_node_matrix[cell_id];
             for (size_t i_node = 0; i_node < cell_nodes.size(); ++i_node) {
               const NodeId node_id = cell_nodes[i_node];
-            std::cout << "node_id = " << node_id << "\n";
+              if (not is_boundary_node[node_id]) {
                 compute_theta(cell_id, i_node) = dot(inverse(node_betar[node_id]) * Cjr(cell_id, i_node), v[node_id]);
               }
             }
-        //          });
+          });
         return compute_theta;
       }();