code: def proximal(self):
          
          if self.pointwise_norm.exponent == np.inf:
              return proximal_convex_conj_l1(space=self.domain)
          elif self.pointwise_norm.exponent == 2:
              return proximal_convex_conj_l1_l2(space=self.domain)
          else:
              raise NotImplementedError('`proximal` only implemented for p = 1 '
                                        'or 2')
docs: Return the `proximal factory` of the functional.
  
          See Also
          --------
          odl.solvers.nonsmooth.proximal_operators.proximal_convex_conj_l1 :
              `proximal factory` for the L1-norms convex conjugate.