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

[Erik Schnetter]

Bruno,

thanks for the pointer.  I wanted to rush the correction it, but this  
seems to have been too fast.

-erik

On May 16, 2008, at 10:52:13, Bruno Giacomazzo wrote:

> 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
> -------------------------------------------------
> _______________________________________________
> developers mailing list
> developers at lists.carpetcode.org
> http://lists.carpetcode.org/listinfo/developers
>


-- 
Erik Schnetter <schnetter at cct.lsu.edu>   http://www.cct.lsu.edu/~eschnett/

My email is as private as my paper mail.  I therefore support encrypting
and signing email messages.  Get my PGP key from www.keyserver.net.



-------------- next part --------------
A non-text attachment was scrubbed...
Name: PGP.sig
Type: application/pgp-signature
Size: 194 bytes
Desc: This is a digitally signed message part
Url : /archives/developers/attachments/20080516/13a7a1c9/attachment-0001.pgp