pretty up normalization constant in GaussianUserObject, closes idahol…

…ab#45
rwcarlsen · Oct 4, 2016 · d7965e8 · d7965e8
1 parent c49158d
commit d7965e8
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/userobjects/GaussianUserObject.C b/src/userobjects/GaussianUserObject.C
@@ -45,7 +45,7 @@ GaussianUserObject::GaussianUserObject(const InputParameters & parameters) :
 {
   // normalization constant so integral over domain (theoretically) equals 1
   unsigned int dim = _mesh.dimension();
-  _norm = 1.0/(pow(_sigma,dim)*pow(pow(2.0*pi,dim),0.5));
+  _norm = 1.0/std::pow(_sigma*std::sqrt(2.0*libMesh::pi), Real(dim));
 }
 
 GaussianUserObject::~GaussianUserObject()