From d9bc797ff9c7a524946082ec3b28ae7c9b66d265 Mon Sep 17 00:00:00 2001 From: wolf Date: Thu, 29 Mar 2001 14:26:06 +0000 Subject: [PATCH] Fix one TODO (how to generate infinity). git-svn-id: https://svn.dealii.org/trunk@4326 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/base/source/function_lib.cc | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) 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; -- 2.39.5