def hfunc_Abel(n) : """ Provides the parameters for an Abel integration of a Sersic function It returns a set of factors which depends on the input Sersic n """ X = get_package_data("Abel_param.dat", asfile=False) ndata = X[:,0] nudata = X[:,1] kdata = X[:,2] a0data = X[:,3] a1data = X[:,4] a2data = X[:,5] a3data = X[:,6] a0 = interpolate.interp1d(ndata, a0data) a1 = interpolate.interp1d(ndata, a1data) a2 = interpolate.interp1d(ndata, a2data) a3 = interpolate.interp1d(ndata, a3data) k = interpolate.interp1d(ndata, kdata) nu = interpolate.interp1d(ndata, nudata) if (n > 10) | (n < 0.5) : print "Error: n is too high or too low here \n" return 0.,0.,0.,0.,0., 0. return nu(n), k(n), a0(n), a1(n), a2(n), a3(n)
def hfunc_Trujillo(n) : from . import __path__ as path X = get_package_data("Trujillo.dat", asfile=False) ndata = X[:,0] nudata = X[:,1] pdata = X[:,2] h1data = X[:,3] h2data = X[:,4] h3data = X[:,5] h1 = interpolate.interp1d(ndata, h1data) h2 = interpolate.interp1d(ndata, h2data) h3 = interpolate.interp1d(ndata, h3data) p = interpolate.interp1d(ndata, pdata) nu = interpolate.interp1d(ndata, nudata) if (n > 10) | (n < 0.5) : print "Error: n is too high or too low here \n" return 0.,0.,0.,0.,0. return nu(n), p(n), h1(n), h2(n), h3(n)