C 原地快速排序算法在一些迭代后崩溃

我已经实现了以下快速排序算法来对成对的点（3D空间）进行排序。 每个对定义一条线：目的是将距离小于或等于powR的所有线放置在包含所有对的阵列内。 包含坐标的数组是一维的，每6个元素定义一对，每3个元素定义一个点

当我运行包含3099642个元素的数组的算法时，在处理2799222尝试进入下一个迭代后，该算法停止。如果我从元素2799228启动算法，它将在3066300处停止。 我想不出问题出在哪里，你有什么建议吗

void QuickSort(float *array, int from, int to, float powR){

float pivot[6];
float temp[6];

float x1;
float y1;
float z1;
float x2;
float y2;
float z2;
float d12;

int i;
int j;

if(from >= to)
    return;

pivot[0] = array[from+0];
pivot[1] = array[from+1];
pivot[2] = array[from+2];
pivot[3] = array[from+3];
pivot[4] = array[from+4];
pivot[5] = array[from+5];

i = from;

for(j = from+6; j <= to; j += 6){

    x1 = pivot[0] - array[j+0];
    y1 = pivot[1] - array[j+1];
    z1 = pivot[2] - array[j+2];
    x2 = pivot[3] - array[j+3];
    y2 = pivot[4] - array[j+4];
    z2 = pivot[5] - array[j+5];
    d12 = (x1*x1 + y1*y1 + z1*z1) + (x2*x2 + y2*y2 + z2*z2);
/*the sorting condition i am using is the regular euclidean norm*/
    if (d12 <= powR){
                i += 6;

                temp[0] = array[i+0];
                temp[1] = array[i+1];
                temp[2] = array[i+2];
                temp[3] = array[i+3];
                temp[4] = array[i+4];
                temp[5] = array[i+5];

                array[i+0] = array[j+0];
                array[i+1] = array[j+1];
                array[i+2] = array[j+2];
                array[i+3] = array[j+3];
                array[i+4] = array[j+4];
                array[i+5] = array[j+5];

                array[j+0] = temp[0];
                array[j+1] = temp[1];
                array[j+2] = temp[2];
                array[j+3] = temp[3];
                array[j+4] = temp[4];
                array[j+5] = temp[5];
    }
}

QuickSort(array, i+6, to, powR);
}

---

用于（j=from+6；j使用调试器或函数内部的一些打印来查看它崩溃的确切位置。我正在打印我能打印的所有内容，我唯一得到的信息是，当尝试使用起始索引2799228进入递归迭代时，它崩溃了。此外，我建议您使用较少的数字作为输入，并说明如何调用函数。但k和till是什么？我猜你忘了一些代码。
k
和
在哪里分配到了
呢？对不起，我的错，我忘了从调试版本中清理这些变量以使其更清晰。如果这是真的，为什么它会工作到某个递归？它不会在第一个递归时停止？因为你很幸运，直到某个时候。但是，你很幸运，贝卡使用“您肯定知道您有逻辑错误”。）
QuickSort(LORs, 0, 6000, powR);

    free(LORs);

for(j = from+6; j <= to; j += 6) {
    array[i+0] = array[j+0];
    array[i+1] = array[j+1];
    array[i+2] = array[j+2];
    array[i+3] = array[j+3];
    array[i+4] = array[j+4];
    array[i+5] = array[j+5];
}