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()