1
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 <doc 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 <doc 80 81 <doc 82 83 84 85 86 87 <doc 88 89 <doc 90 91 92 <doc 93 94 95 96 <doc 97 98 99 <doc 100 101 102 103 104 <doc 105 106 107 <doc 108 109 110 <doc 111 112 113 <doc 114 115 116 <doc 117 118 119 <doc 120 121 122 <doc 123 124 125 <doc 126 127 128 <doc 129 130 131 <doc 132 133 134 <doc 135 136 137 <doc 138 139 140 <doc 141 142 143 <doc 144 145 146 <doc 147 148 149 <doc 150 151 152 153 154 155 <doc 156 157 158 159 <doc 160 161 162 <doc 163 164 165 166 167 168 169 <doc 170 171 172 <doc 173 174 175 176 <doc 177 178 179 180 <doc 181 182 183 184 <doc 185 186 187 <doc 188 189 190 191 192 <doc 193 194 195 <doc 196 197 198 <doc 199 200 201 <doc 202 203 204 205 206 207 208 209 210 211 212 213 214 |
/*
* Copyright (c) 2012, Michael Lehn, Klaus Pototzky * * All rights reserved. * * Redistribution and use in source and binary forms, with or without * modification, are permitted provided that the following conditions * are met: * * 1) Redistributions of source code must retain the above copyright * notice, this list of conditions and the following disclaimer. * 2) Redistributions in binary form must reproduce the above copyright * notice, this list of conditions and the following disclaimer in * the documentation and/or other materials provided with the * distribution. * 3) Neither the name of the FLENS development group nor the names of * its contributors may be used to endorse or promote products derived * from this software without specific prior written permission. * * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT * OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. */ #ifndef FLENS_STORAGE_BANDSTORAGE_BANDSTORAGE_H #define FLENS_STORAGE_BANDSTORAGE_BANDSTORAGE_H 1 #include <cxxblas/typedefs.h> #include <flens/storage/indexoptions.h> #include <flens/typedefs.h> namespace flens { template <typename T, typename I, typename A> class Array; template <typename T, typename I, typename A> class ArrayView; template <typename T, typename I, typename A> class ConstArrayView; template <typename T, StorageOrder Order, typename I, typename A> class BandStorageView; template <typename T, StorageOrder Order, typename I, typename A> class ConstBandStorageView; template <typename T, StorageOrder Order = ColMajor, typename I = IndexOptions<>, typename A = std::allocator<T> > class BandStorage { public: typedef T ElementType; typedef typename I::IndexType IndexType; typedef A Allocator; static const StorageOrder order = Order; static const IndexType defaultIndexBase = I::defaultIndexBase; typedef ConstBandStorageView<T, Order, I, A> ConstView; typedef BandStorageView<T, Order, I, A> View; typedef BandStorage NoView; typedef flens::ConstArrayView<T, I, A> ConstArrayView; typedef flens::ArrayView<T, I, A> ArrayView; typedef flens::Array<T, I, A> Array; BandStorage(); BandStorage(IndexType numRows, IndexType numCols, IndexType numSubDiags, IndexType numSuperDiags, IndexType firstIndex = I::defaultIndexBase, const ElementType &value = ElementType(), const Allocator &allocator = Allocator()); BandStorage(const BandStorage &rhs); template <typename RHS> BandStorage(const RHS &rhs); ~BandStorage(); //-- operators --------------------------------------------------------- const ElementType & operator()(IndexType row, IndexType col) const; ElementType & operator()(IndexType row, IndexType col); //-- methods ----------------------------------------------------------- IndexType firstRow() const; IndexType lastRow() const; IndexType firstCol() const; IndexType lastCol() const; IndexType firstIndex() const; IndexType lastIndex() const; IndexType numRows() const; IndexType numCols() const; IndexType dim() const; IndexType numSubDiags() const; IndexType numSuperDiags() const; IndexType leadingDimension() const; const ElementType * data() const; ElementType * data(); const Allocator & allocator() const; bool resize(IndexType numRows, IndexType numCols, IndexType numSubDiags, IndexType numSuperDiags, IndexType firstIndex = I::defaultIndexBase, const ElementType &value = ElementType()); template <typename FS> bool resize(const FS &rhs, const ElementType &value = ElementType()); bool fill(const ElementType &value = ElementType(0)); bool fillRandom(); void changeIndexBase(IndexType firstIndex); // view of bandstorage scheme as an array const ConstArrayView arrayView(IndexType firstViewIndex = I::defaultIndexBase) const; ArrayView arrayView(IndexType firstViewIndex = I::defaultIndexBase); // view of a diagonal const ConstArrayView viewDiag(IndexType diag, IndexType firstViewIndex = I::defaultIndexBase) const; ArrayView viewDiag(IndexType diag, IndexType firstViewIndex = I::defaultIndexBase); const ConstView viewDiags(IndexType fromDiag, IndexType toDiag) const; View viewDiags(IndexType fromDiag, IndexType toDiag); private: void _setIndexBase(IndexType firstIndex); void _raw_allocate(); void _allocate(const ElementType &value = ElementType()); void _release(); ElementType *_data; Allocator _allocator; IndexType _numRows, _numCols; IndexType _numSubDiags, _numSuperDiags; IndexType _firstIndex; }; } // namespace flens #endif // FLENS_STORAGE_BANDSTORAGE_BANDSTORAGE_H |