def function(self,
x,
y,
n_sersic,
R_sersic,
k_eff,
center_x=0,
center_y=0):
"""
:param x: x-coordinate
:param y: y-coordinate
:param n_sersic: Sersic index
:param R_sersic: half light radius
:param k_eff: convergence at half light radius
:param center_x: x-center
:param center_y: y-center
:return:
"""
n = n_sersic
x_red = self._x_reduced(x, y, n, R_sersic, center_x, center_y)
b = self.b_n(n)
#hyper2f2_b = util.hyper2F2_array(2*n, 2*n, 1+2*n, 1+2*n, -b)
hyper2f2_bx = util.hyper2F2_array(2 * n, 2 * n, 1 + 2 * n, 1 + 2 * n,
-b * x_red)
f_eff = np.exp(b) * R_sersic**2 / 2. * k_eff # * hyper2f2_b
f_ = f_eff * x_red**(2 * n) * hyper2f2_bx # / hyper2f2_b
return f_
def function(self, x, y, n_sersic, r_eff, k_eff, center_x=0, center_y=0):
"""
returns Gaussian
"""
n = n_sersic
x_red = self._x_reduced(x, y, n, r_eff, center_x, center_y)
b = self.b_n(n)
#hyper2f2_b = util.hyper2F2_array(2*n, 2*n, 1+2*n, 1+2*n, -b)
hyper2f2_bx = util.hyper2F2_array(2*n, 2*n, 1+2*n, 1+2*n, -b*x_red)
f_eff = np.exp(b)*r_eff**2/2.*k_eff# * hyper2f2_b
f_ = f_eff * x_red**(2*n) * hyper2f2_bx# / hyper2f2_b
return f_