[cig-commits] (no subject)

[dealii.demon at gmail.com]

-e To: cig-commits at geodynamics.org
Subject: commit by heien to /var/svn/dealii/aspect

Revision 1818

Bug fix

U   trunk/aspect/source/initial_conditions/box.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1818&peg=1818

Diff:
Modified: trunk/aspect/source/initial_conditions/box.cc
===================================================================
--- trunk/aspect/source/initial_conditions/box.cc	2013-07-08 22:15:11 UTC (rev 1817)
+++ trunk/aspect/source/initial_conditions/box.cc	2013-07-16 23:37:10 UTC (rev 1818)
@@ -129,7 +129,7 @@
             {
               perturbation = inclusion_temperature - ambient_temperature;
               for (int d=0; d<dim; ++d)
-                perturbation *= exp(-pow(position[d] - center[d],2) / (2 * pow((radius / 4), 2))) / (2 * radius);
+                perturbation *= exp(-8*pow(position.distance(center)/radius, 2));
             }
           else if (inclusion_gradient == "linear")
             {