Python C/C++ IPMU IRAF

Size: px

Start display at page:

Download "Python C/C++ IPMU IRAF"

ふみなすわ
4 years ago
Views:

1 Python C/C++ IPMU IRAF

2 Python Swig Numpy array Image

3 Python swig numpy pyfits 2.3 pyds9 1.1 svn co hjp://svn.scipy.org/svn/numpy/tags/1.5.0/doc/swig swig/numpy.i /usr/local/share/swig/1.3.40/python swig

4 Python

5 Hello Python Python python/hello.py % python Python (r266:84292, Nov , 16:27:48) [GCC (Red Hat )] on linux2 Type help, copyright, credits or license for more information. >>> print Hello Python! Hello Python! >>> print >>> ^D (Crtl+D) % % python hello.py Hello Python! 8 %

6 python/loop.py a = [ abc, def, ghi, jkl ] for b in a: print b for for i in range(5): print i range(5) [0,1,2,3,4]

7 I/O python/io.py file = open( test.txt ) file = open( out.txt, w ) line = file.readline() file.write(line)

8 python/func.py def times(x, y): return x * y import func from func import * x = func.times(3.14, 4)

9 pyfits FITS hdulist = pyfits.open( test.fits ) hdulist[0].header[ OBJECT ] hdulist[0].data[y, x] pyfits/sample.py numpy array 0- indexed ds9 [y, x]

10 pyds9 ds9 X Public Access (XPA) XPA Access Points d = ds9.ds9() ds9 d.set( file test.fits ) pyds9/sample.py array d.set( array [xdim=2080,ydim=4100,bitpix=-32], data)

11 Python Numpy hjp:// pyfits/stats.py Python Numpy

12 Python C/C++ Python C/C++ Python C/C++

13 Python C/C++ boost.python swig ctypes swig

14 Swig C swig/example.*

15 C /* File : example.c */ swig/example.c #include <time.h> double My_variable = 3.0; int fact(int n) { if (n <= 1) return 1; else return n*fact(n-1); } int my_mod(int x, int y) { return (x%y); } char *get_time() { time_t ltime; time(&ltime); return ctime(&ltime); }

16 Swig swig/example.* C Swig

17 Swig swig/example.i /* example.i */ %module example Python %{ /* Put header files here or function declarations like below */ extern double My_variable; extern int fact(int n); extern int my_mod(int x, int y); extern char *get_time(); %} extern double My_variable; extern int fact(int n); extern int my_mod(int x, int y); extern char *get_time();

18 Swig C Swig swig % swig -python example.i example.py, example_wrap.c % gcc fpic c example.c example_wrap.c I/adc/local/include/ python2.6/ example.o, example_wrap.o shared object % gcc -shared example.o example_wrap.o -o _example.so _example.so swig/example.* compile.sh

19 import example swig/sample.py print example.cvar.my_variable print example.fact(4) print example.my_mod(10,3) print example.get_time() <>.cvar.< > Python

20 Swig output/example.* void add(double a, double b, double *result) %include typemaps.i %apply double *OUTPUT { double *result }; typemaps.i void mod(int a, int b, int *quo, int *rem) %include typemaps.i # %apply int *OUTPUT { int*quo, int *rem };

21 Swig import example output/sample.py a = example.add_0(3,4) print a b = example.add(3,4) print b quo, rem = example.mod(7,3) print quo, rem

22 Swig C++ vector double average(std::vector<int> v) std::vector<double> half(const std::vector<double>& v) %include std_vector.i %template(intvector) std::vector<int>; %template(doublevector) std::vector<double>; vector/example.*

23 Swig from example import * iv = IntVector(4) for i in range(0,4): iv[i] = i print average(iv) vector<int> vector/sample.py print average([0,1,2,3]) print half([1,2,3]) Swig

24 Numpy array C Numpy array rms C double rms(double *seq, int n) Python v = rms(array) numpy/rms.* numpy.itypemaps

25 Numpy array C %module rms %{ #define SWIG_FILE_WITH_INIT #include rms.h %} %include numpy.i %init %{ import_array(); %} numpy/rms.i %apply (double *IN_ARRAY1, int DIM1) {(double *seq, int n)}; %include rms.h 1 array numpy_swig.pdf

26 CNumpy array ezrange/ezrange.i Cnumpy array Python %apply (double *IN_ARRAY1, int DIM1) {(double *seq, int n)}; %apply (int *ARGOUT_ARRAY1, int DIM1) {(int *rangevec, int n)}; array Python

27 array 2 array(=) %apply (double *IN_ARRAY1, int DIM1) %apply (double *IN_ARRAY2, int DIM1, DIM2) C array double mean(double *array, int n, int m) numpy2/stats.* Numpy

28 Python C++

29 Image image0/image.* Numpy array Image Numpy array vector class Image { int npx, npy, npix; double *data; public: Image(); Image(double *array, int n, int m); ~Image(); int getnpx() const; int getnpy() const; int getnpix() const; void getnaxis(int *npx, int *npy); double getvalue(int x, int y) const; };

30 Image image0/image.* Image::Image(double *array, int n, int m) { data = new double[n*m]; // for (int j = 0; j < n; j++) { for (int i = 0; i < m; i++) { data[i+j*m] = array[i+j*m]; // } } npx = m; npy = n; npix = n*m; }

31 Image image0/sample.py pyfits Image img = Image.Image(hdulist[0].data) ds9 image01/image.*

32 median median image1/image.* medianmedian

33 Image::Image(int npx, int npy) void Image::setValue(int x, int y double v) void Image::getData(float *oarray, int n) numpy array Image& Image::operator/=(double a) a double Image::median() median static Image Image::getMedian(std::vector<Image>& v) vmedian image1/image.*

34 ds9 image2/image.* n x n

35 Image& Image::operator/=(const Image& a) image2/image.* Image Image::boxavg(int width) width x width std::vector<point> findpeaks(double thres) thres vector Point class Point { int x, y; public: Point(int x, int y); int getx(); int gety(); }

36 Image median Gaussian

37 Swig C/C++ Python Swig Python C++

Debugging Common Issues in Multithreaded Applications

Debugging Common Issues in Multithreaded Applications Python と C/C++ の混合デバッグ MIXED-LANGUAGE DEBUGGING WITH PYTHON AND C/C++ ( 日本語参考訳 ) PYPL *1 を始めとするプログラミング言語の人気ランキングサイトによると Python の人気は 10 年以上にわたって着実に伸びています使いやすさ明確な構文サードパーティ製の豊富なドキュメント多数の拡張ライブラリがあるため