PythonのNumPyでlaggrid3d()を使い、係数の2次元配列からx・y・zのデカルト積上の3次元ラゲール級数を評価する方法

Pythonでx、y、zのデカルト積（直積）上の点における3次元ラゲール級数を評価するには、numpy.polynomial.laguerre.laggrid3d()メソッドを使用します。このメソッドは、x、y、zのデカルト積に含まれる各点における3次元ラゲール級数の値を返します。

係数配列cが3次元未満の場合は、形状が3次元になるように暗黙的に1が追加されます。結果の形状は c.shape[3:] + x.shape + y.shape + z.shape となります。

パラメータの説明

第1パラメータ（x, y, z）：3次元級数を評価する対象となる点の座標です。x、y、zのデカルト積の各点で評価が実行されます。引数がリストやタプルの場合はまずndarrayに変換され、すでにndarrayであればそのまま使用されます。どちらでもない場合はスカラーとして扱われます。

第2パラメータ（c）：係数を格納した配列です。次数i、jの項に対応する係数が c[i,j] に含まれるように並べられています。cの次元が3より大きい場合、残りのインデックスは複数の係数セットを列挙することになります。

実装の手順

まず、必要なライブラリをインポートします。

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

次に、係数の2次元配列を作成します。

c = np.arange(4).reshape(2,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)

x、y、zのデカルト積上で3次元ラゲール級数を評価するには、polynomial.laguerre.laggrid3d()メソッドを使用します。

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

完全なコード例

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

# 係数の2次元配列を作成
c = np.arange(4).reshape(2,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)

# polynomial.laguerre.laggrid3d()メソッドを使って
# x、y、zのデカルト積上で3次元ラゲール級数を評価
# このメソッドは、デカルト積の各点における3次元ラゲール級数の値を返す

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

実行結果

Our Array...
    [[0 1]
   [2 3]]

Dimensions of our Array...
2

Datatype of our Array object...
int64

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

Result...
    [[ 0. 2.]
   [-1. -1.]]

このように、laggrid3d()メソッドを使えば、複数の座標点の組み合わせに対するラゲール級数の値を一度に効率的に計算できます。科学技術計算や多項式近似の分野で非常に便利な機能です。