deeplearning4j/deeplearning4j

/* ******************************************************************************
 *
 *
 * 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 Yurii Shyrma (iuriish@yahoo.com), created on 05.02.2018
//

#include <system/op_boilerplate.h>
#if NOT_EXCLUDED(OP_meshgrid)

#include <ops/declarable/CustomOperations.h>
#include <ops/declarable/helpers/meshgrid.h>

#include <numeric>

namespace sd {
namespace ops {

CUSTOM_OP_IMPL(meshgrid, -1, -1, false, 0, 0) {
  int rank = block.width();

  if (rank == 1) {
    OUTPUT_VARIABLE(0)->assign(INPUT_VARIABLE(0));
    return sd::Status::OK;
  }

  bool swapFirst2Dims = block.getIArguments()->size() > 0 ? (bool)INT_ARG(0) : true;

  std::vector<NDArray*> inArrs(rank);
  std::vector<NDArray*> outArrs(rank);

  for (int i = 0; i < rank; ++i) {
    inArrs[i] = INPUT_VARIABLE(i);
    outArrs[i] = OUTPUT_VARIABLE(i);
  }

  helpers::meshgrid(block.launchContext(), inArrs, outArrs, swapFirst2Dims);

  return sd::Status::OK;
}

DECLARE_TYPES(meshgrid) {
  getOpDescriptor()->setAllowedInputTypes(DataType::ANY)->setAllowedOutputTypes(DataType::INHERIT)->setSameMode(true);
}

DECLARE_SHAPE_FN(meshgrid) {
  bool swapFirst2Dims = block.getIArguments()->size() > 0 ? (bool)INT_ARG(0) : true;

  int rank = block.width();
  sd::LongType* outShapeInfo = nullptr;
  ALLOCATE(outShapeInfo, block.getWorkspace(), shape::shapeInfoLength(rank), sd::LongType);
  outShapeInfo[0] = rank;
  for (int i = 1; i <= rank; ++i) outShapeInfo[i] = (sd::LongType)shape::length(inputShape->at(i - 1));

  if (swapFirst2Dims && rank > 1) math::sd_swap<sd::LongType>(outShapeInfo[1], outShapeInfo[2]);

  auto in = inputShape->at(0);
  ShapeUtils::updateStridesAndType(outShapeInfo, in, shape::order(in));

  auto shapes = SHAPELIST();
  auto resultShape = CONSTANT(outShapeInfo);
  shapes->push_back(resultShape);

  for (int i = 2; i <= rank; ++i) {
    shapes->push_back(resultShape);
  }

  return shapes;
}

}  // namespace ops
}  // namespace sd

#endif