Draft implementation of Cov NG terms

[carlosggarcia]

I have started the implementation of the missing terms in the Cov NG. So far I have implemented the T^2h_13 term and checked that runs but haven't checked if the result makes sense.

I'm following Eq. 3.25 of https://arxiv.org/pdf/1912.08209.pdf

[damonge]

Good start! I think I'd leave the 2,2 part for now until we have a good idea for how to integrate over the angle...

[c-d-leonard]

FYI @carlosggarcia this is failing the tests due to some flake8 formatting issues. You can see them if you climb on the Details above for the continuous integration tests. Let me know if you need any help with fixing these.

[c-d-leonard]

@carlosggarcia could we get a status update on this? :) anywhere you need help to move forward?

[carlosggarcia]

@c-d-leonard The update is that there's no update. I was focusing on adding NaMaster to TJPCov and now on getting it work with TXPipe. So the problem has been time, not an actual problem with the implementation!

[c-d-leonard]

Hi @carlosggarcia just checking whether the update is still there's no update? no worries if so :) just doing an inventory before next telecon

[carlosggarcia]

@c-d-leonard I'm afraid it is. Sorry! I'll try get back to this in a few weeks.

[c-d-leonard]

Looks like some progress here :) @carlosggarcia any support you need on this at this point?

[carlosggarcia]

@c-d-leonard, thanks for asking! @damonge is already helping me with the maths. I think that will be enough for now.

[c-d-leonard]

Just wondering how things are going with this at this point @carlosggarcia? No pressure, just want to check how it's going and if you need anything.

[carlosggarcia]

@c-d-leonard same as two months ago. I'll try to give it a push in the next two weeks.

[carlosggarcia]

Finally, this is ready for review! @damonge

[damonge]

Thanks a lot @carlosggarcia , this is really great.

See comments below. One overall comment: the new functions in pk_4pt could be optimized by comparing profiles, and avoiding recomputing some integrals if they've been already computed (you can check out what's being done in the 1h trispectrum). If you can see an easy way of doing that now, go for it, but I'm also happy to just leave a comment and optimise this later on (since this is already a huge PR!).

[carlosggarcia]

I have addressed the easy comments. I'll look now how to optimize it.

[carlosggarcia]

back to you @damonge

[damonge]

A few more comments (plus remaining previous conversations)

[carlosggarcia]

I think all it's done now. Back to you @damonge. Thanks for reviewing!

[damonge]

Final comments on the tests. We're almost there!

[damonge]

LGTM. Feel free to merge if the tests pass!