hessian-matrix- все статьи тега

Самый быстрый способ создать разреженную матрицу вида 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 ======== ...