kt                 61 src/secOrd_op_initData_impl.cpp   const double r0 = 1., r_end = 2., kt = .5 * t;
kt                 62 src/secOrd_op_initData_impl.cpp   return r0 / (r_end * exp(-kt) + r0 * (1 - exp(-kt)) );