code: def _call_vecfield_p(self, vf, out):
          
          # Optimization for 1 component - just absolute value (maybe weighted)
          if len(self.domain) == 1:
              vf[0].ufuncs.absolute(out=out)
              if self.is_weighted:
                  out *= self.weights[0] ** (1 / self.exponent)
              return
  
          # Initialize out, avoiding one copy
          self._abs_pow_ufunc(vf[0], out=out, p=self.exponent)
          if self.is_weighted:
              out *= self.weights[0]
  
          tmp = self.range.element()
          for fi, wi in zip(vf[1:], self.weights[1:]):
              self._abs_pow_ufunc(fi, out=tmp, p=self.exponent)
              if self.is_weighted:
                  tmp *= wi
              out += tmp
  
          self._abs_pow_ufunc(out, out=out, p=(1 / self.exponent))
docs: Implement ``self(vf, out)`` for exponent 1 < p < ``inf``.