more stable linearization energy algorithm
work on omp in linearization