I = MultiindexSet.createCompleteOrderSet(2, 4)
print I
from spuq.statistics import *
N = NormalDistribution()
U = UniformDistribution()
print N, U
gpc1 = GPCBasis(I, [N, U])
gpc2 = GPCBasis(I, [N, N])
print gpc1.sample(3)
s1 = gpc1.sample(100)
s2 = gpc2.sample(100)
from spuq.utils.plot.plotter import Plotter
Plotter.figure(1)
Plotter.scatterplot(I.arr[:7, :], s1)
Plotter.figure(2)
Plotter.scatterplot(I.arr[8:15, :], s2)
Plotter.figure(3)
Plotter.histplot(s1[:6, :], bins=50)
# multiplication tensor: first as numpy array 3d
# then:
# class MultiplicationTensor
# def contract( vector ) -> (sparse) matrix
# def contract( vector, vector ) -> vector
# storage: full, sparse matrix, own sparse format?
class MultiplicationTensor:
def __init__(self):
