[논문 리뷰] TaylUR, an arbitrary-order automatic differentiation package for Fortran 95

Peter H. von Hippel|arXiv (Cornell University)|2005. 06. 29.

Numerical Methods and Algorithms인용 수 5

한 줄 요약

TaylUR는 복소수 함수에 대해 혼합 도함수를 제외한 임의의 순서의 자동 미분을 가능하게 하는 Fortran 95 패키지로, 산술 연산자와 내장 함수를 오버로드하여 테일러 급수 객체 위에서 작동하도록 구현한다. 이는 객체 지향 프로그래밍과 연산자 오버로드를 활용하여 라이프니츠의 법칙과 파우 디 브루노의 공식을 이용해 도함수를 효율적으로 계산한다.

ABSTRACT

We present TaylUR, a Fortran 95 module to automatically compute the numerical values of a complex-valued function's derivatives w.r.t. several variables up to an arbitrary order in each variable, but excluding mixed derivatives. Arithmetic operators and Fortran intrinsics are overloaded to act correctly on objects of defined type taylor, which encodes a function along with its first few derivatives w.r.t. the user-defined independent variables. Derivatives of products and composite functions are computed using Leibniz's rule and Faa di Bruno's formula. TaylUR makes heavy use of operator overloading and other object-oriented Fortran 95 features.

연구 동기 및 목표

• Fortran 95에서 복소수 함수의 고차 도함수를 계산하기 위한 유연하고 확장 가능한 프레임워크를 제공하는 것.
• 혼합 도함수를 제외하고 각 독립 변수에 대해 임의의 순서의 도함수를 지원하는 것.
• 산술 연산자와 Fortran 내장 함수의 오버로드를 통해 원활한 도함수 계산을 가능하게 하는 것.
• 라이프니츠의 공식과 파우 디 브루노의 공식과 같은 수학적 규칙을 이용해 복합 함수와 곱셈 함수의 도함수 계산을 구현하는 것.
• 현대적 Fortran 95 객체 지향 기능을 활용하여 함수와 도함수 정보를 타입 안정적인 방식으로 캡슐화하는 것.

제안 방법

• 각 독립 변수에 대해 지정된 순서까지의 함수와 도함수를 표현하기 위해 'taylor'라는 유도형을 정의한다.
• 산술 연산자(+, -, *, / 등)와 Fortran 내장 함수를 'taylor' 객체 위에서 작동하도록 오버로드하여 도함수를 정확하게 전파한다.
• 곱셈 연산에서 함수의 도함수를 유지하기 위해 라이프니츠의 법칙을 적용한다.
• 복합 함수의 고차 도함수를 계산하기 위해 파우 디 브루노의 공식을 사용하여 올바른 테일러 급수 전파를 보장한다.
• 타입 바운드 프로시저와 연산자 오버로드를 구현하여 표준 Fortran 표현식 내에서 투명하게 도함수 계산을 지원한다.
• 혼합 도함수를 제외하고 순수한 도함수 계산에 국한하여 복잡도와 메모리 오버헤드를 줄인다.

실험 결과

연구 질문

• RQ1 혼합 도함수를 제외하고 Fortran 95에서 복소수 함수의 임의의 순서 도함수를 어떻게 효율적으로 계산할 수 있는가?
• RQ2 Fortran 95의 연산자 오버로드와 객체 지향 기능을 얼마나 효과적으로 활용하여 견고한 자동 미분 시스템을 구현할 수 있는가?
• RQ3 라이프니츠의 공식과 파우 디 브루노의 공식을 타입 시스템에 효과적으로 통합하여 고차 도함수 전파를 지원할 수 있는가?
• RQ4 TaylUR의 설계는 컴iles된 언어인 Fortran에서 성능, 정확성, 확장성의 균형을 어떻게 유지하는가?
• RQ5 도함수 모델에서 혼합 도함수를 배제함으로써 발생하는 제한 사항과 상충 관계는 무엇인가?

주요 결과

• TaylUR은 Fortran 95의 타입 안정적이고 객체 지향적인 설계를 통해 각 독립 변수에 대해 임의의 순서까지 도함수를 성공적으로 계산한다.
• 연산자 오버로드의 사용으로 기존 Fortran 코드베이스에 구문적 오버헤드 없이 도함수 계산을 원활하게 통합할 수 있다.
• 파우 디 브루노의 공식을 통해 복합 함수의 도함수를 정확하게 계산하여 중첩 함수 평가에서 정확성을 보장한다.
• 혼합 도함수를 배제함으로써 저장소와 계산 복잡도를 줄여 계산 효율성을 유지한다.
• 현대적 Fortran 95 기능이 강력하고 확장 가능한 자동 미분 시스템을 구축하는 데 충분함을 보여준다.
• 이 방법은 다변수에서 고차 도함수가 필요한 과학 계산 워크로드에 대해 고정밀 도함수 계산을 가능하게 한다.