Mercurial > octave
diff liboctave/CmplxSVD.cc @ 457:3d4b4f0fa5ba
[project @ 1994-06-06 00:33:33 by jwe]
Initial revision
author | jwe |
---|---|
date | Mon, 06 Jun 1994 00:33:51 +0000 |
parents | |
children | 4ecbfd3c3710 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/liboctave/CmplxSVD.cc Mon Jun 06 00:33:51 1994 +0000 @@ -0,0 +1,89 @@ +// -*- C++ -*- +/* + +Copyright (C) 1992, 1993, 1994 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 2, 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, write to the Free +Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. + +*/ + +#ifdef HAVE_CONFIG_H +#include "config.h" +#endif + +#if defined (__GNUG__) +#pragma implementation +#endif + +#include "CmplxSVD.h" +#include "mx-inlines.cc" +#include "f77-uscore.h" + +extern "C" +{ + int F77_FCN (zgesvd) (const char*, const char*, const int*, + const int*, Complex*, const int*, double*, + Complex*, const int*, Complex*, const int*, + Complex*, const int*, double*, int*, long, long); +} + +int +ComplexSVD::init (const ComplexMatrix& a) +{ + int info; + + int m = a.rows (); + int n = a.cols (); + + char jobu = 'A'; + char jobv = 'A'; + + Complex *tmp_data = dup (a.data (), a.length ()); + + int min_mn = m < n ? m : n; + int max_mn = m > n ? m : n; + + Complex *u = new Complex[m*m]; + double *s_vec = new double[min_mn]; + Complex *vt = new Complex[n*n]; + + int lwork = 2*min_mn + max_mn; + Complex *work = new Complex[lwork]; + + int lrwork = 5*max_mn; + double *rwork = new double[lrwork]; + + F77_FCN (zgesvd) (&jobu, &jobv, &m, &n, tmp_data, &m, s_vec, u, &m, + vt, &n, work, &lwork, rwork, &info, 1L, 1L); + + left_sm = ComplexMatrix (u, m, m); + sigma = DiagMatrix (s_vec, m, n); + ComplexMatrix vt_m (vt, n, n); + right_sm = ComplexMatrix (vt_m.hermitian ()); + + delete [] tmp_data; + delete [] work; + + return info; +} + +/* +;;; Local Variables: *** +;;; mode: C++ *** +;;; page-delimiter: "^/\\*" *** +;;; End: *** +*/