[Einstein Toolkit] #851: ReflectionSymmetry Interpolation Errors and WeylScal4

List overview All Threads
Download

newer

older

[Einstein Toolkit] #856:...

[Einstein Toolkit] #686: Handle...

Einstein Toolkit

1 May 2012 1 May '12

4:09 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 --------------------+------------------------------------------------------- Reporter: tbode | Owner: tbode Type: defect | Status: new Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Keywords: | --------------------+------------------------------------------------------- There is an inconsistency in ReflectionSymmetry on whether a scalar can have a tensor parity tag that isn't unity (pseudoscalars). In ReflectionSymmetry's apply.c, a scalar GF's tensorparity flag is taken into account, while in interpolation this is ignored. In the ET Toolkit, this only effects those using the ET version of WeylScal4 since only Psi4i, Psi2i, and Psi0i have all been registered as scalars with tensorparity=-1. When interpolating, e.g. Psi4, onto a sphere for mode decomposition, the interpolated quantities have sign errors in Psi4i and give erroneous gravitational wave modes.

This problem applies to both Maxwell release and current development branch. Symptoms reported to me by Jim Healy. I would suggest backporting the resulting fix to Maxwell.

Either ReflectionSymmetry's interpolation has to allow for tensorparity=-1, or apply.c's acceptance of tensorparity=-1 should be removed and WeylScal4's GFs reverted to manually specified parities.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Show replies by date

Einstein Toolkit

1 May 1 May

4:42 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: new Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by rhaas):

So that's what this was. Jim had reported this to me and I had looked at this but only at the grid function level (where everything was ok, in fact the one with symmetries was better). I think ReflectionSymmetry should take the parity into account. I suspect this would affect eg. the magnetic field as well, yes? Do you have a patch?

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:1 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

4:48 p.m.

Comment (by tbode):

This interpolation problem seems isolated to the cases where the interpolated variables are pseudo-scalars. Pseudovectors aren't affected. I'm looking at possible fixes. No patch quite yet.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:2 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

5:58 p.m.

Comment (by rhaas):

It turns out I had a patch since #254 and had apparently forgotten to include it in the commit. The fact that it fixes pseudoscalars is accidental since #254 was intended to fix pseudovectors. Has been in use at GT for a while.

Ok to apply?

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:3 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

6 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: review Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------ Changes (by eschnett):

* status: new => review

Comment:

Please apply.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:4 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

6 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------ Changes (by eschnett):

* status: review => reviewed_ok

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:5 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

6:05 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by eschnett):

Can you add a test case as well?

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:6 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

6:06 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by rhaas):

yes.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:7 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

7:30 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by rhaas):

I applied the code patch. Tanja will provide a test either a distinct one for ReflectionSymmetry or extend WeylScal4's teukolskyID test which is one of the few thorns actually using pseudo-quantities.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:8 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

2 May 2 May

3:56 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by tbode):

I have an extension for the teukolskyID testsuite which uses InterpToArray to dump a small 2d interpolated slice across the symmetry boundary which I will push to WeylScal4 soon. This clearly shows the bug without the patch. I highly recommend the patch to ReflectionSymmetry be backported to Maxwell so the release version can be used in vanilla form for gravitational wave extraction and pseudo-vector evolution/analysis. I have confirmed that the patch does not make any Maxwell testsuites fail.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:9 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

3:58 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by knarf):

I agree that the backport should be done.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:10 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

4:16 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by rhaas):

Merged into Maxwell. Should we also merge into older releases?

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:11 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

4:42 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: reviewed_ok Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: | Keywords: ---------------------+------------------------------------------------------

Comment (by knarf):

No, I don't think so. We support the current release only, besides trunk/

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:12 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

7:17 p.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: closed Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: fixed | Keywords: ---------------------+------------------------------------------------------ Changes (by rhaas):

* status: reviewed_ok => closed * resolution: => fixed

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:13 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

Einstein Toolkit

3 May 3 May

3:21 a.m.

#851: ReflectionSymmetry Interpolation Errors and WeylScal4 ---------------------+------------------------------------------------------ Reporter: tbode | Owner: tbode Type: defect | Status: closed Priority: major | Milestone: ET_2012_05 Component: Cactus | Version: development version Resolution: fixed | Keywords: ---------------------+------------------------------------------------------

Comment (by hinder):

Would it be possible to make this a separate test? In general, it's best for the tests to be a fine-grained as possible. In this case, the same parameter file could be used apart from the output sections. Don't worry if you don't have time though.

-- Ticket URL: https://trac.einsteintoolkit.org/ticket/851#comment:14 Einstein Toolkit http://einsteintoolkit.org The Einstein Toolkit

5272

Age (days ago)

5274

Last active (days ago)

trac@lists.einsteintoolkit.org

14 comments

1 participants

tags (0)

participants (1)

Einstein Toolkit