* TODO: the intent is for DeviceBuffer to become a class.
*
* \author Artem Zhmurov <zhmurov@gmail.com>
+ * \author Erik Lindahl <erik.lindahl@gmail.com>
+ * \author Andrey Alekseenko <al42and@gmail.com>
*
* \inlibraryapi
*/
#include "gromacs/gpu_utils/device_context.h"
#include "gromacs/gpu_utils/device_stream.h"
#include "gromacs/gpu_utils/devicebuffer_datatype.h"
+#include "gromacs/gpu_utils/gmxsycl.h"
#include "gromacs/gpu_utils/gpu_utils.h" //only for GpuApiCallBehavior
-#include "gromacs/gpu_utils/gputraits_ocl.h"
+#include "gromacs/gpu_utils/gputraits_sycl.h"
#include "gromacs/utility/gmxassert.h"
#include "gromacs/utility/stringutil.h"
+#ifndef DOXYGEN
+template<typename T>
+class DeviceBuffer<T>::ClSyclBufferWrapper : public cl::sycl::buffer<T, 1>
+{
+ using cl::sycl::buffer<T, 1>::buffer; // Get all the constructors
+};
+
+template<typename T>
+using ClSyclBufferWrapper = typename DeviceBuffer<T>::ClSyclBufferWrapper;
+
+//! Constructor.
+template<typename T>
+DeviceBuffer<T>::DeviceBuffer() : buffer_(nullptr)
+{
+}
+
+//! Destructor.
+template<typename T>
+DeviceBuffer<T>::~DeviceBuffer() = default;
+
+//! Copy constructor (references the same underlying SYCL buffer).
+template<typename T>
+DeviceBuffer<T>::DeviceBuffer(DeviceBuffer<T> const& src) :
+ buffer_(new ClSyclBufferWrapper(*src.buffer_))
+{
+}
+
+//! Move constructor.
+template<typename T>
+DeviceBuffer<T>::DeviceBuffer(DeviceBuffer<T>&& src) noexcept = default;
+
+//! Copy assignment (references the same underlying SYCL buffer).
+template<typename T>
+DeviceBuffer<T>& DeviceBuffer<T>::operator=(DeviceBuffer<T> const& src)
+{
+ buffer_.reset(new ClSyclBufferWrapper(*src.buffer_));
+ return *this;
+}
+
+//! Move assignment.
+template<typename T>
+DeviceBuffer<T>& DeviceBuffer<T>::operator=(DeviceBuffer<T>&& src) noexcept = default;
+
+#endif // #ifndef DOXYGEN
+
/*! \libinternal \brief
* Allocates a device-side buffer.
* It is currently a caller's responsibility to call it only on not-yet allocated buffers.
*
* \tparam ValueType Raw value type of the \p buffer.
* \param[in,out] buffer Pointer to the device-side buffer.
- * \param[in] numValues Number of values to accomodate.
+ * \param[in] numValues Number of values to accommodate.
* \param[in] deviceContext The buffer's device context-to-be.
*/
template<typename ValueType>
-void allocateDeviceBuffer(DeviceBuffer<ValueType>* gmx_unused buffer,
- size_t gmx_unused numValues,
- const DeviceContext& gmx_unused deviceContext)
+void allocateDeviceBuffer(DeviceBuffer<ValueType>* buffer, size_t numValues, const DeviceContext& deviceContext)
{
- // SYCL-TODO
+ /* SYCL does not require binding buffer to a specific context or device. The ::context_bound
+ * property only enforces the use of only given context, and possibly offers some optimizations */
+ const cl::sycl::property_list bufferProperties{ cl::sycl::property::buffer::context_bound(
+ deviceContext.context()) };
+ buffer->buffer_.reset(
+ new ClSyclBufferWrapper<ValueType>(cl::sycl::range<1>(numValues), bufferProperties));
}
/*! \brief
*
* \param[in] buffer Pointer to the buffer to free.
*/
-template<typename DeviceBuffer>
-void freeDeviceBuffer(DeviceBuffer* gmx_unused buffer)
+template<typename ValueType>
+void freeDeviceBuffer(DeviceBuffer<ValueType>* buffer)
{
- // SYCL-TODO
+ buffer->buffer_.reset(nullptr);
}
/*! \brief
* Performs the host-to-device data copy, synchronous or asynchronously on request.
*
- * Note that synchronous copy will not synchronize the stream in case of zero \p numValues
- * because of the early return.
+ * Unlike in CUDA and OpenCL, synchronous call does not guarantee that all previously
+ * submitted operations are complete, only the ones that are required for \p buffer consistency.
*
* \tparam ValueType Raw value type of the \p buffer.
- * \param[in,out] buffer Pointer to the device-side buffer
- * \param[in] hostBuffer Pointer to the raw host-side memory, also typed \p ValueType
+ * \param[in,out] buffer Pointer to the device-side buffer.
+ * \param[in] hostBuffer Pointer to the raw host-side memory, also typed \p ValueType.
* \param[in] startingOffset Offset (in values) at the device-side buffer to copy into.
* \param[in] numValues Number of values to copy.
* \param[in] deviceStream GPU stream to perform asynchronous copy in.
* \param[in] transferKind Copy type: synchronous or asynchronous.
* \param[out] timingEvent A pointer to the H2D copy timing event to be filled in.
- * If the pointer is not null, the event can further be used
- * to queue a wait for this operation or to query profiling information.
+ * Ignored in SYCL.
*/
template<typename ValueType>
-void copyToDeviceBuffer(DeviceBuffer<ValueType>* gmx_unused buffer,
- const ValueType* gmx_unused hostBuffer,
- size_t gmx_unused startingOffset,
- size_t gmx_unused numValues,
- const DeviceStream& gmx_unused deviceStream,
- GpuApiCallBehavior gmx_unused transferKind,
+void copyToDeviceBuffer(DeviceBuffer<ValueType>* buffer,
+ const ValueType* hostBuffer,
+ size_t startingOffset,
+ size_t numValues,
+ const DeviceStream& deviceStream,
+ GpuApiCallBehavior transferKind,
CommandEvent* gmx_unused timingEvent)
{
- // SYCL-TODO
+ if (numValues == 0)
+ {
+ return; // such calls are actually made with empty domains
+ }
+ GMX_ASSERT(buffer, "needs a buffer pointer");
+ GMX_ASSERT(buffer->buffer_, "needs an initialized buffer pointer");
+ GMX_ASSERT(hostBuffer, "needs a host buffer pointer");
+
+ cl::sycl::buffer<ValueType>& syclBuffer = *buffer->buffer_;
+
+ cl::sycl::event ev = deviceStream.stream().submit([&](cl::sycl::handler& cgh) {
+ /* Here and elsewhere in this file, accessor constructor is user instead of a more common
+ * buffer::get_access, since the compiler (icpx 2021.1-beta09) occasionally gets confused
+ * by all the overloads */
+ auto d_bufferAccessor = cl::sycl::accessor<ValueType, 1, cl::sycl::access::mode::discard_write>{
+ syclBuffer, cgh, cl::sycl::range(numValues), cl::sycl::id(startingOffset)
+ };
+ cgh.copy(hostBuffer, d_bufferAccessor);
+ });
+ if (transferKind == GpuApiCallBehavior::Sync)
+ {
+ ev.wait_and_throw();
+ }
}
/*! \brief
* Performs the device-to-host data copy, synchronous or asynchronously on request.
*
- * Note that synchronous copy will not synchronize the stream in case of zero \p numValues
- * because of the early return.
+ * Unlike in CUDA and OpenCL, synchronous call does not guarantee that all previously
+ * submitted operations are complete, only the ones that are required for \p buffer consistency.
*
* \tparam ValueType Raw value type of the \p buffer.
* \param[in,out] hostBuffer Pointer to the raw host-side memory, also typed \p ValueType
- * \param[in] buffer Pointer to the device-side buffer
+ * \param[in] buffer Pointer to the device-side buffer.
* \param[in] startingOffset Offset (in values) at the device-side buffer to copy from.
* \param[in] numValues Number of values to copy.
* \param[in] deviceStream GPU stream to perform asynchronous copy in.
* \param[in] transferKind Copy type: synchronous or asynchronous.
* \param[out] timingEvent A pointer to the H2D copy timing event to be filled in.
- * If the pointer is not null, the event can further be used
- * to queue a wait for this operation or to query profiling information.
+ * Ignored in SYCL.
*/
template<typename ValueType>
-void copyFromDeviceBuffer(ValueType* gmx_unused hostBuffer,
- DeviceBuffer<ValueType>* gmx_unused buffer,
- size_t gmx_unused startingOffset,
- size_t gmx_unused numValues,
- const DeviceStream& gmx_unused deviceStream,
- GpuApiCallBehavior gmx_unused transferKind,
+void copyFromDeviceBuffer(ValueType* hostBuffer,
+ DeviceBuffer<ValueType>* buffer,
+ size_t startingOffset,
+ size_t numValues,
+ const DeviceStream& deviceStream,
+ GpuApiCallBehavior transferKind,
CommandEvent* gmx_unused timingEvent)
{
- // SYCL-TODO
+ if (numValues == 0)
+ {
+ return; // such calls are actually made with empty domains
+ }
+ GMX_ASSERT(buffer, "needs a buffer pointer");
+ GMX_ASSERT(hostBuffer, "needs a host buffer pointer");
+
+ cl::sycl::buffer<ValueType>& syclBuffer = *buffer->buffer_;
+
+ cl::sycl::event ev = deviceStream.stream().submit([&](cl::sycl::handler& cgh) {
+ const auto d_bufferAccessor = cl::sycl::accessor<ValueType, 1, cl::sycl::access::mode::read>{
+ syclBuffer, cgh, cl::sycl::range(numValues), cl::sycl::id(startingOffset)
+ };
+ cgh.copy(d_bufferAccessor, hostBuffer);
+ });
+ if (transferKind == GpuApiCallBehavior::Sync)
+ {
+ ev.wait_and_throw();
+ }
}
/*! \brief
* Clears the device buffer asynchronously.
*
* \tparam ValueType Raw value type of the \p buffer.
- * \param[in,out] buffer Pointer to the device-side buffer
+ * \param[in,out] buffer Pointer to the device-side buffer.
* \param[in] startingOffset Offset (in values) at the device-side buffer to start clearing at.
* \param[in] numValues Number of values to clear.
* \param[in] deviceStream GPU stream.
*/
template<typename ValueType>
-void clearDeviceBufferAsync(DeviceBuffer<ValueType>* gmx_unused buffer,
- size_t gmx_unused startingOffset,
- size_t gmx_unused numValues,
- const DeviceStream& gmx_unused deviceStream)
+void clearDeviceBufferAsync(DeviceBuffer<ValueType>* buffer,
+ size_t startingOffset,
+ size_t numValues,
+ const DeviceStream& deviceStream)
{
- // SYCL-TODO
+ if (numValues == 0)
+ {
+ return;
+ }
+ GMX_ASSERT(buffer, "needs a buffer pointer");
+
+ const ValueType pattern{};
+ cl::sycl::buffer<ValueType>& syclBuffer = *(buffer->buffer_);
+
+ cl::sycl::event ev = deviceStream.stream().submit([&](cl::sycl::handler& cgh) {
+ auto d_bufferAccessor = cl::sycl::accessor<ValueType, 1, cl::sycl::access::mode::discard_write>{
+ syclBuffer, cgh, cl::sycl::range(numValues), cl::sycl::id(startingOffset)
+ };
+ cgh.fill(d_bufferAccessor, pattern);
+ });
}
/*! \brief Check the validity of the device buffer.
*
- * Checks if the buffer is not nullptr and if its allocation is big enough.
+ * Checks if the buffer is valid and if its allocation is big enough.
*
* \param[in] buffer Device buffer to be checked.
* \param[in] requiredSize Number of elements that the buffer will have to accommodate.
*
- * \returns Whether the device buffer can be set.
+ * \returns Whether the device buffer exists and has enough capacity.
*/
template<typename T>
-static bool gmx_unused checkDeviceBuffer(DeviceBuffer<T> gmx_unused buffer, int gmx_unused requiredSize)
+static gmx_unused bool checkDeviceBuffer(DeviceBuffer<T> buffer, int requiredSize)
{
- // SYCL-TODO
+ return buffer.buffer_ && (static_cast<int>(buffer.buffer_->get_count()) >= requiredSize);
}
-//! Device texture wrapper.
-using DeviceTexture = void*;
-
/*! \brief Create a texture object for an array of type ValueType.
*
* Creates the device buffer and copies read-only data for an array of type ValueType.
- *
- * \todo Decide if using image2d is most efficient.
+ * Like OpenCL, does not really do anything with textures, simply creates a buffer
+ * and initializes it.
*
* \tparam ValueType Raw data type.
*
* \param[out] deviceBuffer Device buffer to store data in.
- * \param[out] deviceTexture New texture object
* \param[in] hostBuffer Host buffer to get date from.
* \param[in] numValues Number of elements in the buffer.
* \param[in] deviceContext GPU device context.
*/
template<typename ValueType>
-void initParamLookupTable(DeviceBuffer<ValueType>* gmx_unused deviceBuffer,
- DeviceTexture* gmx_unused deviceTexture,
- const ValueType* gmx_unused hostBuffer,
- int gmx_unused numValues,
- const DeviceContext& gmx_unused deviceContext)
+void initParamLookupTable(DeviceBuffer<ValueType>* deviceBuffer,
+ DeviceTexture* /* deviceTexture */,
+ const ValueType* hostBuffer,
+ int numValues,
+ const DeviceContext& deviceContext)
{
- // SYCL-TODO
+ GMX_ASSERT(hostBuffer, "Host buffer should be specified.");
+ GMX_ASSERT(deviceBuffer, "Device buffer should be specified.");
+
+ /* Constructing buffer with cl::sycl::buffer(T* data, size_t size) will take ownership
+ * of this memory region making it unusable, which might lead to side-effects.
+ * On the other hand, cl::sycl::buffer(InputIterator<T> begin, InputIterator<T> end) will
+ * initialize the buffer without affecting ownership of the memory, although
+ * it will consume extra memory on host. */
+ const cl::sycl::property_list bufferProperties{ cl::sycl::property::buffer::context_bound(
+ deviceContext.context()) };
+ deviceBuffer->buffer_.reset(new ClSyclBufferWrapper<ValueType>(
+ hostBuffer, hostBuffer + numValues, bufferProperties));
}
/*! \brief Release the OpenCL device buffer.
* \tparam ValueType Raw data type.
*
* \param[in,out] deviceBuffer Device buffer to store data in.
- * \param[in] deviceTexture Reference to texture object
*/
template<typename ValueType>
-void destroyParamLookupTable(DeviceBuffer<ValueType>* gmx_unused deviceBuffer,
- DeviceTexture& gmx_unused deviceTexture)
+void destroyParamLookupTable(DeviceBuffer<ValueType>* deviceBuffer, DeviceTexture& /* deviceTexture */)
{
- // SYCL-TODO
+ deviceBuffer->buffer_.reset(nullptr);
}
#endif // GMX_GPU_UTILS_DEVICEBUFFER_SYCL_H