_integrate_part#

sunkit_spex.legacy.emission._integrate_part( *, model, photon_energies, maxfcn, rerr, eelow, eebrk, eehigh, p, q, z, a_lg, b_lg, ll, efd, integrator=None, )[source]#

Perform numerical Gaussian-Legendre Quadrature integration for thick- and thin-target models.

This integration is intended to be performed over continuous portions of the electron distribution.

Parameters:

model (str) – Either thick-target or thin-target
maxfcn (int) – Maximum number of points used in Gaussian quadrature integration
rerr (float) – Desired relative error for integral evaluation. For example, rerr = 0.01 indicates that the estimate of the integral is to be correct to one digit, whereas rerr = 0.001 calls for two digits to be correct.
photon_energies (numpp.array) – Photon energies
eelow (float) – Low energy electron cut off
eebrk (float) – Break energy
eehigh (float) – High energy cutoff
p (float) – Slope below the break energy
q (float) – Slope above the break energy
z (float) – Mean atomic number of plasma
a_lg (numpy.array) – Logarithm of lower integration limits
b_lg (numpy.array) – Logarithm of upper integration limit
ll (numpy.array) – Indices for which to carry out integration
efd (boolean) – True (default) electron flux density distribution, False electron density distribution. This input is not used in the main routine, but is passed to thin_target_integrand

Returns:

Array of integrated photon fluxes evaluation and array of integration status (0 converged, 1 not converged)

Return type:

tuple

References

See SSW Brm2_DmlinO_int.pro and brm2_dmlin.pro.