Improved calculation of hyperdual sqrt

2 jobs for feature/localhyperdualstiffness in 4 minutes and 19 seconds (queued for 1 second)