Python NumPy: 4D 계수 배열로 3D 라게르(Laguerre) 급수 평가하기

Python의 NumPy에서 polynomial.laguerre.lagval3d() 메서드를 사용하면 점 (x, y, z)에서 3차원 라게르(Laguerre) 급수를 평가할 수 있습니다. 이 메서드는 x, y, z에서 대응하는 값들의 튜플(triple)로 구성된 점들에 대해 다차원 다항식의 값을 반환합니다.

lagval3d() 메서드의 주요 특징

계수 배열 c가 3차원 미만이면, 자동으로 형태(shape)에 1을 추가하여 3차원으로 만듭니다. 결과의 형태는 c.shape[3:] + x.shape가 됩니다.

첫 번째 매개변수: x, y, z

3차원 급수는 점 (x, y, z)에서 평가되며, 이때 x, y, z는 반드시 동일한 형태를 가져야 합니다. 만약 x, y, z 중 하나라도 리스트(list)나 튜플(tuple)이라면 먼저 ndarray로 변환되고, 그렇지 않으면 그대로 유지됩니다. ndarray가 아닌 경우에는 스칼라(scalar)로 취급됩니다.

두 번째 매개변수: C (계수 배열)

C는 다중 차수(multidegree) i, j, k에 해당하는 항의 계수가 c[i, j, k]에 포함되도록 정렬된 계수 배열입니다. c의 차원이 3보다 크면 나머지 인덱스들은 여러 개의 계수 집합을 열거합니다.

구현 단계

먼저 필요한 라이브러리를 임포트합니다.

import numpy as np
from numpy.polynomial import laguerre as L

4차원 계수 배열을 생성합니다.

c = np.arange(48).reshape(2,2,6,2)

배열을 출력합니다.

print("Our Array...\n",c)

배열의 차원을 확인합니다.

print("\nDimensions of our Array...\n",c.ndim)

배열의 데이터 타입을 확인합니다.

print("\nDatatype of our Array object...\n",c.dtype)

배열의 형태(shape)를 확인합니다.

print("\nShape of our Array object...\n",c.shape)

점 (x, y, z)에서 3D 라게르 급수를 평가하려면 polynomial.laguerre.lagval3d() 메서드를 사용합니다.

print("\nResult...\n",L.lagval3d([1,2],[1,2],[1,2], c))

전체 예제 코드

import numpy as np
from numpy.polynomial import laguerre as L

# 4차원 계수 배열 생성
c = np.arange(48).reshape(2,2,6,2)

# 배열 출력
print("Our Array...\n",c)

# 차원 확인
print("\nDimensions of our Array...\n",c.ndim)

# 데이터 타입 확인
print("\nDatatype of our Array object...\n",c.dtype)

# 형태 확인
print("\nShape of our Array object...\n",c.shape)

# lagval3d() 메서드로 3D 라게르 급수 평가
print("\nResult...\n",L.lagval3d([1,2],[1,2],[1,2], c))

실행 결과

Our Array...
    [[[[ 0  1]
    [  2  3]
    [  4  5]
    [  6  7]
    [  8  9]
    [10 11]]

    [[12 13]
    [14 15]
    [16 17]
    [18 19]
    [20 21]
    [22 23]]]


    [[[24 25]
    [26 27]
    [28 29]
    [30 31]
    [32 33]
    [34 35]]

    [[36 37]
    [38 39]
    [40 41]
    [42 43]
    [44 45]
    [46 47]]]]

Dimensions of our Array...
4

Datatype of our Array object...
int64

Shape of our Array object...
(2, 2, 6, 2)

Result...
    [[-15.66666667  0.  ]
    [-16.925       0.  ]]

정리

이 예제에서는 np.arange(48).reshape(2,2,6,2)를 통해 형태가 (2, 2, 6, 2)인 4차원 int64 배열을 생성했습니다. 앞의 세 차원은 다중 차수(multidegree) i, j, k에 대한 계수를 나타내고, 마지막 네 번째 차원은 여러 개의 계수 집합을 열거합니다. lagval3d() 메서드에 점 [1, 2]를 각 축에 전달하면 각 계수 집합에 대해 두 점에서의 라게르 급수 평가 결과가 2×2 배열로 반환됩니다.