code: def obfn_reg(self):
          
  
          rl1 = np.linalg.norm((self.wl1 * self.obfn_gvar()).ravel(), 1)
          return (self.lmbda*rl1, rl1)
docs: Compute regularisation term and contribution to objective
          function.