//! Debug: print electron repulsion integrals (ERIs) for a single carbon
//! atom in STO-3G.
//!
//! Category: debug
use sci_form::hf::basis::build_sto3g_basis;
use sci_form::hf::integrals::{compute_eris, get_eri};
fn main() {
let elements = [6];
let positions = [[0.0, 0.0, 0.0]];
let basis = build_sto3g_basis(&elements, &positions);
let eris = compute_eris(&basis);
let n = basis.n_basis();
// px is 2, py is 3
println!("ERI[2,3,2,3] : {}", get_eri(&eris, 2, 3, 2, 3, n));
println!("ERI[2,3,3,2] : {}", get_eri(&eris, 2, 3, 3, 2, n));
}