[논문 리뷰] fenics-optim -- Convex optimization interface in FEniCS

Jérémy Bleyer|arXiv (Cornell University)|2020. 01. 10.

Embedded Systems Design Techniques인용 수 6

한 줄 요약

이 논문은 FEniCS에서 볼록 변분 문제를 콘형 프로그래밍과 유한요소 메시화를 활용하여 자동으로 공식화하고 해결하는 파이썬 패키지 fenics-optim을 소개한다. 사용자는 몇 줄의 코드로도 이미지 처리에서 계산 역학에 이르기까지 다양한 분야에서 고정밀도 및 확장 가능한 성능을 보이는 Mosek 내부점수법을 통해 복잡한 PDE 제약 최적화 문제를 효과적으로 표현하고 해결할 수 있다.

ABSTRACT

<strong>fenics-optim</strong> is a Python library that aims at providing a simple user interface in FEniCS for solving convex optimization problems. In particular, it relies on the <strong>Mosek</strong> mathematical programming library. In particular, Mosek is a state-of-the art solver for linear programming (LP), convex quadratic programming (QP), second-order conic programming (SOCP) and semi-definite programming (SDP). The project can be accessed at https://gitlab.enpc.fr/navier-fenics/fenics-optim

연구 동기 및 목표

• 과학 계산 분야에서 비연속 항을 포함한 볼록 변분 문제를 공식화하고 해결하는 데 도전하는 데 목적이 있다. 특히 이미지 처리 및 계산 역학 분야에서 적용된다.
• 최적화 문제의 연속 수학적 공식화를 그대로 반영하는 통합적이고 고수준의 인터페이스를 제공하여 구현 복잡도를 낮추는 것.
• 일관된 프레임워크 내에서 표준 라그랑주, DG, H(div)/H(curl) 요소를 포함한 다양한 요소 유형에 대해 유한요소 메시화를 자동화하는 것.
• 특히 비연속 항이 포함된 문제에 대해 최신 내부점수법을 통해 Mosek를 활용하여 대규모 볼록 문제를 효율적으로 수치적으로 해결하는 것.
• FEniCS의 기호적 PDE 공식화와 현대 콘형 최적화 간 격차를 메우며, 연속 및 비연속 볼록 함수형을 모두 지원하는 것.

제안 방법

• FEniCS의 통합 형식 언어(UFL)를 사용해 볼록 변분 문제를 공식화하고, 비연속 항을 포함한 볼록 함수를 콘형 프로그래밍 프레임워크 내에서 표현한다.
• FEniCS의 유한요소 기능을 활용해 문제를 메시화하며, 혼합 및 혼합-하이브리드 형식을 위한 혼합 유한요소 공간을 지원한다.
• 고정밀도 해를 다항 시간 내에 달성하기 위해 Mosek를 내부점수법 솔버로 활용한다. 특히 복잡한 제약 조건과 비연속 항이 포함된 문제에 매우 효과적이다.
• 변수와 제약 조건을 블록으로 정의하는 블록 구조 최적화 프레임워크를 도입하여 모듈러하고 확장 가능한 문제 구성이 가능하도록 한다.
• 사용자 정의 함수 공간(예: 적분 기반 요소, Quadrature 요소)을 사용해 라그랑주 승수 및 쌍대 변수와 같은 보조 변수를 콘형 공식화 내에서 점별로 표현한다.
• 함수 기반 제약 조건 정의를 통해 등호 및 부등호 제약 조건을 모두 지원하며, UFL과 FEniCS의 형식 컴파일러를 사용해 선형 형식을 자동으로 조립한다.

실험 결과

연구 질문

• RQ1 이미지 처리 및 고체 역학에서 발생하는 볼록 변분 문제를 고수준의 유한요소 프레임워크를 통해 최소한의 사용자 코드로 공식화하고 해결할 수 있는가?
• RQ2 비연속 볼록 항(예: 총 변동량, L1 노름)을 콘형 최적화를 활용해 유한요소 맥락 내에서 효율적으로 표현하고 해결할 수 있는가?
• RQ3 FEniCS의 유한요소 기능과 Mosek의 내부점수법 솔버를 통합함으로써 대규모 볼록 PDE 제약 최적화 문제의 정확도와 효율성이 얼마나 향상되는가?
• RQ4 복잡한 다중장 변분 문제와 혼합 유한요소를 지원하는 모듈러하고 블록 구조 최적화 인터페이스를 설계할 수 있는가?
• RQ5 비연속 목표 함수를 포함한 문제에서, 제안된 프레임워크는 일반적으로 이미지 처리 분야에서 사용되는 1차 방법에 비해 성능과 안정성 면에서 얼마나 뛰어나게 되는가?

주요 결과

• fenics-optim 패키지는 복잡한 볼록 변분 문제를 파이썬 코드 몇 줄로 공식화할 수 있으며, 연속적인 수학적 공식화와 매우 유사하게 표현된다.
• 이 프레임워크는 혼합 및 하이브리드 요소를 포함한 광범위한 유한요소를 지원하며, 비연속 항을 포함한 문제의 자동 메시화를 가능하게 한다.
• Mosek의 내부점수법를 활용해 다항 시간 내에 고정밀도 해를 달성하며, 해의 정밀도 면에서 1차 방법보다 뛰어난 성능을 보인다.
• 비연속 항의 콘형 재구성(예: 장애물 문제에 대한 회전된 2차 콘)을 통해 현대 콘형 솔버를 통한 효율적이고 안정적인 수치적 해법이 가능해진다.
• 블록 구조 최적화 제약 조건 정의 방식을 통해 모듈러하고 확장 가능한 문제 구성이 가능하며, FEniCS의 형식 컴파일러를 통한 희박성 및 행렬 조립의 자동 처리가 이루어진다.
• 적분 기반 변수와 콘형 제약 조건을 사용해 장애물 문제의 명시적 구성은 약한 형태에서의 점별 제약 조건을 효과적으로 처리할 수 있음을 보여주며, 정확하고 안정적인 해를 달성한다.