hydroable version

src/scheme/HyperelasticSolver.cpp

@@ -235,7 +235,7 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS
         for (size_t j = 0; j < node_to_cell.size(); ++j) {
           const CellId J       = node_to_cell[j];
           const unsigned int R = node_local_number_in_its_cells[j];
-          br += Ajr(J, R) * u[J] + sigma[J] * Cjr(J, R);
+          br += Ajr(J, R) * u[J] - sigma[J] * Cjr(J, R);
         }
 
         b[r] = br;
@@ -361,7 +361,7 @@ class HyperelasticSolverHandler::HyperelasticSolver final : public HyperelasticS
         const auto& cell_nodes = cell_to_node_matrix[j];
 
         for (size_t r = 0; r < cell_nodes.size(); ++r) {
-          F(j, r) = Ajr(j, r) * (u[j] - ur[cell_nodes[r]]) + sigma[j] * Cjr(j, r);
+          F(j, r) = -Ajr(j, r) * (u[j] - ur[cell_nodes[r]]) + sigma[j] * Cjr(j, r);
         }
       });