[Carpet] Assertion `leveltimes.at(ml).at(rl) == global_time' failed

[Bruno Giacomazzo]

Erik,

On Fri, 16 May 2008, Erik Schnetter wrote:
> yes, I introduced this check recently to ensure that the times on all 
> refinement levels remain consistent.  I have now made these comparisons less 
> strict.

 	I have still problems. Now it stops at the beginning (iteration 1) 
with the following error:

cactus_WhiskyExp_CarpetGit-5(mpi:15584 at node048): 
/data1/bgiacoma/Cactus/configs/WhiskyExp_CarpetGit/build/Carpet/Evolve.cc:103: 
int Carpet::Evolve(tFleshConfig *): Assertion `abs 
(leveltimes.at(ml).at(rl) - global_time) == eps * global_time' failed

Cheers,
Bruno


> I am currently trying to ensure that the times associated with the individual 
> refinement levels are consistent with the global time in Carpet.  There were 
> some inconsistencies, reported by Ian Hinder and Nils Dorband, which should 
> be corrected now.  This code is an additional check that should ensure 
> consistency.  According to the IEEE floating point standard these conditions 
> should hold, but maybe your compiler is using strong optimisations.  If not, 
> then there is still some inconsistency left.
>
> These inconsistencies are small, of the order of floating point round-off, 
> but they accumulate over time, and I want to prevent this.
>
> -erik
>
>

-- 
Dr. Bruno Giacomazzo
Max Planck Institute for Gravitational Physics
Albert Einstein Institute
Am Muehlenberg 1
D-14476 Potsdam
Germany

Tel.  : +49 331 567 7183
Fax   : +49 331 567 7252
cell. : +49 173 826 4488
email : bgiacoma at aei.mpg.de

-------------------------------------------------
There are only 10 types of people in the world:
Those who understand binary, and those who don't
-------------------------------------------------