u_fun             165 src/secOrd_op_rhs_impl.cpp   auto u_fun = [&tp = tp](auto x, auto y){
u_fun             172 src/secOrd_op_rhs_impl.cpp     factor = a_til + ((a * a_til + e ) * mc - delta *u_fun(d[0], d[1]) )/ norm;