libigl_Martin/include/igl/svd.h

//THIS DEPENDS ON BLAS. WHY? WHY NOT EIGEN?

//#ifndef IGL_SVD_H
//#define IGL_SVD_H
//namespace igl
//{
//  // Compute 3x3 SVD using lapack's dgesdd_ function
//  //
//  // Input:
//  //   a  pointer to 3x3 matrix in COLUMN MAJOR order
//  // Outputs:
//  //   u  pointer to 3x3 matrix in COLUMN MAJOR order
//  //   s  pointer to 3 vector 
//  //   vt  pointer to 3x3 matrix in COLUMN MAJOR order, or think of this as v in
//  //     row-major order
//  // Returns true on success, false on failure
//  // 
//  // Known bugs: This only compiles on Mac and depends on Lapack rather than eigen
//  bool svd3x3(double * a, double * u, double * s, double * vt);
//};
//
//#ifdef IGL_HEADER_ONLY
//#include "svd.cpp"
//#endif
//
//#endif