eigsrt.cpp

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

void  eigsrt(double d[], double v[4][4], int n)
{
	double p;
	int i,j,k;
    for (i = 1; i<=n - 1; i++)
	{
        k = i;
        p = d[i];
        for (j = i + 1; j<=n; j++)
		{
            if (d[j] >= p)
			{
				k = j;
                p = d[j];
			}
           
        }
        if (k != i )
		{
            d[k] = d[i];
            d[i] = p;
            for (j = 1; j<=n; j++)
			{
                p = v[j][ i];
                v[j][i] = v[j][k];
                v[j][k] = p;
            }
		}
    }
}