wehub-resource-sync
92 lines
3.4 KiB
C++
92 lines
3.4 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 George A. Shulinok <sgazeos@gmail.com>, created on 4/18/2019.
|
|
//
|
|
|
|
#include <system/op_boilerplate.h>
|
|
#if NOT_EXCLUDED(OP_barnes_symmetrized)
|
|
|
|
#include <ops/declarable/CustomOperations.h>
|
|
#include <ops/declarable/helpers/BarnesHutTsne.h>
|
|
|
|
namespace sd {
|
|
namespace ops {
|
|
NDArray* rowCountsPtr = nullptr;
|
|
|
|
CUSTOM_OP_IMPL(barnes_symmetrized, 3, 3, false, 0, -1) {
|
|
auto rowP = INPUT_VARIABLE(0);
|
|
auto colP = INPUT_VARIABLE(1);
|
|
auto valP = INPUT_VARIABLE(2);
|
|
auto N = rowP->lengthOf() - 1;
|
|
auto outputRows = OUTPUT_VARIABLE(0);
|
|
auto outputCols = OUTPUT_VARIABLE(1);
|
|
auto outputVals = OUTPUT_VARIABLE(2);
|
|
|
|
if (block.getIArguments()->size() > 0) N = INT_ARG(0);
|
|
|
|
if (rowCountsPtr) {
|
|
helpers::barnes_symmetrize(rowP, colP, valP, N, outputRows, outputCols, outputVals, rowCountsPtr);
|
|
delete rowCountsPtr;
|
|
return Status::OK;
|
|
}
|
|
return Logger::logKernelFailureMsg("barnes_symmetrized: Cannot loop due wrong input data.");
|
|
}
|
|
|
|
DECLARE_TYPES(barnes_symmetrized) {
|
|
getOpDescriptor()
|
|
->setAllowedInputTypes(0, {INT32})
|
|
->setAllowedInputTypes(1, {INT32})
|
|
->setAllowedInputTypes(2, {ALL_INTS, ALL_FLOATS})
|
|
->setAllowedOutputTypes(1, {INT32})
|
|
->setAllowedOutputTypes(1, {INT32})
|
|
->setAllowedOutputTypes(2, {ALL_INTS, ALL_FLOATS})
|
|
->setSameMode(false);
|
|
}
|
|
|
|
DECLARE_SHAPE_FN(barnes_symmetrized) {
|
|
auto valPShapeInfo = inputShape->at(2);
|
|
LongType* outShapeInfo;
|
|
auto rowP = INPUT_VARIABLE(0);
|
|
auto colP = INPUT_VARIABLE(1);
|
|
auto N = rowP->lengthOf() - 1;
|
|
if (block.getIArguments()->size() > 0) N = INT_ARG(0);
|
|
auto dataType = rowP->dataType(); // ArrayOptions::dataType(inputShape->at(0));
|
|
std::vector<sd::LongType> shape = {N};
|
|
NDArray* rowCounts = NDArrayFactory::create_<int>('c',shape, block.launchContext());
|
|
LongType len = helpers::barnes_row_count(rowP, colP, N, *rowCounts);
|
|
rowCounts->syncToHost();
|
|
if (len <= 0) {
|
|
// CRITICAL: Clean up allocated array before throwing exception to prevent memory leak
|
|
delete rowCounts;
|
|
THROW_EXCEPTION("barnes_symmetrized: Cannot allocate shape due non-positive len.");
|
|
}
|
|
rowCountsPtr = rowCounts;
|
|
outShapeInfo =
|
|
ShapeBuilders::createShapeInfo(ArrayOptions::dataType(valPShapeInfo), 'c', {1, len}, block.getWorkspace());
|
|
auto outColsShapeInfo = ShapeBuilders::createShapeInfo(dataType, 'c', {1, len}, block.getWorkspace());
|
|
auto outRowsShapeInfo = ShapeBuilders::createShapeInfo(dataType, 'c', {1, N + 1}, block.getWorkspace());
|
|
return SHAPELIST(CONSTANT(outRowsShapeInfo), CONSTANT(outColsShapeInfo), CONSTANT(outShapeInfo));
|
|
}
|
|
|
|
} // namespace ops
|
|
} // namespace sd
|
|
|
|
#endif
|