[Lapackpp-cvs] CVS: lapackpp/src lasvd.cc,NONE,1.1 Makefile.am,1.5,1.6

[Christian S. <cs...@us...>]

Update of /cvsroot/lapackpp/lapackpp/src
In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv20813/src

Modified Files:
	Makefile.am 
Added Files:
	lasvd.cc 
Log Message:
2004-05-04  Christian Stimming <sti...@tu...>

        * include/lasvd.h, src/lasvd.cc: Initial work for SVD functions.



--- NEW FILE: lasvd.cc ---
//
//              LAPACK++ 1.1 Linear Algebra Package 1.1
//               University of Tennessee, Knoxvilee, TN.
//            Oak Ridge National Laboratory, Oak Ridge, TN.
//        Authors: J. J. Dongarra, E. Greaser, R. Pozo, D. Walker
//                 (C) 1992-1996 All Rights Reserved
//
//                             NOTICE
//
// Permission to use, copy, modify, and distribute this software and
// its documentation for any purpose and without fee is hereby granted
// provided that the above copyright notice appear in all copies and
// that both the copyright notice and this permission notice appear in
// supporting documentation.
//
// Neither the Institutions (University of Tennessee, and Oak Ridge National
// Laboratory) nor the Authors make any representations about the suitability 
// of this software for any purpose.  This software is provided ``as is'' 
// without express or implied warranty.
//
// LAPACK++ was funded in part by the U.S. Department of Energy, the
// National Science Foundation and the State of Tennessee.

#ifdef HAVE_CONFIG_H
# include <config.h>
#endif

#include <iostream>
#include "lapack.h"
#include "lapackc.h"
#include "lafnames.h"
#include LA_GEN_MAT_DOUBLE_H
#include LA_VECTOR_DOUBLE_H
#include LA_VECTOR_LONG_INT_H 
#include LA_GEN_MAT_COMPLEX_H
#include LA_VECTOR_COMPLEX_H
#include LA_EXCEPTION_H

#include LA_SOLVE_DOUBLE_H
#include LA_UTIL_H
#include "lasvd.h"

using std::min;
#if !LAPACK_USE_GCC3
#  ifndef min
#    define min(A,B) ((A) < (B) ? (A) : (B) )
#  endif
#endif


#ifdef LA_COMPLEX_SUPPORT
    
// General QR solver
//
//                          M x N              N x nrhs           M  x nrhs
//
void LaSVD_IP(LaGenMatComplex& A, LaVectorDouble &Sigma, LaGenMatComplex& U, LaGenMatComplex& VT )
{
#ifndef HPPA
    const char fname[] = "LaSVD_IP(LaGenMatComplex &A, &X, &B)";
#else
    char *fname = NULL;  // HP C++ does not support string initalization!
#endif

    // let's not worry about non-unit column strides for the moment
    if ( A.inc(0) != 1 || A.inc(1) != 1)
        throw(LaException(fname, "A is  non-contiguous."));

    char jobz = 'A';
    integer info;
    int M = A.size(0);
    int N = A.size(1);
    integer Ml = M;
    integer Nl = N;
    integer lda = A.inc(0) * A.gdim(0);
//     int nrhs = X.size(1);
//     integer nrhsl = nrhs;

    if (Sigma.size() != min(M,N))
	throw LaException(fname, "Sigma is not of correct size");
    if (U.size(0) != U.size(1) || U.size(0) != M)
	throw LaException(fname, "U is not of correct size");
    if (VT.size(0) != VT.size(1) || VT.size(0) != N)
	throw LaException(fname, "U is not of correct size");
        
    //#define MAX3(A,B,C) ((A)>(B) ? ((A)>(C) ? (A) : (C)) : ((B)>(C) ? (B) : (C)))
    //#define MIN(A,B) (A < B ? A : B )

    //std::cout << "Block size: " << LaEnvBlockSize("DGELSV", A) << std::endl;
    //int nb = 32;
    // int nb = LaEnvBlockSize("ZGESDD", A);
    //     integer W = std::min(M,N) + nb * MAX3(M, N, nrhs);

    //     LaVectorComplex work((int) W);         

    integer ldu = U.inc(0) * U.gdim(0);
    integer ldvt = VT.inc(0) * VT.gdim(0);

    F77NAME(zgesdd)(&jobz, &Ml, &Nl, &A(0,0), &lda, 
		    &Sigma(0), &U(0,0), &ldu, &VT(0,0), &ldvt, &info);

    if (info != 0) {
        throw(LaException(fname, "Internal error in LAPACK: zgesdd()"));
    }

}





#endif // LA_COMPLEX_SUPPORT

Index: Makefile.am
===================================================================
RCS file: /cvsroot/lapackpp/lapackpp/src/Makefile.am,v
retrieving revision 1.5
retrieving revision 1.6
diff -C2 -d --brief -r1.5 -r1.6
Files /tmp/cvsiQqDMN and /tmp/cvsuTxmfa differ