求矩阵行列式的值（基于LU分解法）

算法名称： 求矩阵行列式的值（基于LU 分解法）

 

算法描述：

       对矩阵进行LU 分解，显然，分解后的矩阵对角线上元素的乘积即为原始矩阵行列式的值，

       注： 由于在进行 LU 分解时可能会进行行交换的情况，而每次行交换都会带来行列式的符号变化，所以我们要记录行交换次数的奇偶性，奇数则符号改变，偶数则符号不变，请密切注意下面程序中出现的用来标记变换次数奇偶性的变量 parity 。

大矩阵情况的应对策略：

     对实际出现的大矩阵，很可能出现行列式值上溢或下溢，超出了计算机所能表示的浮点数的范围。在这种情况下，可以修改下述程序中的循环部分，例如除以以十的幂次并单独保存该比例因子，或求出各乘数值的对数和，并单独保存其符号。

 

运行示例：

Origin matrix:

 | 0.0 2.0 0.0 1.0 |

 | 2.0 2.0 3.0 2.0 |

 | 4.0 -3.0 0.0 1.0 |

 | 6.0 1.0 -6.0 -5.0 |

-----------------------------------------------

After LU decomposition:

 | 6.0 1.0 -6.0 -5.0 |

 | 0.0 2.0 0.0 1.0 |

 | 0.3333333333333333 0.8333333333333334 5.0 2.833333333333333 |

 | 0.6666666666666666 -1.8333333333333333 0.8 3.8999999999999995 |

-----------------------------------------------

Its determinate value:-233.99999999999997

示例程序：

package  com.nc4nr.chapter02.det;

public   class  Det  ... {

    double[][] a = ...{
            ...{ 0.0, 2.0, 0.0, 1.0},
            ...{ 2.0, 2.0, 3.0, 2.0},
            ...{ 4.0, -3.0, 0.0, 1.0},