output residuum
added omp parallelization