def luminosity_integral(self, x, axis_ratio):
r = x * axis_ratio<EOL>return <NUM_LIT:2> * np.pi * r * self.intensities_from_grid_radii(x)<EOL>
Routine to integrate the luminosity of an elliptical light profile.

        The axis ratio is set to 1.0 for computing the luminosity within a circle