def testChebcoef(self):
S = PiecewiseFunction(fun=lambda x: sin(4*(x-0.5)), breakPoints=[-1, 1]).toInterpolated()
seg = S.segments[0]
Xs, Ys = seg.f.Xs, seg.f.Ys
print("Xs=", Xs)
print("Ys=", Ys)
print("Cs=", chebt2(Ys))
print("Xs=", ichebt2(chebt2(Ys)))
#figure()
#S.plot(color="r")
D = S.diff()
#D.plot(color="k")
#show()
self.assertTrue(0 < 1)
def testChebcoef(self):
S = PiecewiseFunction(fun=lambda x: sin(4 * (x - 0.5)),
breakPoints=[-1, 1]).toInterpolated()
seg = S.segments[0]
Xs, Ys = seg.f.Xs, seg.f.Ys
print("Xs=", Xs)
print("Ys=", Ys)
print("Cs=", chebt2(Ys))
print("Xs=", ichebt2(chebt2(Ys)))
#figure()
#S.plot(color="r")
D = S.diff()
#D.plot(color="k")
#show()
self.assertTrue(0 < 1)
if __name__ == "__main__":
from pylab import *
from pacal import *
# B= BetaDistr(1,1) * UniformDistr(0,3)
# B =(UniformDistr(0,3)+UniformDistr(0,1)+UniformDistr(0,1)+UniformDistr(0,1)) * (UniformDistr(0,1)+UniformDistr(0,1)+UniformDistr(0,1))
# B = BetaDistr(4,4) * (UniformDistr(-1,1)+UniformDistr(-1,2))
# B.summary(show_moments=True)
# print B.get_piecewise_pdf()
from pacal.segments import PiecewiseFunction
B = PiecewiseFunction(fun=lambda x:sin(3*x), breakPoints=[-1,0,1])
B = B.toInterpolated()
print(B.segments[0].f.__class__)
#B = B.trimInterpolators(abstol=1e-15)
print(B.segments[0].f.Ys, B.segments[0].f.__class__)
D = B.diff()
D2 = D.diff()
D3 = D2.diff()
D4 = D3.diff()
D5 = D4.diff()
print(D.segments[0].f.Ys, D.segments[0].f.__class__)
print(D2.segments[0].f.Ys)
print(D.roots())
figure()
B.plot()
D.plot()
D2.plot()
D3.plot()
D4.plot()
D5.plot()
show()
if __name__ == "__main__":
from pylab import *
from pacal import *
# B= BetaDistr(1,1) * UniformDistr(0,3)
# B =(UniformDistr(0,3)+UniformDistr(0,1)+UniformDistr(0,1)+UniformDistr(0,1)) * (UniformDistr(0,1)+UniformDistr(0,1)+UniformDistr(0,1))
# B = BetaDistr(4,4) * (UniformDistr(-1,1)+UniformDistr(-1,2))
# B.summary(show_moments=True)
# print B.get_piecewise_pdf()
from pacal.segments import PiecewiseFunction
B = PiecewiseFunction(fun=lambda x:sin(3*x), breakPoints=[-1,0,1])
B = B.toInterpolated()
print B.segments[0].f.__class__
#B = B.trimInterpolators(abstol=1e-15)
print B.segments[0].f.Ys, B.segments[0].f.__class__
D = B.diff()
D2 = D.diff()
D3 = D2.diff()
D4 = D3.diff()
D5 = D4.diff()
print D.segments[0].f.Ys, D.segments[0].f.__class__
print D2.segments[0].f.Ys
print D.roots()
figure()
B.plot()
D.plot()
D2.plot()
D3.plot()
D4.plot()
D5.plot()
show()
