naming scheme    所有blas/lapack 函数名字格式都是 xyyzzz, 最后一个z也可以空
 x 代表数据类型   s 实型（real)   d  double precision 双精度浮点   c   complex 复型   z  double complex 双精度复型  对于一类函数统一使用首字母x，  比如 xgesv 就表示 sgesv, cgesv, dgesv, zgesv中的一个或全部
yy 代表矩阵类型，大部分函数适用both 实型 和 复型  bd  bidiagonal 双对角矩阵  di   diagonal  对角矩阵  gd  general band 带状矩阵  ge  general    gg   general matrices ?what's the difference with gg?  gt  general tridiagonal 通用三对角  hb  hermitian 带状矩阵  sb  symmetric band 对称带状矩阵  sp  symmetric, packed storage  对称矩阵，压缩存储  st   symmetric, tridiagonal 对称三对角矩阵  sy  symmetric  一般对称矩阵  tb   triangular band  三对角带状矩阵  tz trapezoidal 梯形矩阵?   more 当指一类函数（只是作用对象（矩阵）格式不同），采用 “xyy”＋zzz
zzz 描述函数的功效，包括（线性代数方程， 正交分解，最小二乘，对称／非对称特征值，奇异值分解，对称／非对称正定特征值问题）
computational routines  need some familiar with matrix computing
linear equations:   A X = B
general matric,  A = PLU
general band  A = LU
symmetric positive definite A = U^T U = L L^T
symmetric positive definite tridiagonal A= UDU^T = L D L^T