nld_T_product

ompy.nld_T_product(double[::1] nld, double[::1] T, double[::1] resolution, double[::1] E_nld, double[::1] Eg, double[::1] Ex)

Computes first generation matrix from nld and gSF

Uses the equation

P(Ex, Eγ) ∝ ρ(Ex-Eγ)·T(Eγ)

The energy resolution is exploited to only perform the product up to the diagonal, and employs multithreading along each row.

Parameters:
  • nld – Nuclear level density

  • T – Transmission coefficients

  • resolution – The resolution at each bin of Eg

  • E_nld – The energy array for the nuclear level density

  • Eg – The gamma energy

  • Ex – The excitation energy

Returns:

The first generation matrix

Todo

  • Is the indexing optimal?