Checks if the entry is neglected due to frozen orbitals being included and absorbs entries into the lower order matrix elements if required
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(inout) | :: | matel | |||
| integer(kind=int64), | intent(inout) | :: | indices(num_inds) |
subroutine freeze_lmat(matel, indices)
HElement_t(dp), intent(inout) :: matel
integer(int64), intent(inout) :: indices(num_inds)
! Offset the orbital indexing
call map_indices(indices)
call add_core_en(matel, indices)
end subroutine freeze_lmat