levmartrif.cpp: a/b==a*(1/b), const