#end dist = pc.Gamma(2) print pc.orth_bert(2, dist) # [1.0, q0-2.0, q0^2-6.0q0+6.0] #end dist = pc.Uniform(-1, 1) print dist.ttr([0, 1, 2, 3]) # [[ 0. 0. 0. 0. ] # [-0. 0.33333333 0.26666667 0.25714286]] #end dist = pc.Lognormal(0.01) orths = pc.orth_ttr(2, dist) print orths # [1.0, q0-1.00501252086, q0^2-2.04042818514q0+0.842739860094] #end dist = pc.Iid(pc.Gamma(1), 2) orths = pc.orth_ttr(2, dist) print orths # [1.0, q1-1.0, q0-1.0, q1^2-4.0q1+2.0, q0q1-q1-q0+1.0, q0^2-4.0q0+2.0] #end q = pc.variable() dist = pc.Normal() print pc.E(q, dist) # 0.0 #end
dist = pc.Gamma(2) print pc.orth_bert(2, dist) # [1.0, q0-2.0, q0^2-6.0q0+6.0] #end dist = pc.Uniform(-1,1) print dist.ttr([0,1,2,3]) # [[ 0. 0. 0. 0. ] # [-0. 0.33333333 0.26666667 0.25714286]] #end dist = pc.Lognormal(0.01) orths = pc.orth_ttr(2, dist) print orths # [1.0, q0-1.00501252086, q0^2-2.04042818514q0+0.842739860094] #end dist = pc.Iid(pc.Gamma(1), 2) orths = pc.orth_ttr(2, dist) print orths # [1.0, q1-1.0, q0-1.0, q1^2-4.0q1+2.0, q0q1-q1-q0+1.0, q0^2-4.0q0+2.0] #end q = pc.variable() dist = pc.Normal() print pc.E(q, dist) # 0.0 #end
from pylab import * import current as pc dist = pc.Iid(pc.Gamma(1), 3) x, y, z = pc.variable(3) poly = x + y + z + x * y * z print poly # q0q1q2+q2+q1+q0 #end print pc.E_cond(poly, z, dist) # 2.0q2+2.0 #end print pc.E_cond(poly, [x, y], dist) # q0q1+q1+q0+1.0 #end
from numpy import * import current as pc x, y = pc.variable(2) print x # q0 #end polys = pc.Poly([1, x, x * y]) print polys # [1, q0, q0q1] #end print pc.basis(4) # [1, q0, q0^2, q0^3, q0^4] #end print pc.basis(1, 2, dim=2) # [q0, q1, q0^2, q0q1, q1^2] #end print pc.basis(1, [1, 2]) # [q0, q1, q0q1, q1^2, q0q1^2] #end print pc.basis(1, 2, dim=2, sort="GRI") # [q0^2, q0q1, q1^2, q0, q1] #end poly = pc.Poly([1, x**2, x * y]) print poly(2, 3)
from numpy import * import current as pc x,y = pc.variable(2) print x # q0 #end polys = pc.Poly([1, x, x*y]) print polys # [1, q0, q0q1] #end print pc.basis(4) # [1, q0, q0^2, q0^3, q0^4] #end print pc.basis(1, 2, dim=2) # [q0, q1, q0^2, q0q1, q1^2] #end print pc.basis(1, [1, 2]) # [q0, q1, q0q1, q1^2, q0q1^2] #end print pc.basis(1, 2, dim=2, sort="GRI") # [q0^2, q0q1, q1^2, q0, q1] #end poly = pc.Poly([1, x**2, x*y]) print poly(2, 3)