adding RTS constraint without recreating matrices and vectors
changed output dea to use multi_glpk (5x speedup)
return thetaOpt and feasible in a vector