code: def _shapeletOutput(self, x, y, beta, shapelets, precalc=True):
          
          n = len(np.atleast_1d(x))
          if n <= 1:
              values = 0.
          else:
              values = np.zeros(len(x[0]))
          n = 0
          k = 0
          i = 0
          num_n = len(shapelets)
          while i < num_n * (num_n+1)/2:
              values += self._function(x, y, shapelets[n-k][k], beta, n-k, k, precalc=precalc)
              k += 1
              if k == n + 1:
                  n += 1
                  k = 0
              i += 1
          return values
docs: returns the the numerical values of a set of shapelets at polar coordinates
          :param shapelets: set of shapelets [l=,r=,a_lr=]
          :type shapelets: array of size (n,3)
          :param coordPolar: set of coordinates in polar units
          :type coordPolar: array of size (n,2)
          :returns:  array of same size with coords [r,phi]
          :raises: AttributeError, KeyError