From: wolf Date: Wed, 12 Jan 2000 10:33:22 +0000 (+0000) Subject: Ack. Wrong equation. Also wrong order for boundary and hanging node constraints. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=45c3c2b3ddf2a76521dbcafaa23575af9486ca6a;p=dealii-svn.git 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 --- 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); };