1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
%module "Math::GSL::Poly"
// this brakes stuff
// %include "typemaps.i"
%include "gsl_typemaps.i"
%include "renames.i"
%{
#include "gsl/gsl_sys.h"
%}
%typemap(in) double * (double dvalue) {
SV* tempsv;
if (!SvROK($input)) {
croak("$input is not a reference!\n");
}
tempsv = SvRV($input);
if ((!SvNOK(tempsv)) && (!SvIOK(tempsv))) {
croak("$input is not a reference to number!\n");
}
dvalue = SvNV(tempsv);
$1 = &dvalue;
}
/* gsl_complex gsl_complex_poly_complex_eval (const gsl_complex c [], const int len, const gsl_complex z); */
%typemap(argout) gsl_complex {
AV* tempav = newAV();
double x,y;
if (argvi >= items) {
EXTEND(sp,1);
}
//fprintf(stderr,"--> %g <--\n", GSL_REAL($1));
//fprintf(stderr,"--> %g <--\n", GSL_IMAG($1));
$result = sv_newmortal();
x = GSL_REAL($1);
y = GSL_IMAG($1);
/* the next 2 lines blow up
sv_setnv($result, x);
argvi++;
*/
};
%typemap(argout) double * {
SV *tempsv;
tempsv = SvRV($input);
sv_setnv(tempsv, *$1);
}
%typemap(in) gsl_complex const [] {
AV *tempav;
I32 len;
int i, magic, stuff;
double x,y;
gsl_complex z;
SV **elem, **helem, **real, **imag;
HV *hash, *htmp;
SV *svtmp, tmp;
double result[2];
printf("gsl_complex typemap\n");
if (!SvROK($input))
croak("Math::GSL : $input is not a reference!");
if (SvTYPE(SvRV($input)) != SVt_PVAV)
croak("Math::GSL : $input is not an array ref!");
z = gsl_complex_rect(0,0);
tempav = (AV*)SvRV($input);
len = av_len(tempav);
$1 = (gsl_complex *) malloc((len+1)*sizeof(gsl_complex));
for (i = 0; i <= len; i++) {
elem = av_fetch(tempav, i, 0);
hash = (HV*) SvRV(*elem);
helem = hv_fetch(hash, "dat", 3, 0);
magic = mg_get(*helem);
if ( magic != 0)
croak("FETCH magic failed!\n");
printf("magic = %d\n", magic);
if( *helem == NULL)
croak("Structure does not contain 'dat' element\n");
printf("helem is:\n");
//Perl_sv_dump(*helem);
if( i == 0){
svtmp = (SV*)SvRV(*helem);
//Perl_sv_dump(svtmp);
}
printf("re z = %f\n", GSL_REAL(z) );
printf("im z = %f\n", GSL_IMAG(z) );
$1[i] = z;
}
}
%{
#include "gsl/gsl_inline.h"
#include "gsl/gsl_nan.h"
#include "gsl/gsl_poly.h"
#include "gsl/gsl_complex.h"
#include "gsl/gsl_complex_math.h"
%}
%include "gsl/gsl_inline.h"
%include "gsl/gsl_nan.h"
%include "gsl/gsl_poly.h"
%include "gsl/gsl_complex.h"
%include "gsl/gsl_complex_math.h"
%include "../pod/Poly.pod"
|
