Merge pull request #15 from gjbex/development

Update for Cython 3.x

Author: gjbex

python_for_hpc.pptx

@@ -0,0 +1 @@
+points.c

source-code/cython/Classes/.gitignore

@@ -0,0 +1,12 @@
+VERSION = cpython-311-x86_64-linux-gnu
+POINTS_LIB = points.$(VERSION).so
+POINTS_PURE_LIB = points_pure.$(VERSION).so
+
+all: $(POINTS_LIB)
+
+$(POINTS_LIB): points.pyx points_pure.py
+	python setup.py build_ext --inplace
+
+clean:
+	python setup.py clean
+	$(RM) points.c points_pure.c $(POINTS_LIB) $(POINTS_PURE_LIB)

source-code/cython/Classes/Makefile

@@ -0,0 +1,22 @@
+# Classes
+
+Illustration of using Cython extension types (aka cdef classes).
+
+## What is it?
+
+1. `points.pyx`: implementation of a cdef class, and a Python
+   child class thereof.
+1. `points_pure.py`: implementation in Cython's pure Python
+   syntax.
+1. `points_python.py`: Python implementation of the class.
+1. `driver.py`: Python script that uses both classes.
+1. `setup.py`: Python installation file to build the Cython
+   extension.
+1. `distances_cython.py`: Python script to compute distances between
+   points using Cython class.
+1. `distances_python.py`: Python script to compute distances between
+   points using Python class.
+1. `distances_internal.py`: Python script to compute distances between
+   points using Cython class, computing the distances using a static
+   class method..
+1. `Makefile`: make file to build the extension.

source-code/cython/Classes/README.md

@@ -0,0 +1,28 @@
+#!/usr/bin/env python
+
+import argparse
+import itertools
+import pyximport
+pyximport.install(pyimport=True, language_level='3str')
+from points_pure import Point
+import random
+
+
+def main():
+    arg_parser = argparse.ArgumentParser(description='compute distances')
+    arg_parser.add_argument('--n', type=int, default=10,
+                            help='number of points')
+    options = arg_parser.parse_args()
+    points = [Point(random.random(), random.random()) for _ in range(options.n)]
+    min_distance = 2.0
+    max_distance = 0.0
+    for i, p1 in enumerate(points):
+        for p2 in points[i+1:]:
+            d = p1.distance(p2)
+            min_distance = min(d, min_distance)
+            max_distance = max(d, max_distance)
+    print(f'min. distance: {min_distance}')
+    print(f'max. distance: {max_distance}')
+
+if __name__ == '__main__':
+    main()

source-code/cython/Classes/distances_cython.py

@@ -0,0 +1,22 @@
+#!/usr/bin/env python
+
+import argparse
+import itertools
+import pyximport
+pyximport.install(pyimport=True, language_level='3str')
+from points_pure import Point
+import random
+
+
+def main():
+    arg_parser = argparse.ArgumentParser(description='compute distances')
+    arg_parser.add_argument('--n', type=int, default=10,
+                            help='number of points')
+    options = arg_parser.parse_args()
+    points = [Point(random.random(), random.random()) for _ in range(options.n)]
+    min_distance, max_distance = Point.min_max_distance(points)
+    print(f'min. distance: {min_distance}')
+    print(f'max. distance: {max_distance}')
+
+if __name__ == '__main__':
+    main()

source-code/cython/Classes/distances_internal.py

@@ -0,0 +1,26 @@
+#!/usr/bin/env python
+
+import argparse
+import itertools
+from points_python import Point
+import random
+
+
+def main():
+    arg_parser = argparse.ArgumentParser(description='compute distances')
+    arg_parser.add_argument('--n', type=int, default=10,
+                            help='number of points')
+    options = arg_parser.parse_args()
+    points = [Point(random.random(), random.random()) for _ in range(options.n)]
+    min_distance = 2.0
+    max_distance = 0.0
+    for i, p1 in enumerate(points):
+        for p2 in points[i+1:]:
+            d = p1.distance(p2)
+            min_distance = min(d, min_distance)
+            max_distance = max(d, max_distance)
+    print(f'min. distance: {min_distance}')
+    print(f'max. distance: {max_distance}')
+
+if __name__ == '__main__':
+    main()

source-code/cython/Classes/distances_python.py

@@ -0,0 +1,11 @@
+#!/usr/bin/env python
+
+from points import Point, ColoredPoint
+
+p = Point(1.0, -2.0)
+print(p)
+print(f'point = {p.x}, {p.y}')
+print(p.distance())
+
+p1 = ColoredPoint(1.0, -2.0, 'blue')
+print(p1.color)

source-code/cython/Classes/driver.py

@@ -0,0 +1,14 @@
+#!/usr/bin/env python
+
+import pyximport
+pyximport.install(pyimport=True, language_level='3str')
+
+from points_pure import Point, ColoredPoint
+
+p = Point(1.0, -2.0)
+print(p)
+print(f'point = {p.x}, {p.y}')
+print(p.distance())
+
+p1 = ColoredPoint(1.0, -2.0, 'blue')
+print(p1.color)

source-code/cython/Classes/driver_pure.py

@@ -0,0 +1,39 @@
+from libc.math cimport sqrt
+
+cdef class Point:
+    
+    cdef double _x, _y
+    
+    def __init__(self, x, y):
+        self._x = x
+        self._y = y
+        
+    cpdef distance(self, other):
+        return sqrt((self._x - other._x)**2 + (self._y - other._y)**2)
+    
+    property x:
+        def __get__(self):
+            return self._x
+        def __set__(self, value):
+            self._x = float(value)
+            
+    property y:
+        def __get__(self):
+            return self._y
+        def __set__(self, value):
+            self._y = float(value)            
+
+
+class ColoredPoint(Point):
+
+    def __init__(self, x, y, color):
+        super().__init__(x, y)
+        self._color = color
+
+    @property
+    def color(self):
+        return self._color
+
+    @color.setter
+    def color(self, value):
+        self._color = str(value)