similar to: optim "CG" bug w/patch proposal (PR#8786)

Dear R team, when using optim with method "CG" I got the wrong $value for the reported $par. Example: f<-function(p) { if (!all(p>-.7)) return(2) if (!all(p<.7)) return(2) sin((p[1])^2)*sin(p[2]) } optim(c(0.1,-0.1),f,method="CG",control=list(trace=0,type=1)) $par 19280.68 -10622.32 $value -0.2346207 # should be 2!

Hi, I get a segmentation fault in eval.c: do_set (call=0x8ce166c, op=0x82b43b8, args=0x8ce1650, rho=0x8acf218) at eval.c:1309 1309 switch (NAMED(s)) { The instruction is 0x080becf1 <do_set+177>: movzbl (%eax),%edx and %eax contains a non-accessible address (gdb) x/x $eax 0x4a5bf008: Cannot access memory at address 0x4a5bf008 in another try it was (gdb) x/x

Hi, I have a problem with the queue. My system is 'Asterisk 1.8.20.0 built by mockbuild @ buildvm-24.phx2.fedoraproject.org on a x86_64 running Linux' (core show version) and my OS is Red Hat Enterprise Linux Server release 6.3 (Santiago). I have six queues...into the top five I run the application background with a choice (1-5) then redirect to a queue, but here the caller doesn't

I am trying to get something working that is just not doing quite what I want. It may not be possible, but I figured it was worth asking about. The details: Asterisk 11.6.0 Polycom SoundPoint IP650 phones running 4.03 firmware. We have a queue with 4 phones in it. ringinuse is set to yes and the stategy is ringall. In sip.conf, we have notifyringing set to yes as well. Asterisk is sending

Version 0.9.0, a first release of cg, is published on CRAN and its mirrors. The "cg" name stands for "compare groups." Its genesis and evolution are driven by common "in-practice" needs to compare samples, treatments, administrations, conditions, etc. in medicine research & development. The current version provides comparisons of unpaired samples, i.e. a linear

Version 0.9.0, a first release of cg, is published on CRAN and its mirrors. The "cg" name stands for "compare groups." Its genesis and evolution are driven by common "in-practice" needs to compare samples, treatments, administrations, conditions, etc. in medicine research & development. The current version provides comparisons of unpaired samples, i.e. a linear

Does anybody know if NIC Realtek RTL8151GH-CG is supported on FreeBSD 9.2 or 10.0? It is included in HP ProDesk 490 MT which I may buy, but I cannot find what hardware parts are supported by FreeBSD and what not. I am running FreeBSD for about 13 years on servers. I didn't run it on desktop grade HW for more than 8 years. So any tips for some solid home PC with i7-4770 (or i7-3770?) is

Hi Daniel, On Wed, Aug 12, 2015 at 04:16:25PM +0200, Daniel Vetter wrote: > > * Reprobing if the inactive GPU initializes before the apple-gmux module: > > v1 used Matthew Garrett's approach of adding a driver callback. > > v2 simply generates a hotplug event instead. nouveau polls its outputs > > every 10 seconds so we want it to poll immediately once apple-gmux

The R-1.2.1 builds have gone fairly smoothly on most of my UNIX architectures, but there were many warnings reported by the picky SGI compiler that could be eliminated in future releases. Generally, I have found the SGI compilers very helpful in ferreting out portability problems, unused and/or obsolete variables, unexpected datatype coercions, etc., and they compile faster than any other

… actually, now that I’m able double-check this, I’m quite surprised to find that we didn’t define fmax(+0,–0) in IEEE–754, which says [paraphrased]: minNum(x,y) is x if x < y, y if y < x, and the number if one is a number and the other is NaN. Otherwise, it is either x or y (this means results might differ among implementations). So I think your proposed semantics are perfectly

I was trying to find source for optimize and I ran across function (f, interval, ..., lower = min(interval), upper = max(interval), maximum = FALSE, tol = .Machine$double.eps^0.25) { if (maximum) { val <- .Internal(fmin(function(arg) -f(arg, ...), lower, upper, tol)) list(maximum = val, objective = f(val, ...)) } else { val <-

Hi Carter, I would strongly advise you against this direction. I’m aware of two directions that existing languages go in defining min/max operations: - IEEE 754, C, Fortran, Matlab, OpenCL, and HLSL all define it not to propagate NaNs - C++ (std::min/std::max) and OpenGL define it in the trinary operator manner: (a < b) ? a : b What you’re proposing does not match any existing languages

This is a problem with all floating point folding, not just with these operations. What Matt is proposing is consistent with how we fold other libm intrinsics. —Owen > On Aug 18, 2014, at 1:22 AM, Mueller-Roemer, Johannes Sebastian <Johannes.Sebastian.Mueller-Roemer at igd.fraunhofer.de> wrote: > > Wouldn’t it be better to use the target’s implementation (if there is one)

Hi, I’d like to re-propose adding intrinsics for fmin / fmax. These can be used to implement the equivalent libm functions as defined in C99 and OpenCL, which R600 and AArch64 at least have instructions with the same semantics. This is not equivalent to a simple fcmp + select due to its handling of NaNs. This has been proposed before, but never delivered

Hi, The specification for the llvm.minnum/llvm.maxnum intrinsics is too unclear right now to usefully optimize. There are two problems. First the expected behavior for signaling NaNs needs to be clarified. Second, whether the returned value is expected to be canonicalized (as if by llvm.canonicalize). Currently according to the LangRef: Follows the IEEE-754 semantics for minNum, which also

Hi, Is there something like the Matlab function FMINS in any R packages? Here is the help file of fmins in case something similar exist under R with a different name. Thanks in advance. Yves Gauvreau *********************** function [x, options] = fmins(funfcn,x,options,grad,varargin) %FMINS Minimize function of several variables. % X = FMINS('F',X0) attempts to return a vector X

Dear R Users/Experts, I am using a function called logitreg() originally described in MASS (the book 4th Ed.) by Venebles & Ripley, p445. I used the code as provided but made couple of changes to run a 'constrained' logistic regression, I set the method = "L-BFGS-B", set lower/upper values for the variables. Here is the function, logitregVR <- function(x, y, wt =

Sorry, I deleted my old post. Pasting the new query below. Dear R Users/Experts, I am using a function called logitreg() originally described in MASS (the book 4th Ed.) by Venebles & Ripley, p445. I used the code as provided but made couple of changes to run a 'constrained' logistic regression, I set the method = "L-BFGS-B", set lower/upper values for the variables. Here

Hi, I wanted to use R's lbfgsb method for minimization from C. Unfortunately, my toy examples always crashes (segmentation fault). What's wrong with it? double eval(int n, double* par, void *ex) { double result = 0; for (int i=0; i<n; ++i) { result += par[i]*par[i]; } printf("result=%.2f\n", result); return result; } void grad(int n, double *par, double *gr,

This kind of error seems to surprise R users. It surprises me that it doesn't happen much more frequently. The "BFGS" method of optim() from the 1990 Pascal version of my book was called the Variable Metric method as per Fletcher's 1970 paper it was drawn from. It really works much better with analytic gradients, and the Rvmmin package which is an all-R version that adds bounds