1// This file is part of Eigen, a lightweight C++ template library
2// for linear algebra.
3//
4// Copyright (C) 2008-2009 Gael Guennebaud <gael.guennebaud@inria.fr>
5//
6// This Source Code Form is subject to the terms of the Mozilla
7// Public License v. 2.0. If a copy of the MPL was not distributed
8// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
9
10#ifndef EIGEN_SPARSETRANSPOSE_H
11#define EIGEN_SPARSETRANSPOSE_H
12
13namespace Eigen {
14
15template<typename MatrixType> class TransposeImpl<MatrixType,Sparse>
16  : public SparseMatrixBase<Transpose<MatrixType> >
17{
18    typedef typename internal::remove_all<typename MatrixType::Nested>::type _MatrixTypeNested;
19  public:
20
21    EIGEN_SPARSE_PUBLIC_INTERFACE(Transpose<MatrixType> )
22
23    class InnerIterator;
24    class ReverseInnerIterator;
25
26    inline Index nonZeros() const { return derived().nestedExpression().nonZeros(); }
27};
28
29// NOTE: VC10 and VC11 trigger an ICE if don't put typename TransposeImpl<MatrixType,Sparse>:: in front of Index,
30// a typedef typename TransposeImpl<MatrixType,Sparse>::Index Index;
31// does not fix the issue.
32// An alternative is to define the nested class in the parent class itself.
33template<typename MatrixType> class TransposeImpl<MatrixType,Sparse>::InnerIterator
34  : public _MatrixTypeNested::InnerIterator
35{
36    typedef typename _MatrixTypeNested::InnerIterator Base;
37    typedef typename TransposeImpl::Index Index;
38  public:
39
40    EIGEN_STRONG_INLINE InnerIterator(const TransposeImpl& trans, typename TransposeImpl<MatrixType,Sparse>::Index outer)
41      : Base(trans.derived().nestedExpression(), outer)
42    {}
43    typename TransposeImpl<MatrixType,Sparse>::Index row() const { return Base::col(); }
44    typename TransposeImpl<MatrixType,Sparse>::Index col() const { return Base::row(); }
45};
46
47template<typename MatrixType> class TransposeImpl<MatrixType,Sparse>::ReverseInnerIterator
48  : public _MatrixTypeNested::ReverseInnerIterator
49{
50    typedef typename _MatrixTypeNested::ReverseInnerIterator Base;
51    typedef typename TransposeImpl::Index Index;
52  public:
53
54    EIGEN_STRONG_INLINE ReverseInnerIterator(const TransposeImpl& xpr, typename TransposeImpl<MatrixType,Sparse>::Index outer)
55      : Base(xpr.derived().nestedExpression(), outer)
56    {}
57    typename TransposeImpl<MatrixType,Sparse>::Index row() const { return Base::col(); }
58    typename TransposeImpl<MatrixType,Sparse>::Index col() const { return Base::row(); }
59};
60
61} // end namespace Eigen
62
63#endif // EIGEN_SPARSETRANSPOSE_H
64