Rajat, The first equation is the usual equilibrium equation. It is solved to give > displacements. The stiffness matrix is then multiplied by this global > displacement vector to give the global traction vector which forms the rhs > for the second equation. > The easiest solution to deal with constraints to deal with them while assembling the matrix. Therefore, I would not compute the right-hand side for the second equation beforehand. Instead just do it while assembling the second equation and use constraints.distribute_local_to_global().
Best, Daniel -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
