A wide range of cost functionals that describe the criteria for designing optimal pulses can be reduced to two basic functionals by the introduction of product spaces. We extend previous monotonically convergent algorithms to solve the generalized pulse design equations derived from those basic functionals. The new algorithms are proved to exhibit monotonic convergence. Numerical tests are implemented in four-level model systems employing stationary and/or nonstationary targets in the absence and/or presence of relaxation. Trajectory plots that conveniently present the global nature of the convergence behavior show that slow convergence may often be attributed to "trapping" and that relaxation processes may remove such unfavorable behavior. (C) 2004 American Institute of Physics.