[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