[논문 리뷰] The Dune Python Module

Andreas Dedner, Martin Nolte|arXiv (Cornell University)|2018. 07. 13.

Computational Physics and Python Applications참고 문헌 7인용 수 13

한 줄 요약

이 논문은 Dune C++ 유한요소 프레임워크를 위한 Just-in-Time 컴파일링된 Python 인터페이스인 Dune-Python을 소개한다. 이는 Dune의 정적 다형성 기반 그리드 인터페이스에 직접적이고 효율적인 접근을 가능하게 한다. Pybind11와 유형 레지스트리 기반으로, Python 프로토타입을 고성능 C++ 코드로 원활하게 변환하면서도 원래 Dune 생태계의 완전한 유연성과 모듈성을 유지한다.

ABSTRACT

In this paper we present the new Dune-Python module which provides Python bindings for the Dune core, which is a C++ environment for solving partial differential equations. The aim of this new module is to firstly provide the general infrastructure for exporting realizations of statically polymorphic interfaces based on just-in-time compilation and secondly to provide bindings for the central interfaces of the dune core modules. In the first release we focus on the grid interface. Our aim is to only introduce a thin layer when passing objects into Python which can be removed when the object is passed back into a C++ algorithm. Thus no efficiency is lost and little additional code maintenance cost is incurred. To make the transition for Dune users to the Python environment straightforward the Python classes provide a very similar interface to their C++ counterparts. In addition, vectorized versions of many interfaces allow for more efficient code on the Python side. The infrastructure for exporting these interfaces and the resulting bindings for a Dune grid are explained in detail in this paper for both experienced Dune users and others interested in a flexible Python environment for implementing grid based schemes for solving partial differential equations.

연구 동기 및 목표

• 성능이나 모듈성을 희생시키지 않고 Dune C++ 코어에 얇고 효율적인 Python 인터페이스를 제공하기 위해.
• 유한요소 스킴을 Python에서 원활하게 프로토타이핑하고, 이를 직접 고성능 C++ 코드로 번역할 수 있도록 하기 위해.
• 기존 접근 방식이 동적 다형성 또는 유형 추상화된 인터페이스를 사용해 성능와 유연성을 저해하는 문제를 해결하기 위해.
• 정적 다형성 기반 C++ 템플릿을 직접 바인딩하여 Dune의 그리드 기반 수치 해법 전반을 지원하기 위해.
• Dune의 핵심 그리드 인터페이스를 Python에서 노출시켜, 유한요소 방법에 대한 교육과 실험의 장벽을 낮추기 위해.

제안 방법

• Pybind11를 사용하여 Dune의 C++ 인터페이스, 특히 그리드 모듈에 대한 Python 바인딩을 생성하기 위해.
• Dune의 정적 다형성 설계에서 발생하는 수많은 가능한 타입을 처리하기 위해, C++ 템플릿 인스턴시에이션을 실시간으로 Just-in-Time(JIT) 컴파일하기 위해.
• 파이썬 클래스 인스턴스를 해당하는 C++ 템플릿 타입과 필요한 포함 파일로 매핑하는 유형 레지스트리 구현하기 위해.
• C++ 인터페이스를 최소한으로 반영하는 파이썬 래퍼 레이어를 설계하여, 파이썬과 C++ 간의 직접적인 코드 이행을 가능하게 하기 위해.
• 인터언어 호출 오버헤드를 줄이고 파이썬 측 성능을 향상시키기 위해 벡터화된 인터페이스 도입하기 위해.
• C++ 함수 템플릿 컴파일을 지원하도록 바인딩 인fra스트럭처 확장하기 위해, 하이브리드 C++/파이썬 워크플로우 가능하게 하기 위해.

실험 결과

연구 질문

• RQ1 Dune 프레임워크의 정적 다형성 기반 C++ 인터페이스를 런타임 성능 저하 없이 파이썬에 효율적으로 노출할 수 있는가?
• RQ2 JIT 컴파일링 접근 방식을 통해 Dune의 그리드 인터페이스를 파이썬에서 완전히 접근하면서도 원래 C++ 코드의 모듈성과 성능을 유지할 수 있는가?
• RQ3 파이썬에서 작성한 유한요소 스킴의 프로토타입이 이 바인딩 레이어를 통해 최적화된 C++ 코드로 직접 번역될 수 있는 정도는 어느 정도인가?
• RQ4 Dune 기반 수치 해법의 파이썬 구현에서 인터언어 호출 오버헤드를 어떻게 최소화할 수 있는가?
• RQ5 이 바인딩 시스템이 복잡한 PDE 솔버의 교육 및 빠른 프로토타이핑에 실현 가능할까?

주요 결과

• Dune-Python 모듈은 C++ 템플릿 인스턴시에이션의 JIT 컴파일링을 통해 파이썬에서 Dune의 그리드 인터페이스에 직접적이고 효율적인 접근을 성공적으로 가능하게 하였다.
• 바인딩 레이어는 최소한의 런타임 오버헤드를 유발하며, Dune의 정적 다형성 설계와 완전히 호환되어 성능과 모듈성을 유지한다.
• Dune-Python 인터페이스를 사용해 작성한 파이썬 코드는 C++ 대응 코드와 매우 유사하여 고성능 C++ 구현으로의 직접적인 번역이 가능하다.
• 벡터화된 인터페이스는 인터언어 호출 횟수를 크게 줄여 소형에서 중형 문제에 대해 성능을 향상시켰다.
• 유형 레지스트리 메커니즘은 파이썬 클래스 인스턴스를 해당하는 C++ 템플릿 타입과 포함 파일로 신뢰성 있고 확장 가능한 방식으로 매핑할 수 있게 하였다.
• 이 접근 방식은 하이브리드 개발 워크플로우를 지원하여, 성능에 민감한 컴포넌트는 C++로, 고수준 로직은 파이썬으로 유지할 수 있게 하였다.