  The time-independent nonlinear Schr\"odinger equation is solved for two
attractive delta-function shaped potential wells where an imaginary loss term
is added in one well, and a gain term of the same size but with opposite sign
in the other. We show that for vanishing nonlinearity the model captures all
the features known from studies of PT symmetric optical wave guides, e.g., the
coalescence of modes in an exceptional point at a critical value of the
loss/gain parameter, and the breaking of PT symmetry beyond. With the
nonlinearity present, the equation is a model for a Bose-Einstein condensate
with loss and gain in a double well potential. We find that the nonlinear
Hamiltonian picks as stationary eigenstates exactly such solutions which render
the nonlinear Hamiltonian itself PT symmetric, but observe coalescence and
bifurcation scenarios different from those known from linear PT symmetric
Hamiltonians.
