IndexReduction.cpp

Go to the documentation of this file.

// ----------------------------------------------------------------------------
// -                        CloudViewer: www.cloudViewer.org                  -
// ----------------------------------------------------------------------------
// Copyright (c) 2018-2024 www.cloudViewer.org
// SPDX-License-Identifier: MIT
// ----------------------------------------------------------------------------
 
#include "cloudViewer/core/kernel/IndexReduction.h"
 
#include <Logging.h>
 
namespace cloudViewer {
namespace core {
namespace kernel {
 
void IndexAdd_(int64_t dim,
               const Tensor& index,
               const Tensor& src,
               Tensor& dst) {
    // Permute the reduction dimension to the first.
    SizeVector permute = {};
    for (int64_t d = 0; d <= dim; ++d) {
        if (d == 0) {
            permute.push_back(dim);
        } else {
            permute.push_back(d - 1);
        }
    }
    for (int64_t d = dim + 1; d < src.NumDims(); ++d) {
        permute.push_back(d);
    }
 
    auto src_permute = src.Permute(permute);
    auto dst_permute = dst.Permute(permute);
 
    if (dst.IsCPU()) {
        IndexAddCPU_(dim, index, src_permute, dst_permute);
    } else if (dst.IsSYCL()) {
#ifdef BUILD_SYCL_MODULE
        IndexAddSYCL_(dim, index, src_permute, dst_permute);
#endif
    } else if (dst.IsCUDA()) {
#ifdef BUILD_CUDA_MODULE
        IndexAddCUDA_(dim, index, src_permute, dst_permute);
#endif
    } else {
        utility::LogError("IndexAdd_: Unimplemented device");
    }
}
 
}  // namespace kernel
}  // namespace core
}  // namespace cloudViewer