openmp do loops

[Stefan Behnel]

Pascal, 04.11.2010 11:13:
> I would like to parallelize this loop:
> do i=1,hklsize
> fcalctable(i)=structfact(hkltable(1,i),hkltable(2,i),hkltable(3,i))
> end do
>
> I thought I would do this:
> !$OMP PARALLEL DO default(private) shared(hkltable, fcalctable,hklsize)
> do i=1,hklsize
> fcalctable(i)=structfact(hkltable(1,i),hkltable(2,i),hkltable(3,i))
> end do
> !$OMP END PARALLEL DO

Seeing this makes me seriously happy that I can write my code in Python (or 
Cython, for that purpose).

Stefan