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

[Bruno Giacomazzo]

On Fri, 16 May 2008, Bruno Giacomazzo wrote:
>
> 	            assert (abs (leveltimes.at(ml).at(rl) - global_time) <=
>                       eps * global_time);

 	with this change it is now running and it has also passed the 
point where it was crashing before.

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