changed to use mkl for second order tensor mode product
CULA inclusion and preallocated memory for tensor mode product
bug fix