这是用于线性方程组求解的ILUK预处理算法的实现。在VC++编译通过。矩阵采用压缩稀疏行格式存储（CSR），采用如下结构存储：struct Distmatrix {double **ma int **ja,dimension,*nnzrow } 很容易移植到自己定义数值计算软件包中。经本人测试计算效率比Fortran写的高很多（比如与Sparskit2比较）。
