Running HF with anisotropic contracted basis set

I don’t believe this is possible without rewriting much of our integral and AO code.