Частичная оболочка для задачи оптимизации с невыпуклыми ограничениями

Основная информация

• ID статьи: 2510.22223
• Название: Partial Envelope for Optimization Problem with Nonconvex Constraints
• Авторы: Xiaoyin Hu, Xin Liu, Kim-Chuan Toh, Nachuan Xiao
• Классификация: math.OC (Математическая оптимизация и управление)
• Дата подачи: 25 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.22223v1

Аннотация

В данной работе исследуется задача нелинейной оптимизации с ограничениями (NCP) вида $\{x \in \mathcal{X}: c(x) = 0\}$, где $\mathcal{X}$ — замкнутое выпуклое подмножество $\mathbb{R}^n$. На основе прямо-обратной оболочечной схемы на $\mathcal{X}$ авторы предлагают метод прямо-обратной частичной оболочки (FBSE). Данный метод устраняет ограничение $x \in \mathcal{X}$ посредством специально разработанной оболочечной схемы, сохраняя при этом ограничение $x \in \mathcal{M} := \{x \in \mathbb{R}^n: c(x) = 0\}$. Доказано, что FBSE хорошо определена и локально липшицева гладкая в окрестности $\mathcal{M}$, и что NCP и FBSE имеют одинаковые точки первого порядка стационарности в окрестности $\mathcal{X} \cap \mathcal{M}$. Кроме того, авторы разработали неточный метод проективного градиентного спуска и установили его глобальную сходимость и сложность итераций $O(\varepsilon^{-2})$.

Исследовательский контекст и мотивация

Решаемая проблема

В работе исследуется задача оптимизации с ограничениями следующего вида: $\min_{x \in \mathbb{R}^n} f(x) + I_{\mathcal{X}}(x) \quad \text{при условии} \quad x \in \mathcal{M} := \{x \in \mathbb{R}^n: c(x) = 0\}$

где $I_{\mathcal{X}}(x)$ — индикаторная функция множества $\mathcal{X}$, а $\mathcal{X}$ — компактное выпуклое подмножество с легко вычисляемым оператором проекции. Данная задача эквивалентна минимизации $f(x)$ на множестве $\{x \in \mathcal{X}: c(x) = 0\}$.

Важность проблемы

Данный класс задач оптимизации охватывает несколько важных моделей оптимизации:

1. Оптимизация с равенствами и неравенствами
2. Задачи конического программирования (например, полуопределённое программирование)
3. Оптимизация на многообразиях

Области применения включают:

• Задачи машинного обучения
• Обработку сигналов
• Проектирование механизмов и другие

Ограничения существующих методов

Ограничения традиционных оболочечных методов:

1. Прямо-обратная оболочка (Forward-Backward Envelope) и оболочка Моро зависят от выпуклости множества ограничений
2. При рассмотрении NCP как задачи без ограничений с индикаторной функцией $I_{\mathcal{X} \cap \mathcal{M}}$ оболочечная функция становится негладкой из-за невыпуклости $\mathcal{M} \cap \mathcal{X}$
3. Проекция на $\mathcal{X} \cap \mathcal{M}$ вычислительно дорогостояща, даже если $\Pi_{\mathcal{M}}$ и $\Pi_{\mathcal{X}}$ легко вычисляются

Ограничения методов растворения ограничений: Недавно предложенные методы растворения ограничений (constraint dissolving approach) разделяют ограничения через точную штрафную функцию: $\min_{x \in \mathcal{X}} h_{cdf}(x) := f(A(x)) + \frac{\beta}{2}\|c(x)\|^2$

но требуют выбора штрафного параметра $\beta$, что на практике является сложной задачей.

Исследовательская мотивация

Авторы ставят центральный вопрос:

Можно ли разработать оболочечный метод для задач оптимизации с ограничениями вида NCP, который не вводит никаких штрафных параметров?

Основные вклады

