I am developing functional domain specific embedded language within C++ to translate formulas into working code as concisely and accurately as possible.
I posted a proto
T = T - T'/e;
or, if you're not operating on all of T
T(0:i,0:j) = T(0:i,0:j) - T(0:i,0:j)'/e(0:i,0:j)