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