EOF的C＃代码

Mc.Fish

把EOF翻译为了C＃

Fortran实在不怎么用，也不太方便，
以前都是把数据处理好，手工修改Fortran代码运行，
数据多了不好办，改成了C#方便一点。

调用：
                //Raw[nT,mSp]
                double[,] V1 = new double[mSpa, mSpa];//排序后的特征向量－－空间分布

                double[,] T = new double[nT, mSpa];      //             时间系数

                double[] H1 = new double[mSpa];          //             累计解释方差

                double[] D = new double[mSpa];            //             特征值

                WSAA.Base.EOF.SEOF(mSpa, nT, WSAA.Base.EOF.ContralParameter.NORMALIZED, Raw, ref V1, ref T, ref H1, ref D);

有个问题：

有时nT和 mSpa差别大，进行时空交换

https://wenku.baidu.com/view/e64ba4691eb91a37f1115c0b.html        P48

里提到 这样不改变特征根

可是一个特征值是  nT个， 一个特征值是 mSpa个，怎么会不改变呢

建议 大家把常用的统计程序写成C＃ 共享，建立一个库，方便大家使用。

C＃代码在附件

funvers

楼主的C#源程序还在吗？能提供一份吗？490110000@qq.coom

tulalang

楼主比较过C#和fortran运行速度么

funvers

tulalang 发表于 2022-6-15 22:46
楼主比较过C#和fortran运行速度么

C#1秒内完成分解输出最后的矩阵

funvers

楼主，我测试了这个EOF的C#代码。在雅克比方法里面，有一段内容可能翻译错了。
if (KM1 < 0)
                            {
                                for (int j = 0; j < KM1; j++)
                                {
                                    G = A[j, k];
                                    H = A[j, L];
                                    A[j, k] = c * G - s * H;
                                    A[j, L] = s * G + c * H;
                                }
                            }

                            int L1 = L - 1;
                            if (L1 - k1 <= 0)
                            {
                                for (int j = k1; j < L1; j++)
                                {
                                    G = A[k, L];
                                    H = A[j, L];
                                    A[k, j] = c * G - s * H;
                                    A[j, L] = s * G + c * H;
                                }
                            }

                            L1 = L + 1;
if(KM1<0)后面永远不会执行，同时，这份代码分解的结果还会出现贡献率为负数的问题。