wehub-resource-sync
113 lines
3.9 KiB
C++
113 lines
3.9 KiB
C++
/* ******************************************************************************
|
|
*
|
|
*
|
|
* This program and the accompanying materials are made available under the
|
|
* terms of the Apache License, Version 2.0 which is available at
|
|
* https://www.apache.org/licenses/LICENSE-2.0.
|
|
*
|
|
* See the NOTICE file distributed with this work for additional
|
|
* information regarding copyright ownership.
|
|
* Unless required by applicable law or agreed to in writing, software
|
|
* distributed under the License is distributed on an "AS IS" BASIS, WITHOUT
|
|
* WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the
|
|
* License for the specific language governing permissions and limitations
|
|
* under the License.
|
|
*
|
|
* SPDX-License-Identifier: Apache-2.0
|
|
******************************************************************************/
|
|
|
|
//
|
|
// @author sgazeos@gmail.com
|
|
//
|
|
#include <ops/declarable/helpers/roll.h>
|
|
#if NOT_EXCLUDED(OP_roll)
|
|
namespace sd {
|
|
namespace ops {
|
|
namespace helpers {
|
|
|
|
template <typename T>
|
|
static void rollFunctorLinear_(NDArray* input, NDArray* output, int shift, bool inplace) {
|
|
auto source = input;
|
|
if (!inplace) output->assign(input);
|
|
|
|
int fullLen = source->lengthOf();
|
|
int actualShift = shift; // % fullLen; // shift already non-negative then
|
|
if (actualShift < 0) {
|
|
actualShift -= fullLen * (actualShift / fullLen - 1);
|
|
} else
|
|
actualShift %= fullLen;
|
|
|
|
if (actualShift) {
|
|
int shiftCount = fullLen / actualShift - 1;
|
|
int remainShift = fullLen % actualShift;
|
|
|
|
// stage 1) swap last actualShift elements with first ones.
|
|
for (int e = 0; e < actualShift; ++e) {
|
|
int sourceIndex = fullLen - actualShift + e;
|
|
|
|
auto _e0 = output->e<T>(e);
|
|
auto _e1 = output->e<T>(sourceIndex);
|
|
|
|
output->p<T>(e, _e1);
|
|
output->p<T>(sourceIndex, _e0);
|
|
}
|
|
|
|
// stage 2) swap swapped actualShift elements with rest remainShiftCount times.
|
|
for (int count = 1; count < shiftCount; ++count) {
|
|
for (int e = 0; e < actualShift; ++e) {
|
|
int destinationIndex = fullLen - (count + 1) * actualShift + e;
|
|
int sourceIndex = fullLen - count * actualShift + e;
|
|
|
|
auto _e0 = output->e<T>(destinationIndex);
|
|
auto _e1 = output->e<T>(sourceIndex);
|
|
|
|
output->p<T>(destinationIndex, _e1);
|
|
output->p<T>(sourceIndex, _e0);
|
|
}
|
|
}
|
|
|
|
// stage 3) swap remainder of items.
|
|
if (remainShift && shiftCount)
|
|
for (int i = actualShift; i < 2 * actualShift; ++i) {
|
|
auto _e0 = output->e<T>(i);
|
|
auto _e1 = output->e<T>(i + remainShift);
|
|
output->p<T>(i, _e1);
|
|
output->p<T>(i + remainShift, _e0);
|
|
}
|
|
}
|
|
}
|
|
|
|
void rollFunctorFull(sd::LaunchContext* context, NDArray* input, NDArray* output, const std::vector<LongType>& shifts,
|
|
const std::vector<LongType>& axes, bool inplace) {
|
|
if (!inplace) output->assign(input);
|
|
|
|
auto source = output; // input;
|
|
for (size_t i = 0; i < axes.size(); i++) {
|
|
int axe = axes[i];
|
|
ResultSet listOfTensors = source->allTensorsAlongDimension({axe});
|
|
ResultSet listOfOutTensors = output->allTensorsAlongDimension({axe});
|
|
int fullLen = listOfTensors.size();
|
|
sd_debug("Roll: fullLen at last dimension is %d\n", fullLen);
|
|
int theShift = shifts[i];
|
|
if (theShift > 0) {
|
|
theShift %= fullLen;
|
|
} else {
|
|
theShift -= fullLen * (theShift / fullLen - 1);
|
|
}
|
|
for (int k = 0; k < fullLen; k++) {
|
|
rollFunctorLinear(context, listOfTensors.at(k), listOfOutTensors.at(k), theShift, true);
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
void rollFunctorLinear(sd::LaunchContext* context, NDArray* input, NDArray* output, int shift, bool inplace) {
|
|
BUILD_SINGLE_SELECTOR(input->dataType(), rollFunctorLinear_, (input, output, shift, inplace), SD_COMMON_TYPES);
|
|
}
|
|
|
|
BUILD_SINGLE_TEMPLATE( void rollFunctorLinear_, (NDArray * input, NDArray* output, int shift, bool inplace),
|
|
SD_COMMON_TYPES);
|
|
} // namespace helpers
|
|
} // namespace ops
|
|
} // namespace sd
|
|
#endif |