GSL中的稀疏线性代数

稀疏线性代数

    本章描述求解稀疏线性方程组的函数。本库提供了直接操作gsl_spmatrix和gsl_vector对象的线性代数函数。

本章描述的函数声明在头文件gsl_splinalgh中。

45.1 概述

    本章主要讨论如下线性系统的解，

                                                       Ax=b

    其中A是一个一般矩形n乘n非奇异稀疏矩阵，x是一个未知的n乘1向量，b是一个给定的n乘1右边的向量。求解稀疏线性系统的方法有很多，大体上可分为直接或迭代两类。直接方法包括LU分解和QR分解，迭代方法从向量x的初始猜测开始，通过迭代更新猜测直到收敛。GSL目前不提供任何直接稀疏解算器。