Four-index integral transformation exploiting symmetry

Method of solution: The sparsity of molecular integrals is exploited. The program treats only nonzero integrals. The length of running indices in DO loops is reduced using the block-diagonal form of the MO coefficient matrix. In the present program, the point group is limited to D2h and its subgroups.