[Rd] system() in packages
Prof Brian Ripley
ripley at stats.ox.ac.uk
Fri Apr 27 16:20:24 CEST 2007
The first comment is that will only work on Unix-alikes, since '|' needs a
shell.
So, if this is on a Unix-alike you need to establish if the program is in
the path at run time and cache the result. I have no idea if this would
actually work, but for example system('gp --version') might provide a
suitable test.
On Fri, 27 Apr 2007, Robin Hankin wrote:
> Hello
>
> Quite often I need higher precision, or larger numbers,
> than IEEE double precision allows. One strategy I sometimes
> use is to use system() to call pari/gp, which is not constrained by
> IEEE.
>
> [pari/gp is a GPL high-level mathematical programming language
> geared towards pure mathematics]
>
> Two of my packages contain lines like the following:
>
> > system(" echo '1.12^66' | gp -q f",intern=TRUE)
> [1]
> "1771.697189476241729649767636564084681203806302318041262248838950177194
> 116346432205160921568393661760"
>
> Note the high precision of the answer.
>
> My question is, how to deal with the possibility that pari/gp is not
> installed?
>
> If the system cannot find gp for some reason, I get:
>
> > system(" echo '1.12^66' | gp -q f",intern=TRUE)
> sh: line 1: gp: command not found
> character(0)
> >
>
> What's the recommended way to handle this eventuality gracefully? The
> functions that do use pari/gp have "pure" R equivalents (but much
> slower and less accurate) so I want users to be able to install the
> package without pari/gp.
>
>
> --
> Robin Hankin
> Uncertainty Analyst
> National Oceanography Centre, Southampton
> European Way, Southampton SO14 3ZH, UK
> tel 023-8059-7743
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
--
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-devel
mailing list