def TD_to_1D_All_Point(pk_or_kr):
x = []
y = []
for i in range(number):
for j in range(number):
for k in range(number):
x.append(
math.sqrt((i - number / 2)**2 + (j - number / 2)**2 +
(k - number / 2)**2))
y.append(pk_or_kr[i][j][k])
xx = list(set(x))
xx.sort()
stop = ml.Fib_Search(number / 2, xx)
del xx[stop + 1:]
del yy[stop + 1:]
block = len(xx)
yy_n = [0] * block
yy = [0] * block
for i in range(len(x)):
j = ml.Fib_Search(x[i], xx)
if j <> -1:
yy_n[j] += 1
yy[j] += y[i]
for j in range(block):
yy[j] = yy[j] / yy_n[j]
return xx, yy
