lgf               202 src/secOrd_op_initData_impl.cpp     lgf = rE*rA/(rE * ekt + rA * (1 - ekt));
lgf               204 src/secOrd_op_initData_impl.cpp   y *= lgf / norm;