1. Предложение метода прямо-обратной частичной оболочки (FBSE): новая оболочечная схема, которая устраняет только выпуклое ограничение $x \in \mathcal{X}$, сохраняя невыпуклое равенство $c(x) = 0$, без введения штрафных параметров
2. Установление теоретической эквивалентности: доказано, что в окрестности $\mathcal{X} \cap \mathcal{M}$ NCP и FBSE имеют одинаковые точки первого порядка стационарности (для достаточно малого параметра оболочки $\mu$)
3. Доказательство хороших свойств гладкости: показано, что FBSE хорошо определена в окрестности $\mathcal{M}$, непрерывно дифференцируема, и её градиент локально липшицев непрерывен
4. Разработка эффективного алгоритма: предложен неточный метод проективного градиентного спуска, избегающий вычисления членов Гессиана в полном градиенте $H(x)$, с доказанными:
   • Глобальной сходимостью
   • Сложностью итераций $O(\varepsilon^{-2})$
5. Численная верификация: эксперименты на задачах оптимизации с ограничениями конуса полуопределённости показывают, что метод превосходит существующие решатели по точности и эффективности

Подробное описание метода

Определение задачи

Исходная задача (NCP): $\min_{x \in \mathbb{R}^n} f(x) + I_{\mathcal{X}}(x) \quad \text{при условии} \quad c(x) = 0$

Ключевые предположения (Assumption 1.1):

1. $f: \mathbb{R}^n \to \mathbb{R}$ дважды дифференцируема на $\mathbb{R}^n$
2. $c: \mathbb{R}^n \to \mathbb{R}^p$ дважды дифференцируема с локально липшицевой второй производной
3. Условие невырожденности ограничений: для всех $x \in \mathcal{K} := \mathcal{X} \cap \mathcal{M}$, $\nabla c(x)^\top \text{lin}(T_{\mathcal{X}}(x)) = \mathbb{R}^p$

Архитектура основного метода

1. Проективное отображение (Projective Mapping)

Определяется отображение $Q: \mathbb{R}^n \to \mathbb{S}^{n \times n}_+$, удовлетворяющее:

• $Q(x)$ локально липшицева гладкая
• Для всех $x \in \mathcal{X}$, $\text{null}(Q(x)) = \text{range}(N_{\mathcal{X}}(x))$

Отображение растворения ограничений: $A(x) = x - Q(x)\nabla c(x)(\nabla c(x)^\top Q(x)\nabla c(x) + \tau(x)I_p)^{-1}c(x)$

где $\tau(x) := L_\tau(\|c(x)\|^2 + \text{dist}(x, \mathcal{X})^2)$, $L_\tau > 0$ — предустановленный параметр.

2. Прямо-обратная частичная оболочка (FBSE)

Задача FBSE: $\min_{x \in \mathbb{R}^n} \psi_\mu(x) \quad \text{при условии} \quad x \in \mathcal{M}$

где функция частичной оболочки определяется как: $\psi_\mu(x) := \min_{w \in \mathcal{X}} f(x) + \langle J(x)\nabla f(x), w - x \rangle + \frac{1}{2\mu}\|w - x\|^2$

Ключевое отображение: $J(x) := I_n - \nabla c(x)(\nabla c(x)^\top Q(x)\nabla c(x) + \tau(x)I_p)^{-1}\nabla c(x)^\top Q(x)$

Оптимальное решение: $T_\mu(x) := \arg\min_{w \in \mathcal{X}} f(x) + \langle J(x)\nabla f(x), w - x \rangle + \frac{1}{2\mu}\|w - x\|^2 = \Pi_{\mathcal{X}}(x - \mu J(x)\nabla f(x))$

3. Выражение для градиента

Согласно Lemma 3.7, градиент $\psi_\mu$ имеет вид: $\nabla \psi_\mu(x) = \frac{1}{\mu}(I_n - \mu H(x))(x - T_\mu(x)) + (I_n - J(x))\nabla f(x)$

где $H(x) = J(x)\nabla^2 f(x) + \nabla J(x)[\nabla f(x)]$.

