Add general frozen velocity operator/form

This is mainly to add the frozen velocity operator/form, which can be used when the analytical form of density is not known.

Major change is the extra form.