BelosMultiVecTraits.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//                 Belos: Block Linear Solvers Package
//
// Copyright 2004-2016 NTESS and the Belos contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
//
#ifndef BELOS_MULTI_VEC_TRAITS_HPP
#define BELOS_MULTI_VEC_TRAITS_HPP
 
 
#include "BelosTypes.hpp"
#include "BelosStubTsqrAdapter.hpp"
#include "Teuchos_Range1D.hpp"
#include "Teuchos_RCP.hpp"
 
// This is included for backwards compatibility
// for all codes using the default template parameter
// DM = DefaultDenseMatrix<int,ScalarType>
#include "BelosDenseMatTraits.hpp" 
 
namespace Belos {
 
  template<class ScalarType, class MV, class DM>
  struct UndefinedMultiVecTraits
  {
    static inline ScalarType notDefined() {
      return MV::this_type_is_missing_a_specialization();
    }
  };
 
  template<class ScalarType, class MV, class DM = DefaultDenseMatrix<int, ScalarType>>
  class MultiVecTraits {
  public:
 
 
    static Teuchos::RCP<MV> Clone( const MV& mv, const int numvecs )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneCopy( const MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneCopy( const MV& mv, const std::vector<int>& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneCopy( const MV& mv, const Teuchos::Range1D& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneViewNonConst( MV& mv, const std::vector<int>& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneViewNonConst( MV& mv, const Teuchos::Range1D& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<const MV> CloneView( const MV& mv, const std::vector<int>& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
    static Teuchos::RCP<MV> CloneView( MV& mv, const Teuchos::Range1D& index )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return Teuchos::null; }
 
 
 
 
    static ptrdiff_t GetGlobalLength( const MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return 0; }
 
    static int GetNumberVecs( const MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return 0; }
 
    static bool HasConstantStride( const MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); return false; }
 
 
 
 
    static void MvTimesMatAddMv( const ScalarType alpha, const MV& A,
                                 const DM& B,
                                 const ScalarType beta, MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvAddMv( const ScalarType alpha, const MV& A, const ScalarType beta, const MV& B, MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvScale ( MV& mv, const ScalarType alpha )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvScale ( MV& mv, const std::vector<ScalarType>& alpha )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvTransMv( const ScalarType alpha, const MV& A, const MV& mv, DM& B)
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvDot ( const MV& mv, const MV& A, std::vector<ScalarType> &b)
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
 
 
    static void MvNorm( const MV& mv, std::vector<typename Teuchos::ScalarTraits<ScalarType>::magnitudeType>& normvec, NormType type = TwoNorm )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
 
 
 
    static void SetBlock( const MV& A, const std::vector<int>& index, MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void SetBlock( const MV& A, const Teuchos::Range1D& index, MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void Assign( const MV& A, MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvRandom( MV& mv )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
    static void MvInit( MV& mv, const ScalarType alpha = Teuchos::ScalarTraits<ScalarType>::zero() )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
 
 
 
    static void MvPrint( const MV& mv, std::ostream& os )
    { UndefinedMultiVecTraits<ScalarType, MV, DM>::notDefined(); }
 
 
#ifdef HAVE_BELOS_TSQR
    typedef Belos::details::StubTsqrAdapter<MV, DM> tsqr_adaptor_type;
#endif // HAVE_BELOS_TSQR
  };
 
} // namespace Belos
 
#endif // BELOS_MULTI_VEC_TRAITS_HPP