ReadFPLOBasisFunctions(orbitals,file1,file2) reads the files specified by file1 and file2 that contain the wavefunctions obtained from the FPLO DFT calculation and returns interpolating functions corresponding to the names given by orbitals.
For non-relativistic wave functions the orbitals are labeled “nl” where n is the principle quantum number and l is the letter corresponding to the angular momentum quantum number. For relativistic wave functions the orbitals are labeled “nlj” where j is the total angular momentum quantum number.
orbs = {"3s","3p","3d","4d"} f = ReadFPLOBasisFunctions(orbs,"+fval.001.1") orbs = {"3d3/2","3d5/2"} G, F = ReadFPLOBasisFunctions(orbs,"+fval.001.1","+fval.001.2")