From 45c3c2b3ddf2a76521dbcafaa23575af9486ca6a Mon Sep 17 00:00:00 2001 From: wolf Date: Wed, 12 Jan 2000 10:33:22 +0000 Subject: [PATCH] Ack. Wrong equation. Also wrong order for boundary and hanging node constraints. git-svn-id: https://svn.dealii.org/trunk@2208 0785d39b-7218-0410-832d-ea1e28bc413d --- .../Attic/examples/step-by-step/step-8/step-8.cc | 10 +++------- deal.II/examples/step-8/step-8.cc | 10 +++------- 2 files changed, 6 insertions(+), 14 deletions(-) diff --git a/deal.II/deal.II/Attic/examples/step-by-step/step-8/step-8.cc b/deal.II/deal.II/Attic/examples/step-by-step/step-8/step-8.cc index e46a322f01..75a1b9bd2d 100644 --- a/deal.II/deal.II/Attic/examples/step-by-step/step-8/step-8.cc +++ b/deal.II/deal.II/Attic/examples/step-by-step/step-8/step-8.cc @@ -543,10 +543,6 @@ void ElasticProblem::assemble_system () shape_grads[j][q_point][component_j] * (lambda_values[q_point] + mu_values[q_point])) - + // (mu d_i u_j, d_i v_j) - (shape_grads[i][q_point][component_j] * - shape_grads[j][q_point][component_i] * - mu_values[q_point]) + // (mu d_i v_j, d_i v_j) ((component_i == component_j) ? (shape_grads[i][q_point] * @@ -592,6 +588,9 @@ void ElasticProblem::assemble_system () }; }; + hanging_node_constraints.condense (system_matrix); + hanging_node_constraints.condense (system_rhs); + // The interpolation of the // boundary values needs a small // modification: since the solution @@ -621,9 +620,6 @@ void ElasticProblem::assemble_system () system_matrix, solution, system_rhs); - - hanging_node_constraints.condense (system_matrix); - hanging_node_constraints.condense (system_rhs); }; diff --git a/deal.II/examples/step-8/step-8.cc b/deal.II/examples/step-8/step-8.cc index e46a322f01..75a1b9bd2d 100644 --- a/deal.II/examples/step-8/step-8.cc +++ b/deal.II/examples/step-8/step-8.cc @@ -543,10 +543,6 @@ void ElasticProblem::assemble_system () shape_grads[j][q_point][component_j] * (lambda_values[q_point] + mu_values[q_point])) - + // (mu d_i u_j, d_i v_j) - (shape_grads[i][q_point][component_j] * - shape_grads[j][q_point][component_i] * - mu_values[q_point]) + // (mu d_i v_j, d_i v_j) ((component_i == component_j) ? (shape_grads[i][q_point] * @@ -592,6 +588,9 @@ void ElasticProblem::assemble_system () }; }; + hanging_node_constraints.condense (system_matrix); + hanging_node_constraints.condense (system_rhs); + // The interpolation of the // boundary values needs a small // modification: since the solution @@ -621,9 +620,6 @@ void ElasticProblem::assemble_system () system_matrix, solution, system_rhs); - - hanging_node_constraints.condense (system_matrix); - hanging_node_constraints.condense (system_rhs); }; -- 2.39.5