Технические инновации

1. Стратегия частичной оболочки

Ключевая инновация: в отличие от традиционных оболочечных методов, обрабатывающих всё множество ограничений $\mathcal{X} \cap \mathcal{M}$, FBSE использует стратегию "частичной оболочки":

• Устраняет выпуклое ограничение $x \in \mathcal{X}$ через оболочечную технику
• Сохраняет невыпуклое равенство $c(x) = 0$
• Избегает вычислительных трудностей проекции на невыпуклое множество

2. Специальные свойства отображения $J(x)$

Lemma 3.2: для всех $x \in \mathcal{X} \cap \mathcal{M}$, $J(x)\nabla c(x) = 0$

Lemma 3.3: для всех $d \in \text{range}(N_{\mathcal{X}}(x))$, $J(x)d = d$

Эти свойства гарантируют:

• В допустимых точках $J(x)$ проецирует градиент в касательное пространство
• Сохраняется информация о направлениях нормального конуса

3. Теория эквивалентности

Proposition 3.9: если $x \in \mathcal{X} \cap \mathcal{M}$ — точка первого порядка стационарности NCP, то $x$ — точка первого порядка стационарности FBSE.

Theorem 3.10 (основной теоретический результат): для достаточно малого $\mu \leq \mu_{\max}$, если $x \in \mathcal{K}_\rho$ — точка первого порядка стационарности FBSE, то $x$ — точка первого порядка стационарности NCP.

Ключ доказательства: через доказательство $\|T_\mu(x) - x\| = 0$, в сочетании с положительной определённостью нижней границы $\nabla c(x)^\top Q(T_\mu(x))\nabla c(x)$ ($\geq \sigma_Q/4$).

4. Неточный метод градиента

Конструкция алгоритма (уравнение 3.20): $g_k = \frac{1}{\mu}(I_n - \nabla c(x_k)\nabla c(x_k)^\dagger)(x_k - T_\mu(x_k))$$x_{k+1} = \Pi_{\mathcal{M}}(x_k - \eta_k g_k)$

Преимущества:

• Использует $\frac{1}{\mu}(x - T_\mu(x))$ как неточную оценку $\nabla \psi_\mu$
• Избегает вычисления $H(x)$ (содержит Гессиан)
• Проекция на $\text{null}(\nabla c(x_k)^\top)$ (касательное пространство $\mathcal{M}$)

Proposition 3.13: свойство достаточного убывания $\langle (I_n - \nabla c(x)\nabla c(x)^\dagger)\nabla \psi_\mu(x), T_\mu(x) - x \rangle \leq -\frac{1}{2\mu}\left(\frac{\sigma_Q}{8M_QM_c^2 + 2\sigma_Q}\right)^2\|x - T_\mu(x)\|^2$

Экспериментальная установка

Наборы данных

Эксперимент 1: Конус полуопределённости и сферическое ограничение

Задача оптимизации: $\min_{X \in \mathbb{S}^{n \times n}} \langle B, X \rangle + \frac{1}{2}\langle X, H(X) \rangle + \frac{\nu}{6}\|X\|_F^3$$\text{при условии} \quad \|X\|_F^2 = 1, \quad X \succeq 0, \quad \|X\|_2 \leq M$

• Тестируемые размеры: $n \in \{10, 20, 30, 50\}$
• $B \in \mathbb{S}^{n \times n}$ генерируется случайно (стандартное нормальное распределение)
• $H: \mathbb{S}^{n \times n} \to \mathbb{S}^{n \times n}$ — самосопряжённое линейное отображение
• Параметры: $\nu = 1.0$, $M = 10^6$, $\mu = 0.01$

Эксперимент 2: Конус полуопределённости и линейные ограничения

