Add the "Eckart bottleneck potential" to the library
Fix computation of gradient packets in multidimensional case

Note: The computed gradients (basis shape and coefficients) were NOT
wrong. But the option 'as_packet = True' failed for multidimensional
packets.
Fix kinetic overlap matrix computation for multidimensional packets