feat($Zhuang09-OpenMP): Parallel version

    Parallel version. Seems to work.
    A lot of fixes  including:
	- lr/lw_iter needs to be of size nnz since they are private to each thread.
	- p has a offset of 1.
		For example first iteration is for p=0,
		DDA(0) should be 1 and wave(0) should be superior to -1
		so every occurence of p is including the offset of 1
	- the arrays wavestart, wavefronts are done by only one thread as it needs the number of waves.