def scale_nodes(L,N,delta=0.5,s=1.,t=1.): from spyctral.mapjpoly.quad import gq from spyctral.common.scaling import scale_factor [x,w] = gq(N,s=s,t=t) return scale_factor(L,x,delta=delta)
def scale_nodes(L,N,delta=0.5,mu=0.): from spyctral.hermite.quad import pgq from spyctral.common.scaling import scale_factor [x,w] = pgq(N,mu=0.) return scale_factor(L,x,delta=delta)
def scale_nodes(self,L,physical_scale_ratio): """ Sets the affine scaling factor self.scale so that (physical_scale_ratio x N) of the canonical nodes lie inside [-L,L]. """ from spyctral.common.scaling import scale_factor if self.canonical_quadrature is None: print "Error: cannot scale...I don't have a canonical set of nodes" else: x = self.canonical_quadrature().nodes scale = scale_factor(L,x,delta=physical_scale_ratio) self.rehash_parameters(scale=scale) self.initialize_quadrature(None,None) self.make_nodal_differentiation_matrix()