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

[Bruno Giacomazzo]

Erik,
 	btw, why did you write:

             assert (abs (leveltimes.at(ml).at(rl) - global_time) ==
                       eps * global_time);

 	shouldn't be :

 	            assert (abs (leveltimes.at(ml).at(rl) - global_time) <=
                       eps * global_time);


Cheers,
Bruno

On Fri, 16 May 2008, Bruno Giacomazzo wrote:

> 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
-------------------------------------------------