similar to: new package: elliptic

Hi List I'm thinking of putting together some R functionality (maybe a package) of elliptic functions. An elliptic function is one that is doubly periodic on the complex plane and I think they are cool. One objective of mine is to reproduce Abramowitz and Stegun's beautiful diagram on p643 (figure 18.5). The gsl package does sn, cn, dn et seq, but before I put a huge amount of

Hi has anyone coded up domain colouring for visualizing complex analytic functions (such as elliptic functions)? [ the idea is to depict a complex function f(z) using a filled.contour() variant in which the hue is given by Arg(f(z)), and the saturation by Mod(f(z)). ] -- Robin Hankin Uncertainty Analyst Southampton Oceanography Centre European Way, Southampton SO14 3ZH, UK tel

Dear List I have just uploaded a new package, "onion", to CRAN. It provides some functionality for manipulating and visualizing quaternions and octonions. All comments welcome! enjoy rksh -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743 _______________________________________________

Hi try this: p3dpairs <- function(x,x1, xlim=NULL,ylim=NULL,zlim=NULL,col=par("col"), pch=par("pch"), cex=par("cex"), ...){ if(is.matrix(x)){ z <- x[,3] y <- x[,2] x <- x[,1] } if(is.matrix(x1)){ z1 <- x1[,3] y1 <- x1[,2] x1 <- x1[,1] } if(missing(zlim)) { z.grid <-

Hi [MacOsX 10.4.6; R-2.2.1] I have a bundle that comprises three packages. I want to run R CMD check on each one individually, as it takes a long time to run on all three. I am having problems. The bundle as a whole passes R CMD check, but fails when I cd to the bundle directory and run R CMD check on a package directory. The whole bundle passes: octopus:~/scratch% R CMD check

Jacobi's theta functions crop up here and there in various branches of mathematics, such as unimodular and elliptic functions. They have Taylor expansions, but the powers increase as the square of n, as in 1 + z + z^4 + z^9 + z^16 + z^25 + . . . so they converge very quickly, provided |z|<1 The following toy function shows how I'm implementing these objects. I just add terms

Hello everybody. I am trying to upgrade from R-1.9.0 to R-1.9.1 on a RedHat linux 2.4.18 system. I get the following error after "tar -xvzf R-1.9.1.tgz ; cd ./R-1.9.1/ ; ./configure" and "make" : [make works for 10 minutes ... snip ...] varExp text html latex varFixed text html latex varFunc

Hi I have a package that I'm testing. It seems to install fine and it works, as far as I can tell. For example, I can install the package, and use it, and source the test suite with no errors. My problem is with R CMD check. It passes on R-2.2-0: Robin-Hankins-Computer:~/scratch% R CMD check ./partitions_1.1-0.tar.gz * checking for working latex ... OK [snip] make[1]: Leaving directory

Hi R version 2.4.0 alpha (2006-09-06 r39158) MacOSX 10.4.7 There was a thread some time ago as to whether the structure created by package.skeleton() would pass R CMD check. I have an example where package.skeleton() gives an R file that gives an error when sourced. If I type setClass("brob", representation = representation

Dear List please find on CRAN a new R bundle, BACCO, for Bayesian analysis of random functions, comprising two packages: emulator and calibrator. Package calibrator implements: ?Bayesian calibration of computer models?, M. C. Kennedy and A. O'Hagan 2001. Journal of the Royal Statistical Society B, 63(3) pp425-464 and package emulator implements: J. Oakley 2004. ?Estimating percentiles

Hi I'm developing an R package that needs to execute some code written in pari/gp. I've used this before from an R package (elliptic) but the interface is very basic: the R function creates a string such as the following: string <- echo ' ellwp ([ 2+0*I , 0+2*I ], 1+0*I )' | gp -q And then system(string) returns the output from gp which then needs to be text processed

Hello everyone I am finding complex numbers to be printed oddly in R-alpha. > 1e100+1e100i [1] 1e+100+1e+100i > 1e100 + 1e44i [1] 1e+100+100000000000000008821361405306422640701865984i > dput(1e100+1e44i) 1e+100+100000000000000008821361405306422640701865984i > [same thing at home on a linux box] The first line is what I would expect, the second and third lines aren't: I

Dear developers Neither Math.Rd nor Log.Rd mention the branch cuts that appear for complex arguments. I think it's important to include such information. Please find following two context diffs for Log.Rd and Math.Rd. [The pedants amongst us will observe that both sqrt() and log() have a branch point at complex infinity, which is not mentioned in the patch. Comments anyone?] rksh

Professor Ripley thanks for this. Very much appreciated. The original subject line reflected my late-night conviction that the answer might involve passing a strange list to do.call(). Anyway, package magic is broken (only in R-devel, I might add) because I have a function called "%eq%". R-2.0.0 CMD check is stopping (I think) because it interprets the "%" as a

Dear list R package "gsl" is now on CRAN. This is a wrapper for the special functions of the Gnu Scientific Library (GSL). Functions include elliptic integrals, Airy functions, hypergeometric functions, and so on. Most functions optionally return error estimates. This library is a little odd in that documentation is limited to a pointer to the GSL reference manual, and occasional

Dear list R package "gsl" is now on CRAN. This is a wrapper for the special functions of the Gnu Scientific Library (GSL). Functions include elliptic integrals, Airy functions, hypergeometric functions, and so on. Most functions optionally return error estimates. This library is a little odd in that documentation is limited to a pointer to the GSL reference manual, and occasional

Hi If string <- "xyz" f <- function(x){1 + sin(cos(x)) + exp(x^2)} How do I manipulate "string" and f() to give the string "1 + sin(cos(xyz)) + exp(xyz^2)" ? -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743

(this after asking the package author) Hi I cannot install the rmvnorm package under R-2.3.0, or R-2.3.1 beta. It installs fine under R-2.2.1. transcript for installation under R-2.3.0 follows. Robin-Hankins-Computer:~/scratch% R --version R version 2.3.0 (2006-04-24) Copyright (C) 2006 R Development Core Team R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to

Hi I have several methods for evaluating a function. The methods are algebraically identical but use different numerical techniques. The different methods work better (converge faster, etc) in different parts of the function's domain. I am compiling a test suite for a package, and would like to verify that the different methods return approximately identical results. Toy example

Hi I have been looking at complex.c and want to access z_cos() and z_sin () from C in one of my packages. There doesn't seem to be a corresponding header file: there is no complex.h file. Where are the prototypes of z_sin() and z_cos() for these functions? grepping didn't help me: find ~/downloads/R-2.1.1/ -name "*.h " | xargs egrep "z_cos"