最近仿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混用的方法。