xref: /petsc/src/sys/objects/device/interface/mark_dcontext.cxx (revision 271378f25b3157c7b80658e8fc508cbdf0c8645c)
10e6b6b59SJacob Faibussowitsch #include "petscdevice_interface_internal.hpp" /*I <petscdevice.h> I*/
20e6b6b59SJacob Faibussowitsch 
30e6b6b59SJacob Faibussowitsch #include <petsc/private/cpp/object_pool.hpp>
40e6b6b59SJacob Faibussowitsch #include <petsc/private/cpp/utility.hpp>
50e6b6b59SJacob Faibussowitsch 
60e6b6b59SJacob Faibussowitsch #include <unordered_map>
70e6b6b59SJacob Faibussowitsch #include <algorithm> // std::remove_if(), std::find_if()
80e6b6b59SJacob Faibussowitsch #include <vector>
90e6b6b59SJacob Faibussowitsch #include <string>
100e6b6b59SJacob Faibussowitsch #include <sstream> // std::ostringstream
110e6b6b59SJacob Faibussowitsch 
12a966a306SPierre Jolivet #if defined(__clang__)
13a966a306SPierre Jolivet   #pragma clang diagnostic push
14a966a306SPierre Jolivet   #pragma clang diagnostic ignored "-Wgnu-zero-variadic-macro-arguments"
15a966a306SPierre Jolivet #endif
16a966a306SPierre Jolivet 
170e6b6b59SJacob Faibussowitsch // ==========================================================================================
180e6b6b59SJacob Faibussowitsch // PetscEvent
190e6b6b59SJacob Faibussowitsch // ==========================================================================================
200e6b6b59SJacob Faibussowitsch 
210e6b6b59SJacob Faibussowitsch struct PetscEventAllocator : public Petsc::AllocatorBase<PetscEvent> {
22d71ae5a4SJacob Faibussowitsch   PETSC_NODISCARD static PetscErrorCode create(PetscEvent *event) noexcept
23d71ae5a4SJacob Faibussowitsch   {
240e6b6b59SJacob Faibussowitsch     PetscFunctionBegin;
250e6b6b59SJacob Faibussowitsch     PetscCall(PetscNew(event));
260e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
270e6b6b59SJacob Faibussowitsch   }
280e6b6b59SJacob Faibussowitsch 
29d71ae5a4SJacob Faibussowitsch   PETSC_NODISCARD static PetscErrorCode destroy(PetscEvent event) noexcept
30d71ae5a4SJacob Faibussowitsch   {
310e6b6b59SJacob Faibussowitsch     PetscFunctionBegin;
320e6b6b59SJacob Faibussowitsch     PetscCall(reset(event));
330e6b6b59SJacob Faibussowitsch     PetscCall(PetscFree(event));
340e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
350e6b6b59SJacob Faibussowitsch   }
360e6b6b59SJacob Faibussowitsch 
37*271378f2SJacob Faibussowitsch   PETSC_NODISCARD static PetscErrorCode reset(PetscEvent event) noexcept
38d71ae5a4SJacob Faibussowitsch   {
390e6b6b59SJacob Faibussowitsch     PetscFunctionBegin;
400e6b6b59SJacob Faibussowitsch     if (auto &destroy = event->destroy) {
410e6b6b59SJacob Faibussowitsch       PetscCall((*destroy)(event));
420e6b6b59SJacob Faibussowitsch       destroy = nullptr;
430e6b6b59SJacob Faibussowitsch     }
44*271378f2SJacob Faibussowitsch     PetscAssert(!event->data, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Event failed to destroy its data member: %p", event->data);
450e6b6b59SJacob Faibussowitsch     event->dctx_id    = 0;
460e6b6b59SJacob Faibussowitsch     event->dctx_state = 0;
470e6b6b59SJacob Faibussowitsch     event->dtype      = PETSC_DEVICE_DEFAULT();
480e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
490e6b6b59SJacob Faibussowitsch   }
500e6b6b59SJacob Faibussowitsch };
510e6b6b59SJacob Faibussowitsch 
520e6b6b59SJacob Faibussowitsch static Petsc::ObjectPool<PetscEvent, PetscEventAllocator> event_pool;
530e6b6b59SJacob Faibussowitsch 
54d71ae5a4SJacob Faibussowitsch static PetscErrorCode PetscDeviceContextCreateEvent_Private(PetscDeviceContext dctx, PetscEvent *event)
55d71ae5a4SJacob Faibussowitsch {
560e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
570e6b6b59SJacob Faibussowitsch   PetscValidDeviceContext(dctx, 1);
580e6b6b59SJacob Faibussowitsch   PetscValidPointer(event, 2);
590e6b6b59SJacob Faibussowitsch   PetscCall(event_pool.allocate(event));
600e6b6b59SJacob Faibussowitsch   PetscCall(PetscDeviceContextGetDeviceType(dctx, &(*event)->dtype));
610e6b6b59SJacob Faibussowitsch   PetscTryTypeMethod(dctx, createevent, *event);
620e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
630e6b6b59SJacob Faibussowitsch }
640e6b6b59SJacob Faibussowitsch 
65d71ae5a4SJacob Faibussowitsch static PetscErrorCode PetscEventDestroy_Private(PetscEvent *event)
66d71ae5a4SJacob Faibussowitsch {
670e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
680e6b6b59SJacob Faibussowitsch   PetscValidPointer(event, 1);
690e6b6b59SJacob Faibussowitsch   if (*event) PetscCall(event_pool.deallocate(Petsc::util::exchange(*event, nullptr)));
700e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
710e6b6b59SJacob Faibussowitsch }
720e6b6b59SJacob Faibussowitsch 
73d71ae5a4SJacob Faibussowitsch static PetscErrorCode PetscDeviceContextRecordEvent_Private(PetscDeviceContext dctx, PetscEvent event)
74d71ae5a4SJacob Faibussowitsch {
750e6b6b59SJacob Faibussowitsch   PetscObjectId    id;
760e6b6b59SJacob Faibussowitsch   PetscObjectState state;
770e6b6b59SJacob Faibussowitsch 
780e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
790e6b6b59SJacob Faibussowitsch   PetscValidDeviceContext(dctx, 1);
800e6b6b59SJacob Faibussowitsch   PetscValidPointer(event, 2);
810e6b6b59SJacob Faibussowitsch   id    = PetscObjectCast(dctx)->id;
820e6b6b59SJacob Faibussowitsch   state = PetscObjectCast(dctx)->state;
830e6b6b59SJacob Faibussowitsch   // technically state can never be less than event->dctx_state (only equal) but we include
840e6b6b59SJacob Faibussowitsch   // it in the check just in case
850e6b6b59SJacob Faibussowitsch   if ((id == event->dctx_id) && (state <= event->dctx_state)) PetscFunctionReturn(0);
860e6b6b59SJacob Faibussowitsch   if (dctx->ops->recordevent) {
870e6b6b59SJacob Faibussowitsch     // REVIEW ME:
880e6b6b59SJacob Faibussowitsch     // TODO maybe move this to impls, as they can determine whether they can interoperate with
890e6b6b59SJacob Faibussowitsch     // other device types more readily
900e6b6b59SJacob Faibussowitsch     if (PetscDefined(USE_DEBUG) && (event->dtype != PETSC_DEVICE_HOST)) {
910e6b6b59SJacob Faibussowitsch       PetscDeviceType dtype;
920e6b6b59SJacob Faibussowitsch 
930e6b6b59SJacob Faibussowitsch       PetscCall(PetscDeviceContextGetDeviceType(dctx, &dtype));
940e6b6b59SJacob Faibussowitsch       PetscCheck(event->dtype == dtype, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Event type %s does not match device context type %s", PetscDeviceTypes[event->dtype], PetscDeviceTypes[dtype]);
950e6b6b59SJacob Faibussowitsch     }
960e6b6b59SJacob Faibussowitsch     PetscUseTypeMethod(dctx, recordevent, event);
970e6b6b59SJacob Faibussowitsch   }
980e6b6b59SJacob Faibussowitsch   event->dctx_id    = id;
990e6b6b59SJacob Faibussowitsch   event->dctx_state = state;
1000e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
1010e6b6b59SJacob Faibussowitsch }
1020e6b6b59SJacob Faibussowitsch 
103d71ae5a4SJacob Faibussowitsch static PetscErrorCode PetscDeviceContextWaitForEvent_Private(PetscDeviceContext dctx, PetscEvent event)
104d71ae5a4SJacob Faibussowitsch {
1050e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
1060e6b6b59SJacob Faibussowitsch   PetscValidDeviceContext(dctx, 1);
1070e6b6b59SJacob Faibussowitsch   PetscValidPointer(event, 2);
1080e6b6b59SJacob Faibussowitsch   // empty data implies you cannot wait on this event
1090e6b6b59SJacob Faibussowitsch   if (!event->data) PetscFunctionReturn(0);
1100e6b6b59SJacob Faibussowitsch   if (PetscDefined(USE_DEBUG)) {
1110e6b6b59SJacob Faibussowitsch     const auto      etype = event->dtype;
1120e6b6b59SJacob Faibussowitsch     PetscDeviceType dtype;
1130e6b6b59SJacob Faibussowitsch 
1140e6b6b59SJacob Faibussowitsch     PetscCall(PetscDeviceContextGetDeviceType(dctx, &dtype));
1150e6b6b59SJacob Faibussowitsch     PetscCheck(etype == dtype, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Event type %s does not match device context type %s", PetscDeviceTypes[etype], PetscDeviceTypes[dtype]);
1160e6b6b59SJacob Faibussowitsch   }
1170e6b6b59SJacob Faibussowitsch   if (PetscObjectCast(dctx)->id == event->dctx_id) PetscFunctionReturn(0);
1180e6b6b59SJacob Faibussowitsch   PetscTryTypeMethod(dctx, waitforevent, event);
1190e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
1200e6b6b59SJacob Faibussowitsch }
1210e6b6b59SJacob Faibussowitsch 
1220e6b6b59SJacob Faibussowitsch // ==========================================================================================
1230e6b6b59SJacob Faibussowitsch // PetscStackFrame
1240e6b6b59SJacob Faibussowitsch //
1250e6b6b59SJacob Faibussowitsch // A helper class that (when debugging is enabled) contains the stack frame from which
1260e6b6b59SJacob Faibussowitsch // PetscDeviceContextMakrIntentFromID(). It is intended to be derived from, since this enables
1270e6b6b59SJacob Faibussowitsch // empty-base-class optimization to kick in when debugging is disabled.
1280e6b6b59SJacob Faibussowitsch // ==========================================================================================
1290e6b6b59SJacob Faibussowitsch 
1300e6b6b59SJacob Faibussowitsch template <bool use_debug>
1310e6b6b59SJacob Faibussowitsch struct PetscStackFrame;
1320e6b6b59SJacob Faibussowitsch 
1330e6b6b59SJacob Faibussowitsch template <>
1340e6b6b59SJacob Faibussowitsch struct PetscStackFrame</* use_debug = */ true> {
1350e6b6b59SJacob Faibussowitsch   std::string file{};
1360e6b6b59SJacob Faibussowitsch   std::string function{};
1370e6b6b59SJacob Faibussowitsch   int         line{};
1380e6b6b59SJacob Faibussowitsch 
1390e6b6b59SJacob Faibussowitsch   PetscStackFrame() = default;
1400e6b6b59SJacob Faibussowitsch 
1410e6b6b59SJacob Faibussowitsch   PetscStackFrame(const char *file_, const char *func_, int line_) noexcept : file(split_on_petsc_path_(file_)), function(func_), line(line_) { }
1420e6b6b59SJacob Faibussowitsch 
1430e6b6b59SJacob Faibussowitsch   bool operator==(const PetscStackFrame &other) const noexcept { return line == other.line && file == other.file && function == other.function; }
1440e6b6b59SJacob Faibussowitsch 
1457c441f3aSJacob Faibussowitsch   PETSC_NODISCARD std::string to_string() const noexcept
1467c441f3aSJacob Faibussowitsch   {
1477c441f3aSJacob Faibussowitsch     std::string ret;
1487c441f3aSJacob Faibussowitsch 
1497c441f3aSJacob Faibussowitsch     ret = '(' + function + "() at " + file + ':' + std::to_string(line) + ')';
1507c441f3aSJacob Faibussowitsch     return ret;
1517c441f3aSJacob Faibussowitsch   }
1527c441f3aSJacob Faibussowitsch 
1530e6b6b59SJacob Faibussowitsch private:
154d71ae5a4SJacob Faibussowitsch   static std::string split_on_petsc_path_(std::string &&in) noexcept
155d71ae5a4SJacob Faibussowitsch   {
1560e6b6b59SJacob Faibussowitsch     auto pos = in.find("petsc/src");
1570e6b6b59SJacob Faibussowitsch 
1580e6b6b59SJacob Faibussowitsch     if (pos == std::string::npos) pos = in.find("petsc/include");
1590e6b6b59SJacob Faibussowitsch     if (pos == std::string::npos) pos = 0;
1600e6b6b59SJacob Faibussowitsch     return in.substr(pos);
1610e6b6b59SJacob Faibussowitsch   }
1620e6b6b59SJacob Faibussowitsch 
163d71ae5a4SJacob Faibussowitsch   friend std::ostream &operator<<(std::ostream &os, const PetscStackFrame &frame)
164d71ae5a4SJacob Faibussowitsch   {
1657c441f3aSJacob Faibussowitsch     os << frame.to_string();
1660e6b6b59SJacob Faibussowitsch     return os;
1670e6b6b59SJacob Faibussowitsch   }
1687c441f3aSJacob Faibussowitsch 
1697c441f3aSJacob Faibussowitsch   friend void swap(PetscStackFrame &lhs, PetscStackFrame &rhs) noexcept
1707c441f3aSJacob Faibussowitsch   {
1717c441f3aSJacob Faibussowitsch     using std::swap;
1727c441f3aSJacob Faibussowitsch 
1737c441f3aSJacob Faibussowitsch     swap(lhs.file, rhs.file);
1747c441f3aSJacob Faibussowitsch     swap(lhs.function, rhs.function);
1757c441f3aSJacob Faibussowitsch     swap(lhs.line, rhs.line);
1767c441f3aSJacob Faibussowitsch   }
1770e6b6b59SJacob Faibussowitsch };
1780e6b6b59SJacob Faibussowitsch 
1790e6b6b59SJacob Faibussowitsch template <>
1800e6b6b59SJacob Faibussowitsch struct PetscStackFrame</* use_debug = */ false> {
1810e6b6b59SJacob Faibussowitsch   template <typename... T>
182d71ae5a4SJacob Faibussowitsch   constexpr PetscStackFrame(T &&...) noexcept
183d71ae5a4SJacob Faibussowitsch   {
184d71ae5a4SJacob Faibussowitsch   }
1850e6b6b59SJacob Faibussowitsch 
1860e6b6b59SJacob Faibussowitsch   constexpr bool operator==(const PetscStackFrame &) const noexcept { return true; }
1870e6b6b59SJacob Faibussowitsch 
1887c441f3aSJacob Faibussowitsch   PETSC_NODISCARD static std::string to_string() noexcept { return "(unknown)"; }
1897c441f3aSJacob Faibussowitsch 
190d71ae5a4SJacob Faibussowitsch   friend std::ostream &operator<<(std::ostream &os, const PetscStackFrame &) noexcept
191d71ae5a4SJacob Faibussowitsch   {
1920e6b6b59SJacob Faibussowitsch     os << "(unknown)";
1930e6b6b59SJacob Faibussowitsch     return os;
1940e6b6b59SJacob Faibussowitsch   }
1950e6b6b59SJacob Faibussowitsch };
1960e6b6b59SJacob Faibussowitsch 
1970e6b6b59SJacob Faibussowitsch // ==========================================================================================
1980e6b6b59SJacob Faibussowitsch // MarkedObjectMap
1990e6b6b59SJacob Faibussowitsch //
2000e6b6b59SJacob Faibussowitsch // A mapping from a PetscObjectId to a PetscEvent and (if debugging is enabled) a
2010e6b6b59SJacob Faibussowitsch // PetscStackFrame containing the location where PetscDeviceContextMarkIntentFromID was called
2020e6b6b59SJacob Faibussowitsch // ==========================================================================================
2030e6b6b59SJacob Faibussowitsch 
2040e6b6b59SJacob Faibussowitsch class MarkedObjectMap : public Petsc::RegisterFinalizeable<MarkedObjectMap> {
2050e6b6b59SJacob Faibussowitsch public:
2060e6b6b59SJacob Faibussowitsch   // Note we derive from PetscStackFrame so that the empty base class optimization can kick
2070e6b6b59SJacob Faibussowitsch   // in. If it were just a member it would still take up storage in optimized builds
2080e6b6b59SJacob Faibussowitsch   class snapshot_type : private PetscStackFrame<PetscDefined(USE_DEBUG)> {
2090e6b6b59SJacob Faibussowitsch   public:
2100e6b6b59SJacob Faibussowitsch     using frame_type = PetscStackFrame<PetscDefined(USE_DEBUG)>;
2110e6b6b59SJacob Faibussowitsch 
2120e6b6b59SJacob Faibussowitsch     snapshot_type() = default;
2130e6b6b59SJacob Faibussowitsch     snapshot_type(PetscDeviceContext, frame_type) noexcept;
2140e6b6b59SJacob Faibussowitsch 
2150e6b6b59SJacob Faibussowitsch     ~snapshot_type() noexcept;
2160e6b6b59SJacob Faibussowitsch 
2170e6b6b59SJacob Faibussowitsch     // movable
2180e6b6b59SJacob Faibussowitsch     snapshot_type(snapshot_type &&) noexcept;
2190e6b6b59SJacob Faibussowitsch     snapshot_type &operator=(snapshot_type &&) noexcept;
2200e6b6b59SJacob Faibussowitsch 
2210e6b6b59SJacob Faibussowitsch     // not copyable
2220e6b6b59SJacob Faibussowitsch     snapshot_type(const snapshot_type &) noexcept            = delete;
2230e6b6b59SJacob Faibussowitsch     snapshot_type &operator=(const snapshot_type &) noexcept = delete;
2240e6b6b59SJacob Faibussowitsch 
2250e6b6b59SJacob Faibussowitsch     PETSC_NODISCARD PetscEvent        event() const noexcept { return event_; }
2260e6b6b59SJacob Faibussowitsch     PETSC_NODISCARD const frame_type &frame() const noexcept { return *this; }
2270e6b6b59SJacob Faibussowitsch     PETSC_NODISCARD frame_type       &frame() noexcept { return *this; }
2287c441f3aSJacob Faibussowitsch 
2297c441f3aSJacob Faibussowitsch     PETSC_NODISCARD PetscObjectId dctx_id() const noexcept
2307c441f3aSJacob Faibussowitsch     {
2317c441f3aSJacob Faibussowitsch       PetscFunctionBegin;
2327c441f3aSJacob Faibussowitsch       PetscAssertAbort(event(), PETSC_COMM_SELF, PETSC_ERR_ARG_WRONGSTATE, "Snapshot %s does not contain an event!", frame().to_string().c_str());
2337c441f3aSJacob Faibussowitsch       PetscFunctionReturn(event()->dctx_id);
2347c441f3aSJacob Faibussowitsch     }
2357c441f3aSJacob Faibussowitsch 
2367c441f3aSJacob Faibussowitsch     PETSC_NODISCARD PetscErrorCode ensure_event(PetscDeviceContext) noexcept;
2377c441f3aSJacob Faibussowitsch 
2387c441f3aSJacob Faibussowitsch     friend void swap(snapshot_type &, snapshot_type &) noexcept;
2390e6b6b59SJacob Faibussowitsch 
2400e6b6b59SJacob Faibussowitsch   private:
2410e6b6b59SJacob Faibussowitsch     PetscEvent event_{}; // the state of device context when this snapshot was recorded
2420e6b6b59SJacob Faibussowitsch 
2430e6b6b59SJacob Faibussowitsch     PETSC_NODISCARD static PetscEvent init_event_(PetscDeviceContext) noexcept;
2440e6b6b59SJacob Faibussowitsch   };
2450e6b6b59SJacob Faibussowitsch 
2460e6b6b59SJacob Faibussowitsch   // the "value" each key maps to
2470e6b6b59SJacob Faibussowitsch   struct mapped_type {
2480e6b6b59SJacob Faibussowitsch     using dependency_type = std::vector<snapshot_type>;
2490e6b6b59SJacob Faibussowitsch 
2500e6b6b59SJacob Faibussowitsch     PetscMemoryAccessMode mode = PETSC_MEMORY_ACCESS_READ;
2510e6b6b59SJacob Faibussowitsch     snapshot_type         last_write{};
2520e6b6b59SJacob Faibussowitsch     dependency_type       dependencies{};
2530e6b6b59SJacob Faibussowitsch   };
2540e6b6b59SJacob Faibussowitsch 
2550e6b6b59SJacob Faibussowitsch   using map_type = std::unordered_map<PetscObjectId, mapped_type>;
2560e6b6b59SJacob Faibussowitsch 
2570e6b6b59SJacob Faibussowitsch   map_type map;
2580e6b6b59SJacob Faibussowitsch 
2590e6b6b59SJacob Faibussowitsch private:
260*271378f2SJacob Faibussowitsch   friend RegisterFinalizeable;
2610e6b6b59SJacob Faibussowitsch 
2620e6b6b59SJacob Faibussowitsch   PETSC_NODISCARD PetscErrorCode finalize_() noexcept;
2630e6b6b59SJacob Faibussowitsch };
2640e6b6b59SJacob Faibussowitsch 
2650e6b6b59SJacob Faibussowitsch // ==========================================================================================
2660e6b6b59SJacob Faibussowitsch // MarkedObejctMap Private API
2670e6b6b59SJacob Faibussowitsch // ==========================================================================================
2680e6b6b59SJacob Faibussowitsch 
269d71ae5a4SJacob Faibussowitsch inline PetscErrorCode MarkedObjectMap::finalize_() noexcept
270d71ae5a4SJacob Faibussowitsch {
2710e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
2720e6b6b59SJacob Faibussowitsch   PetscCall(PetscInfo(nullptr, "Finalizing marked object map\n"));
273*271378f2SJacob Faibussowitsch   map.clear();
2740e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
2750e6b6b59SJacob Faibussowitsch }
2760e6b6b59SJacob Faibussowitsch 
2770e6b6b59SJacob Faibussowitsch // ==========================================================================================
2780e6b6b59SJacob Faibussowitsch // MarkedObejctMap::snapshot_type Private API
2790e6b6b59SJacob Faibussowitsch // ==========================================================================================
2800e6b6b59SJacob Faibussowitsch 
281d71ae5a4SJacob Faibussowitsch inline PetscEvent MarkedObjectMap::snapshot_type::init_event_(PetscDeviceContext dctx) noexcept
282d71ae5a4SJacob Faibussowitsch {
2830e6b6b59SJacob Faibussowitsch   PetscEvent event = nullptr;
2840e6b6b59SJacob Faibussowitsch 
2850e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
2860e6b6b59SJacob Faibussowitsch   PetscCallAbort(PETSC_COMM_SELF, PetscDeviceContextCreateEvent_Private(dctx, &event));
2870e6b6b59SJacob Faibussowitsch   PetscCallAbort(PETSC_COMM_SELF, PetscDeviceContextRecordEvent_Private(dctx, event));
2880e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(event);
2890e6b6b59SJacob Faibussowitsch }
2900e6b6b59SJacob Faibussowitsch 
2910e6b6b59SJacob Faibussowitsch // ==========================================================================================
2920e6b6b59SJacob Faibussowitsch // MarkedObejctMap::snapshot_type Public API
2930e6b6b59SJacob Faibussowitsch // ==========================================================================================
2940e6b6b59SJacob Faibussowitsch 
2950e6b6b59SJacob Faibussowitsch MarkedObjectMap::snapshot_type::snapshot_type(PetscDeviceContext dctx, frame_type frame) noexcept : frame_type(std::move(frame)), event_(init_event_(dctx)) { }
2960e6b6b59SJacob Faibussowitsch 
297d71ae5a4SJacob Faibussowitsch MarkedObjectMap::snapshot_type::~snapshot_type() noexcept
298d71ae5a4SJacob Faibussowitsch {
2990e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
3000e6b6b59SJacob Faibussowitsch   PetscCallAbort(PETSC_COMM_SELF, PetscEventDestroy_Private(&event_));
3010e6b6b59SJacob Faibussowitsch   PetscFunctionReturnVoid();
3020e6b6b59SJacob Faibussowitsch }
3030e6b6b59SJacob Faibussowitsch 
3040e6b6b59SJacob Faibussowitsch // movable
3050e6b6b59SJacob Faibussowitsch MarkedObjectMap::snapshot_type::snapshot_type(snapshot_type &&other) noexcept : frame_type(std::move(other)), event_(Petsc::util::exchange(other.event_, nullptr)) { }
3060e6b6b59SJacob Faibussowitsch 
307d71ae5a4SJacob Faibussowitsch MarkedObjectMap::snapshot_type &MarkedObjectMap::snapshot_type::operator=(snapshot_type &&other) noexcept
308d71ae5a4SJacob Faibussowitsch {
3090e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
3100e6b6b59SJacob Faibussowitsch   if (this != &other) {
3110e6b6b59SJacob Faibussowitsch     frame_type::operator=(std::move(other));
3120e6b6b59SJacob Faibussowitsch     PetscCallAbort(PETSC_COMM_SELF, PetscEventDestroy_Private(&event_));
3130e6b6b59SJacob Faibussowitsch     event_ = Petsc::util::exchange(other.event_, nullptr);
3140e6b6b59SJacob Faibussowitsch   }
3150e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(*this);
3160e6b6b59SJacob Faibussowitsch }
3170e6b6b59SJacob Faibussowitsch 
3187c441f3aSJacob Faibussowitsch PetscErrorCode MarkedObjectMap::snapshot_type::ensure_event(PetscDeviceContext dctx) noexcept
3197c441f3aSJacob Faibussowitsch {
3207c441f3aSJacob Faibussowitsch   PetscFunctionBegin;
3217c441f3aSJacob Faibussowitsch   if (PetscUnlikely(!event_)) PetscCall(PetscDeviceContextCreateEvent_Private(dctx, &event_));
3227c441f3aSJacob Faibussowitsch   PetscFunctionReturn(0);
3237c441f3aSJacob Faibussowitsch }
3247c441f3aSJacob Faibussowitsch 
3257c441f3aSJacob Faibussowitsch void swap(MarkedObjectMap::snapshot_type &lhs, MarkedObjectMap::snapshot_type &rhs) noexcept
3267c441f3aSJacob Faibussowitsch {
3277c441f3aSJacob Faibussowitsch   using std::swap;
3287c441f3aSJacob Faibussowitsch 
3297c441f3aSJacob Faibussowitsch   swap(lhs.frame(), rhs.frame());
3307c441f3aSJacob Faibussowitsch   swap(lhs.event_, rhs.event_);
3317c441f3aSJacob Faibussowitsch }
3327c441f3aSJacob Faibussowitsch 
3330e6b6b59SJacob Faibussowitsch // A mapping between PetscObjectId (i.e. some PetscObject) to the list of PetscEvent's encoding
3340e6b6b59SJacob Faibussowitsch // the last time the PetscObject was accessed
3350e6b6b59SJacob Faibussowitsch static MarkedObjectMap marked_object_map;
3360e6b6b59SJacob Faibussowitsch 
3370e6b6b59SJacob Faibussowitsch // ==========================================================================================
3380e6b6b59SJacob Faibussowitsch // Utility Functions
3390e6b6b59SJacob Faibussowitsch // ==========================================================================================
3400e6b6b59SJacob Faibussowitsch 
3410e6b6b59SJacob Faibussowitsch template <typename T>
342d71ae5a4SJacob Faibussowitsch static PetscErrorCode PetscDeviceContextMapIterVisitor(PetscDeviceContext dctx, T &&callback) noexcept
343d71ae5a4SJacob Faibussowitsch {
3440e6b6b59SJacob Faibussowitsch   const auto dctx_id    = PetscObjectCast(dctx)->id;
3450e6b6b59SJacob Faibussowitsch   auto      &dctx_deps  = CxxDataCast(dctx)->deps;
3460e6b6b59SJacob Faibussowitsch   auto      &object_map = marked_object_map.map;
3470e6b6b59SJacob Faibussowitsch 
3480e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
3490e6b6b59SJacob Faibussowitsch   for (auto &&dep : dctx_deps) {
3500e6b6b59SJacob Faibussowitsch     const auto mapit = object_map.find(dep);
3510e6b6b59SJacob Faibussowitsch 
3520e6b6b59SJacob Faibussowitsch     // Need this check since the final PetscDeviceContext may run through this *after* the map
3530e6b6b59SJacob Faibussowitsch     // has been finalized (and cleared), and hence might fail to find its dependencies. This is
3540e6b6b59SJacob Faibussowitsch     // perfectly valid since the user no longer cares about dangling dependencies after PETSc
3550e6b6b59SJacob Faibussowitsch     // is finalized
3560e6b6b59SJacob Faibussowitsch     if (PetscLikely(mapit != object_map.end())) {
3570e6b6b59SJacob Faibussowitsch       auto      &deps = mapit->second.dependencies;
3580e6b6b59SJacob Faibussowitsch       const auto end  = deps.end();
3590e6b6b59SJacob Faibussowitsch       const auto it   = std::remove_if(deps.begin(), end, [&](const MarkedObjectMap::snapshot_type &obj) { return obj.dctx_id() == dctx_id; });
3600e6b6b59SJacob Faibussowitsch 
3610e6b6b59SJacob Faibussowitsch       PetscCall(callback(mapit, deps.cbegin(), static_cast<decltype(deps.cend())>(it)));
3620e6b6b59SJacob Faibussowitsch       // remove ourselves
3630e6b6b59SJacob Faibussowitsch       PetscCallCXX(deps.erase(it, end));
3640e6b6b59SJacob Faibussowitsch       // continue to next object, but erase this one if it has no more dependencies
3650e6b6b59SJacob Faibussowitsch       if (deps.empty()) PetscCallCXX(object_map.erase(mapit));
3660e6b6b59SJacob Faibussowitsch     }
3670e6b6b59SJacob Faibussowitsch   }
3680e6b6b59SJacob Faibussowitsch   PetscCallCXX(dctx_deps.clear());
3690e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
3700e6b6b59SJacob Faibussowitsch }
3710e6b6b59SJacob Faibussowitsch 
372d71ae5a4SJacob Faibussowitsch PetscErrorCode PetscDeviceContextSyncClearMap_Internal(PetscDeviceContext dctx)
373d71ae5a4SJacob Faibussowitsch {
3740e6b6b59SJacob Faibussowitsch   using map_iterator = MarkedObjectMap::map_type::const_iterator;
3750e6b6b59SJacob Faibussowitsch   using dep_iterator = MarkedObjectMap::mapped_type::dependency_type::const_iterator;
3760e6b6b59SJacob Faibussowitsch 
3770e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
3780e6b6b59SJacob Faibussowitsch   PetscCall(PetscDeviceContextMapIterVisitor(dctx, [&](map_iterator mapit, dep_iterator it, dep_iterator end) {
3790e6b6b59SJacob Faibussowitsch     PetscFunctionBegin;
3800e6b6b59SJacob Faibussowitsch     if (PetscDefined(USE_DEBUG_AND_INFO)) {
3810e6b6b59SJacob Faibussowitsch       std::ostringstream oss;
3820e6b6b59SJacob Faibussowitsch       const auto         mode = PetscMemoryAccessModeToString(mapit->second.mode);
3830e6b6b59SJacob Faibussowitsch 
3840e6b6b59SJacob Faibussowitsch       oss << "synced dctx " << PetscObjectCast(dctx)->id << ", remaining leaves for obj " << mapit->first << ": {";
3850e6b6b59SJacob Faibussowitsch       while (it != end) {
3860e6b6b59SJacob Faibussowitsch         oss << "[dctx " << it->dctx_id() << ", " << mode << ' ' << it->frame() << ']';
3870e6b6b59SJacob Faibussowitsch         if (++it != end) oss << ", ";
3880e6b6b59SJacob Faibussowitsch       }
3890e6b6b59SJacob Faibussowitsch       oss << '}';
3900e6b6b59SJacob Faibussowitsch       PetscCall(PetscInfo(nullptr, "%s\n", oss.str().c_str()));
3910e6b6b59SJacob Faibussowitsch     }
3920e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
3930e6b6b59SJacob Faibussowitsch   }));
3940e6b6b59SJacob Faibussowitsch   {
3950e6b6b59SJacob Faibussowitsch     // the recursive sync clear map call is unbounded in case of a dependenct loop so we make a
3960e6b6b59SJacob Faibussowitsch     // copy
3970e6b6b59SJacob Faibussowitsch     // clang-format off
3980e6b6b59SJacob Faibussowitsch     const std::vector<CxxData::upstream_type::value_type> upstream_copy(
3990e6b6b59SJacob Faibussowitsch       std::make_move_iterator(CxxDataCast(dctx)->upstream.begin()),
4000e6b6b59SJacob Faibussowitsch       std::make_move_iterator(CxxDataCast(dctx)->upstream.end())
4010e6b6b59SJacob Faibussowitsch     );
4020e6b6b59SJacob Faibussowitsch     // clang-format on
4030e6b6b59SJacob Faibussowitsch 
4040e6b6b59SJacob Faibussowitsch     // aftermath, clear our set of parents (to avoid infinite recursion) and mark ourselves as no
4050e6b6b59SJacob Faibussowitsch     // longer contained (while the empty graph technically *is* always contained, it is not what
4060e6b6b59SJacob Faibussowitsch     // we mean by it)
4070e6b6b59SJacob Faibussowitsch     PetscCall(CxxDataCast(dctx)->clear());
4080e6b6b59SJacob Faibussowitsch     //dctx->contained = PETSC_FALSE;
4090e6b6b59SJacob Faibussowitsch     for (auto &&upstrm : upstream_copy) {
4100e6b6b59SJacob Faibussowitsch       // check that this parent still points to what we originally thought it was
4110e6b6b59SJacob Faibussowitsch       PetscCheck(upstrm.second.id == PetscObjectCast(upstrm.first)->id, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Upstream dctx %" PetscInt64_FMT " no longer exists, now has id %" PetscInt64_FMT, upstrm.second.id, PetscObjectCast(upstrm.first)->id);
4120e6b6b59SJacob Faibussowitsch       PetscCall(PetscDeviceContextSyncClearMap_Internal(upstrm.first));
4130e6b6b59SJacob Faibussowitsch     }
4140e6b6b59SJacob Faibussowitsch   }
4150e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
4160e6b6b59SJacob Faibussowitsch }
4170e6b6b59SJacob Faibussowitsch 
418d71ae5a4SJacob Faibussowitsch PetscErrorCode PetscDeviceContextCheckNotOrphaned_Internal(PetscDeviceContext dctx)
419d71ae5a4SJacob Faibussowitsch {
4200e6b6b59SJacob Faibussowitsch   std::ostringstream oss;
4210e6b6b59SJacob Faibussowitsch   //const auto         allow = dctx->options.allow_orphans, contained = dctx->contained;
4220e6b6b59SJacob Faibussowitsch   const auto allow = true, contained = true;
4230e6b6b59SJacob Faibussowitsch   auto       wrote_to_oss = false;
4240e6b6b59SJacob Faibussowitsch   using map_iterator      = MarkedObjectMap::map_type::const_iterator;
4250e6b6b59SJacob Faibussowitsch   using dep_iterator      = MarkedObjectMap::mapped_type::dependency_type::const_iterator;
4260e6b6b59SJacob Faibussowitsch 
4270e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
4280e6b6b59SJacob Faibussowitsch   PetscCall(PetscDeviceContextMapIterVisitor(dctx, [&](map_iterator mapit, dep_iterator it, dep_iterator end) {
4290e6b6b59SJacob Faibussowitsch     PetscFunctionBegin;
4300e6b6b59SJacob Faibussowitsch     if (allow || contained) PetscFunctionReturn(0);
4310e6b6b59SJacob Faibussowitsch     wrote_to_oss = true;
4320e6b6b59SJacob Faibussowitsch     oss << "- PetscObject (id " << mapit->first << "), intent " << PetscMemoryAccessModeToString(mapit->second.mode) << ' ' << it->frame();
4330e6b6b59SJacob Faibussowitsch     if (std::distance(it, end) == 0) oss << " (orphaned)"; // we were the only dependency
4340e6b6b59SJacob Faibussowitsch     oss << '\n';
4350e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
4360e6b6b59SJacob Faibussowitsch   }));
4370e6b6b59SJacob Faibussowitsch   PetscCheck(!wrote_to_oss, PETSC_COMM_SELF, PETSC_ERR_ORDER, "Destroying PetscDeviceContext ('%s', id %" PetscInt64_FMT ") would leave the following dangling (possibly orphaned) dependants:\n%s\nMust synchronize before destroying it, or allow it to be destroyed with orphans",
4380e6b6b59SJacob Faibussowitsch              PetscObjectCast(dctx)->name ? PetscObjectCast(dctx)->name : "unnamed", PetscObjectCast(dctx)->id, oss.str().c_str());
4390e6b6b59SJacob Faibussowitsch   PetscCall(CxxDataCast(dctx)->clear());
4400e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
4410e6b6b59SJacob Faibussowitsch }
4420e6b6b59SJacob Faibussowitsch 
4437c441f3aSJacob Faibussowitsch #define DEBUG_INFO(mess, ...) PetscDebugInfo(dctx, "dctx %" PetscInt64_FMT " (%s) - obj %" PetscInt64_FMT " (%s): " mess, PetscObjectCast(dctx)->id, PetscObjectCast(dctx)->name ? PetscObjectCast(dctx)->name : "unnamed", id, name, ##__VA_ARGS__)
4440e6b6b59SJacob Faibussowitsch 
4457c441f3aSJacob Faibussowitsch // The current mode is compatible with the previous mode (i.e. read-read) so we need only
4467c441f3aSJacob Faibussowitsch // update the existing version and possibly appeand ourselves to the dependency list
4477c441f3aSJacob Faibussowitsch 
4487c441f3aSJacob Faibussowitsch template <bool use_debug>
4497c441f3aSJacob Faibussowitsch static PetscErrorCode MarkFromID_CompatibleModes(MarkedObjectMap::mapped_type &marked, PetscDeviceContext dctx, PetscObjectId id, PetscMemoryAccessMode mode, PetscStackFrame<use_debug> &frame, const char *PETSC_UNUSED name, bool *update_object_dependencies)
4507c441f3aSJacob Faibussowitsch {
4517c441f3aSJacob Faibussowitsch   const auto dctx_id             = PetscObjectCast(dctx)->id;
4527c441f3aSJacob Faibussowitsch   auto      &object_dependencies = marked.dependencies;
4530e6b6b59SJacob Faibussowitsch   const auto end                 = object_dependencies.end();
4540e6b6b59SJacob Faibussowitsch   const auto it                  = std::find_if(object_dependencies.begin(), end, [&](const MarkedObjectMap::snapshot_type &obj) { return obj.dctx_id() == dctx_id; });
4550e6b6b59SJacob Faibussowitsch 
4567c441f3aSJacob Faibussowitsch   PetscFunctionBegin;
4577c441f3aSJacob Faibussowitsch   PetscCall(DEBUG_INFO("new mode (%s) COMPATIBLE with %s mode (%s), no need to serialize\n", PetscMemoryAccessModeToString(mode), object_dependencies.empty() ? "default" : "old", PetscMemoryAccessModeToString(marked.mode)));
4580e6b6b59SJacob Faibussowitsch   if (it != end) {
4597c441f3aSJacob Faibussowitsch     using std::swap;
4607c441f3aSJacob Faibussowitsch 
4610e6b6b59SJacob Faibussowitsch     // we have been here before, all we must do is update our entry then we can bail
4620e6b6b59SJacob Faibussowitsch     PetscCall(DEBUG_INFO("found old self as dependency, updating\n"));
4630e6b6b59SJacob Faibussowitsch     PetscAssert(CxxDataCast(dctx)->deps.find(id) != CxxDataCast(dctx)->deps.end(), PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscDeviceContext %" PetscInt64_FMT " listed as dependency for object %" PetscInt64_FMT " (%s), but does not have the object in private dependency list!", dctx_id, id, name);
4647c441f3aSJacob Faibussowitsch     swap(it->frame(), frame);
4650e6b6b59SJacob Faibussowitsch     PetscCall(PetscDeviceContextRecordEvent_Private(dctx, it->event()));
4667c441f3aSJacob Faibussowitsch     *update_object_dependencies = false;
4670e6b6b59SJacob Faibussowitsch     PetscFunctionReturn(0);
4680e6b6b59SJacob Faibussowitsch   }
4690e6b6b59SJacob Faibussowitsch 
4700e6b6b59SJacob Faibussowitsch   // we have not been here before, need to serialize with the last write event (if it exists)
4710e6b6b59SJacob Faibussowitsch   // and add ourselves to the dependency list
4720e6b6b59SJacob Faibussowitsch   if (const auto event = marked.last_write.event()) PetscCall(PetscDeviceContextWaitForEvent_Private(dctx, event));
4737c441f3aSJacob Faibussowitsch   PetscFunctionReturn(0);
4747c441f3aSJacob Faibussowitsch }
4757c441f3aSJacob Faibussowitsch 
4767c441f3aSJacob Faibussowitsch template <bool use_debug>
4777c441f3aSJacob Faibussowitsch static PetscErrorCode MarkFromID_IncompatibleModes_UpdateLastWrite(MarkedObjectMap::mapped_type &marked, PetscDeviceContext dctx, PetscObjectId id, PetscMemoryAccessMode mode, PetscStackFrame<use_debug> &frame, const char *PETSC_UNUSED name, bool *update_object_dependencies)
4787c441f3aSJacob Faibussowitsch {
4797c441f3aSJacob Faibussowitsch   const auto      dctx_id    = PetscObjectCast(dctx)->id;
4807c441f3aSJacob Faibussowitsch   auto           &last_write = marked.last_write;
4817c441f3aSJacob Faibussowitsch   auto           &last_dep   = marked.dependencies.back();
4827c441f3aSJacob Faibussowitsch   PetscDeviceType dtype;
4837c441f3aSJacob Faibussowitsch 
4847c441f3aSJacob Faibussowitsch   PetscFunctionBegin;
4857c441f3aSJacob Faibussowitsch   PetscAssert(marked.dependencies.size() == 1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Can only have a single writer as dependency, have %zu!", marked.dependencies.size());
4867c441f3aSJacob Faibussowitsch   PetscCall(PetscDeviceContextGetDeviceType(dctx, &dtype));
4877c441f3aSJacob Faibussowitsch   if (last_dep.event()->dtype != dtype) {
4887c441f3aSJacob Faibussowitsch     PetscCall(DEBUG_INFO("moving last write dependency (intent %s)\n", PetscMemoryAccessModeToString(marked.mode)));
4897c441f3aSJacob Faibussowitsch     last_write = std::move(last_dep);
4907c441f3aSJacob Faibussowitsch     PetscFunctionReturn(0);
4917c441f3aSJacob Faibussowitsch   }
4927c441f3aSJacob Faibussowitsch 
4937c441f3aSJacob Faibussowitsch   // we match the device type of the dependency, we can reuse its event!
4947c441f3aSJacob Faibussowitsch   auto      &dctx_upstream_deps     = CxxDataCast(dctx)->deps;
4957c441f3aSJacob Faibussowitsch   const auto last_write_was_also_us = last_write.event() && (last_write.dctx_id() == dctx_id);
4967c441f3aSJacob Faibussowitsch   using std::swap;
4977c441f3aSJacob Faibussowitsch 
4987c441f3aSJacob Faibussowitsch   PetscCall(DEBUG_INFO("we matched the previous write dependency's (intent %s) device type (%s), swapping last dependency with last write\n", PetscMemoryAccessModeToString(marked.mode), PetscDeviceTypes[dtype]));
4997c441f3aSJacob Faibussowitsch   if (last_dep.event()->dctx_id != dctx_id) dctx_upstream_deps.emplace(id);
5007c441f3aSJacob Faibussowitsch   PetscAssert(dctx_upstream_deps.find(id) != dctx_upstream_deps.end(), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Did not find id %" PetscInt64_FMT "in object dependencies, but we have apparently recorded the last dependency %s!", id,
5017c441f3aSJacob Faibussowitsch               last_write.frame().to_string().c_str());
5027c441f3aSJacob Faibussowitsch   swap(last_write, last_dep);
5037c441f3aSJacob Faibussowitsch   if (last_write_was_also_us) {
5047c441f3aSJacob Faibussowitsch     PetscCall(DEBUG_INFO("we were also the last write event (intent %s), updating\n", PetscMemoryAccessModeToString(mode)));
5057c441f3aSJacob Faibussowitsch     // we are both the last to write *and* the last to leave a write event. This is the
5067c441f3aSJacob Faibussowitsch     // fast path, we only need to update the frame and update the recorded event
5077c441f3aSJacob Faibussowitsch     swap(last_dep.frame(), frame);
5087c441f3aSJacob Faibussowitsch     // last used to be last_write which is not guaranteed to have an event, so must
5097c441f3aSJacob Faibussowitsch     // create it now
5107c441f3aSJacob Faibussowitsch     PetscCall(last_dep.ensure_event(dctx));
5117c441f3aSJacob Faibussowitsch     PetscCall(PetscDeviceContextRecordEvent_Private(dctx, last_dep.event()));
5127c441f3aSJacob Faibussowitsch     *update_object_dependencies = false;
5137c441f3aSJacob Faibussowitsch   }
5147c441f3aSJacob Faibussowitsch   PetscFunctionReturn(0);
5157c441f3aSJacob Faibussowitsch }
5167c441f3aSJacob Faibussowitsch 
5177c441f3aSJacob Faibussowitsch // The current mode is NOT compatible with the previous mode. We must serialize with all events
5187c441f3aSJacob Faibussowitsch // in the dependency list, possibly clear it, and update the previous write event
5197c441f3aSJacob Faibussowitsch 
5207c441f3aSJacob Faibussowitsch template <bool use_debug>
5217c441f3aSJacob Faibussowitsch static PetscErrorCode MarkFromID_IncompatibleModes(MarkedObjectMap::mapped_type &marked, PetscDeviceContext dctx, PetscObjectId id, PetscMemoryAccessMode mode, PetscStackFrame<use_debug> &frame, const char *name, bool *update_object_dependencies)
5227c441f3aSJacob Faibussowitsch {
5237c441f3aSJacob Faibussowitsch   auto &old_mode            = marked.mode;
5247c441f3aSJacob Faibussowitsch   auto &object_dependencies = marked.dependencies;
5257c441f3aSJacob Faibussowitsch 
5267c441f3aSJacob Faibussowitsch   PetscFunctionBegin;
5277c441f3aSJacob Faibussowitsch   // we are NOT compatible  with the previous mode
5280e6b6b59SJacob Faibussowitsch   PetscCall(DEBUG_INFO("new mode (%s) NOT COMPATIBLE with %s mode (%s), serializing then clearing (%zu) %s\n", PetscMemoryAccessModeToString(mode), object_dependencies.empty() ? "default" : "old", PetscMemoryAccessModeToString(old_mode),
5290e6b6b59SJacob Faibussowitsch                        object_dependencies.size(), object_dependencies.size() == 1 ? "dependency" : "dependencies"));
5300e6b6b59SJacob Faibussowitsch 
5317c441f3aSJacob Faibussowitsch   for (const auto &dep : object_dependencies) PetscCall(PetscDeviceContextWaitForEvent_Private(dctx, dep.event()));
5327c441f3aSJacob Faibussowitsch   // if the previous mode wrote, update the last write node with it
5337c441f3aSJacob Faibussowitsch   if (PetscMemoryAccessWrite(old_mode)) PetscCall(MarkFromID_IncompatibleModes_UpdateLastWrite(marked, dctx, id, mode, frame, name, update_object_dependencies));
5340e6b6b59SJacob Faibussowitsch 
5350e6b6b59SJacob Faibussowitsch   old_mode = mode;
5367c441f3aSJacob Faibussowitsch   // clear out the old dependencies if are about to append ourselves
5377c441f3aSJacob Faibussowitsch   if (*update_object_dependencies) object_dependencies.clear();
5387c441f3aSJacob Faibussowitsch   PetscFunctionReturn(0);
5390e6b6b59SJacob Faibussowitsch }
5407c441f3aSJacob Faibussowitsch 
5417c441f3aSJacob Faibussowitsch template <bool use_debug>
5427c441f3aSJacob Faibussowitsch static PetscErrorCode PetscDeviceContextMarkIntentFromID_Private(PetscDeviceContext dctx, PetscObjectId id, PetscMemoryAccessMode mode, PetscStackFrame<use_debug> frame, const char *name)
5437c441f3aSJacob Faibussowitsch {
5447c441f3aSJacob Faibussowitsch   auto &marked                     = marked_object_map.map[id];
5457c441f3aSJacob Faibussowitsch   auto &object_dependencies        = marked.dependencies;
5467c441f3aSJacob Faibussowitsch   auto  update_object_dependencies = true;
5477c441f3aSJacob Faibussowitsch 
5487c441f3aSJacob Faibussowitsch   PetscFunctionBegin;
5497c441f3aSJacob Faibussowitsch   if ((marked.mode == PETSC_MEMORY_ACCESS_READ) && (mode == PETSC_MEMORY_ACCESS_READ)) {
5507c441f3aSJacob Faibussowitsch     PetscCall(MarkFromID_CompatibleModes(marked, dctx, id, mode, frame, name, &update_object_dependencies));
5517c441f3aSJacob Faibussowitsch   } else {
5527c441f3aSJacob Faibussowitsch     PetscCall(MarkFromID_IncompatibleModes(marked, dctx, id, mode, frame, name, &update_object_dependencies));
5537c441f3aSJacob Faibussowitsch   }
5547c441f3aSJacob Faibussowitsch   if (update_object_dependencies) {
5550e6b6b59SJacob Faibussowitsch     // become the new leaf by appending ourselves
5560e6b6b59SJacob Faibussowitsch     PetscCall(DEBUG_INFO("%s with intent %s\n", object_dependencies.empty() ? "dependency list is empty, creating new leaf" : "appending to existing leaves", PetscMemoryAccessModeToString(mode)));
5570e6b6b59SJacob Faibussowitsch     PetscCallCXX(object_dependencies.emplace_back(dctx, std::move(frame)));
5580e6b6b59SJacob Faibussowitsch     PetscCallCXX(CxxDataCast(dctx)->deps.emplace(id));
5590e6b6b59SJacob Faibussowitsch   }
5607c441f3aSJacob Faibussowitsch   PetscFunctionReturn(0);
5617c441f3aSJacob Faibussowitsch }
5627c441f3aSJacob Faibussowitsch 
5637c441f3aSJacob Faibussowitsch #undef DEBUG_INFO
5640e6b6b59SJacob Faibussowitsch 
5650e6b6b59SJacob Faibussowitsch /*@C
5660e6b6b59SJacob Faibussowitsch   PetscDeviceContextMarkIntentFromID - Indicate a `PetscDeviceContext`s access intent to the
5670e6b6b59SJacob Faibussowitsch   auto-dependency system
5680e6b6b59SJacob Faibussowitsch 
5690e6b6b59SJacob Faibussowitsch   Not Collective
5700e6b6b59SJacob Faibussowitsch 
5710e6b6b59SJacob Faibussowitsch   Input Parameters:
5720e6b6b59SJacob Faibussowitsch + dctx - The `PetscDeviceContext`
5730e6b6b59SJacob Faibussowitsch . id   - The `PetscObjectId` to mark
5740e6b6b59SJacob Faibussowitsch . mode - The desired access intent
5750e6b6b59SJacob Faibussowitsch - name - The object name (for debug purposes, ignored in optimized builds)
5760e6b6b59SJacob Faibussowitsch 
5770e6b6b59SJacob Faibussowitsch   Notes:
5780e6b6b59SJacob Faibussowitsch   This routine formally informs the dependency system that `dctx` will access the object
5790e6b6b59SJacob Faibussowitsch   represented by `id` with `mode` and adds `dctx` to `id`'s list of dependencies (termed
5800e6b6b59SJacob Faibussowitsch   "leaves").
5810e6b6b59SJacob Faibussowitsch 
5820e6b6b59SJacob Faibussowitsch   If the existing set of leaves have an incompatible `PetscMemoryAccessMode` to `mode`, `dctx`
5830e6b6b59SJacob Faibussowitsch   will be serialized against them.
5840e6b6b59SJacob Faibussowitsch 
5850e6b6b59SJacob Faibussowitsch   Level: intermediate
5860e6b6b59SJacob Faibussowitsch 
5870e6b6b59SJacob Faibussowitsch .seealso: `PetscDeviceContextWaitForContext()`, `PetscDeviceContextSynchronize()`,
5880e6b6b59SJacob Faibussowitsch `PetscObjectGetId()`, `PetscMemoryAccessMode`
5890e6b6b59SJacob Faibussowitsch @*/
590d71ae5a4SJacob Faibussowitsch PetscErrorCode PetscDeviceContextMarkIntentFromID(PetscDeviceContext dctx, PetscObjectId id, PetscMemoryAccessMode mode, const char name[])
591d71ae5a4SJacob Faibussowitsch {
5920e6b6b59SJacob Faibussowitsch #if PetscDefined(USE_DEBUG)
5930e6b6b59SJacob Faibussowitsch   const auto index    = petscstack.currentsize > 2 ? petscstack.currentsize - 2 : 0;
5940e6b6b59SJacob Faibussowitsch   const auto file     = petscstack.file[index];
5950e6b6b59SJacob Faibussowitsch   const auto function = petscstack.function[index];
5960e6b6b59SJacob Faibussowitsch   const auto line     = petscstack.line[index];
5970e6b6b59SJacob Faibussowitsch #else
5980e6b6b59SJacob Faibussowitsch   constexpr const char *file     = nullptr;
5990e6b6b59SJacob Faibussowitsch   constexpr const char *function = nullptr;
6000e6b6b59SJacob Faibussowitsch   constexpr auto        line     = 0;
6010e6b6b59SJacob Faibussowitsch #endif
6020e6b6b59SJacob Faibussowitsch 
6030e6b6b59SJacob Faibussowitsch   PetscFunctionBegin;
6040e6b6b59SJacob Faibussowitsch   PetscCall(PetscDeviceContextGetOptionalNullContext_Internal(&dctx));
6050e6b6b59SJacob Faibussowitsch   if (name) PetscValidCharPointer(name, 4);
6060e6b6b59SJacob Faibussowitsch   PetscCall(marked_object_map.register_finalize());
6076a4a1270SPierre Jolivet   PetscCall(PetscLogEventBegin(DCONTEXT_Mark, dctx, nullptr, nullptr, nullptr));
6080e6b6b59SJacob Faibussowitsch   PetscCall(PetscDeviceContextMarkIntentFromID_Private(dctx, id, mode, MarkedObjectMap::snapshot_type::frame_type{file, function, line}, name ? name : "unknown object"));
6096a4a1270SPierre Jolivet   PetscCall(PetscLogEventEnd(DCONTEXT_Mark, dctx, nullptr, nullptr, nullptr));
6100e6b6b59SJacob Faibussowitsch   PetscFunctionReturn(0);
6110e6b6b59SJacob Faibussowitsch }
612a966a306SPierre Jolivet 
613a966a306SPierre Jolivet #if defined(__clang__)
614a966a306SPierre Jolivet   #pragma clang diagnostic pop
615a966a306SPierre Jolivet #endif
616