  We construct two strong-field path integral methods with full Coulomb
distortion, in which the quantum pathways are mimicked by interfering electron
orbits: the rate-based CQSFA (R-CQSFA) and the hybrid forward-boundary CQSFA
(H-CQSFA). The methods have the same starting point as the standard Coulomb
quantum-orbit strong-field approximation (CQSFA), but their implementation does
not require pre-knowledge of the orbits' dynamics. These methods are applied to
ultrafast photoelectron holography. In the rate-based method, electron orbits
are forward propagated and we derive a non-adiabatic ionization rate from the
CQSFA, which includes sub-barrier Coulomb corrections and is used to weight the
initial orbit ensemble. In the H-CQSFA, the initial ensemble provides initial
guesses for a subsequent boundary problem and serves to include or exclude
specific momentum regions, but the ionization probabilities associated with
individual trajectories are computed from sub-barrier complex integrals. We
perform comparisons with the standard CQSFA and \textit{ab-initio} methods,
which show that the standard, purely boundary-type implementation of the CQSFA
leaves out whole sets of trajectories. We show that the sub-barrier Coulomb
corrections broaden the resulting photoelectron momentum distributions (PMDs)
and improve the agreement of the R-CQSFA with the H-CQSFA and other approaches.
We probe different initial sampling distributions, uniform and otherwise, and
their influence on the PMDs. We find that initial biased sampling emphasizes
rescattering ridges and interference patterns in high-energy ranges, while an
initial uniform sampling guarantees accurate modeling of the holographic
patterns near the ionization threshold or polarization axis. Our results are
explained using the initial to final momentum mapping for different types of
interfering trajectories.
