]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Ack. Wrong equation. Also wrong order for boundary and hanging node constraints.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 12 Jan 2000 10:33:22 +0000 (10:33 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 12 Jan 2000 10:33:22 +0000 (10:33 +0000)
git-svn-id: https://svn.dealii.org/trunk@2208 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/Attic/examples/step-by-step/step-8/step-8.cc
deal.II/examples/step-8/step-8.cc

index e46a322f0151428778304b975dfc9d0a29c35197..75a1b9bd2da559039f39216e6e0757867bc6e66e 100644 (file)
@@ -543,10 +543,6 @@ void ElasticProblem<dim>::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<dim>::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<dim>::assemble_system ()
                                           system_matrix,
                                           solution,
                                           system_rhs);
-
-  hanging_node_constraints.condense (system_matrix);
-  hanging_node_constraints.condense (system_rhs);
 };
 
 
index e46a322f0151428778304b975dfc9d0a29c35197..75a1b9bd2da559039f39216e6e0757867bc6e66e 100644 (file)
@@ -543,10 +543,6 @@ void ElasticProblem<dim>::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<dim>::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<dim>::assemble_system ()
                                           system_matrix,
                                           solution,
                                           system_rhs);
-
-  hanging_node_constraints.condense (system_matrix);
-  hanging_node_constraints.condense (system_rhs);
 };
 
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.