From: wolf Date: Tue, 3 May 2005 15:10:23 +0000 (+0000) Subject: Again fix up one more TODO. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ce2fc0fd23b508764d31d12dda154f75456ab1d7;p=dealii-svn.git Again fix up one more TODO. git-svn-id: https://svn.dealii.org/trunk@10620 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-18/step-18.cc b/deal.II/examples/step-18/step-18.cc index 22d697c3ee..90a5558607 100644 --- a/deal.II/examples/step-18/step-18.cc +++ b/deal.II/examples/step-18/step-18.cc @@ -1531,22 +1531,32 @@ namespace QuasiStaticElasticity // @sect4{TopLevel::solve_linear_problem} - // Solving the linear system again works - // mostly as before. The only difference is - // that we want to only keep a complete - // local copy of the solution vector - // instead of the distributed one that we - // get as output from PETSc's solver - // routines. To this end, we declare a - // local temporary variable for the - // distributed vector, solve with it, and - // at the end of the function copy it again - // into the complete local vector that we - // declared as a member variable. Hanging - // node constraints are then distributed + // Solving the linear system again + // works mostly as before. The only + // difference is that we want to + // only keep a complete local copy + // of the solution vector instead + // of the distributed one that we + // get as output from PETSc's + // solver routines. To this end, we + // declare a local temporary + // variable for the distributed + // vector and initialize it with + // the contents of the local + // variable (remember that the + // ``apply_boundary_values'' + // function called in + // ``assemble_system'' preset the + // values of boundary nodes in this + // vector), solve with it, and at + // the end of the function copy it + // again into the complete local + // vector that we declared as a + // member variable. Hanging node + // constraints are then distributed // only on the local copy, - // i.e. independently of each other on each - // of the processors: + // i.e. independently of each other + // on each of the processors: template unsigned int TopLevel::solve_linear_problem () { @@ -1554,7 +1564,6 @@ namespace QuasiStaticElasticity distributed_incremental_displacement (mpi_communicator, dof_handler.n_dofs(), n_local_dofs); -//TODO document distributed_incremental_displacement = incremental_displacement; //TODO: make more robust against changes in the size of the domain!