Beefy Boxes and Bandwidth Generously Provided by pair Networks
go ahead... be a heretic
 
PerlMonks  

Using GP with Math::Pari

by hv (Parson)
on Dec 18, 2007 at 14:00 UTC ( #657677=perlquestion: print w/ replies, xml ) Need Help??
hv has asked for the wisdom of the Perl Monks concerning the following question:

I was asking this the other day, here's what I found out.

I wanted to know how to use Math::Pari to declare and then use a new function written in GP. After a brief exchange with the author, I discovered that the only problem was that Math::Pari hiccups on braces and spaces - and these are in any case not needed (at least for dealing with numbers), they are pure "syntactic sugar".

So I knocked up a local function:

use Math::Pari qw/ PARI /; sub PARIdo { my $s = shift; $s =~ s/[{}\s]+//g; PARI($s); }
.. and was then able to use it to write nice clean GP function declarations:
PARIdo(q{ mygcd(a, b) = { if (a > b, { mygcd(b, a); }, { if (a == 0, { b; }, { mygcd(b % a, a); }); }); } }); print PARI("mygcd(12,8)");

I don't yet have a good solution for passing arguments in. You can interpolate them:

my $v1 = PARI(12); my $v2 = PARI(8); print PARI(qq{mygcd($v1,$v2)});
.. but it is distressingly inefficient to convert large numbers to string and back. I've done some experiments with setting up a PARIvar, but haven't worked out yet if it is possible to assign to such a variable from perl space.

I'd welcome any other tips and tricks for getting the most out of perl/PARI hybrid code.

Hugo

Comment on Using GP with Math::Pari
Select or Download Code

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: perlquestion [id://657677]
Approved by Corion
Front-paged by Corion
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others lurking in the Monastery: (7)
As of 2014-09-17 03:47 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    How do you remember the number of days in each month?











    Results (57 votes), past polls