| 169 | |
| 170 | template <typename T> |
| 171 | static T |
| 172 | getElementPacked( |
| 173 | clblasOrder order, |
| 174 | clblasTranspose trans, |
| 175 | clblasUplo uplo, |
| 176 | size_t row, |
| 177 | size_t column, |
| 178 | T *A, |
| 179 | size_t rows) |
| 180 | { |
| 181 | // Needed for Macro : testDG.h |
| 182 | int vectorLength = 1; |
| 183 | const T* data = A; |
| 184 | clblasUplo fUplo = (trans == clblasNoTrans) ? uplo : ((uplo == clblasUpper) ? clblasLower : clblasUpper); |
| 185 | |
| 186 | if(fUplo == clblasLower) //Should not access elements out of bounds. |
| 187 | { |
| 188 | if (column > row) |
| 189 | return FNAN<T>(); |
| 190 | } |
| 191 | else |
| 192 | { |
| 193 | if (column < row) |
| 194 | return FNAN<T>(); |
| 195 | } |
| 196 | switch (order) { |
| 197 | case clblasRowMajor: |
| 198 | if (fUplo == clblasLower) |
| 199 | { |
| 200 | return *RMLPacked(row, column); |
| 201 | } |
| 202 | else { |
| 203 | return *RMUPacked(row, column); |
| 204 | } |
| 205 | break; |
| 206 | case clblasColumnMajor: |
| 207 | if (fUplo == clblasLower) |
| 208 | { |
| 209 | return *CMLPacked(row, column); |
| 210 | } |
| 211 | else { |
| 212 | return *CMUPacked(row, column); |
| 213 | } |
| 214 | break; |
| 215 | default: return FNAN<T>(); |
| 216 | } |
| 217 | } |
| 218 | |
| 219 | |
| 220 | template <typename T> |