From 42ddea690b3bc3b36f1f317709366e4808a578a7 Mon Sep 17 00:00:00 2001
From: schnetter <>
Date: Tue, 3 Aug 2004 17:39:00 +0000
Subject: 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
---
 Carpet/CarpetRegrid/src/manualcoordinatelist.cc | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'Carpet/CarpetRegrid')

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],
-- 
cgit v1.2.3