We developed a new numerical simulation method to calculate the amplification of arbitrary laser pulses in double-pass amplifiers in which two polarization-orthogonal pulses inside the gain medium overlap. The proposed method agrees very well with the experimental results given in this work for the pulse energy and temporal shape.