search for: ifail

...ction nlm. Inside this function I have to integrate, hence I am using the function adapt. I don't understand why it is giving the following warnings: At the beginning: Warning: a final empty element has been omitted the part of the args list of 'c' being evaluated was: (paste("Ifail=2, lenwrk was too small. -- fix adapt() !\n", "Check the returned relerr!"), "Ifail=3: ndim > 20 -- rewrite the fortran code ;-) !", "Ifail=4, minpts > maxpts; should not happen!", "Ifail=5, internal non-convergence; should not happen!", ) When it...

R-listers, In using the adapt function, I am getting the following warning: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(ndim = 2, lower = lower.limit, upper = upper.limit, functn = pr.set, Would someone explain what the 'lenwrk' value indicates in order to help diagnose this issue. Also, what are the possible codes for Ifail,...

Dear People, I'm trying to use the function adapt, from the adapt library package, which does multidimensional numerical integration. I think I must be using the wrong syntax or something, because even a simple example does not work. Consider foo <- function(x){x[1]*x[2]} and adapt(2, lo = c(-1,-1), up = c(1,1), functn = foo) This simply hangs. A more complicated example crashes R,

...kc <- x[2] + f(a) * g(kc) * H(a,kc) } + adapt(2,lower.limits,upper.limits,minpts=200,maxpts=20000,prob.of.fracture,e ps=1e-3) } > risk(eifs.at.time.using.density(5000,eifs,1e-4,2),g,H,c(0.001,20),c(.5,60)) [1] 1.999471 value relerr minpts lenwrk ifail 28.39852 0.0003455776 627 2803 0 > > eifs.5000 <- eifs.at.time.using.density(5000,eifs,1e-4,2) [1] 1.999471 > risk(eifs.5000,g,H,c(0.001,20),c(.5,60)) value relerr minpts lenwrk ifail 5.461088e-13 0.0006449058...

...gt; dummy(-100) [1] 0 > dummy(-1000) [1] NaN > dummy(-10000) [1] NaN If I choose the lower boundary to be too small integrate causes a segmentation fault: > library(integrate) > integrate(dummy, -100, 0)$value [1] 1.387143 > integrate(dummy, -1000, 0)$value [1] NaN Warning message: Ifail=1, maxpts was too small. Check the returned relerr! in: adapt(1, lower, upper, minpts, maxpts, functn, eps, ...) > integrate(dummy, -10000, 0)$value Segmentation fault I'm running R.1.2.0 on Debian Linux. Regards Achim > version _ platform i686-pc-linux-gnu ar...

Hi all, I need to calculate a multidimensional integration on R. I am using the command "adapt" (from library adapt), although sometimes I get the following error message: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(3, linf, lsup, functn = Integrando1) I guess it happens because the domain of integration is too small, although I tried a change of variables to avoid this problem and it didn?t help. The command adapt calls a fortran...

...use 'memory not mapped' Traceback: 1: .C("cadapt", as.integer(ndim), as.double(lower), as.double(upper), minpts = as.integer(minpts), maxpts = as.integer(maxpts), ff, rho = environment(), as.double(eps), relerr = double(1), lenwrk = as.integer(lenwrk), value = double(1), ifail = integer(1), PACKAGE = "adapt") 2: adapt(2, c(-5, 0, -0.5), c(5, 10, 0.5), f = unsc.post, max = 5e+05, val = data) On a perhaps related note, if I increase the length of data to 1000+10 in the original example, I get a garbage answer and: Warning message: Ifail=2, lenwrk was too s...

Dear R-helpers! What it's means "Rank failure in Choleski decomposition" by using of function surf.gls {spatial}? Sincerely Yours Branimir K. Hackenberger

how does one numerically intergate the following: A=function(x,y) { xy } over the range: 2<x<0 4<y<10 say. ie how would one set up the integrate function? i forgot!

Dear r-users, Is there any command in R allowing to evaluate a double integral? for instance let say I want to evaluate the following integral: integrate[lo=(0,1),up=(2,3)] f(x,y)=x^2+y^2 where lo is the vector of lower bounds and up that of upper bounds. I thaught the function "adapt" would work but it did not. Many thanks, Dominique K.

Hi, I need to integrate a 2D function over range where the limits depend on the other e.g integrate f(x,y)=x*y over {x,0,1} and {y,x,1}. i.e \int_0^1 \int_x^1 xy dydx I checked adapt but it doesn't seem to help here. Are they any packages for this sort of thing? I tried RSitesearch but couldn't find the answer to this. Many thanks for you help. Regards Saptarshi Saptarshi Guha

Hi all, Running R 0.90.1 on a RH 6.1 system. Installation of the integrate_2.1-2 package went smoothly. My code contains a loop in which integrate() is called several times in each pass. I get a segmentation fault after what seems to be a random number of calls to integrate(). Debug output shows: Program received signal SIGSEGV, Segmentation fault. promiseArgs (el=0x40276414,

Hello, there, I got this error when i tried to run " data.kr <- surf.gls(2, expcov, data, d=0.7);" "Error: cannot allocate vector of size 382890 Kb Execution halted" My data is 100x100 grid. the following is the summary of "data": > summary(data); x y z Min. : 1.00 Min. : 1.00 Min. :-1.0172 1st Qu.: 26.00