From: wolf Date: Mon, 27 May 2002 13:19:40 +0000 (+0000) Subject: Revert the sign of the jump residual and again when summing up, to X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=42743b245daebb254f9db318c308b2902cc0125e;p=dealii-svn.git Revert the sign of the jump residual and again when summing up, to keep with common notation. git-svn-id: https://svn.dealii.org/trunk@5879 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-14/step-14.cc b/deal.II/examples/step-14/step-14.cc index f1d4eeac1b..be09ab869f 100644 --- a/deal.II/examples/step-14/step-14.cc +++ b/deal.II/examples/step-14/step-14.cc @@ -2991,9 +2991,10 @@ namespace LaplaceSolver // sure that the contributions of // each of the faces are there, // and add them up. Only take - // half of the jump term, since - // the other half will be taken - // by the neighboring cell. + // minus one half of the jump + // term, since the other half + // will be taken by the + // neighboring cell. unsigned int present_cell=0; for (active_cell_iterator cell=DualSolver::dof_handler.begin_active(); cell!=DualSolver::dof_handler.end(); @@ -3005,7 +3006,7 @@ namespace LaplaceSolver face_integrals.end(), ExcInternalError()); error_indicators(present_cell) - += 0.5*face_integrals[cell->face(face_no)]; + -= 0.5*face_integrals[cell->face(face_no)]; }; std::cout << " Estimated error=" << std::accumulate (error_indicators.begin(), @@ -3368,7 +3369,7 @@ namespace LaplaceSolver // normal vector: for (unsigned int p=0; pface(face_no)] = face_integral; @@ -3523,7 +3524,7 @@ namespace LaplaceSolver // the other function: for (unsigned int p=0; p