Correct for squaring the eigenvalues of the gyration tensor and simplifying the expressions for the shape parameters
- use MathSpecial::square(x) instead of pow(x,2) for improved precision and handling of small and negative numbers - remove unused include statements - no need to refetch the compute in every step. during init() is sufficient
The right file with the cpp code
