tdvp_order throws an error for order 4

TDVP_order which determines the number of forward and backward sweeps does not work for order=4.
This is because the function in tdvporder.jl

function orderings(::TDVPOrder{4,direction}) where {direction}
  return [direction, Base.ReverseOrdering(direction)]
end

returns an array of length 2. However, the for loop in the file tdvp_step.jl runs over the length of sub_time_steps, which is an array of length 6 according to the definition given in tdvporder.jl.

  for substep in 1:length(sub_time_steps)
    psi, PH, info = tdvp_sweep(
      orderings[substep], solver, PH, sub_time_steps[substep], psi; current_time, kwargs...
    )

This causes an error as attached in the image. One can overcome this error by changing the function as

function orderings(::TDVPOrder{4,direction}) where {direction}
  return [direction, Base.ReverseOrdering(direction),direction, Base.ReverseOrdering(direction),direction, Base.ReverseOrdering(direction)]
end

But not sure if this a correct approach. Please check this.
Thank you

Start by reading the orderings definition in tdvporder.jl alongside the sub_time_steps handling in tdvp_step.jl. Trace the order=4 sweep sequence and verify that every substep has a corresponding ordering; done means order=4 completes without an indexing error and preserves the intended forward/backward sweep pattern.