Если мы выделили память в динамической библиотеке и попытаемся освободить из другого модуля могут возникнуть проблемы. Кажется это не характерно для линукса, но под виндой это обычное правило.

Ну так этим занимается код из инициализации/финализации dll-ки, так что никаких проблем с разными аллокаторами

Там вроде дело в разных рантаймах С/С++. Если добиться, чтоб он был один, то всё можно

так у меня owner'ов множество

Так если у нас остались живые shared_ptr они не дадут удалить объект при выгрузке dll. То есть освобождение памяти произойдет из другого модуля. Если живых shared_ptr не было, то и проблемы нет.

Помогите определиться. Задача - найти индекс любого ненулевого элемента вектора. (Весь код запускаем на гпу, по факту это часть поиска обратной матрицы). Можно просто идти последовательно, пока не выполнится условие. Я хочу завести по потоку на каждый элемент, поток проверяет условие и записывает по указателю индекс, если условие выполнено. Синхронизации нет, результат недетерминирован. Но цель найти любой элемент, так что по хорошему все будет ОК. Вправе ли я так сделать, или все таки есть вероятность разбудить Ктулху?

Можно наверное что-то типа parallel_for использовать? Он сам параллелит, я так понял элементы не связаны. Или на гпу оно не взлетит?

Это вообще мелкая  функция в коде поиска обратной матрицы методом Гаусса-Жордана. Ее по любому выполнять на гпу надо, просто есть вариант запустить в одном потоке, что для гпу странно, либо как выше описал. Второй вариант предпочтительнее, просто я не уверен, могу ли я так делать

Потому что с точки зрения стандарта это же, наверное, UB

Просто руками запускать потоки, а ещё следить чтобы не вышло рейсов итд это как-то неправильно. Я так понял нужно записать один индекс

С помощью чего на GPU запускаете?

cuda

Если в коде, то по факту выбор между двумя вариантами

 template <typename T>
    __global__ void first_not_null_row(T* A, int* idx, int col, int order)
    {
        for (int i = 0; i < order; ++i)
        {
            if (A[i * order + col] != Zero_v<T>)
            {
                *idx = i;
                return;
            }
        }
    }

    template <typename T>
    __global__ void first_not_null_row(T* A, int* idx, int col, int order)
    {
        int i = blockIdx.x * blockDim.x + threadIdx.x;
        if (A[i * order + col] != Zero_v<T>)
        {
            *idx = i;
        }
    }

с какого C++ можно юзать циклы в constexpr?

С 14-го

то есть я могу compile-time хеширование строк сделать?

Можешь

Да и в С++11 можешь при желании

Рекурсия наше всё

Дарова