4
/* SORTI(): Place all the components of the Lagrangian into a large
5
** matrix, I (moinfo.I), which we also symmetrize by computing Ipq =
6
** 1/2 (Ipq + Iqp). This matrix is later written to disk in dump()
7
** for subsequent backtransformation. Note that some of the
8
** components of the Lagrangian computed into the IIJ, Iij, IIA, and
9
** Iia matrices remain non-symmetric (e.g., IIJ neq IJI). I re-used
10
** my sortone.c code here, so don't let some of the variable names
13
void sortI_ROHF(void);
18
if(params.ref == 0 || params.ref == 1) sortI_ROHF();
19
else if(params.ref == 2) sortI_UHF();