c-api vs. ctypes

[HerrHagen]

Wäre nett wenn du etwas Code posten würdest wenn du ein Beispiel mit OpenMP zum laufen gekriegt hast. Die Thematik interessiert mich schon seit längeren...

[HerrHagen]

Vieleicht interessant zu diesem Thema:
http://conference.scipy.org/proceedings ... l_text.pdf

[HWK]

Hier einmal eine Lösung mit f2py:
File test.f95

Code: Alles auswählen

subroutine test(arrayin, arrayout, m, n)
  implicit none
  REAL, dimension(m, n), intent(in) :: arrayin
  REAL, dimension(m, n), intent(out) :: arrayout
  integer, intent(in) :: m, n
  arrayout = arrayin + 3
end subroutine test

File testpy.py

Code: Alles auswählen

import numpy
import test

a = numpy.arange(9).reshape(3, 3)
b = test.test(a)
print a
print b

Ausgabe

Code: Alles auswählen

>C:\Programme\Python26\pythonw -u "testpy.py"
[[0 1 2]
 [3 4 5]
 [6 7 8]]
[[  3.   4.   5.]
 [  6.   7.   8.]
 [  9.  10.  11.]]
>Exit code: 0

MfG
HWK

[brennsuppa]

@openMP: ich konnte noch keine Performancesteigerungen feststellen, auch wenn beide Cores jetzt ausgelastet sind...muss den Code wahrsch. noch optimieren!

wenn ich den Code beim compilieren optimiere (-O2) rennt es schneller als auf meinem Dual Core!