Faster gradients

Improved performance in geometry optimizations and numerical frequency runs due to a new implementation of the calculation of the gradients that now use linear scaling techniques. Especially important for large molecules.