Alexandria  2.25.0
SDC-CH common library for the Euclid project
Operations.cpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2012-2021 Euclid Science Ground Segment
3  *
4  * This library is free software; you can redistribute it and/or modify it under
5  * the terms of the GNU Lesser General Public License as published by the Free
6  * Software Foundation; either version 3.0 of the License, or (at your option)
7  * any later version.
8  *
9  * This library is distributed in the hope that it will be useful, but WITHOUT
10  * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11  * FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
12  * details.
13  *
14  * You should have received a copy of the GNU Lesser General Public License
15  * along with this library; if not, write to the Free Software Foundation, Inc.,
16  * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
17  */
18 
19 #include "NdArray/Operations.h"
20 
21 namespace Euclid {
22 namespace NdArray {
23 
25  std::vector<std::size_t> coords(shape.size());
26  std::size_t i = coords.size();
27  do {
28  --i;
29  std::size_t stride = shape[i];
30  coords[i] = index % stride;
31  index /= stride;
32  } while (i > 0);
33  if (index > 0) {
34  throw std::out_of_range("Index out of bounds");
35  }
36  return coords;
37 }
38 
39 } // namespace NdArray
40 } // namespace Euclid
std::vector< std::size_t > unravel_index(std::size_t index, const std::vector< std::size_t > &shape)
Definition: Operations.cpp:24
T size(T... args)