Parallel computing and some changes

- add openmp for parallel computing
- rb_gauss_seidel_relaxation takes omega as relaxation parameter
- Wcycle() is added
- interpolate() doesn't create an extra array now