AbstractWe solve the Nonlinear Schrödinger Equation (NLSE) in 1D in presence of one, two and several Dirac delta potentials. With the help of an equivalent central force problem we obtain the analytical solutions in terms of a biparametric family containing the Jacobi functions. Elliptic Jacobi functions are already reported in the literature but they have not been used in the context of a scattering problem under causal boundary conditions. In the simplest examples of one or two Dirac deltas we analyze how the nonlinear term of the equation affects the modulus and phase profiles of the wave function. We also study the transmission curves under the nonlinear modification of the tunneling behavior for the first time. For a Fabri-Perot configuration made of two deltas, we obtain the effect of nonlinear coupling in the positions of the local maxima (resonances). We lay the foundations for nonlinear Anderson localization of 1D BECs in a speckle field. Upon redefinition of parameters these novel results describe the dynamics of a stationary Higgs field in 1D. Finally, we discuss the conditions for soliton formation under the influence of a Dirac comb potential, giving rise to fully correlated locations and intensities of the defects.