Esempio n. 1
0
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)
Esempio n. 2
0
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)
Esempio n. 3
0
    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()