t3evol - Numerical Solution of Twist-three Evolution Equations

A program dedicated to the numerical solution of the evolution equations for twist-three multiparton correlation functions is presented. The solutions are obtained by direct integration on a discretized momentum fraction grid. Both flavor nonsinglet and flavor singlet evolution (in both $C$-parity sectors) can be addressed. Physical applications include single spin asymmetries and the subleading twist contribution to the polarized structure function $g_2$. An arbitrary input is accepted for the initial distributions.