def vmodel_dataset(pars, i, xy_mesh, r_0 = None, r_space = ring_space): parvals = pars.valuesdict() pa = parvals['pa'] inc = parvals['inc'] Vsys = parvals['Vsys'] Vrot = parvals['Vrot_%i'% i] x0,y0 = parvals['x0'],parvals['y0'] modl = (CIRC_MODEL(xy_mesh,Vrot,pa,inc,x0,y0))*weigths_w(xy_mesh,shape,pa,inc,x0,y0,r_0,r_space,pixel_scale=pixel_scale) return modl,Vsys
def evaluate_model(vmode,
xy_mesh,
Vrot_i,
r_0,
r_space=ring_space,
Vrad_i=None,
Vtan_i=None):
if vmode == "circular":
modl = (CIRC_MODEL(xy_mesh, Vrot_i, pa, inc, x0, y0)) * weigths_w(
xy_mesh, shape, pa, inc, x0, y0, r_0, r_space, pixel_scale)
return modl, Vsys
if vmode == "radial":
modl = (RADIAL_MODEL(
xy_mesh, Vrot_i, Vrad_i, pa, inc, x0, y0)) * weigths_w(
xy_mesh, shape, pa, inc, x0, y0, r_0, r_space, pixel_scale)
return modl, Vsys
if vmode == "bisymmetric":
modl = (BISYM_MODEL(xy_mesh, Vrot_i, Vrad_i, pa, inc, x0, y0,
Vtan_i, phi_b)) * weigths_w(
xy_mesh, shape, pa, inc, x0, y0, r_0,
r_space, pixel_scale)
return modl, Vsys
def vmodel_dataset(pars, i, xy_mesh, r_0 = None, r_space = ring_space): parvals = pars.valuesdict() pa = parvals['pa'] inc = parvals['inc'] Vsys = parvals['Vsys'] Vrot = parvals['Vrot_%i'% i] Vrad = parvals['Vrad_%i'% i] Vtan = parvals['Vtan_%i'% i] x0,y0 = parvals['x0'],parvals['y0'] phi_b = parvals['phi_b'] modl = (BISYM_MODEL(xy_mesh,Vrot,Vrad,pa,inc,x0,y0,Vtan,phi_b))*weigths_w(xy_mesh,shape,pa,inc,x0,y0,r_0,r_space,pixel_scale=pixel_scale) return modl,Vsys
def evaluate_model(vmode,xy_mesh, Vrot_i, r_0, r_space = ring_space): modl = Vrot_i*np.ones(len(xy_mesh[0]),)*weigths_w(xy_mesh,shape,pa,inc,x0,y0,r_0,r_space,pixel_scale) return modl,Vsys