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