最近仿08年TOG上一篇骨架提取的文章Skeleton extraction by mesh contraction，其中涉及到线性方程组的最小二乘解问题，即Ax = b。
最开始使用了Armadillo库进行求解，程序写完后发现矩阵A的规模与顶点数的平方成正比，不使用稀疏矩阵的话只能计算很小的模型，但Armadillo没有提供稀疏矩阵模块。听说Eigen库有稀疏矩阵模块，又查了下Eigen库，但是发现Eigen库的稀疏矩阵求解线性方程组的功能只能用于A为方阵的情况。最后考虑用matlab结合vs2010的方式实现。
Armadillo和Eigen都是很绿色的线性代数库，都是以泛型编程的方式实现。Armadillo非常简练，文档也小巧精悍，上手很快，底层依赖lapack和blas库，环境配置方式写在了之前的一篇文档里http://www.cnblogs.com/youthlion/archive/2012/05/15/2501465.html。Eigen库更加重量级一些，功能更加全面，文档详细，不依赖于任何其他底层库。环境配置也很简单，不多说了。下面主要记录一下vs2010和matlab混用的方法。