A multistage algorithm is presented to allow for the convenient and accurate computation of the finite amplitude-distorted acoustic fields produced by phased arrays and nonplanar sources. The algorithm also allows for non-normal (angled) beam computations.