From: wolf Date: Thu, 29 Mar 2001 14:26:06 +0000 (+0000) Subject: Fix one TODO (how to generate infinity). X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d9bc797ff9c7a524946082ec3b28ae7c9b66d265;p=dealii-svn.git Fix one TODO (how to generate infinity). git-svn-id: https://svn.dealii.org/trunk@4326 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/base/source/function_lib.cc b/deal.II/base/source/function_lib.cc index bf5221eb31..f859343934 100644 --- a/deal.II/base/source/function_lib.cc +++ b/deal.II/base/source/function_lib.cc @@ -859,11 +859,13 @@ LSingularityFunction::gradient (const Point<2> &p, if ((x>=0) && (y>=0)) { -//TODO: should return infinity, but how to do that? - static const double infty[2] = {0., 0.}; + // in this region, the gradient + // does not exist, so return + // infinity + static const double infty[2] = { atof("Inf"), atof("Inf") }; return Tensor<1,2>(infty); } - + // double phi = atan2(y,-x)+M_PI; // double r2 = x*x+y*y;