[Carpet] IOBasic output issues

[Erik Schnetter]

On Sep 19, 2006, at 13:07:47, Wolfgang Kastaun wrote:

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

Ah, the confusion due to terminology.  In Cactus, "ghost zones" are  
always points used for inter-processor communication.  Near the outer  
boundary, there are either "boundary zones" or "symmetry zones".

> I attach the par file.

Thanks.

At first glance, this assertion seems to be over-zealous; it should  
not be tested near outer boundaries.  As a short-term work-around,  
you could try commenting out this line.  I haven't looked whether  
later code makes the same assumption, so you may get into trouble for  
that.

You can also set different numbers of ghost zones (Cactus ghost  
zones, not boundary zones) in different directions.  This may be the  
more appropriate thing to do.

-erik

-- 
Erik Schnetter <schnetter at cct.lsu.edu>

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: 186 bytes
Desc: This is a digitally signed message part
Url : /archives/developers/attachments/20060919/0073fc31/attachment.pgp