malloc.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 <cstdlib>
 
namespace nda::mem {


  template <AddressSpace AdrSp>
  void *malloc(size_t size) {
    check_adr_sp_valid<AdrSp>();
    static_assert(nda::have_device == nda::have_cuda, "Adjust function for new device types");
 
    void *ptr = nullptr;
    if constexpr (AdrSp == Host) {
      ptr = std::malloc(size); // NOLINT (we want to return a void*)
    } else if constexpr (AdrSp == Device) { // NOLINT (branch is not repeated)
      device_error_check(cudaMalloc((void **)&ptr, size), "cudaMalloc");
    } else {
      device_error_check(cudaMallocManaged((void **)&ptr, size), "cudaMallocManaged");
    }
    return ptr;
  }

  template <AddressSpace AdrSp>
  void free(void *p) {
    check_adr_sp_valid<AdrSp>();
    static_assert(nda::have_device == nda::have_cuda, "Adjust function for new device types");
 
    if constexpr (AdrSp == Host) {
      std::free(p); // NOLINT (we want to call free with a void*)
    } else {
      device_error_check(cudaFree(p), "cudaFree");
    }
  }

 
} // namespace nda::mem