Added new methods to NmssmSoftsusy class:  P1SfSfCouplings (+ 1-->2, 1-->3), getP1HiggsTriCoup (+1-->2,3), getP1NeutralinoCoup (+1-->2,3) and the pseudoscalar self energies pip1p1, pip1p2 etc in the basis (Im(Hd), Im(Hu), Im(S)).  Tested in two independent versions and against Flexi SUSY.  Matches Flexi expect for b0 gauge boson parts where there is discrepancy between results from SARAH and from Degrassi and Slavich.  We think Degrassi, Slavich are correct but hopefully will get confirmation from Florian Staub soon.