pyopensn.aquad.GLCTriangularQuadrature3DXYZ

class pyopensn.aquad.GLCTriangularQuadrature3DXYZ

Triangular Gauss-Legendre-Chebyshev quadrature for 3D, XYZ geometry.

For each polar level away from the equator, there is 1 less azimuthal angle per octant. The maximum number of azimuthal angles (at the equator) is automatically computed as 2 * n_polar.

Wrapper of opensn::GLCTriangularQuadrature3DXYZ.

GetDiscreteToMomentOperator(self: pyopensn.aquad.AngularQuadrature) → numpy.ndarray[numpy.float64]

Get the discrete-to-moment operator as a numpy array.

GetMomentToDiscreteOperator(self: pyopensn.aquad.AngularQuadrature) → numpy.ndarray[numpy.float64]

Get the moment-to-discrete operator as a numpy array.

GetMomentToHarmonicsIndexMap(self: pyopensn.aquad.AngularQuadrature) → list[pyopensn.aquad.HarmonicIndices]

__init__(self: pyopensn.aquad.GLCTriangularQuadrature3DXYZ, **kwargs) → None

Construct a Triangular Gauss-Legendre-Chebyshev quadrature for 3D, XYZ geometry.

Parameters:

• n_polar (int) – Number of polar angles. The maximum number of azimuthal angles (at the equator) is automatically computed as 2 * n_polar.

• scattering_order (int) – Maximum scattering order supported by the angular quadrature. Optional when operator_method='galerkin_one', in which case the scattering order is automatically determined so that the number of moments equals the number of angles.

• operator_method ({'standard', 'galerkin_one', 'galerkin_three'}, default='standard') – Method used to construct the discrete-to-moment and moment-to-discrete operators.

• verbose (bool, default=False) – Verbosity.

property abscissae

Vector of polar and azimuthal angles.

property omegas

Vector of direction vectors.

property weights

Quadrature weights.