Задача оптимизации: $\min_{X \in \mathbb{R}^{n \times n}} \langle B_0, X \rangle + \frac{1}{2}\langle X, H(X) \rangle + \frac{\nu}{6}\|X\|_F^3$$\text{при условии} \quad \mathcal{B}(X) = b, \quad X \succeq 0, \quad \|X\|_2 \leq M$

• Тестируемые размеры: $n \in \{10, 20, 30, 50\}$
• $\mathcal{B}: \mathbb{S}^{n \times n} \to \mathbb{R}^m$ — линейное отображение
• Параметры: $\nu = 1.0$, $\mu = 0.001$

Метрики оценки

1. Стационарность: $\text{dist}(0, \nabla f(y) + N_{\mathcal{X}}(y) + \text{range}(\nabla c(y)))$, где $y = \Pi_{\mathcal{X}}(x)$
2. Нарушение допустимости: $\|c(\Pi_{\mathcal{X}}(x))\|$
3. Значение целевой функции
4. Количество итераций и количество вычислений функции
5. Время CPU (в секундах)

Методы сравнения

1. PGD: предложенный в работе метод проективного градиентного спуска (с адаптивным шагом Барзилаи-Борвейна и немонотонным поиском линии)
2. TRCON: решатель оптимизации с доверительной областью из SciPy
3. SLSQP: последовательное программирование наименьших квадратов из SciPy
4. RGD: риманов градиентный спуск из PyManopt
5. RCG: риманов сопряжённый градиент из PyManopt

Детали реализации

• Среда программирования: Python 3.12.2
• Оборудование: AMD Ryzen 7 5700 CPU, 16 GB RAM
• Допуск: $10^{-5}$
• Максимальное время выполнения: 300 секунд
• Оператор проекции (Эксперимент 1): $Q(X): Y \mapsto \Phi(X^2\Theta_M(X)^2 Y)$ где $\Phi(M) = (M + M^\top)/2$ — оператор симметризации

Результаты экспериментов

Основные результаты

Эксперимент 1: Конус полуопределённости и сферическое ограничение (Таблица 4)

Ключевые выводы:

1. Высокая точность: PGD и TRCON достигают допуска $10^{-5}$, значения целей совпадают
2. Эффективность: PGD при $n=50$ работает в 28.7 раз быстрее TRCON (1.082s vs 31.039s)
3. Отказ риманова метода: показатели стационарности RGD и RCG на уровне $10^{-1}$, далеко от сходимости
4. Отказ SLSQP: превышает время ожидания при $n \geq 30$

Эксперимент 2: Конус полуопределённости и линейные ограничения (Таблица 5)

Ключевые выводы:

1. Масштабируемость: PGD решает при $n=50$, когда TRCON превышает время ожидания
2. Преимущество в скорости: при $n=30$ PGD работает в 5.7 раз быстрее TRCON
3. Полный отказ SLSQP: все тестовые примеры не сходятся или численно нестабильны

Экспериментальные выводы

1. Верификация эквивалентности: эксперименты подтверждают теоретическую эквивалентность NCP и FBSE в точках первого порядка стационарности (PGD и TRCON получают одинаковые значения целей)
2. Эффективность неточного градиента: использование $\frac{1}{\mu}(x - T_\mu(x))$ как приближённого градиента, избегая вычисления $H(x)$, всё ещё гарантирует сходимость
3. Ограничения риманова метода:
   • RGD/RCG оптимизируют на сферическом многообразии, но не учитывают PSD ограничение
   • Плохие показатели стационарности указывают на неполную сходимость к стационарной точке NCP
4. Вызовы универсальных решателей:
   • SLSQP чувствителен к невыпуклым ограничениям, численно нестабилен
   • TRCON надёжен, но вычислительно дорогостоящ
5. Преимущества FBSE:
   • Преобразует задачу с невыпуклыми ограничениями в задачу с равенствами
   • Сохраняет структуру задачи
   • Позволяет разработку эффективных алгоритмов

Связанные работы

Оболочечные методы

1. Прямо-обратная оболочка (Forward-Backward Envelope)

