Sparse matrix solver

为写Photography的作业，要用到sparse matrix solver。原以为OpenCV中有这个solver，毕竟里面定义了CvSparseMat这个结构，还有cvSolver()等函数，我还真的相信了。花了些时间把OpenCV里的Matrix运算的看了一遍，然后也写好了代码，结果cvSolver()一直报错，很认真地检查code，改了不少matrix运算的策略还是不行。上yahoo 的group搜了一把，才发现opencv里的CvSparseMat就是个噱头，只是用来存数据的，只有histogram这个函数能用sparse matrix，其他函数都不能用。我晕。。。

想着IPP应该有sparse matrix solver的，想着自己对IPP还比较熟，OpenCV不行还有IPP。结果down下来也没有，这样就只有MKL这唯一的希望了。不过真地不想用Intel的这些库，虽然效率高，但是用起来真地不太方便，而且把数据从opencv的格式导成MKL/IPP的格式也挺麻烦，我就放弃了。最后还是选用了TAUCS这个库，现在用起来还不错，它可以同时解多个方程，如果Ax=b中A不变。这样我就可以对3个通道同时求解了，cool。