direct_solve, matrix trans

Improved performance of triangular solver, created benchmark direct_solve,at this time works on Cuda and OpenMP backends
B=trans(B) resize problem solved, improved B=trans(B) performance