a_til             170 src/secOrd_op_rhs_impl.cpp     a_til = k * (1 - norm/rE),
a_til             172 src/secOrd_op_rhs_impl.cpp     factor = a_til + ((a * a_til + e ) * mc - delta *u_fun(d[0], d[1]) )/ norm;