Add ODE_DEFAULT_NORM overload for GTPSA

[mattsignorelli]

Checklist

• Appropriate tests were added
• Any code changes were done in a way that does not break public API
• All documentation related to code changes were updated
• The new code follows the
  contributor guidelines, in particular the SciML Style Guide and
  COLPRAC.
• Any new documentation only uses public API

Additional context

I added an overload for ODE_DEFAULT_NORM with TPS types where it calculates the 2-norm of the scalar (0th order) part of the TPS, and is necessary because the GTPSA C library does not allow sqrt(TPS(0)) as the derivatives are undefined.

The calculation of the norm of the scalar part follows from what I see in the Tracker extension. If instead it should calculate a norm including all monomial coefficients/partial derivatives, I can quickly change that to the function GTPSA.normTPS.

[ChrisRackauckas]

If instead it should calculate a norm including all monomial coefficients/partial derivatives, I can quickly change that to the function GTPSA.normTPS

That would be better because then it would be guaranteed convergent.

[mattsignorelli]

Done.