[Carpet] poison in Rotatingsymmetry90

corvino at fis.unipr.it corvino at fis.unipr.it
Tue Jun 26 10:19:46 CEST 2007


Thanks Erik,

I will check the output with the right grid parameters and find out what
to do.

Giovanni



> Hi Giovanni.
>
> On Jun 19, 2007, at 04:38:53, corvino at fis.unipr.it wrote:
>
>> Hi, I'm running a simulation in bitant symmetry with reflection
>> along the
>> z axis. I poisoned the symmetry boundary using the parameter in
>> RotatingSymmetry90 and I don't understand the output. In attach
>> there are
>> the parameter file, the std output and error of the run, I'm using
>> carpet
>> stable 3.
>
> You are really only simulating an octant, since you are also using
> RotatingSymmetry90.  Your grid setup is inconsistent with that; the
> lower boundary in the x and y directions should then be at 0, not at
> -55.
>
>> My questions are:
>> Does this output mean that the poison is the cause of the
>> simulation crash?
>
> Yes.
>
>> How do I interpret that list of triplets of number: are they the
>> coordinates or indexes of the point where poison was found?
>
> Yes, they are the indices.
>
>> These numbers refer to points where poison was found or not?
>
> Where poison was found.
>
>> In this case how many poisoned points does the code print to screen?
>
> All poisoned points.  You can also look at the source code, file
> rotatingsymmetry90.c.  Look for the place where the string "ijk"
> appears -- this is where the output is generated.
>
>
>
> The error that RotatingSymmetry90 reports means that it could not
> apply the symmetry condition.  The most likely cause is that the grid
> structure itself does not have the necessary symmetry.
>
> Do you want to use a rotating symmetry?  Depending on this, you
> should either deactivate the thorn or change the grid structure.
>
> -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.
>
>
>
> _______________________________________________
> developers mailing list
> developers at lists.carpetcode.org
> http://lists.carpetcode.org/listinfo/developers
>




More information about the developers mailing list