| 82 | } |
| 83 | |
| 84 | void PermuationMatrixGeneration(gnsstk::Matrix<double>& Pno, gnsstk::LUDecomp<double>& LUno) |
| 85 | { |
| 86 | //Identity Matrix Construction |
| 87 | int iCount = 0; |
| 88 | for (int i = 0; i < Pno.rows(); i++) |
| 89 | for (int j = 0; j < Pno.cols(); j++) |
| 90 | if ((j == iCount) && (i == iCount)) { |
| 91 | Pno(i,j) = 1; |
| 92 | iCount += 1;} |
| 93 | else |
| 94 | Pno(i,j) = 0; |
| 95 | |
| 96 | //Make Identity Martix the Permuation Matrix |
| 97 | for (int i = 0; i < LUno.Pivot.size(); i++) |
| 98 | Pno.swapRows(i, LUno.Pivot[i]); |
| 99 | } |
| 100 | |
| 101 | void LUDecompInitializer(void) |
| 102 | { |
no test coverage detected