Karim Belabas on Mon, 17 Jun 2013 11:24:00 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: Ploth function |
* Ha Tran [2013-06-13 09:42]: > I have a problem in using the function "ploth" in pari gp. > > Actually, I need to compute a infinite sums: > > ff(x,vp,re,wp)=log > (sum(s1=-\infinity,\infinity,sum(s2=-\infinity,\infinity, > exp(g(x,vp,re,wp)))) > for a<x<b; > > and with the accuracy epsilon =10^(-5), by using a function to limit this > sum, let call limitshortvec(x,vp,re), I have the bounds of s1 and s2 (lis1 > and lis2), then I define: > > ff(x,vp,re,wp)=limitshortvec(x,vp,re);log > (sum(s1=-lis1,lis1,sum(s2=-lis2,lis2, exp(g(x,vp,re,wp)))) for > a<x<b; > > Finally, I have to "link" these functions "ff"" together, for example: > > h(x)=if(a<=x &&x<=b,ff(x,vp1,re1,wpr1),if(b<x &&x<=c,ff(x,vp2,re2,wpr2))) > > then plot h: > > ploth(x=a,c,h(x)) > > then Pari gp gives an announcement: > > â1347953757 [main] gp 6648 child_copy: linked dll data write copy failed, > 0x2CC000..0x2CF1B0, done 0, windows pid 6648, Win32 error 487â I have never seen this message in relation to PARI/GP. :-( (in other contexts, it's related to the --enable-auto-image-base gcc flag, which I thought was a thing of the past) 1) What exact version of Pari are you using. I.e. what is the exact name of the binary you installed ? 2) In case this is not the one you used, please try to reproduce the problem with http://pari.math.u-bordeaux1.fr/pub/pari/windows/Pari-2-6-0.exe 3) Can you send me the complete GP script in a separate file, including the limitshortvec() function ? Please make the example minimal: it should contain no code that's not absolutely necessary for you to reproduce the problem. Just reading in the file in the gp interpreter should trigger the bug. Cheers, K.B. -- Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17 Universite Bordeaux 1 Fax: (+33) (0)5 40 00 69 50 351, cours de la Liberation http://www.math.u-bordeaux1.fr/~kbelabas/ F-33405 Talence (France) http://pari.math.u-bordeaux1.fr/ [PARI/GP] `