diff options
| author | schnetter <> | 2004-08-03 17:39:00 +0000 |
|---|---|---|
| committer | schnetter <> | 2004-08-03 17:39:00 +0000 |
| commit | 42ddea690b3bc3b36f1f317709366e4808a578a7 (patch) | |
| tree | 22a8aa6917927cdd5cdc3dcc9bcf06ca96349ad5 /Carpet | |
| parent | 97ebb5c43caf9d42225a368f52e9b5dab5c97efc (diff) | |
Use pow() instead of ipow() because the exponent can be negative.
Use pow() instead of ipow() because the exponent can be negative.
Fixes a bug with negative convergence levels.
darcs-hash:20040803173945-07bb3-cd6785a069013a9ae3370ac476f02f24e6769bb6.gz
Diffstat (limited to 'Carpet')
| -rw-r--r-- | Carpet/CarpetRegrid/src/manualcoordinatelist.cc | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc index acb2683ec..d41855b23 100644 --- a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc +++ b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc @@ -13,7 +13,7 @@ #include "regrid.hh" extern "C" { - static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetRegrid/src/manualcoordinatelist.cc,v 1.10 2004/08/02 11:41:05 schnetter Exp $"; + static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetRegrid/src/manualcoordinatelist.cc,v 1.11 2004/08/03 19:39:45 schnetter Exp $"; CCTK_FILEVERSION(Carpet_CarpetRegrid_manualcoordinatelist_cc); } @@ -82,7 +82,7 @@ namespace CarpetRegrid { for (size_t c=0; c<newobss.at(rl).size(); ++c) { for (int d=0; d<dim; ++d) { assert (mglevel==0); - rvect const spacing = base_spacing * ipow(mgfact, basemglevel) / ipow(reffact, rl+1); + rvect const spacing = base_spacing * pow(CCTK_REAL(mgfact), basemglevel) / ipow(reffact, rl+1); ierr = ConvertFromPhysicalBoundary (dim, &physical_min[0], &physical_max[0], &interior_min[0], &interior_max[0], |