BLAS_Cpp_Types.hpp

// @HEADER
// ***********************************************************************
// 
// Moocho: Multi-functional Object-Oriented arCHitecture for Optimization
//                  Copyright (2003) Sandia Corporation
// 
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
// 
// 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.1 of the
// License, 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; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 
// 
// ***********************************************************************
// @HEADER
//
// Types for BLAS_Cpp namespace
//

#ifndef BLAS_CPP_TYPES_H
#define BLAS_CPP_TYPES_H

#include "Moocho_ConfigDefs.hpp"

namespace BLAS_Cpp {


typedef size_t size_type;


enum Side{
  left
  ,right
};
enum Transp{
  no_trans     
  ,trans       
  ,conj_trans  
};
enum Uplo {upper, lower};
inline Uplo operator!(Uplo uplo)
{ return uplo == upper ? lower : upper; }
enum Diag {unit, nonunit};



inline Transp bool_to_trans(bool return_trans) {
  return return_trans ? trans : no_trans;
}

inline bool trans_to_bool(Transp _trans) {
  return ( _trans == trans );
}

inline Transp operator!(Transp _trans)  // does not work with conj_trans
{ return _trans == no_trans ? trans : no_trans; }

inline Transp trans_not(Transp _trans)  // does not work with conj_trans
{ return _trans == no_trans ? trans : no_trans; }

inline Transp trans_trans(Transp _trans1, Transp _trans2) // does not work with conj_trans
{ return _trans1 == _trans2 ? no_trans : trans; }

inline const char*  trans_to_string(Transp _trans)
{
  return _trans == no_trans ? "no_trans" : "trans";
}

inline size_type rows(size_type rows, size_type cols
  , BLAS_Cpp::Transp _trans)
{
  return _trans == BLAS_Cpp::no_trans ? rows : cols;
}

inline size_type cols(size_type rows, size_type cols
  , BLAS_Cpp::Transp _trans)
{
  return _trans == BLAS_Cpp::no_trans ? cols : rows;
}



} // end namespace BLAS_Cpp

#endif // BLAS_CPP_TYPES_H