vector

have made .size() out of for loop
have tried change vector to flat array, not big difference, so I keep vector
have tried less private variables in omp for, got slower
have tried pass by reference function call, got slower, might try this again
might try BLAS for matrix/ vector operation