From a13f7d36d7d60172604247dd4dd0ab4d01ca9ee9 Mon Sep 17 00:00:00 2001
From: Drouard <axelle.drouard@orange.fr>
Date: Wed, 25 May 2022 11:25:34 +0200
Subject: [PATCH] fix theta_jr
---
src/scheme/ScalarNodalScheme.cpp | 21 ++++++++++-----------
1 file changed, 10 insertions(+), 11 deletions(-)
diff --git a/src/scheme/ScalarNodalScheme.cpp b/src/scheme/ScalarNodalScheme.cpp
index 90bcff7d4..4b65e08eb 100644
--- a/src/scheme/ScalarNodalScheme.cpp
+++ b/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) {
- 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";
- compute_theta(cell_id, i_node) = dot(inverse(node_betar[node_id]) * Cjr(cell_id, i_node), v[node_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];
+ 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;
}();
--
GitLab