From 38317c31398b21de9c782f30da233612e0583767 Mon Sep 17 00:00:00 2001
From: "Ralf W. Grosse-Kunstleve" <rwgk@google.com>
Date: Sun, 8 Oct 2023 22:22:37 -0700
Subject: [PATCH] LazyInitializeAtLeastOnceDestroyNever v1

---
 include/pybind11/numpy.h | 39 +++++++++++++++++++++++++++++++++------
 1 file changed, 33 insertions(+), 6 deletions(-)

diff --git a/include/pybind11/numpy.h b/include/pybind11/numpy.h
index 23c38660e..5a60f18b2 100644
--- a/include/pybind11/numpy.h
+++ b/include/pybind11/numpy.h
@@ -14,12 +14,14 @@
 
 #include <algorithm>
 #include <array>
+#include <cassert>
 #include <cstdint>
 #include <cstdlib>
 #include <cstring>
 #include <functional>
 #include <numeric>
 #include <sstream>
+#include <stdalign.h>
 #include <string>
 #include <type_traits>
 #include <typeindex>
@@ -42,6 +44,30 @@ class array; // Forward declaration
 
 PYBIND11_NAMESPACE_BEGIN(detail)
 
+// Main author of this class: jbms@
+template <typename T>
+class LazyInitializeAtLeastOnceDestroyNever {
+public:
+    template <typename Initialize>
+    T &Get(Initialize &&initialize) {
+        if (!initialized_) {
+            assert(PyGILState_Check());
+            // Multiple threads may run this concurrently, but that is fine.
+            auto value = initialize(); // May release and re-acquire the GIL.
+            if (!initialized_) {       // This runs with the GIL held,
+                new                    // therefore this is reached only once.
+                    (reinterpret_cast<T *>(value_storage_)) T(std::move(value));
+                initialized_ = true;
+            }
+        }
+        return *reinterpret_cast<T *>(value_storage_);
+    }
+
+private:
+    alignas(T) char value_storage_[sizeof(T)];
+    bool initialized_ = false;
+};
+
 template <>
 struct handle_type_name<array> {
     static constexpr auto name = const_name("numpy.ndarray");
@@ -206,8 +232,8 @@ struct npy_api {
     };
 
     static npy_api &get() {
-        static npy_api api = lookup();
-        return api;
+        static LazyInitializeAtLeastOnceDestroyNever<npy_api> api_init;
+        return api_init.Get(lookup);
     }
 
     bool PyArray_Check_(PyObject *obj) const {
@@ -643,10 +669,11 @@ public:
     char flags() const { return detail::array_descriptor_proxy(m_ptr)->flags; }
 
 private:
-    static object _dtype_from_pep3118() {
-        module_ m = detail::import_numpy_core_submodule("_internal");
-        static PyObject *obj = m.attr("_dtype_from_pep3118").cast<object>().release().ptr();
-        return reinterpret_borrow<object>(obj);
+    static object &_dtype_from_pep3118() {
+        static detail::LazyInitializeAtLeastOnceDestroyNever<object> imported_obj;
+        return imported_obj.Get([]() {
+            return detail::import_numpy_core_submodule("_internal").attr("_dtype_from_pep3118");
+        });
     }
 
     dtype strip_padding(ssize_t itemsize) {