code: def decomposition(self, image, x, y, n_max, beta, deltaPix, center_x=0, center_y=0):
          
          num_param = self.shapelets.num_param(n_max)
          param_list = np.zeros(num_param)
          amp_norm = 1. * deltaPix**2
          L_list = self._pre_calc(x, y, beta, n_max, center_x, center_y)
          for i in range(num_param):
              base = self._pre_calc_function(L_list, i) * amp_norm
              param = np.sum(image*base)
              n, m, complex_bool = self.shapelets.index2poly(i)
              if m != 0:
                  param *= 2
              param_list[i] = param
          return param_list
docs: decomposes an image into the shapelet coefficients in same order as for the function call
          :param image:
          :param x:
          :param y:
          :param n_max:
          :param beta:
          :param center_x:
          :param center_y:
          :return: