linear-algebra- все статьи тега


Очень медленная производительность анализа cusparse csrsv

Я написал решатель сопряженного градиента (для линейной системы уравнений) с предварительным условием LU, я использовал статьи доктора Максима Наумова на исследовательском сообществе nvidia в качестве руководства, шаг обновления остатков, который требует решения нижней треугольной матричной системы, а затем решение верхней треугольной матричной системы делится на две фазы: ]} Фаза анализа (которая использует модель разреженности и определяет уровень распараллеливания). фаза решения сам. Сог ...

Самый быстрый способ создать разреженную матрицу вида A. T * diag(b) * A + C?

Я пытаюсь оптимизировать фрагмент кода, который решает большую разреженную нелинейную систему, используя метод внутренней точки. На этапе обновления это включает вычисление матрицы Гессена H, градиента g, а затем решение для d в H * d = -g, чтобы получить новое направление поиска. Матрица Гессена имеет симметричную трехдиагональную структуру вида: A. T * diag (b) * A + C Я убежал.line_profiler о конкретной функции в вопросе: Line # Hits Time Per Hit % Time Line Contents ======== ...

Каковы наиболее широко используемые библиотеки векторной / матричной математики / линейной алгебры C++ и их компромиссы затрат и выгод? [закрытый]

похоже, что многие проекты медленно приходят к необходимости делать матричную математику и попадают в ловушку первого построения некоторых векторных классов и медленного добавления функциональности, пока их не поймают, создавая наполовину заданную пользовательскую библиотеку линейной алгебры и в зависимости от нее. Я хотел бы избежать этого, не создавая зависимости от некоторой тангенциально связанной библиотеки (например, OpenCV, OpenSceneGraph). каковы обычно используемые матричная математ ...