>Enter ATOM, relativistic (Y/N) with mass correction (Y/N) A,y,y >Gradient or Slater form? (G/S): g > Enter Maximum and minimum values of 2*J 4,2 >Enter eigenvalues: one line per term, eigenvalues separated by commas 2*J = 4 2,3,4 2*J = 2 1,2 .......
To compute the number of matrix elements, and respectively to allocate memory, alcmat() and alctab() rely on a call to jmem(), which in effect reads hnr.lst, hzeta.lst, hspin.lst and computes the number of nonzero elements in the interaction matrix. If the requested memory is larger than the what is available, the program continues the calculation, by allocating memory for the largest column. In this case the dvdson() performs matrix-vector multiplication (in dssbmv), keeping in memory only a single column.
The matrix elements are assembled in (hmx_lsj()), then the interaction matrix is diagonalized (lsjmat()). Next, dvdson() finds the requested eigenvalues. Finally, brci_zeeman(), computes the Zeman factor. The requested eigenvectors are saved in a <name>.j file.