Non thermal models clean up and refactor

[samaloney]

PR Description

Less code, more readable code, and some performance improvements

• Remove some unused and unnecessary thin wrappers
• The old gauss_legendre port and the fixed_quad_batch are almost identical which should now be clearer the later is numerically better as the scaling operations are not compounded in the integration as they are done at the end - changed default to fixed_quad_batch
• Figure out source of some of the numeric constants and added comment with units
• Pull out a bunch of static constants to module level vars
• Modernised some of the code structure and eliminated unnecessary checks
• More / improved comments
• More edge case tests
• ThickTarget ~ 1.5x and ThinTarget ~2x faster

• Still importing constants from legacy should we just switch to the astropy values at this stage or also pull out the constant function from legacy

AI Assistance Disclosure

AI tools were used for:

• Code generation (e.g., when writing an implementation or fixing a bug)
• [x ] Test/benchmark generation
• Documentation (including examples)
• [x ] Research and understanding
• No AI tools were used

Regardless of AI use, the human contributor remains fully responsible for correctness, design choices, licensing compatibility, and long-term maintainability.

[samaloney]

It much easier to see that gauss_legendre and gauss_legendre do almost the same thing the only difference is when the mapping to the real limits is done before gauss_legendre or after in gauss_legendre which mean the numerical error aren't folded in to the integration