Automatic implementation of finite strain anisotropic hyperelastic models using hyper-dual numbers