]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
scalar boundary residuals
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 8 Aug 2012 14:36:32 +0000 (14:36 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 8 Aug 2012 14:36:32 +0000 (14:36 +0000)
git-svn-id: https://svn.dealii.org/trunk@25773 0785d39b-7218-0410-832d-ea1e28bc413d

tests/integrators/laplacian_01.cc

index 6ce7df48931ebd248672751e8b91ecbe5d1f2837..80b13e99b300aedf246d42b91daaf40bfa03cfd3 100644 (file)
@@ -63,12 +63,20 @@ void test_boundary(const FEValuesBase<dim>& fev)
     {
       u = 0.;
       u(i) = 1.;
+      w = 0.;
       fev.get_function_values(u, indices, uval, true);
       fev.get_function_gradients(u, indices, ugrad, true);
       nitsche_residual(w, fev, make_slice(uval), make_slice(ugrad), make_slice(null_val), 17);
       M.vmult(v,u);
       w.add(-1., v);
       deallog << ' ' << w.l2_norm();
+      if (d==1)
+       {
+         nitsche_residual(w, fev, uval[0], ugrad[0], null_val[0], 17);
+         M.vmult(v,u);
+         w.add(-1., v);
+         deallog << " e" << w.l2_norm();         
+       }
     }
   deallog << std::endl;
 }

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.