From 58c6adf58f21f11ab79477b6f4718bbf1ea64796 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 18 Aug 2009 21:24:31 +0000 Subject: [PATCH] Remove assertions. It is true that they are satisfied all the time but just not on very coarse meshes where sometimes we have quadrature points inside the inner ring because we only use a Q1 mapping. git-svn-id: https://svn.dealii.org/trunk@19307 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index e757edb002..f1f46a3a61 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -138,9 +138,6 @@ namespace EquationData const double rho = (r-R0)/h; - Assert (rho >= -.001, ExcInternalError()); - Assert (rho <= 1.001, ExcInternalError()); - return T1+(T0-T1)*((1-rho)*(1-rho)); } @@ -1306,7 +1303,7 @@ compute_viscosity (const std::vector &old_temperature, // side is cheap anyway so we won't even // notice that this part is not parallized // by threads. - // + // // Regarding the implementation of // inhomogeneous Dirichlet boundary // conditions: Since we use the temperature @@ -2462,7 +2459,7 @@ local_assemble_temperature_rhs (const std::pair global_T_range, time_step * nu * ext_grad_T * scratch.grad_phi_T[i] + - time_step * + time_step * scratch.gamma_values[q] * scratch.phi_T[i]) * scratch.temperature_fe_values.JxW(q); @@ -2470,15 +2467,15 @@ local_assemble_temperature_rhs (const std::pair global_T_range, if (temperature_constraints.is_inhomogeneously_constrained(data.local_dof_indices[i])) { for (unsigned int j=0; j