]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Don't forget to also multiply the preconditioner term with eta.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 23 Oct 2008 18:37:54 +0000 (18:37 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 23 Oct 2008 18:37:54 +0000 (18:37 +0000)
git-svn-id: https://svn.dealii.org/trunk@17326 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-31/step-31.cc
deal.II/examples/step-32/step-32.cc

index 80fd1c401662c62074541cb0a3d1d0fe749b9857..e1ca873c0a176b10bc1dd14f4013bb6e0ea6ad9a 100644 (file)
@@ -1277,7 +1277,8 @@ BoussinesqFlowProblem<dim>::assemble_stokes_preconditioner ()
          
          for (unsigned int i=0; i<dofs_per_cell; ++i)
            for (unsigned int j=0; j<dofs_per_cell; ++j)
-             local_matrix(i,j) += (scalar_product (phi_grad_u[i], phi_grad_u[j])
+             local_matrix(i,j) += (EquationData::eta *
+                                   scalar_product (phi_grad_u[i], phi_grad_u[j])
                                    +
                                    phi_p[i] * phi_p[j])
                                   * stokes_fe_values.JxW(q);
index 0c7e52e46ef0e177d735ac7ea508a6219c6d45df..55426bf18a3105b1abda272416fd6e4eb93b4b2d 100644 (file)
@@ -766,7 +766,8 @@ BoussinesqFlowProblem<dim>::assemble_stokes_preconditioner ()
 
            for (unsigned int i=0; i<dofs_per_cell; ++i)
              for (unsigned int j=0; j<dofs_per_cell; ++j)
-               local_matrix(i,j) += (scalar_product (phi_grad_u[i], phi_grad_u[j])
+               local_matrix(i,j) += (EquationData::eta *
+                                     scalar_product (phi_grad_u[i], phi_grad_u[j])
                                      +
                                      phi_p[i] * phi_p[j])
                                    * stokes_fe_values.JxW(q);

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.