The class Lapack_svd provides an algorithm to solve in the least square sense a linear system with a singular value decomposition. The field type is double.
#include <CGAL/Lapack/Linear_algebra_lapack.h>