def plot(self):
l = self.tree.leaves()
for leaf in l:
print leaf.boundary(), leaf.error()
plotfunc(self.f, *self.tree.boundary())
plottree(self.tree)
plt.show()
TreeGrapher([self.tree], "tree.pdf", unique=False, debug=False).graph()
def plot(self):
l = self.tree_hp.leaves()
dof = 0
for leaf in l:
dof = dof + leaf._p
print leaf._p, leaf.boundary(), leaf.error()
print dof, self.tree_hp.sum_of_leaves()
plotfunc(self.f, *self.tree.boundary())
plottree(self.tree_hp)
plt.show()
TreeGrapher([self.tree_hp], "tree_hp.pdf", unique=False, debug=True).graph()
TreeGrapher([self.tree], "tree.pdf", unique=False, debug=True).graph()
from error_functionals import TwoNormOrth as no from error_functionals import TwoNorm as n from plotter import plotfunc import matplotlib.pyplot as plt import numpy as np f = lambda x: x ** 3 node = [0, 10] e = n(f) e2 = no(f) pol, coeffs = e.bestpoly(node, d=17) pol2, coeffs2 = e2.bestpoly(node, d=17) plotfunc(lambda x: pol(x) - pol2(x), *node) plt.show()
