blas1pp.h

Go to the documentation of this file.

// -*-C++-*- 

// Copyright (C) 2004 
// Christian Stimming <stimming@tuhh.de>

// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2, or (at
// your option) any later version.

// This library 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 Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public License along
// with this library; see the file COPYING.  If not, write to the Free
// Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
// USA.

//      LAPACK++ (V. 1.1)
//      (C) 1992-1996 All Rights Reserved.

#ifndef _BLAS1_PP_H_
#define _BLAS1_PP_H_

#include "blas1.h"
#include "lafnames.h"
#include "arch.h" // needed for DLLIMPORT
#include LA_VECTOR_DOUBLE_H
#ifdef LA_COMPLEX_SUPPORT
# include LA_VECTOR_COMPLEX_H
#endif
#include <cmath>

#ifdef LA_COMPLEX_SUPPORT

DLLIMPORT
void Blas_Add_Mult(LaVectorComplex &dy, COMPLEX da, const LaVectorComplex &dx);

DLLIMPORT
void Blas_Mult(LaVectorComplex &dy, COMPLEX da, const LaVectorComplex &dx);

DLLIMPORT
COMPLEX Blas_U_Dot_Prod(const LaVectorComplex &cx, const LaVectorComplex &cy);

DLLIMPORT
COMPLEX Blas_H_Dot_Prod(const LaVectorComplex &cx, const LaVectorComplex &cy);

DLLIMPORT
void Blas_Copy(const LaVectorComplex &dx, LaVectorComplex &dy);

DLLIMPORT
void Blas_Scale(COMPLEX da, LaVectorComplex &dx);

DLLIMPORT
void Blas_Swap(LaVectorComplex &dx, LaVectorComplex &dy);
#endif // LA_COMPLEX_SUPPORT





DLLIMPORT
void Blas_Add_Mult(LaVectorDouble &dy, double da, const LaVectorDouble &dx);

DLLIMPORT
void Blas_Mult(LaVectorDouble &dy, double da, const LaVectorDouble &dx);

DLLIMPORT
double Blas_Dot_Prod(const LaVectorDouble &dx, const LaVectorDouble &dy);

DLLIMPORT
void Blas_Apply_Plane_Rot(LaVectorDouble &dx, LaVectorDouble &dy, 
                          double &c, double &s);

DLLIMPORT
void Blas_Gen_Plane_Rot(double &da, double &db, double &c, double &s);

DLLIMPORT
void Blas_Copy(const LaVectorDouble &dx, LaVectorDouble &dy);

DLLIMPORT
void Blas_Scale(double da, LaVectorDouble &dx);

DLLIMPORT
void Blas_Swap(LaVectorDouble &dx, LaVectorDouble &dy);





#ifdef LA_COMPLEX_SUPPORT

DLLIMPORT
double Blas_Norm1(const LaVectorComplex &dx);

DLLIMPORT
double Blas_Norm2(const LaVectorComplex &dx);

DLLIMPORT
int Blas_Index_Max(const LaVectorComplex &dx);

inline double Blas_Norm_Inf(const LaVectorComplex &x)
{   
    integer index = Blas_Index_Max(x);
    return la::abs(LaComplex(x(index)));
}
#endif // LA_COMPLEX_SUPPORT

DLLIMPORT
double Blas_Norm1(const LaVectorDouble &dx);

DLLIMPORT
double Blas_Norm2(const LaVectorDouble &dx);

DLLIMPORT
int Blas_Index_Max(const LaVectorDouble &dx);

inline double Blas_Norm_Inf(const LaVectorDouble &x)
{   
    integer index = Blas_Index_Max(x);
    return std::fabs(x(index));
}

#ifndef DOXYGEN_IGNORE
# ifdef LA_COMPLEX_SUPPORT

inline double Norm_Inf(const LaVectorComplex &A) { return Blas_Norm_Inf(A); }
# endif // LA_COMPLEX_SUPPORT

inline double Norm_Inf(const LaVectorDouble &A) { return Blas_Norm_Inf(A); }
#endif // DOXYGEN_IGNORE



#endif  
    /* _BLAS1_PP_H_ */

---