• Patrinos & Bemporad (2013): первое предложение для выпуклой композитной оптимизации
• Stella et al. (2017): квазиньютоновские методы
• Themelis et al. (2018): алгоритмы с немонотонным поиском линии
• Ограничение: требует выпуклость $\mathcal{X}$, неприменимо к $\mathcal{X} \cap \mathcal{M}$

2. Оболочка Моро

• Moreau (1965): классическая техника сглаживания
• Davis & Drusvyatskiy (2019): случайный субградиентный метод для слабо выпуклых функций
• Ограничение: подзадачи обычно не имеют замкнутого решения, практически невычислимы

Методы оптимизации с ограничениями

1. Методы растворения ограничений

• Xiao et al. (2025): предложение отображения растворения ограничений $A(x)$ и точной штрафной функции
• Отличие данной работы: FBSE избегает введения штрафного параметра, напрямую обрабатывает равенства

2. Традиционные методы

• Последовательное квадратичное программирование (SQP): требует информацию второго порядка
• Увеличенный метод Лагранжа: требует настройку штрафного параметра и множителя Лагранжа
• Преимущество данной работы: требует только информацию первого порядка, простой выбор параметров

Оптимизация на многообразиях

• Absil et al. (2008): алгоритмы оптимизации на многообразиях
• Связь с данной работой: когда $\mathcal{M}$ — многообразие, FBSE можно рассматривать как частный случай оптимизации на многообразиях
• Расширение данной работы: обработка более общих нелинейных равенств

Заключение и обсуждение

Основные выводы

1. Теоретические вклады:
   • Установление эквивалентности NCP и FBSE в точках первого порядка стационарности (Theorem 3.10)
   • Доказательство липшицевой гладкости $\psi_\mu$ (Lemma 3.7)
   • Соотношение между $\varepsilon$-стационарными точками (Theorem 3.12)
2. Вклады в алгоритмы:
   • Предложение неточного метода проективного градиентного спуска, избегающего вычисления Гессиана
   • Доказательство сложности итераций $O(\varepsilon^{-2})$ (Theorem 3.17)
   • Экспериментальная верификация эффективности алгоритма
3. Методологические вклады:
   • Стратегия "частичной оболочки": избирательная обработка ограничений
   • Отсутствие штрафного параметра: избегает трудностей настройки параметров
   • Модульная конструкция: может комбинироваться с существующими решателями равенств

Ограничения

1. Теоретические предположения

• Условие невырожденности ограничений (Assumption 1.1(3)): требует $\nabla c(x)^\top \text{lin}(T_{\mathcal{X}}(x)) = \mathbb{R}^p$, может не выполняться в некоторых приложениях
• Локальные свойства: эквивалентность справедлива только в окрестности $\mathcal{K}_\rho$, где $\rho$ зависит от нескольких констант

2. Выбор параметров

• Параметр оболочки $\mu$: требует $\mu \leq \mu_{\max}$, где $\mu_{\max}$ вычисляется через несколько трудно оцениваемых констант (Таблицы 1-2)
• На практике: статья рекомендует адаптивную оценку или методы Монте-Карло, но не обсуждает подробно

3. Конструкция проективного отображения

• Зависит от структуры задачи: требует конструкции $Q(x)$, удовлетворяющей Assumption 1.2 для конкретного $\mathcal{X}$
• Таблица 3 охватывает только типичные случаи: для сложных ограничений конструкция $Q(x)$ может быть нетривиальной

4. Численные эксперименты

• Ограниченные размеры тестов: максимум $n=50$, большие задачи не тестировались
• Узкий класс задач: тестировались только SDP задачи, другие приложения не проверены

Направления будущих исследований

1. Теоретические расширения:
   • Ослабление условия невырожденности ограничений
   • Анализ глобальной сходимости (вместо локальной эквивалентности)
   • Исследование свойств сходимости второго порядка
