44 #ifndef RTOPPACK_TYPES_HPP 45 #define RTOPPACK_TYPES_HPP 101 {
public:
UnknownError(
const std::string& what_arg) : std::logic_error(what_arg) {}};
104 {
public:
InvalidUsage(
const std::string& what_arg) : std::logic_error(what_arg) {}};
107 {
public:
InvalidNumVecs(
const std::string& what_arg) : std::logic_error(what_arg) {}};
143 template<
class Scalar>
156 {
initialize(globalOffset_in, subDim_in, values_in, stride_in); }
172 subDim_in*std::abs(Teuchos::as<int>(stride_in)) - 1 <= values_in.
upperOffset());
246 template<
class Scalar>
287 template<
class Scalar>
294 for(
int i = 0; i < sv.
subDim(); ++i ) {
304 template<
class Scalar>
305 std::ostream& operator<<(std::ostream &out, const ConstSubVectorView<Scalar> &sv)
309 <<
"globalOffset="<<sv.globalOffset()
310 <<
",subDim="<<sv.subDim()
311 <<
",values="<<sv.values()
312 <<
",stride="<<sv.stride()
344 template<
class Scalar>
361 initialize(globalOffset_in, subDim_in, colOffset_in, numSubCols_in, values_in,
472 template<
class Scalar>
482 Teuchos::arcp<Scalar>(numRows_in*numCols_in), numRows_in)
491 colOffset_in, numSubCols_in, values_in, leadingDim_in)
505 subDim_in, colOffset_in, numSubCols_in, values_in, leadingDim_in);
510 return Teuchos::arcp_const_cast<Scalar>(
534 template<
class Scalar>
544 (*msmv)(i,j) = smv(i,j);
565 template <
class Scalar,
class ConcreteObj>
572 {
return Scalar::this_type_is_missing_a_specialization(); }
575 {
return Scalar::this_type_is_missing_a_specialization(); }
578 {
return Scalar::this_type_is_missing_a_specialization(); }
587 Scalar::this_type_is_missing_a_specialization(obj);
597 *obj = Scalar::this_type_is_missing_a_specialization();
606 template <
class Scalar>
625 assertInput(primitiveObjs, indexObjs, charObjs);
626 primitiveObjs[0] = obj;
636 assertInput(primitiveObjs, indexObjs, charObjs);
637 *obj = primitiveObjs[0];
648 || charObjs.
size()!=0 );
655 template <
class Scalar>
676 assertInput(primitiveObjs, indexObjs, charObjs);
687 assertInput(primitiveObjs, indexObjs, charObjs);
699 || charObjs.
size()!=0 );
705 #if defined(HAVE_COMPLEX) && defined(HAVE_TEUCHOS_COMPLEX) 711 template <
class Scalar>
712 class PrimitiveTypeTraits<
std::complex<Scalar>, std::complex<Scalar> > {
715 typedef PrimitiveTypeTraits<Scalar,Scalar> ScalarPrimitiveTypeTraits;
717 typedef typename ScalarPrimitiveTypeTraits::primitiveType
primitiveType;
720 {
return 2*ScalarPrimitiveTypeTraits::numPrimitiveObjs(); }
727 const std::complex<Scalar> &obj,
728 const ArrayView<primitiveType> &primitiveObjs,
729 const ArrayView<index_type> &indexObjs,
730 const ArrayView<char> &charObjs
734 const int numScalarPrimitiveObjs =
735 ScalarPrimitiveTypeTraits::numPrimitiveObjs();
736 assertInput(primitiveObjs, indexObjs, charObjs);
737 ScalarPrimitiveTypeTraits::extractPrimitiveObjs(
738 obj.real(), primitiveObjs(0,numScalarPrimitiveObjs), null, null );
739 ScalarPrimitiveTypeTraits::extractPrimitiveObjs(
740 obj.imag(), primitiveObjs(numScalarPrimitiveObjs,numScalarPrimitiveObjs), null, null );
744 const ArrayView<const primitiveType> &primitiveObjs,
745 const ArrayView<const index_type> &indexObjs,
746 const ArrayView<const char> &charObjs,
747 const Ptr<std::complex<Scalar> > &obj
751 using Teuchos::outArg;
752 assertInput(primitiveObjs, indexObjs, charObjs);
753 const int numScalarPrimitiveObjs =
754 ScalarPrimitiveTypeTraits::numPrimitiveObjs();
756 ScalarPrimitiveTypeTraits::loadPrimitiveObjs(
757 primitiveObjs(0,numScalarPrimitiveObjs), null, null,
759 ScalarPrimitiveTypeTraits::loadPrimitiveObjs(
760 primitiveObjs(numScalarPrimitiveObjs,numScalarPrimitiveObjs), null, null,
762 *obj = std::complex<Scalar>( real, imag );
765 static void assertInput(
766 const ArrayView<const primitiveType> &primitiveObjs,
767 const ArrayView<const index_type> &indexObjs,
768 const ArrayView<const char> &charObjs
773 primitiveObjs.size()!=2*ScalarPrimitiveTypeTraits::numPrimitiveObjs()
774 || indexObjs.size()!=0
775 || charObjs.size()!=0 );
784 template <
class Scalar>
785 class PrimitiveTypeTraits<
std::complex<Scalar>, Scalar> {
788 typedef PrimitiveTypeTraits<Scalar,Scalar> ScalarPrimitiveTypeTraits;
790 typedef typename ScalarPrimitiveTypeTraits::primitiveType
primitiveType;
793 {
return ScalarPrimitiveTypeTraits::numPrimitiveObjs(); }
801 const ArrayView<primitiveType> &primitiveObjs,
802 const ArrayView<index_type> &indexObjs,
803 const ArrayView<char> &charObjs
807 assertInput(primitiveObjs, indexObjs, charObjs);
808 ScalarPrimitiveTypeTraits::extractPrimitiveObjs(
809 obj, primitiveObjs, null, null );
813 const ArrayView<const primitiveType> &primitiveObjs,
814 const ArrayView<const index_type> &indexObjs,
815 const ArrayView<const char> &charObjs,
816 const Ptr<Scalar > &obj
820 assertInput(primitiveObjs, indexObjs, charObjs);
821 ScalarPrimitiveTypeTraits::loadPrimitiveObjs(
822 primitiveObjs, null, null, obj );
825 static void assertInput(
826 const ArrayView<const primitiveType> &primitiveObjs,
827 const ArrayView<const index_type> &indexObjs,
828 const ArrayView<const char> &charObjs
833 primitiveObjs.size()!=ScalarPrimitiveTypeTraits::numPrimitiveObjs()
834 || indexObjs.size()!=0
835 || charObjs.size()!=0 );
841 #endif // defined(HAVE_COMPLEX) && defined(HAVE_TEUCHOS_COMPLEX) 851 template<
class Scalar>
class RTOpT;
857 #endif // RTOPPACK_TYPES_HPP ConstSubMultiVectorView(Ordinal globalOffset_in, Ordinal subDim_in, Ordinal colOffset_in, Ordinal numSubCols_in, const ArrayRCP< const Scalar > &values_in, Ordinal leadingDim_in)
ArrayRCP< const Scalar > values_
bool is_null(const boost::shared_ptr< T > &p)
ConstSubMultiVectorView()
static int numPrimitiveObjs()
ConstSubVectorView< Scalar > col(const Ordinal j) const
Return a ConstSubVectorView view of the jth sub-column (Preconditions: values()!=NULL && (0<=j<numSub...
const Scalar & operator()(Ordinal i) const
Zero-based indexing (Preconditions: values()!=NULL && (0 <= i < subDim())).
static void extractPrimitiveObjs(const Scalar &obj, const ArrayView< primitiveType > &primitiveObjs, const ArrayView< index_type > &indexObjs, const ArrayView< char > &charObjs)
Class for a changeable sub-vector.
SubMultiVectorView(const SubMultiVectorView< Scalar > &smv)
void assign_entries(const Ptr< const SubVectorView< Scalar > > &msv, const ConstSubVectorView< Scalar > &sv)
static void loadPrimitiveObjs(const ArrayView< const primitiveType > &primitiveObjs, const ArrayView< const index_type > &indexObjs, const ArrayView< const char > &charObjs, const Ptr< Scalar > &obj)
Teuchos_Ordinal index_type
IncompatibleVecs(const std::string &what_arg)
ConstSubVectorView(Ordinal globalOffset_in, Ordinal subDim_in, const ArrayRCP< const Scalar > &values_in, ptrdiff_t stride_in)
Ordinal globalOffset() const
static int numIndexObjs()
SubMultiVectorView(Ordinal globalOffset_in, Ordinal subDim_in, Ordinal colOffset_in, Ordinal numSubCols_in, const ArrayRCP< Scalar > &values_in, Ordinal leadingDim_in)
Class for a non-changeable sub-vector.
SubVectorView(const ArrayRCP< Scalar > &values_in)
InvalidNumVecs(const std::string &what_arg)
SubMultiVectorView(Ordinal numRows_in, Ordinal numCols_in)
#define TEUCHOS_ASSERT_IN_RANGE_UPPER_EXCLUSIVE(index, lower_inclusive, upper_exclusive)
const ArrayRCP< const Scalar > values() const
static int numIndexObjs()
SubVectorView(Ordinal subDim_in)
TEUCHOS_ORDINAL_TYPE Teuchos_Ordinal
Specialization where the scalar type is the same as the concrete object type.
static void assertInput(const ArrayView< const primitiveType > &primitiveObjs, const ArrayView< const index_type > &indexObjs, const ArrayView< const char > &charObjs)
A templated traits class for decomposing object into an array of primitive objects.
UnknownError(const std::string &what_arg)
static int numPrimitiveObjs()
SubVectorView(const SubVectorView< Scalar > &sv)
ArrayRCP< T > persistingView(size_type lowerOffset, size_type size) const
Class for a changeable sub-vector.
static void assertInput(const ArrayView< const primitiveType > &primitiveObjs, const ArrayView< const index_type > &indexObjs, const ArrayView< const char > &charObjs)
SubVectorView(Ordinal globalOffset_in, Ordinal subDim_in, const ArrayRCP< Scalar > &values_in, ptrdiff_t stride_in)
void initialize(Ordinal globalOffset_in, Ordinal subDim_in, Ordinal colOffset_in, Ordinal numSubCols_in, const ArrayRCP< Scalar > &values_in, Ordinal leadingDim_in)
void initialize(Ordinal globalOffset_in, Ordinal subDim_in, const ArrayRCP< Scalar > &values_in, ptrdiff_t stride_in)
SubVectorView< Scalar > col(const Ordinal j) const
Return a SubVectorView view of the jth sub-column (Preconditions: values()!=NULL && && (0<=j<numSubCo...
size_type upperOffset() const
const Scalar & operator[](Ordinal i) const
Zero-based indexing (Preconditions: values()!=NULL && (0 <= i < subDim())).
static void extractPrimitiveObjs(const Scalar &obj, const ArrayView< primitiveType > &primitiveObjs, const ArrayView< index_type > &indexObjs, const ArrayView< char > &charObjs)
void setGlobalOffset(Ordinal globalOffset_in)
ArrayRCP< const Scalar > values_
InvalidNumTargVecs(const std::string &what_arg)
ConstSubVectorView(const ConstSubVectorView< Scalar > &sv)
const ArrayRCP< const Scalar >::iterator valuesBegin() const
static int numIndexObjs()
void setGlobalOffset(Ordinal globalOffset_in)
ConstSubVectorView(const ArrayRCP< const Scalar > &values_in)
const Scalar & operator()(Ordinal i, Ordinal j) const
Zero-based indexing (Preconditions: values()!=NULL && (0<=i<subDim()) && (0<=j< numSubCols()).
basic_FancyOStream< char > FancyOStream
Ordinal numSubCols() const
const ArrayRCP< Scalar > values() const
static int numPrimitiveObjs()
Scalar & operator()(Ordinal i, Ordinal j) const
Zero-based indexing (Preconditions: values()!=NULL && (0<=i< subDim()) && (0<=j<numSubCols()).
InvalidUsage(const std::string &what_arg)
Class for a non-changeable sub-multi-vector (submatrix).
Scalar & operator()(Ordinal i) const
Zero-based indexing (Preconditions: values()!=NULL && (0 <= i < subDim())).
ScalarPrimitiveTypeTraits::primitiveType primitiveType
void initialize(Ordinal globalOffset_in, Ordinal subDim_in, Ordinal colOffset_in, Ordinal numSubCols_in, const ArrayRCP< const Scalar > &values_in, Ordinal leadingDim_in)
static void loadPrimitiveObjs(const ArrayView< const primitiveType > &primitiveObjs, const ArrayView< const index_type > &indexObjs, const ArrayView< const char > &charObjs, const Ptr< index_type > &obj)
static void loadPrimitiveObjs(const ArrayView< const primitiveType > &primitiveObjs, const ArrayView< const index_type > &indexObjs, const ArrayView< const char > &charObjs, const Ptr< Scalar > &obj)
const ArrayRCP< Scalar > values() const
Ordinal globalOffset() const
PrimitiveTypeTraits< Scalar, Scalar > ScalarPrimitiveTypeTraits
Scalar & operator[](Ordinal i) const
Zero-based indexing (Preconditions: values()!=NULL && (0 <= i < subDim())).
#define TEUCHOS_ASSERT(assertion_test)
Ordinal colOffset() const
#define TEUCHOS_ASSERT_EQUALITY(val1, val2)
static void extractPrimitiveObjs(const index_type &obj, const ArrayView< primitiveType > &primitiveObjs, const ArrayView< index_type > &indexObjs, const ArrayView< char > &charObjs)
size_type lowerOffset() const
IncompatibleReductObj(const std::string &what_arg)
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
ConstSubMultiVectorView(const ConstSubMultiVectorView< Scalar > &smv)
Ordinal leadingDim() const
const ArrayRCP< const Scalar > values() const
void initialize(Ordinal globalOffset_in, Ordinal subDim_in, const ArrayRCP< const Scalar > &values_in, ptrdiff_t stride_in)