[Carpet] IOBasic output issues

Wolfgang Kastaun kastaun at tat.physik.uni-tuebingen.de
Tue Sep 19 20:07:47 CEST 2006


Erik Schnetter wrote:

> I think the "triggers" mechanism is too simplistic for mesh  refinement
> with subcycling in time.  It may need to be extended.
> 
>> -If I perform 1D/2D runs with no ghostzones in  one or two directions,
>>  and demand reductions of variables which are only computed when
>>  triggered by an IO routine, I get:
>>  build/CarpetReduce/reduce.cc:1103: int CarpetReduce::ReduceGVs (const
>> cGH*, int, int, int, void*, int,
>>  const int*, const CarpetReduce::reduction*):
>>  Assertion `nghostzones[d]>=0 && 2*nghostzones[d]<=lsh[d]' failed.
>>  Maybe this is connected to the problem above.
> 
> 
> This seems unrelated to me.  If you have zero ghost zones, then the 
> assertion should succeed.  What are the values of nghostzones and of 
> lsh at these lines?
> 
lsh=1, nghostzones=2
I was talking of outer ghost zones (boundary points). Carpet
seems to check if internal ghostzones fit into the domain ?
I attach the par file.


Wolfgang.
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: igtoy_2d40_20ms_pert.par
Url: /archives/developers/attachments/20060919/1b222f36/attachment.txt 


More information about the developers mailing list