From: Wolfgang Bangerth Date: Tue, 18 Aug 2009 21:24:31 +0000 (+0000) Subject: Remove assertions. It is true that they are satisfied all the time but X-Git-Tag: v8.0.0~7261 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=58c6adf58f21f11ab79477b6f4718bbf1ea64796;p=dealii.git 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 --- 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