On Tuesday, June 7, 2016 at 7:11:08 AM UTC-7, William wrote:
>
> Hi,
>
> Long ago in 2006, I put Mark Watkins amazing C program "sympow" in Sage:
>
>
> http://doc.sagemath.org/html/en/reference/lfunctions/sage/lfunctions/sympow.html
>
>
> In case you're not a number theorist, this program computes values of
> generating functions
> attached to symmetric power representations attached to elliptic
> curves, and I think it's the only open source
> program that does it. Symmetric power representations turn out to be
> extremely important in number
> theory, e.g., they are used to prove the Sato-Tate conjecture, which
> is a major result.
>
> The documentation cited above begins
>
> sage: sympow('-new_data 2') # not tested
>
> So.. it's not tested, and with Sage, as we all know, "if it isn't
> tested, then it's completely horribly broken".
>
> I just tried to use Sympow, since it would be extremely useful for
> some research I'm doing with Barry
> Mazur during my visit to Harvard this week. However, it just
> segfaults in both sage-6.10 **and** sage-7.3.beta.
> It silently fails when doing the above "not tested" thing, but if you
> try it at the command line it segfaults:
>
> -----
> (sage-sh) 95d92fa7cb50414ea35d9897eabe44de@compute4-us:sympow$ sympow
> -new_data 2
> Make data for symmetric power 2
> Running the new_data script for -sp 2
> Making the datafiles for -sp 2
>
> Rewarping the param_data file
> Left with 13 entries in param_data
> Segmentation fault (core dumped)
> -----
>
It works for me on OS X, both within Sage and from the command line:
┌────────────────────────────────────────────────────────────────────┐
│ SageMath version 7.3.beta2, Release Date: 2016-05-28 │
│ Type "notebook()" for the browser-based notebook interface. │
│ Type "help()" for help. │
└────────────────────────────────────────────────────────────────────┘
┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓
┃ Warning: this is a prerelease version, and it may be unstable. ┃
┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛
sage: sympow('-new_data 2')
'Running the new_data script for -sp 2\nMaking the datafiles for -sp
2\n\nRewarping the param_data file\nLeft with 13 entries in
param_data\necho \'Removing any old data files\'\ncd datafiles\nrm -f
P02HM.txt P02HS.txt P02HM.bin\nrm -f P02LM.txt P02LS.txt P02LM.bin\ncd
..\nRemoving any old data files\nRunning the gp script\n\nN=600; dv=0;
mx=1;\n\\p 250\nSTR="P02H";\n\\r
standard1.gp\nF(k)=if(k%2==0,J(k-2,X)/1!*J(k/2-1,X/2)*sinv(k,X),sqrt(Pi)/2*J(k-1,X)/0!*J(k-2,X)/1!*1/J((k-1)/2,X/2)*two1ms(k,X)*sinv(k,X))\n\\r
standard2.gp\n\\l datafiles/P02HM.txt\n\\r standard3.gp\n\\l
datafiles/P02HS.txt\ncoeffs(0);\ncoeffE(1);\nSTR="P02L";\n\\r
standard1.gp\nF(k)=if(k%2==1,J(k-1,X)/1!*J((k-1)/2,X/2)*sinv(k,X),sqrt(Pi)/2*J(k-1,X)/1!*J(k-1,X)/J(k/2-1,X/2)*two1ms(k,X)*sinv(k,X))\n\\r
standard2.gp\n\\l datafiles/P02LM.txt\n\\r standard3.gp\n\\l
datafiles/P02LS.txt\ncoeffs(0);\ncoeffO(1);\n\\q\n\necho \'Trimming the
data files\'\ncd datafiles\ngrep -v \'^?\' P02HM.txt | sed \'s/ E/e/\' >
.tempfile.123\\\n && echo \'END\' >> .tempfile.123 && mv .tempfile.123
P02HM.txt\ngrep -v \'^?\' P02HS.txt | sed \'s/ E/e/\' > .tempfile.123\\\n
&& echo \'END\' >> .tempfile.123 && mv .tempfile.123 P02HS.txt\ngrep -v
\'^?\' P02LM.txt | sed \'s/ E/e/\' > .tempfile.123\\\n && echo \'END\' >>
.tempfile.123 && mv .tempfile.123 P02LM.txt\ngrep -v \'^?\' P02LS.txt | sed
\'s/ E/e/\' > .tempfile.123\\\n && echo \'END\' >> .tempfile.123 && mv
.tempfile.123 P02LS.txt\necho \'Turning the meshes into
binaries\'\nNUM=`grep -c AT P02HM.txt`\n../sympow -txt2bin $NUM P02HM.bin <
P02HM.txt\nNUM=`grep -c AT P02LM.txt`\n../sympow -txt2bin $NUM P02LM.bin <
P02LM.txt\ncd ..\nTrimming the data files\nTurning the meshes into
binaries\nRewarping the param_data file\nLeft with 15 entries in
param_data\nFinished with -sp 2'
(At least, it doesn't segfault or fail silently. I don't know what it's
supposed to return.)
Does the installation log file for sympow say anything interesting?
--
John
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
