==3469338== Memcheck, a memory error detector ==3469338== Copyright (C) 2002-2024, and GNU GPL'd, by Julian Seward et al. ==3469338== Using Valgrind-3.24.0 and LibVEX; rerun with -h for copyright info ==3469338== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R -f boundary.R --restore --save --no-readline --vanilla ==3469338== R Under development (unstable) (2026-07-02 r90204) -- "Unsuffered Consequences" Copyright (C) 2026 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ## In both of these cases boundary fit (i.e. estimate of zero RE > ## variance) is *incorrect*. (Nelder_Mead, restart_edge=FALSE) is the > ## only case where we get stuck; either optimizer=bobyqa or > ## restart_edge=TRUE (default) works > if (.Platform$OS.type != "windows") { + + library(lme4) + library(testthat) + + if(!dev.interactive(orNone=TRUE)) pdf("boundary_plots.pdf") + + ## Stephane Laurent: + dat <- read.csv(system.file("testdata","dat20101314.csv", package="lme4")) + + fit <- lmer(y ~ (1|Operator)+(1|Part)+(1|Part:Operator), data=dat, + control= lmerControl(optimizer="Nelder_Mead")) + fit_b <- lmer(y ~ (1|Operator)+(1|Part)+(1|Part:Operator), data=dat, + control= lmerControl(optimizer="bobyqa", restart_edge=FALSE)) + fit_c <- lmer(y ~ (1|Operator)+(1|Part)+(1|Part:Operator), data=dat, + control= lmerControl(optimizer="Nelder_Mead", restart_edge=FALSE, + check.conv.hess="ignore")) + ## final fit gives degenerate-Hessian warning + ## FIXME: use fit_c with expect_warning() as a check on convergence tests + ## tolerance=1e-5 seems OK in interactive use but not in R CMD check ... ?? + stopifnot(all.equal(getME(fit, "theta") -> th.f, + getME(fit_b,"theta"), tolerance=5e-5), + all(th.f > 0)) + + ## Manuel Koller + + source(system.file("testdata", "koller-data.R", package="lme4")) + ldata <- getData(13) + ## old (backward compatible/buggy) + fm4 <- lmer(y ~ (1|Var2), ldata, control=lmerControl(optimizer="Nelder_Mead", + use.last.params=TRUE), + start=list(theta=1)) + + fm4b <- lmer(y ~ (1|Var2), ldata, + control = lmerControl(optimizer="Nelder_Mead", use.last.params=TRUE, + restart_edge = FALSE, + check.conv.hess="ignore", check.conv.grad="ignore"), + start = list(theta=1)) + ## FIXME: use as convergence test check + stopifnot(getME(fm4b,"theta") == 0) + fm4c <- lmer(y ~ (1|Var2), ldata, control=lmerControl(optimizer="bobyqa", + use.last.params=TRUE), + start=list(theta=1)) + stopifnot(all.equal(getME(fm4, "theta") -> th4, + getME(fm4c,"theta"), tolerance=1e-4), + th4 > 0) + + + ## new: doesn't get stuck at edge any more, but gets stuck somewhere else ... + fm5 <- lmer(y ~ (1|Var2), ldata, control=lmerControl(optimizer="Nelder_Mead", + check.conv.hess="ignore", + check.conv.grad="ignore"), + start=list(theta=1)) + fm5b <- lmer(y ~ (1|Var2), ldata, control=lmerControl(optimizer="Nelder_Mead", + restart_edge=FALSE, + check.conv.hess="ignore", + check.conv.grad="ignore"), + start = list(theta = 1)) + fm5c <- lmer(y ~ (1|Var2), ldata, control=lmerControl(optimizer="bobyqa"), + start = list(theta = 1)) + stopifnot(all.equal(unname(getME(fm5c,"theta")), 0.21067645, tolerance = 1e-7)) + # 0.21067644264 [64-bit, lynne] + + if (require("optimx")) { + ## additional stuff for diagnosing Nelder-Mead problems. + + fm5d <- update(fm5,control=lmerControl(optimizer="optimx", + optCtrl=list(method="L-BFGS-B"))) + + fm5e <- update(fm5, control=lmerControl(optimizer="nloptwrap")) + + mList <- setNames(list(fm4,fm4b,fm4c,fm5,fm5b,fm5c,fm5d,fm5e), + c("NM/uselast","NM/uselast/norestart","bobyqa/uselast", + "NM","NM/norestart","bobyqa","LBFGSB","nloptr/bobyqa")) + pp <- profile(fm5c,which=1) + dd <- as.data.frame(pp) + par(las=1,bty="l") + v <- sapply(mList, + function(x) sqrt(VarCorr(x)[[1]])) + plot(.zeta^2~.sig01, data=dd, type="b") + abline(v=v) + + res <- cbind(VCorr = sapply(mList, function(x) sqrt(VarCorr(x)[[1]])), + theta = sapply(mList, getME,"theta"), + loglik = sapply(mList, logLik)) + res + print(sessionInfo(), locale=FALSE) + } + + ###################### + library(lattice) + ## testing boundary and near-boundary cases + + tmpf <- function(i,...) { + set.seed(i) + d <- data.frame(x=rnorm(60),f=factor(rep(1:6,each=10))) + d$y <- simulate(~x+(1|f),family=gaussian,newdata=d, + newparams=list(theta=0.01,beta=c(1,1),sigma=5))[[1]] + lmer(y~x+(1|f),data=d,...) + } + sumf <- function(m) { + unlist(VarCorr(m))[1] + } + if (FALSE) { + ## figuring out which seeds will give boundary and + ## near-boundary solutions + mList <- lapply(1:201,tmpf) # [FIXME tons of messages "theta parameters vector not named"] + ss <- sapply(mList,sumf)+1e-50 + par(las=1,bty="l") + hist(log(ss),col="gray",breaks=50) + ## values lying on boundary + which(log(ss)<(-40)) ## 5, 7-13, 15, 21, ... + ## values close to boundary (if check.edge not set) + which(log(ss)>(-40) & log(ss) <(-20)) ## 16, 44, 80, 86, 116, ... + } + ## diagnostic plot + tmpplot <- function(i, FUN=tmpf) { + dd <- FUN(i, devFunOnly=TRUE) + x <- 10^seq(-10,-6.5,length=201) + dvec <- sapply(x,dd) + op <- par(las=1,bty="l"); on.exit(par(op)) + plot(x,dvec-min(dvec)+1e-16, log="xy", type="b") + r <- FUN(i) + abline(v = getME(r,"theta"), col=2) + invisible(r) + } + + ## Case #1: boundary estimate with or without boundary.tol + m5 <- tmpf(5) + m5B <- tmpf(5,control=lmerControl(boundary.tol=0)) + stopifnot(getME(m5, "theta")==0, + getME(m5B,"theta")==0) + p5 <- profile(m5) ## bobyqa warnings but results look reasonable + xyplot(p5) + ## reveals slight glitch (bottom row of plots doesn't look right) + expect_warning(splom(p5),"unreliable for singular fits") + p5B <- profile(m5, signames=FALSE) # -> bobyqa convergence warning (code 3) + expect_warning(splom(p5B), "unreliable for singular fits") + + if(lme4:::testLevel() >= 2) { ## avoid failure to warn + ## Case #2: near-boundary estimate, but boundary.tol can't fix it + m16 <- tmpplot(16) + ## sometimes[2014-11-11] fails (??) : + p16 <- profile(m16) ## warning message*s* (non-monotonic profile and more) + plotOb <- xyplot(p16) + ## NB: It's the print()ing of 'plotOb' which warns ==> need to do this explicitly: + expect_warning(print(plotOb), ## warns about linear interpolation in profile for variable 1 + "using linear interpolation") + d16 <- as.data.frame(p16) + xyplot(.zeta ~ .focal|.par, data=d16, type=c("p","l"), + scales = list(x=list(relation="free"))) + try(splom(p16)) ## breaks when calling predict(.) + } + + ## bottom line: + ## * xyplot.thpr could still be improved + ## * most of the near-boundary cases are noisy and can't easily be + ## fixed + + tmpf2 <- function(i,...) { + set.seed(i) + d <- data.frame(x=rnorm(60),f=factor(rep(1:6,each=10)), + w=rep(10,60)) + d$y <- simulate(~x+(1|f),family=binomial, + weights=d$w,newdata=d, + newparams=list(theta=0.01,beta=c(1,1)))[[1]] + glmer(y~x+(1|f),data=d,family=binomial,weights=w,...) + } + + if (FALSE) { + ## figuring out which seeds will give boundary and + ## near-boundary solutions + mList <- lapply(1:201,tmpf2) + ss <- sapply(mList,sumf)+1e-50 + par(las=1,bty="l") + hist(log(ss),col="gray",breaks=50) + ## values lying on boundary + head(which(log(ss)<(-50))) ## 1-5, 7 ... + ## values close to boundary (if check.edge not set) + which(log(ss)>(-50) & log(ss) <(-20)) ## 44, 46, 52, ... + } + + ## m1 <- tmpf2(1) + + ## FIXME: doesn't work if we generate m1 via tmpf2(1) -- + ## some environment lookup problem ... + + set.seed(1) + d <- data.frame(x=rnorm(60),f=factor(rep(1:6,each=10)), + w=rep(10,60)) + d$y <- simulate(~x+(1|f),family=binomial, + weights=d$w,newdata=d, + newparams=list(theta=0.01,beta=c(1,1)))[[1]] + m1 <- glmer(y~x+(1|f),data=d,family=binomial,weights=w) + + p1 <- profile(m1) + xyplot(p1) + expect_warning(splom(p1),"splom is unreliable") + + } ## skip on windows (for speed) Loading required package: Matrix boundary (singular) fit: see help('isSingular') Loading required package: optimx R Under development (unstable) (2026-07-02 r90204) Platform: x86_64-pc-linux-gnu Running under: Fedora Linux 40 (Workstation Edition) Matrix products: default BLAS: /data/blackswan/ripley/R/R-devel-vg/lib/libRblas.so LAPACK: /usr/lib64/liblapack.so.3.12.0 LAPACK version 3.12.0 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] optimx_2025-4.9 testthat_3.3.2 lme4_2.0-1 Matrix_1.7-5 loaded via a namespace (and not attached): [1] R6_2.6.1 numDeriv_2016.8-1.1 lattice_0.22-9 [4] magrittr_2.0.5 splines_4.7.0 lifecycle_1.0.5 [7] Rdpack_2.6.6 cli_3.6.6 nloptr_2.2.1 [10] grid_4.7.0 reformulas_0.4.4 compiler_4.7.0 [13] boot_1.3-32 rbibutils_2.4.1 tools_4.7.0 [16] pracma_2.4.6 brio_1.1.5 nlme_3.1-169 [19] minqa_1.2.8 Rcpp_1.1.2 otel_0.2.0 [22] rlang_1.3.0 MASS_7.3-65 boundary (singular) fit: see help('isSingular') boundary (singular) fit: see help('isSingular') ==3469338== Invalid read of size 8 ==3469338== at 0x247DA545: coeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925) ==3469338== by 0x247DA545: assignCoeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654) ==3469338== by 0x247DA545: assignCoeffByOuterInner (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668) ==3469338== by 0x247DA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347) ==3469338== by 0x247DA545: call_dense_assignment_loop >, Eigen::Product > >, Eigen::Map >, 1>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785) ==3469338== by 0x247DA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954) ==3469338== by 0x247DA545: call_assignment_no_alias >, Eigen::Product > >, Eigen::Map >, 1>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DA545: call_assignment >, Eigen::Product > >, Eigen::Map >, 1>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858) ==3469338== by 0x247DA545: call_assignment >, Eigen::Product > >, Eigen::Map >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836) ==3469338== by 0x247DA545: operator= > >, Eigen::Map >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Assign.h:66) ==3469338== by 0x247DA545: lme4::merPredD::updateXwts(Eigen::Array const&) (packages/tests-vg/lme4/src/predModule.cpp:221) ==3469338== by 0x247B2A2F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:274) ==3469338== by 0x247B2A2F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Address 0xfdcad60 is 48 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D31F4: ploadu (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603) ==3469338== by 0x247D31F4: ploadt (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969) ==3469338== by 0x247D31F4: load (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203) ==3469338== by 0x247D31F4: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcad60 is 48 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D3204: ploadu (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603) ==3469338== by 0x247D3204: ploadt (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969) ==3469338== by 0x247D3204: load (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203) ==3469338== by 0x247D3204: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcad70 is 64 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D3211: ploadu (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603) ==3469338== by 0x247D3211: ploadt (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969) ==3469338== by 0x247D3211: load (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203) ==3469338== by 0x247D3211: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcad80 is 80 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D321E: pmul (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237) ==3469338== by 0x247D321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959) ==3469338== by 0x247D321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95) ==3469338== by 0x247D321E: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcad90 is 96 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D360C: ploadu (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603) ==3469338== by 0x247D360C: ploadt (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969) ==3469338== by 0x247D360C: load (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203) ==3469338== by 0x247D360C: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:201) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcaf20 is 496 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== ==3469338== Invalid read of size 16 ==3469338== at 0x247D361C: pmul (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237) ==3469338== by 0x247D361C: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959) ==3469338== by 0x247D361C: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95) ==3469338== by 0x247D361C: Eigen::internal::general_matrix_vector_product, 0, false, double, Eigen::internal::const_blas_data_mapper, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper const&, Eigen::internal::const_blas_data_mapper const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:201) ==3469338== by 0x247E7D5B: run >, Eigen::Matrix, Eigen::Matrix > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253) ==3469338== by 0x247E7D5B: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388) ==3469338== by 0x247E7D5B: void Eigen::internal::generic_product_impl, 0, Eigen::Stride<0, 0> >, Eigen::Matrix, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo >(Eigen::Matrix&, Eigen::Map, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376) ==3469338== by 0x247DD3BE: scaleAndAddTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361) ==3469338== by 0x247DD3BE: evalTo > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349) ==3469338== by 0x247DD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::Product >, Eigen::Matrix, 0>, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: run, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223) ==3469338== by 0x247DD3BE: call_assignment_no_alias, Eigen::CwiseBinaryOp, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> >, Eigen::internal::assign_op > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890) ==3469338== by 0x247DD3BE: _set_noalias, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797) ==3469338== by 0x247DD3BE: PlainObjectBase, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594) ==3469338== by 0x247DD3BE: Matrix, const Eigen::Product >, Eigen::Matrix, 0>, const Eigen::Product > >, Eigen::Matrix, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423) ==3469338== by 0x247DD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95) ==3469338== by 0x247B286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290) ==3469338== by 0x247B286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330) ==3469338== by 0x247B8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386) ==3469338== by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770) ==3469338== by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278) ==3469338== Address 0xfdcaf30 is 512 bytes inside a block of size 1,008 free'd ==3469338== at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989) ==3469338== by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167) ==3469338== by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951) ==3469338== by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237) ==3469338== by 0x537ADE: Rf_cons (svn/R-devel/src/main/memory.c:2488) ==3469338== by 0x4C042F: Rf_lcons (svn/R-devel/src/include/Rinlinedfuns.h:717) ==3469338== by 0x4C042F: getActiveValue (svn/R-devel/src/main/envir.c:166) ==3469338== by 0x4C54AF: R_findVarInFrame (svn/R-devel/src/main/envir.c:1249) ==3469338== by 0x4C54AF: R_findVar (svn/R-devel/src/main/envir.c:1387) ==3469338== by 0x4F4844: Rf_eval (svn/R-devel/src/main/eval.c:1175) ==3469338== by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683) ==3469338== by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249) ==3469338== by 0x4F7D61: do_if (svn/R-devel/src/main/eval.c:2709) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001) ==3469338== by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230) ==3469338== Block was alloc'd at ==3469338== at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446) ==3469338== by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894) ==3469338== by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609) ==3469338== by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236) ==3469338== by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676) ==3469338== by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573) ==3469338== by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533) ==3469338== by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518) ==3469338== by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167) ==3469338== by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398) ==3469338== by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314) ==3469338== by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333) ==3469338== by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473) ==3469338== boundary (singular) fit: see help('isSingular') > > proc.time() user system elapsed 680.665 3.394 686.598 ==3469338== ==3469338== HEAP SUMMARY: ==3469338== in use at exit: 210,803,133 bytes in 47,286 blocks ==3469338== total heap usage: 1,272,039 allocs, 1,224,753 frees, 779,235,625 bytes allocated ==3469338== ==3469338== LEAK SUMMARY: ==3469338== definitely lost: 0 bytes in 0 blocks ==3469338== indirectly lost: 0 bytes in 0 blocks ==3469338== possibly lost: 0 bytes in 0 blocks ==3469338== still reachable: 210,800,585 bytes in 47,285 blocks ==3469338== suppressed: 2,548 bytes in 1 blocks ==3469338== Reachable blocks (those to which a pointer was found) are not shown. ==3469338== To see them, rerun with: --leak-check=full --show-leak-kinds=all ==3469338== ==3469338== For lists of detected and suppressed errors, rerun with: -s ==3469338== ERROR SUMMARY: 21960 errors from 7 contexts (suppressed: 0 from 0)