code: def _shapeletOutput(self, r, phi, beta, shapelets):
          
          if type(r) == float or type(r) == int or type(r) == type(np.float64(1)) or len(r) <= 1:
              values = 0.
          else:
              values = np.zeros(len(r), 'complex')
          for nl in range(0,len(shapelets)): #sum over different shapelets
              for nr in range(0,len(shapelets)):
                  value = shapelets[nl][nr]*self._chi_lr(r, phi, nl, nr, beta)
                  values += value
          return values.real
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