Currently, the LHFCALC docs (https://www.vasp.at/wiki/LHFCALC) state the following:
It is recommended to select a "direct optimization" algorithm for HF/DFT hybrid functonal type calculations, i.e., ALGO=Damped (IALGO=53) or ALGO=All (IALGO=58) in the INCAR file.
The blocked-Davidson algorithm ALGO=Normal is, with certain caveat, also supported, whereas calculations for the other algorithms (ALGO=Fast) are not properly supported (note: no warning is printed). The blocked-Davidson algorithm (ALGO=Normal) is generally rather slow, and in many cases the Pulay mixer will be unable to determine the proper ground-state.
However, as of VASP 6, LFOCKACE (https://www.vasp.at/wiki/index.php/LFOCKACE) enables a 3x speedup when ALGO = Normal. As such, the recommendation to users to use ALGO = Damped or ALGO=All for hybrid calculations seems not ideal. Since few readers probably find the LFOCKACE docs, they may not be taking advantage of the significant speedup enabled by using ALGO = Normal.
