Evaluate a 3-D Laguerre series on the Cartesian product of x, y and z with 4d array of coefficient in Python

To evaluate a 3-D Laguerre series on the Cartesian product of x, y and z, use the polynomial.laguerre.laggrid3d() method in Python. The method returns the values of the three-dimensional Laguerre series at points in the Cartesian product of x, y and z.

If the coefficient array c has fewer than three dimensions, ones are implicitly appended to its shape to make it 3-D. The shape of the result will be c.shape[3:] + x.shape + y.shape + z.shape.

Syntax

numpy.polynomial.laguerre.laggrid3d(x, y, z, c)

Parameters

x, y, z ? The three-dimensional series is evaluated at the points in the Cartesian product of x, y, and z. If x, y, or z is a list or tuple, it is first converted to an ndarray, otherwise it is left unchanged. If it isn't an ndarray, it is treated as a scalar.

c ? An array of coefficients ordered so that the coefficients for terms of degree i,j are contained in c[i,j]. If c has dimension greater than two, the remaining indices enumerate multiple sets of coefficients.

Example

Let's create a 4D array of coefficients and evaluate the 3-D Laguerre series ?

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

# Create a 4d array of coefficients
c = np.arange(48).reshape(2,2,6,2)

# Display the array
print("Our Array...")
print(c)

# Check the Dimensions
print("\nDimensions of our Array...")
print(c.ndim)

# Get the Datatype
print("\nDatatype of our Array object...")
print(c.dtype)

# Get the Shape
print("\nShape of our Array object...")
print(c.shape)

# Evaluate 3-D Laguerre series on the Cartesian product
print("\nResult...")
print(L.laggrid3d([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   2.        ]
   [ 15.1          3.2       ]]

  [[ 30.2          6.4       ]
   [  0.           0.        ]]]


 [[[-16.925        1.73333333]
   [ 15.1          3.2       ]]

  [[ 30.2          6.4       ]
   [  0.           0.        ]]]]

How It Works

The laggrid3d() method evaluates the 3-D Laguerre series by computing values at each point in the Cartesian product of the input arrays. The coefficient array structure determines how the polynomial terms are combined, with the first three dimensions representing the polynomial degrees and any additional dimensions creating multiple result sets.

Conclusion

The laggrid3d() method efficiently evaluates 3-D Laguerre series on Cartesian products. Use 4D coefficient arrays when you need to evaluate multiple polynomial sets simultaneously across the same coordinate grid.