We present a study of lattice-QCD methods to determine the relevant hadronic form factors for radiative-leptonic decays of pseudoscalar mesons. We provide numerical results for D+s -> e+ nu gamma. Our calculation is performed using a domain-wall action for all quark flavors and on a single RBC/UKQCD lattice gauge -field ensemble. The first part of the study is how to best control two sources of systematic error inherent in the calculation, specifically the unwanted excited states created by the meson interpolating field and unwanted exponentials in the sum over intermediate states. Using a 3D sequential propagator allows for better control over unwanted exponentials from intermediate states, while using a 4D sequential propagator allows for better control over excited states. We perform individual analyses of the 3D and 4D methods, as well as a combined analysis using both methods, and find that the 3D sequential propagator offers good control over both sources of systematic uncertainties for the smallest number of propagator solves. From there, we further improve the use of a 3D sequential propagator by employing an infinite-volume approximation method, which allows us to calculate the relevant form factors over the entire allowed range of photon energies. We then study improvements gained by performing the calculation using a different three-point function, using ratios of three-point functions, averaging over positive and negative photon momentum, and using an improved method for extracting the structure-dependent part of the axial form factor. The optimal combination of methods yields results for the D+s -> e+ nu gamma structure-dependent vector and axial form factors in the entire kinematic range with statistical plus fitting uncertainties of order 5%, using 25 gauge configurations with 64 samples per configuration.