2. Улучшения алгоритма:
   • Разработка стратегий адаптивного выбора $\mu$
   • Комбинирование со второй информацией (например, BFGS) для ускорения
   • Проектирование специализированных алгоритмов для конкретных структур
3. Расширение приложений:
   • Тестирование на большем числе приложений (машинное обучение, обработка сигналов)
   • Обработка больших задач
   • Расширение на неравенства
4. Оболочка Моро (Moreau half-envelope):
   • Статья упоминает, но не обсуждает подробно $\psi_{M,\mu}(x) := \arg\min_{y \in \mathcal{X}} f(y) + \frac{1}{2\mu}\|y - x\|^2$
   • Может быть применима к негладким целевым функциям

Глубокая оценка

Преимущества

1. Теоретическая строгость

• Полная теоретическая схема: от хорошей определённости (Lemma 3.1) к эквивалентности (Theorem 3.10) и сходимости (Theorem 3.17), логика безупречна
• Богатые технические леммы: Lemma 3.2-3.8 обеспечивают прочную основу для основных теорем
• Явные константы: Таблицы 1-2 детально перечисляют все релевантные константы, облегчая теоретический анализ

2. Методологическая инновативность

• Идея частичной оболочки: впервые предложена стратегия избирательной обработки ограничений, преодолевает ограничения традиционных оболочечных методов
• Конструкция без штрафного параметра: в сравнении с методами растворения ограничений избегает трудностей настройки штрафного параметра
• Техника неточного градиента: умело использует $\frac{1}{\mu}(x - T_\mu(x))$, снижает вычислительную сложность

3. Практическая применимость алгоритма

• Лёгкая реализация: проекции на $\mathcal{M}$ и $\mathcal{X}$ имеют известные методы
• Численная стабильность: в экспериментах показатели стационарности достигают $10^{-6}$
• Вычислительная эффективность: значительное ускорение по сравнению с TRCON (максимум 28.7 раз)

4. Ясность изложения

• Логичная структура: от мотивации к теории и экспериментам, иерархия ясна
• Стандартная нотация: Section 2.1 специально определяет символы, избегает путаницы
• Подробные доказательства: ключевые теоремы имеют чёткие шаги доказательства

Недостатки

1. Теоретический разрыв

• Практичность $\mu_{\max}$: определение в Таблице 2 включает $\sup$ и $\inf$, практическое вычисление затруднено
• Отсутствие глобальных свойств: не обсуждается, как алгоритм входит в окрестность $\mathcal{K}_\rho$
• Зависимость от констант: $\rho$ и $\mu_{\max}$ зависят от нескольких трудно оцениваемых констант, может привести к консервативным оценкам

2. Ограничения экспериментов

• Неполное сравнение:
  • Не сравнивается со специализированными SDP решателями (SDPT3, MOSEK)
  • Не тестируется увеличенный метод Лагранжа
• Недостаточное разнообразие задач: тестировались только SDP задачи, не охватываются другие приложения (оптимизация на многообразиях, машинное обучение)
• Неизвестная масштабируемость: максимум $n=50$, производительность на больших задачах неизвестна

3. Применимость метода

• Конструкция проективного отображения:
  • Таблица 3 предоставляет $Q(x)$ для 4 типичных ограничений
  • Для сложных ограничений (например, пересечение нескольких ограничений) конструкция $Q(x)$ может быть сложной
• Ограничения предположений: условие невырожденности ограничений может не выполняться в некоторых задачах

4. Технические детали

• Выбор размера шага: уравнение (3.22) даёт $\eta_{\max}$, но реальный алгоритм использует шаг Барзилаи-Борвейна, связь не ясна
• Требование начальной точки: алгоритм требует $x_0 \in \mathcal{X} \cap \mathcal{M}$, но как получить допустимую начальную точку не обсуждается
• Оболочка Моро: упоминается, но не анализируется подробно, это упущение

Влияние

1. Вклад в область

• Теоретическое значение:
  • Расширяет применимость оболочечных методов (от выпуклых к смешанным ограничениям)
  • Предоставляет новый теоретический инструмент (схема частичной оболочки)
