& , objf_smrsss
_RL objf_smrsss (nsx,nsy)
& objf_smrsss, num_smrsss,
& objf_smrsst,objf_smrsss
& + mult_smrsss * objf_smrsss(bi,bj)
f_smrsss = f_smrsss + objf_smrsss(bi,bj)
objf_smrsss(bi,bj) = 0. _d 0