def METHOD_NAME(self,vec):
    cntr = np.zeros(3)
    for i in range(3):
        cntr[i] = vec[3+i] + (1.0-self.fixpt[0])*(vec[9+i]-vec[6+i])
        cntr[i] += self.fixpt[1]*(vec[i]-vec[3+i])
        cntr[i] += self.fixpt[2]*(vec[6+i]-vec[3+i])
    range = 0.14
    if ((abs(cntr[0]) < range) and (abs(cntr[1]) < range) and (abs(cntr[2]) < range)):
        self.qua.append([cntr[0],cntr[1],cntr[2]])
    radius = 0.0
    for i in range(3):
        radius += cntr[i]*cntr[i]
    radius = sqrt(radius)
    if (radius < range):
        self.sph.append([cntr[0],cntr[1],cntr[2]])
    return