sort3.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 30 行

void sort3(int n, double ra[], double rb[], double rc[], double wksp[], int iwksp[])
{
	int j;
    indexx(n, ra, iwksp);
    for (j=1; j<=n; j++)
	{
        wksp[j] = ra[j];
    }
    for (j=1; j<=n; j++)
	{
        ra[j] = wksp[iwksp[j]];
    }
    for (j=1; j<=n; j++)
	{
        wksp[j] = rb[j];
    }
    for (j=1; j<=n; j++)
	{
        rb[j] = wksp[iwksp[j]];
    }
    for (j=1; j<=n; j++)
	{
        wksp[j] = rc[j];
    }
    for (j=1; j<=n; j++)
	{
        rc[j] = wksp[iwksp[j]];
    }
}