code: def obfn_reg(self):
          r
          l1 = np.sum(mp_wl1*np.abs(self.obfn_gvar()))
          return (self.lmbda*l1, l1)
docs: r"""Compute regularisation term, :math:`\| x \|_1`, and
          contribution to objective function.