//! Debug: compute and print overlap and kinetic integrals for a single
//! carbon atom in STO-3G basis for manual verification.
//!
//! Category: debug
use sci_form::hf::basis::build_sto3g_basis;
use sci_form::hf::overlap_kin::{compute_kinetic_matrix, compute_overlap_matrix};
fn main() {
let elements = [6];
let positions = [[0.0, 0.0, 0.0]];
let basis = build_sto3g_basis(&elements, &positions);
let s = compute_overlap_matrix(&basis);
let t = compute_kinetic_matrix(&basis);
for i in 0..s.nrows() {
println!("S[{},{}] = {}", i, i, s[(i, i)]);
println!("T[{},{}] = {}", i, i, t[(i, i)]);
}
}