快速排序无限循环数值配方

我试图从MATLAB导出一些代码到MEX文件(a.la C)。对于我试图导出的代码段，我需要实现一种可以对2D数组进行排序的算法。打开一份数字食谱，我发现了基于索引的排序。我已经用零索引实现了这个函数，以一种简单的方式，稍后我将对其进行简化。目前的问题是，当我运行我的函数时，它将在传递某些随机数集时进入无限循环。例如:

0.8147,
0.9058,
0.127,
0.9134,
0.6324,
0.0975,
0.2785,
0.5469,
0.9575,
0.9649,
0.1576,
0.9706,
0.9572,
0.4854,
0.8003,
0.1419,
0.4218,
0.9157,
0.7922,
0.9595

这些数字都不重复。下面是我的C代码版本的数值配方。有什么问题吗?我在这上面花了3个小时。该算法为表现良好的随机数集提供适当的索引。

#define SWAP(a,b) itemp=(a);(a)=(b);(b)=itemp;
#define M 7
#define NSTACK 50
long* IndexSort(unsigned long int vectorLength, double* column)
{
unsigned long i, indxt, ir = vectorLength, itemp,j, k, l = 1;
int jstack = 0, *istack;
long *indx;
double a;
istack = (int*)malloc(NSTACK*sizeof(int));
//initalize output
indx = (long*)malloc(vectorLength*sizeof(long));
for (j = 0; j < vectorLength; j++)
{
    indx[j] = j;
}
// 
while (true)
{
    if (ir - l < M)
    {   
        for (j = l+1; j <= ir;j++)
        {
            indxt = indx[j-1];
            a = column[indxt];
            for (i = j-1; i >= 1; i--)
            {
                if (column[indx[i - 1]] <= a)
                {
                    break;
                }
                indx[i + 1 - 1] = indx[i - 1];
            }
            indx[i + 1 - 1] = indxt;
        }
        if (jstack == 0)
        {
            break;
        }
        ir = istack[jstack--];
        l = istack[jstack--];
    }
    else
    {   
        k = (l + ir) >> 1;
        SWAP(indx[k - 1], indx[l + 1 - 1])
        if (column[indx[l + 1 - 1]] > column[indx[ir - 1]]){
            SWAP(indx[l + 1 - 1], indx[ir - 1])
        }
        if (column[indx[l - 1]] > column[indx[ir - 1]]){
            SWAP(indx[l - 1], indx[ir - 1])
        }
        if (column[indx[l + 1 - 1]] > column[indx[l - 1]]){
            SWAP(indx[l + 1 - 1], indx[l - 1])
        }
        i = l + 1;
        j = ir;
        indxt = indx[l - 1];
        a = column[indxt];
        while (true)
        {
            do i++; while (column[indx[i - 1]] < a);
            do j--; while (column[indx[j - 1]] > a);
            if (j < i) break;
            SWAP(indx[i - 1], indx[j - 1])
        }
        indx[l - 1] = indx[j - 1];
        indx[j - 1] = indxt;
        jstack += 2;
        if (jstack > NSTACK) error("NSTACK too small for IndexSort");
        if (ir - i + 1 > j - 1){
            istack[jstack] = ir;
            istack[jstack - 1] = l;
            ir = j - 1;
        }
        else{
            istack[jstack] = j - 1;
            istack[jstack - 1] = l;
            l = i;
        }

    }
}
free(istack);
return indx;
}

编辑1:我在顶部添加了#define。作为对第一波评论的回应:

• 这个函数调用产生一个排序索引数组。我将使用这些值对2D数组进行排序。对不起，我应该说得更清楚一点。
• 至于使用调试器，目前是不可能的。该函数是从一个MEX包装器调用的。我将编写一个main函数来研究这个问题。

编辑2:

• 按照建议移动free命令，仍然没有变化。
• 更新问题说明我在论坛上看到问题，大多数人都崩溃了。我的问题实际上是一个无限循环，我不知道它在哪里开始。

提前感谢你的帮助，