def METHOD_NAME():
  # generate test map
  d = 10
  gridding = flex.grid(d,d,d)
  test_map = flex.double(gridding)
  parameters = (1.5,1.5,1.5,-2.0,-2.0,-2.0,1.25,1.25,1.25,24.0)
  p = fit_peak.gaussian(parameters=parameters)
  for i in range(d):
    for j in range(d):
      for k in range(d):
        test_map[(i,j,k)] = p.get_height(r=(float(i)/d,float(j)/d,float(k)/d))
  # check fitting
  neighbors = fit_peak.pick_map_neighbors(site=(4,4,4),map_in=test_map)
  height_list,xyz_list = neighbors.get_26_nearest_neighbors()
  fp = fit_peak.fit_peak(height_list=height_list,xyz_list=xyz_list,
                         shape="gaussian")
  assert approx_equal(fp.x, parameters)
  assert approx_equal(fp.vertex, p.vertex)