Contents

TENSOR_EMP compute inner products of tensor product of MD polynomial basis functions

A = TENSOR_EMP(name,N,I,PHI,mom) return rows of N-fold tensor indices
with basis functions in (PHI,I).
Computed by moments substitution and stored in name.
REQUIRE: deg(P_i)*N=P*N 1D moments per dimension

INDEX_TENSOR Recursive index generation for PCE-tensor

A = INDEXTENSOR(N,K) return sorted tensor indices for K elements of order N.
Fast & low memory storag property: Reduce total number of K^N permutations
(==permutations of elements 1,..,K w/o repetition of elements)
using ortogonality.