BlackHolePerturbationToolkit / PerturbationEquations

MIT License
5 stars 1 forks source link

SchwarzschildLinearOperator trTensor basis components mixing parity #21

Closed thompsonphys closed 1 year ago

thompsonphys commented 1 year ago

Evaluating SchwarzschildLinearOperator["dG", "Generic", "trTensor", "trTensor"] for the spin-1 components (e.g. "t+", "r-" etc.) produces output containing both even- and odd-parity pieces. An example output is given below.

Screenshot 2023-09-14 at 15 09 06

As we can see, these opposite-parity pieces always show up with an explicit factor of i in front. I'm wondering if one of the vector harmonic projection rules might be broken.

thompsonphys commented 1 year ago

Some minor updates/comments about some strange versioning confusion that I'm encountering (as we discussed a bit on a call). Perhaps @barrywardell knows some of the nuances here.

I've cloned PerturbationEquations and symlinked the repository into my $UserBaseDirectory/Applications folder (equivalently inside the xAct folder). When calling <<xAct`PerturbationEquations`, I recover version 0.0.1 {2023, 6, 30} which does not include some of the recent changes that @adamrpound has implemented.

On the other hand, if I explicitly use SetDirectory to the kernel path of PerturbationEquations and load <<PerturbationEquations.wl, then I load version 0.1.0 {2023, 6, 30} which does include the changes recently committed. Is there any advice on how to potentially avoid this confusion in the future while loading through xAct, or is loading PerturbationEquations.wl directly the recommended approach when working with unreleased development code in Mathematica?

Just as a note, the parity issue above still remains in the more recent version of the code.

barrywardell commented 1 year ago

I recommend only every installing using the Paclet system. If you need to try out a development version from git, then you should follow the developer instructions to temporarily use a version checked out from git. You should avoid putting things in $UserBaseDirectory/Applications or similar as they could cause confusion.

adamrpound commented 1 year ago

@barrywardell Can you reproduce Jonathan's issue with the mixed parity? I don't get the problematic terms when I run the same command as Jonathan.

barrywardell commented 1 year ago

With version 0.1.0 I get

Screenshot 2023-09-16 at 13 57 12

whereas with the current main branch (4705f25) I get

Screenshot 2023-09-16 at 13 58 00
adamrpound commented 1 year ago

I realized that the output is actually correct. The terms with the wrong parity are multiplied by zero. It's just not obvious until you use mutolrule. I've updated the package to set these terms to zero.