Separating each operation for adjoint testing

Noting kdvTLMPropagator L, it becomes in operator
notation:
    L=R.[opSn(nDt+1).opPn(nDt+1) ...
            ... opSn(2).opPn(2)].opE1.I.F

    Each of theses validate adjoint test.
    But the global subroutine does not (it
    reaches half of precision type; it could be
    a wrong implicit type casting to simple
    precision somewhere, a constant
    multiplication maybe.