Self-contained Cartesian-GTO evaluator for OQP basis sets.
Builds AO values on arbitrary points from the basis exposed by
``mol.data.get_basis()``, using OQP's normalization/ordering convention
(verified by reproducing OQP's overlap matrix ``OQP::SM`` to ~1e-6 with an
analytic McMurchie-Davidson overlap, and on a grid to <=1e-3).
Convention (from the OQP molden writer):
* shells store OQP-internal contraction coefficients; the molden/standard
normalized-primitive coefficient is ``d = coef * NORMS[L] * (2a)^-(L/2+3/4)``.
* Cartesian component order: L=0 [s]; L=1 [x,y,z]; L=2 [xx,yy,zz,xy,xz,yz];
L=3 [xxx,yyy,zzz,xyy,xxy,xxz,xzz,yzz,yyz,xyz].