n_p 141 src/secOrd_op_solutionDriven_impl.cpp auto n_p = surface_normal(g); n_p 144 src/secOrd_op_solutionDriven_impl.cpp n_p *= u_p * tp.deltaT() * delta * integral_factor; n_p 148 src/secOrd_op_solutionDriven_impl.cpp X_p += n_p;