mò «jwEc@s>dkZdklZdkZd„ZeƒZd„ZdS(N(s ext_toolsc Cstidddƒ}tidddƒitiƒ}|id}ti ddgtiƒ}|id}d }d}d}d }d }d }tid |ddddddgƒ}|i|ƒ|ii|ƒd} tiii| ƒ\}}d} }dhdg}ti"|ƒD] }||q~ dd!<dg}|D] }||qH~ dd!<d| <}d} tid| ddddgƒ}|ii|ƒ|i|ƒti'd&tiƒ} d!} d"}tid#|d$gƒ}|ii| ƒ|i|ƒ|S('sT Builds an extension module with Bspline basis calculators using weave. t_bsplinetcompilertgcciii f0.40000000000000002f0.5ii s5 double *bspline(double **output, double *x, int nx, double *knots, int nknots, int m, int d, int lower, int upper) { int nbasis; int index, i, j, k; double *result, *b, *b0, *b1; double *f0, *f1; double denom; nbasis = upper - lower; result = *((double **) output); f0 = (double *) malloc(sizeof(*f0) * nx); f1 = (double *) malloc(sizeof(*f1) * nx); if (m == 1) { for(i=0; i= knots[index]) * (x[k] < knots[index+1]); result++; } } else { for (k=0; kdata; bspline(&data, x, Nx[0], knots, Nknots[0], m, d, lower, upper); return_val = (PyObject *) basis; tevaluatetxtknotstmtdtlowertupperiisƒ double *bspline_prod(double *x, int nx, double *knots, int nknots, int m, int l, int r, int dl, int dr) { double *result, *bl, *br; int k; if (fabs(r - l) <= m) { result = (double *) malloc(sizeof(*result) * nx); bl = (double *) malloc(sizeof(*bl) * nx); br = (double *) malloc(sizeof(*br) * nx); bl = bspline(&bl, x, nx, knots, nknots, m, dl, l, l+1); br = bspline(&br, x, nx, knots, nknots, m, dr, r, r+1); for (k=0; k nknots - 1) {upper = nknots-1;} /* upper = nknots - m; */ for (k=lower; kdata; bspline_gram(&data, knots, Nknots[0], m, dl, dr); return_val = (PyObject *) gram; tgramtdltdrii sÁ void invband_compute(double **dataptr, double *L, int n, int m) { /* Note: m is number of bands not including the diagonal so L is of size (m+1)xn */ int i,j,k; int idx, idy; double *data, *odata; double diag; data = *((double **) dataptr); for (i=0; i 0) { data[j*n+i] = 0;} } } for (i=n-1; i>=0; i--) { for (j=1; j <= (mdata; invband_compute(&data, L, NL[1], NL[0]-1); return_val = (PyObject *) invband; tinvbandtLN(ii (,t ext_toolst ext_moduletmodtNtlinspacetastypetfloat64RtshapetnknotstarrayRtnxRRRR t eval_codet eval_ext_codet ext_functiont bspline_evalt add_functiont customizetadd_support_codeR tscipytspecialt orthogonaltp_rootsR R RRt_[1]trealtqt gram_codet gram_ext_codet bspline_gramtzerosRtinvband_support_codetinvband_ext_codeR(RR RR-RR*RR+RR RR/RR,RRRR0RR R RR(RRR ((t$Lib/sandbox/models/bspline_module.pytsetup_bspline_modulesH!  k   tc  (   cCstiƒdS(N(Rtcompile(((R1tbuild_bspline_moduleUs( tnumpyRt scipy.weaveRtscipy.special.orthogonalR$R2RR4(R2RR$RR4R((R1t?s    ÿO