0

以下代码是计算向量的绑定秩的函数(性能关键):

//The function here is to compute tied-ranks: answers.com/topic/tied-rank
mergeSort(x,inds,ci);
//mergeSort(): to sort vector x of length ci, also returns keys (inds) of x.

int tj=0;
double xi=x[0];

for (int j = 1; j < ci; ++j)
{
    if (x[j] > xi)
    {
        double rankvalue = 0.5 * (j - 1 + tj);

        for (int k = tj; k < j; ++k)
        {
            ranks[inds[k]] = rankvalue;
        };

        tj = j;
        xi = x[j];
    };      
};

double rankvalue = 0.5 * (ci - 1 + tj);

for (int k = tj; k < ci; ++k)
{
    ranks[inds[k]] = rankvalue;
};

问题是,假定的性能瓶颈 mergeSort(),即 O(NlogN) 比代码的其他部分(即 O(N))快几倍,这表明其他部分还有很大的改进空间代码,有什么建议吗?

4

1 回答 1

0

该算法似乎具有二次行为:如果x[0]是序列中的最大值,则保持tj内部迭代。您的意思是使用and吗?0cix[inds[0]]x[inds[j]]

于 2012-12-12T22:40:35.700 回答