RcppSimpleTensor: no se puede usar tensorFunction para calcular un objeto multivariante

I want to use the tensorFunction in package RcppSimpleTensor (https://github.com/tlamadon/RcppSimpleTensor) to compute a 3 dimensional function like this one:

DGP <- function( x,y,z ) {
    return((x + y - 5)^2 + (z-5)^2)
}

I create 3 vectors of values for x,y and z:

xdata <- seq(0,10,length=10)
ydata <- seq(0,10,length=10)
zdata <- seq(0,5,length=15)

I define the tensorFunction:

library(RcppSimpleTensor)
RcppVals <- tensorFunction( R[i,j,k] ~ (X[i] + Y[j] - 5)^2 + (Z[k] - 5)^2 )
RcppArray <- RcppVals(xdata,ydata,zdata)

But running this function gives an error "Error in tmpfun(X, Y, Z, k, j, i) : expecting a single value"

¿Qué estoy haciendo mal aquí?

preguntado el 09 de marzo de 12 a las 16:03

You may want to ask the author of that RcppSimpleTensor package. -

yes, I filed a bug on github. it's the seq() in xdata. for example as.array(seq(0,10,le-10)) works. -

1 Respuestas

This is now fixed in the main branch of the code,

mejor,

t.

contestado el 03 de mayo de 12 a las 19:05

No es la respuesta que estás buscando? Examinar otras preguntas etiquetadas or haz tu propia pregunta.