def setUp(self):
#linear SAMPLE objects
ct = 4
lim_t = 31
lbpar_t = cd.LbFixParams([ct], [lim_t])
shf_t = 5.0
band_t = 40.0
dt = band_t / ct
diff = [dt]
shf = [shf_t]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
lb_t = cd.LbGen(lbpar_t)
pt_t = cd.PtGen(ptpar)
dim_lb = (0, 0)
obj_lb = [lb_t]
dims_lb = [dim_lb]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [(0, 0)]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genS = gn.Generator(mst_lb, mst_pt)
self.gdS = gd.Grid(genS)
#BOUNDS
diff = [dt]
shf = [shf_t - 5.0]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
pt_t = cd.PtGen(ptpar)
dims_pt = [(0, 0)]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genB = gn.Generator(mst_lb, mst_pt)
self.gdB = gd.Grid(genB)
return
def adapt_scale(self, band_old, n_old, band_new, n_new, did):
dims = (did, did)
X = band_new[1] - band_new[0]
dx = X / n_new
shf = band_new[0]
par_lb = cd.LbFixParams( [n_new], [n_new * 4 - 1] )
par_pt = cd.PtFixParams( [dx], [shf], [X] )
lb_obj = cd.LbGen(par_lb)
pt_obj = cd.PtGen(par_pt)
mst_lb = mst.MdimStruct( [lb_obj], [(0,0)] )
mst_pt = mst.MdimStruct( [pt_obj], [(0,0)] )
gen = gn.Generator(mst_lb, mst_pt)
rgrid = gd.Grid(gen)
self.samps.replace(rgrid, dims)
self.bounds.replace(rgrid, dims)
return
def factory(bandlims, ns, shf = [0.0], gdtype = 'cartes'):
res = None
if 'cartes' == gdtype:
ct = ns[0]
lim_t = ct * 4 - 1
lbpar_t = cd.LbFixParams( [ct], [lim_t] )
shf_t = shf[0]
T = bandlims[0][1]
dt = T / ct
diff = [dt]
shf = [shf_t]
band = [T]
ptpar = cd.PtFixParams(diff, shf, band)
lb_t = cd.LbGen(lbpar_t)
pt_t = cd.PtGen(ptpar)
dim_lb = (0, 0)
obj_lb = [lb_t]
dims_lb = [dim_lb]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
samps = gd.Grid(gen)
shf = [0.0]
ptpar = cd.PtFixParams(diff, shf, band)
pt_t = cd.PtGen(ptpar)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
bounds = gd.Grid(gen)
bandstr = mst.MdimStruct( [ [0, T] ], [ (0,0) ] )
res = Domain(samps, bounds, None, None, True, bandstr)
elif 'polar' == gdtype:
R = bandlims[0][1]
A = bandlims[1][1]
cr = int( ns[0] )
ca = int( ns[1] )
lim_r = cr * 4 - 1
lim_a = ca * 2 - 1
lbpar_r = pl.LbFixParams( [cr], [lim_r] )
lbpar_a = pl.LbFixParams( [ca], [lim_a] )
shf_r = shf[0]
shf_a = shf[1]
dr = R / cr
da = A / ca
diff = [dr, da]
shf = [shf_r, shf_a]
band = [R, A]
ptpar = pl.PtFixParams(diff, shf, band)
pol_lb_r = pl.PolarLbGen(lbpar_r)
pol_lb_a = pl.PolarLbGen(lbpar_a)
pol_pt = pl.PolarPtGen(ptpar)
dim_lb_r = (0, 0)
dim_lb_a = (1, 1)
obj_lb = [pol_lb_r, pol_lb_a]
dims_lb = [dim_lb_r, dim_lb_a]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [ (0, 1) ]
obj_pt = [pol_pt]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
samps = gd.Grid(gen)
shf_r = 0.0
shf_a = 0.0
ptpar = pl.PtFixParams(diff, shf, band)
pol_pt = pl.PolarPtGen(ptpar)
dims_pt = [ (0, 1) ]
obj_pt = [pol_pt]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
bounds = gd.Grid(gen)
volume = vol.Volume( mst.MdimStruct( [vol.sector], [ (0,1) ] ) )
meas = ms.Measure( mst.MdimStruct( [ms.sector], [ (0,1) ] ) )
bandstr = mst.MdimStruct( [ [0, R], [0, A] ], [ (0,0), (1,1) ] )
res = Domain(samps, bounds, volume, meas, True, bandstr)
else:
pass
return res
def setUp(self):
#linear SAMPLE objects
ct = 4
lim_t = 31
lbpar_t = cd.LbFixParams( [ct], [lim_t] )
shf_t = 0.0
band_t = 1.0
dt = band_t / ct
diff = [dt]
shf = [shf_t]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
lb_t = cd.LbGen(lbpar_t)
pt_t = cd.PtGen(ptpar)
dim_lb = (0, 0)
obj_lb = [lb_t]
dims_lb = [dim_lb]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genS = gn.Generator(mst_lb, mst_pt)
gdS = gd.Grid(genS)
#BOUNDS
diff = [dt]
shf = [shf_t]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
pt_t = cd.PtGen(ptpar)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genB = gn.Generator(mst_lb, mst_pt)
gdB = gd.Grid(genB)
bandstr = mst.MdimStruct( [ [0.0, 1.0] ], [ (0,0) ] )
symb_t = mst.MdimStruct( ['t'], [ (0,0) ] )
dm_t = dm.Domain(gdS, gdB, None, None, True, bandstr)
self.table = tb.Table(dm_t, symb_t, None, None, None, None, True)
return
def setUp(self):
cr = 8
ca = 8
lim_r = 31
lim_a = 7
lbpar_r = pl.LbFixParams([cr], [lim_r])
lbpar_a = pl.LbFixParams([ca], [lim_a])
shf_r = 0.0
shf_a = 0.0
band_r = 8.0
band_a = 2.0 * np.pi
dr = band_r / cr
da = band_a / ca
diff = [dr, da]
shf = [shf_r, shf_a]
band = [band_r, band_a]
ptpar = pl.PtFixParams(diff, shf, band)
pol_lb_r = pl.PolarLbGen(lbpar_r)
pol_lb_a = pl.PolarLbGen(lbpar_a)
pol_pt = pl.PolarPtGen(ptpar)
dim_lb_r = (0, 0)
dim_lb_a = (1, 1)
obj_lb = [pol_lb_r, pol_lb_a]
dims_lb = [dim_lb_r, dim_lb_a]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [(0, 1)]
obj_pt = [pol_pt]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
self.gdPol = gd.Grid(gen)
#linear objects
ct = 4
lim_t = 31
lbpar_t = cd.LbFixParams([ct], [lim_t])
shf_t = 0.0
band_t = 4.0
dt = band_t / ct
diff = [dt]
shf = [shf_t]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
lb_t = cd.LbGen(lbpar_t)
pt_t = cd.PtGen(ptpar)
dim_lb = (0, 0)
obj_lb = [lb_t]
dims_lb = [dim_lb]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [(0, 0)]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
gen = gn.Generator(mst_lb, mst_pt)
self.gdA = gd.Grid(gen)
self.gdB = gd.Grid(gen)
return
def setUp(self):
samp_shf = 0.00
#linear SAMPLE objects
ct = 4
lim_t = 31
lbpar_t = cd.LbFixParams( [ct], [lim_t] )
shf_t = samp_shf
band_t = 1.0
dt = band_t / ct
diff = [dt]
shf = [shf_t]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
lb_t = cd.LbGen(lbpar_t)
pt_t = cd.PtGen(ptpar)
dim_lb = (0, 0)
obj_lb = [lb_t]
dims_lb = [dim_lb]
mst_lb = mst.MdimStruct(obj_lb, dims_lb)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genS = gn.Generator(mst_lb, mst_pt)
self.gdS = gd.Grid(genS)
#BOUNDS
diff = [dt]
shf = [shf_t - samp_shf]
band = [band_t]
ptpar = cd.PtFixParams(diff, shf, band)
pt_t = cd.PtGen(ptpar)
dims_pt = [ (0, 0) ]
obj_pt = [pt_t]
mst_pt = mst.MdimStruct(obj_pt, dims_pt)
genB = gn.Generator(mst_lb, mst_pt)
self.gdB = gd.Grid(genB)
bandstr = mst.MdimStruct([ [0.0, 1.0] ], [(0,0)])
dmA = dm.Domain(self.gdS, self.gdB, None, None, True, bandstr)
self.dmP2d = dmA.prod(dmA)
self.dmP3d = self.dmP2d.prod(dmA)
#set up codomains
T = dmA.bandlims.get_obj((0,0))[1]
n = dmA.ns.get_obj((0,0))
dt = T/n
shf = 2
phi = 1.0 + 0.0j
ampl = 1.0
symb = 't'
fcstr = sgl.cyc_stepfc(dt, T, shf, phi, ampl, symb)
dims = (0,0)
codmstr = mst.MdimStruct( [fcstr], [dims] )
codmA = codm.Codomain(codmstr)
shf = 1
fcstr = sgl.cyc_stepfc(dt, T, shf, phi, ampl, symb)
dims = (0,0)
codmstr = mst.MdimStruct( [fcstr], [dims] )
codmB = codm.Codomain(codmstr)
self.codmP2d = codmA.prod(codmB)
return