new file mode 100644
--- /dev/null
+++ b/liboctave/floatSVD.h
@@ -0,0 +1,92 @@
+/*
+
+Copyright (C) 1994, 1995, 1996, 1997, 2000, 2002, 2003, 2004, 2005,
+              2006, 2007 John W. Eaton
+
+This file is part of Octave.
+
+Octave is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 3 of the License, or (at your
+option) any later version.
+
+Octave is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with Octave; see the file COPYING.  If not, see
+<http://www.gnu.org/licenses/>.
+
+*/
+
+#if !defined (octave_FloatSVD_h)
+#define octave_FloatSVD_h 1
+
+#include <iostream>
+
+#include "fDiagMatrix.h"
+#include "fMatrix.h"
+#include "dbleSVD.h"
+
+class
+OCTAVE_API
+FloatSVD
+{
+public:
+
+  FloatSVD (void) : sigma (), left_sm (), right_sm () { }
+
+  FloatSVD (const FloatMatrix& a, SVD::type svd_type = SVD::std) { init (a, svd_type); }
+
+  FloatSVD (const FloatMatrix& a, octave_idx_type& info, SVD::type svd_type = SVD::std)
+    {
+      info = init (a, svd_type);
+    }
+
+  FloatSVD (const FloatSVD& a)
+    : type_computed (a.type_computed),
+      sigma (a.sigma), left_sm (a.left_sm), right_sm (a.right_sm) { }
+
+  FloatSVD& operator = (const FloatSVD& a)
+    {
+      if (this != &a)
+	{
+	  type_computed = a.type_computed;
+	  sigma = a.sigma;
+	  left_sm = a.left_sm;
+	  right_sm = a.right_sm;
+	}
+
+      return *this;
+    }
+
+  ~FloatSVD (void) { }
+
+  FloatDiagMatrix singular_values (void) const { return sigma; }
+
+  FloatMatrix left_singular_matrix (void) const;
+
+  FloatMatrix right_singular_matrix (void) const;
+
+  friend std::ostream&  operator << (std::ostream& os, const FloatSVD& a);
+
+private:
+
+  SVD::type type_computed;
+
+  FloatDiagMatrix sigma;
+  FloatMatrix left_sm;
+  FloatMatrix right_sm;
+
+  octave_idx_type init (const FloatMatrix& a, SVD::type svd_type = SVD::std);
+};
+
+#endif
+
+/*
+;;; Local Variables: ***
+;;; mode: C++ ***
+;;; End: ***
+*/