def make_poly_regressors(n_samples, order=2):
# mean
X = np.ones((n_samples, 1))
for d in range(order):
poly = Legendre.basis(d + 1)
poly_trend = poly(np.linspace(-1, 1, n_samples))
X = np.hstack((X, poly_trend[:, None]))
return X
def generate_gauss_lobatto_nodes(self, num_points):
nodes = Legendre.basis(num_points - 1).deriv().roots()
return np.concatenate(([-1], nodes, [1]))
def generate_gauss_nodes(self, num_points):
return Legendre.basis(num_points).roots()
