smpsi function

Auxiliary functions for smoothm