memcpy.hpp

Go to the documentation of this file.

// Copyright (c) 2022--present, The Simons Foundation
// This file is part of TRIQS/nda and is licensed under the Apache License, Version 2.0.
// SPDX-License-Identifier: Apache-2.0
// See LICENSE in the root of this distribution for details.

 
#pragma once
 
#include "./address_space.hpp"
#include "../device.hpp"
#include "../macros.hpp"
 
#include <cstring>
 
namespace nda::mem {


  template <AddressSpace DestAdrSp, AddressSpace SrcAdrSp>
  void memcpy(void *dest, void const *src, size_t count) {
    check_adr_sp_valid<DestAdrSp, SrcAdrSp>();
    static_assert(nda::have_device == nda::have_cuda, "Adjust function for new device types");
 
    if constexpr (DestAdrSp == Host && SrcAdrSp == Host) {
      std::memcpy(dest, src, count);
    } else {
      device_error_check(cudaMemcpy(dest, src, count, cudaMemcpyDefault), "cudaMemcpy");
    }
  }

  template <AddressSpace DestAdrSp, AddressSpace SrcAdrSp>
  void memcpy2D(void *dest, size_t dpitch, const void *src, size_t spitch, size_t width, size_t height) {
    EXPECTS(width <= dpitch && width <= spitch);
    check_adr_sp_valid<DestAdrSp, SrcAdrSp>();
    static_assert(nda::have_device == nda::have_cuda, "Adjust function for new device types");
 
    if constexpr (DestAdrSp == Host && SrcAdrSp == Host) {
      auto *desti = static_cast<unsigned char *>(dest);
      auto *srci  = static_cast<const unsigned char *>(src);
      for (size_t i = 0; i < height; ++i, desti += dpitch, srci += spitch) std::memcpy(desti, srci, width);
    } else if (nda::have_device) {
      device_error_check(cudaMemcpy2D(dest, dpitch, src, spitch, width, height, cudaMemcpyDefault), "cudaMemcpy2D");
    }
  }

 
} // namespace nda::mem