• Методологическое значение:
  • Вдохновляет идею "избирательной обработки ограничений"
  • Предлагает новую перспективу для оптимизации с невыпуклыми ограничениями

2. Практическая ценность

• Непосредственное применение: может использоваться для решения SDP, оптимизации на многообразиях и других задач
• Потенциальные приложения: машинное обучение с ограничениями справедливости, ограничения разреженности
• Реализация ПО: команда авторов имеет опыт разработки пакета CDOpt, возможен выпуск инструментария

3. Воспроизводимость

• Преимущества:
  • Описание алгоритма ясно (уравнение 3.20)
  • Экспериментальная установка подробна
  • Проективные отображения имеют явные формулы (Таблица 3)
• Недостатки:
  • Код не опубликован
  • Некоторые детали реализации (параметры немонотонного поиска линии) не даны

4. Направления дальнейших исследований

• Краткосрочные:
  • Ослабление теоретических предположений
  • Расширение на неравенства
  • Больше тестов приложений
• Долгосрочные:
  • Развитие общей теории "частичной оболочки"
  • Комбинирование с другими техниками оптимизации (ADMM, проксимальные методы)
  • Распределённые/стохастические версии

Сценарии применения

1. Идеальные сценарии

• Структура ограничений:
  • $\mathcal{X}$ — простое выпуклое множество (проекция легко вычисляется)
  • $c(x) = 0$ — гладкое равенство
  • Выполняется условие невырожденности ограничений
• Размер задачи: средний ($n \sim 10^2$)
• Требуемая точность: средняя ($\varepsilon \sim 10^{-5}$)

2. Конкретные приложения

• Полуопределённое программирование: эксперименты уже проверены
• Оптимизация на многообразиях: например, оптимизация на многообразии Штифеля
• Машинное обучение:
  • Обучение нейронных сетей с равенствами
  • Классификация с ограничениями справедливости
• Обработка сигналов: задачи восстановления с ограничениями норм

3. Неподходящие сценарии

• Доминирующие неравенства: FBSE обрабатывает только равенства
• Сложная проекция на $\mathcal{X}$: если $\Pi_{\mathcal{X}}$ вычислительно дорога
• Требуемая очень высокая точность: сложность $O(\varepsilon^{-2})$ может быть недостаточна
• Сверхбольшие задачи: проекция и вычисление градиента могут стать узким местом

Избранные ссылки

1. Stella et al. (2017): Forward–backward quasi-newton methods for nonsmooth optimization problems. Computational Optimization and Applications
   • Расширение прямо-обратной оболочки квазиньютоновскими методами
2. Xiao et al. (2023): Dissolving constraints for Riemannian optimization. Mathematics of Operations Research
   • Теоретическая основа методов растворения ограничений
3. Xiao et al. (2025): An exact penalty approach for equality constrained optimization over a convex set. arXiv preprint
   • Предшествующая работа, предложившая отображение растворения ограничений
4. Absil et al. (2008): Optimization algorithms on matrix manifolds. Princeton University Press
   • Классический учебник по оптимизации на многообразиях
5. Rockafellar & Wets (2009): Variational analysis. Springer
   • Теоретическая основа вариационного анализа, используется для анализа проекций и нормальных конусов

Общая оценка: это отличная статья с строгой теорией и инновационным методом. Идея "частичной оболочки" предлагает новый взгляд на решение задач оптимизации со смешанными ограничениями, теоретический анализ полон, численные эксперименты предварительно верифицируют эффективность метода. Основные недостатки — практичность теоретических констант, полнота экспериментов и проверка масштабируемости. Данная работа вносит важный вклад в область невыпуклой оптимизации с ограничениями и имеет высокую академическую ценность и прикладной потенциал. Рекомендуется дальнейшая работа по ослаблению теоретических предположений, расширению тестов приложений и обработке больших задач.