TDHF for even operators about ampsci HOT 4 CLOSED

Also: doesn't work for E1 v_form (but diagram method does).
Not sure if this is the same issue

from ampsci.

Note: It doesn't seem like this is actually a problem with solving TDHF.
Rather, this implies the issue is in the dV equations?
See commit below; introduced a new tests for hyperfine, testing equations for \chi and \eta - and they work extremely well (implying we solve TDHF equations perfectly) - note: the \delte E is included in this, so we're not just skipping the hard part.

1ced43b

from ampsci.

Note: Seems to be different if I swap:

Xx = HF::solveMixedState(Xx.k, Fc, omega, m_vl, m_alpha, *p_core, rhs, eps_ms, nullptr, p_VBr, m_Hmag);
HF::solveMixedState(Xx, Fc, omega, m_vl, m_alpha, *p_core, rhs, eps_ms, nullptr, p_VBr, m_Hmag);

This gives opposite conclusion to last comment; implying issue is with solving equations.
Difference between these two functions is one starts from scratch, other starts from existing solution

from ampsci.

TDHF - physics

• Mixed-states work with non-local DiracODE? (#11)
• Issue for even-parity operators (from de?) (#3)
• Perhaps linked to (#11)
• dV conj - sometimes causes issues; + not consistent
• PNC: TDHF vs Diagram?
• Note: TDHF doesn't work for E1v, but diagram does!

TDHF - performance

• Parallelised inefficiently (#6)
• Diagram RPA: use Qk table?
• Write out to disk?
• Fewer allocations?
• Use Qk,Pk functions?
• RPAD: minor "eps" race cond?

from ampsci.