lacomplex.h

Go to the documentation of this file.

/*-*-c++-*-****************************************************************
 *                     lacomplex.h Helper file for complex numbers 
                       -------------------
 begin                : 2004-01-14
 copyright            : (C) 2004 by Christian Stimming
 email                : 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.

#ifndef _LACOMPLEX_H
#define _LACOMPLEX_H

// Include the version number defines
#include <laversion.h>

// Include the FORTRAN definitions from LAPACK++
#include <f2c.h>


// ////////////////////////////////////////////////////////////

#ifdef LA_COMPLEX_SUPPORT

typedef doublecomplex COMPLEX;

// As opposed to the FORTRAN "COMPLEX", include now the
// std::complex<double> type.
#include <complex>
// And finally include the la::complex<double> which is a copy of
// std::complex<double> with additional type conversions.
#include <lacomplex>

typedef la::complex<double> LaComplex;

inline std::ostream&
operator<<(std::ostream& __os, const COMPLEX& __x)
{
  return __os << LaComplex(__x);
}

inline bool operator==(const COMPLEX& _a, const COMPLEX& _b)
{
    return _a.r == _b.r && _a.i == _b.i;
}
inline bool operator!=(const COMPLEX& _a, const COMPLEX& _b)
{
    return _a.r != _b.r || _a.i != _b.i;
}

#endif /* LA_COMPLEX_SUPPORT */

// ////////////////////////////////////////////////////////////

#endif // _LACOMPLEX_H

---