==3459103== Memcheck, a memory error detector
==3459103== Copyright (C) 2002-2024, and GNU GPL'd, by Julian Seward et al.
==3459103== Using Valgrind-3.24.0 and LibVEX; rerun with -h for copyright info
==3459103== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R -f agridat_gotway.R --restore --save --no-readline --vanilla
==3459103== 

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.

> ## require(agridat)
> ## dat <- gotway.hessianfly
> 
> if (.Platform$OS.type != "windows") {
+ ## don't actually use gotway_hessianfly_fit or gotway_hessianfly_prof,
+ ## so we should be OK even with R< 3.0.1
+ load(system.file("testdata","gotway_hessianfly.rda",package="lme4"))
+ # Block random.  See Glimmix manual, output 1.18.
+ # Note: (Different parameterization)
+ 
+ ## require("lme4.0")
+ ## fit2 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=dat, family=binomial)
+ ## params <- list(fixef=fixef(fit2),theta=getME(fit2,"theta"))
+ ## detach("package:lme4.0")
+ lme4.0fit <- structure(list(fixef = structure(c(1.50345713031203, -0.193853259383803, 
+ -0.540808391060274, -1.43419379979154, -0.203701042949808, -0.978322555343941, 
+ -0.604078624475678, -1.67742449813309, -1.39842466673692, -0.681709344788684, 
+ -1.46295367186169, -1.45908310198959, -3.55285756517073, -2.50731975980307, 
+ -2.08716296677356, -2.96974270029992), .Names = c("(Intercept)", 
+ "genG02", "genG03", "genG04", "genG05", "genG06", "genG07", "genG08", 
+ "genG09", "genG10", "genG11", "genG12", "genG13", "genG14", "genG15", 
+ "genG16")), theta = structure(0.0319087494293615, .Names = "block.(Intercept)")), .Names = c("fixef", 
+ "theta"))
+ 
+ ## start doesn't work because we don't get there
+ library(lme4)
+ m1 <- glmer(cbind(y, n-y) ~ gen + (1|block), data=gotway.hessianfly,
+             family=binomial)
+ m1B <- update(m1,control=glmerControl(optimizer="bobyqa"))
+ max(abs(m1@optinfo$derivs$gradient))  ## 0.0012
+ max(abs(m1B@optinfo$derivs$gradient)) ## 2.03e-5
+ abs(m1@optinfo$derivs$gradient)/abs(m1B@optinfo$derivs$gradient)
+ ## bobyqa gets gradients *at least* 1.64* lower
+ 
+ lme4fit <- list(fixef=fixef(m1),theta=getME(m1,"theta"))
+ 
+ ## hack around slight naming differences
+ lme4fit$theta <- unname(lme4fit$theta)
+ lme4.0fit$theta <- unname(lme4.0fit$theta)
+ ## difference in theta on x86_64-w64-mingw32 (64-bit) with r-devel is  0.000469576
+ stopifnot(all.equal(lme4fit, lme4.0fit, tolerance = 5e-4))
+ 
+ ## Fun stuff: visualize and alternative model
+ 
+ ## library(ggplot2)
+ ## dat$prop <- dat$y/dat$n
+ ## theme_set(theme_bw())
+ ## ggplot(dat,aes(x=gen,y=prop,colour=block))+geom_point(aes(size=n))+
+ ##     geom_line(aes(group=block,colour=block))+
+ ##     geom_smooth(family=binomial,aes(weight=n,colour=block,group=block),method="glm",
+ ##                 alpha=0.1)
+ 
+ ## dat$obs <- factor(seq(nrow(dat)))
+ ## m2 <- glmer(cbind(y, n-y) ~ block+ (1|gen) + (1|obs), data=dat, family=binomial)
+ 
+ }  ## not on windows/CRAN
Loading required package: Matrix
==3459103== Invalid read of size 8
==3459103==    at 0x24FDA545: coeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:925)
==3459103==    by 0x24FDA545: assignCoeff (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:654)
==3459103==    by 0x24FDA545: assignCoeffByOuterInner (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:668)
==3459103==    by 0x24FDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:347)
==3459103==    by 0x24FDA545: call_dense_assignment_loop<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:785)
==3459103==    by 0x24FDA545: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:954)
==3459103==    by 0x24FDA545: call_assignment_no_alias<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:858)
==3459103==    by 0x24FDA545: call_assignment<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:836)
==3459103==    by 0x24FDA545: operator=<Eigen::Product<Eigen::DiagonalWrapper<const Eigen::Map<Eigen::Matrix<double, -1, 1> > >, Eigen::Map<Eigen::Matrix<double, -1, -1> >, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Assign.h:66)
==3459103==    by 0x24FDA545: lme4::merPredD::updateXwts(Eigen::Array<double, -1, 1, 0, -1, 1> const&) (packages/tests-vg/lme4/src/predModule.cpp:221)
==3459103==    by 0x24FB2A2F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:274)
==3459103==    by 0x24FB2A2F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3459103==    by 0x24FB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3459103==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3459103==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==    by 0x4F96BC: do_set (svn/R-devel/src/main/eval.c:3585)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==  Address 0x1c44b290 is 48 bytes inside a block of size 8,240 free'd
==3459103==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3459103==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3459103==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3459103==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3459103==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3459103==    by 0x4F9870: Rf_evalList (svn/R-devel/src/main/eval.c:3685)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Block was alloc'd at
==3459103==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3459103==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3459103==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3459103==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3459103==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3459103==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3459103==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3459103== 
==3459103== Invalid read of size 16
==3459103==    at 0x24FD31F4: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3459103==    by 0x24FD31F4: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3459103==    by 0x24FD31F4: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3459103==    by 0x24FD31F4: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3459103==    by 0x24FE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3459103==    by 0x24FE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3459103==    by 0x24FE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3459103==    by 0x24FDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3459103==    by 0x24FDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3459103==    by 0x24FDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3459103==    by 0x24FDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3459103==    by 0x24FDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3459103==    by 0x24FDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3459103==    by 0x24FB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3459103==    by 0x24FB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3459103==    by 0x24FB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3459103==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3459103==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Address 0x1c44b290 is 48 bytes inside a block of size 8,240 free'd
==3459103==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3459103==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3459103==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3459103==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3459103==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3459103==    by 0x4F9870: Rf_evalList (svn/R-devel/src/main/eval.c:3685)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Block was alloc'd at
==3459103==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3459103==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3459103==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3459103==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3459103==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3459103==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3459103==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3459103== 
==3459103== Invalid read of size 16
==3459103==    at 0x24FD3204: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3459103==    by 0x24FD3204: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3459103==    by 0x24FD3204: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3459103==    by 0x24FD3204: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3459103==    by 0x24FE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3459103==    by 0x24FE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3459103==    by 0x24FE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3459103==    by 0x24FDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3459103==    by 0x24FDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3459103==    by 0x24FDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3459103==    by 0x24FDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3459103==    by 0x24FDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3459103==    by 0x24FDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3459103==    by 0x24FB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3459103==    by 0x24FB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3459103==    by 0x24FB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3459103==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3459103==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Address 0x1c44b2a0 is 64 bytes inside a block of size 8,240 free'd
==3459103==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3459103==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3459103==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3459103==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3459103==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3459103==    by 0x4F9870: Rf_evalList (svn/R-devel/src/main/eval.c:3685)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Block was alloc'd at
==3459103==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3459103==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3459103==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3459103==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3459103==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3459103==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3459103==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3459103== 
==3459103== Invalid read of size 16
==3459103==    at 0x24FD3211: ploadu<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:603)
==3459103==    by 0x24FD3211: ploadt<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:969)
==3459103==    by 0x24FD3211: load<double, 0> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/util/BlasUtil.h:203)
==3459103==    by 0x24FD3211: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3459103==    by 0x24FE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3459103==    by 0x24FE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3459103==    by 0x24FE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3459103==    by 0x24FDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3459103==    by 0x24FDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3459103==    by 0x24FDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3459103==    by 0x24FDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3459103==    by 0x24FDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3459103==    by 0x24FDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3459103==    by 0x24FB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3459103==    by 0x24FB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3459103==    by 0x24FB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3459103==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3459103==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Address 0x1c44b2b0 is 80 bytes inside a block of size 8,240 free'd
==3459103==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3459103==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3459103==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3459103==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3459103==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3459103==    by 0x4F9870: Rf_evalList (svn/R-devel/src/main/eval.c:3685)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Block was alloc'd at
==3459103==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3459103==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3459103==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3459103==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3459103==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3459103==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3459103==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3459103== 
==3459103== Invalid read of size 16
==3459103==    at 0x24FD321E: pmul<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:237)
==3459103==    by 0x24FD321E: pmadd<double> (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GenericPacketMath.h:959)
==3459103==    by 0x24FD321E: pmadd (R-devel/site-library/RcppEigen/include/Eigen/src/Core/arch/Default/ConjHelper.h:95)
==3459103==    by 0x24FD321E: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double) [clone .isra.0] (R-devel/site-library/RcppEigen/include/Eigen/src/Core/products/GeneralMatrixVector.h:177)
==3459103==    by 0x24FE7D5B: run<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/GeneralProduct.h:253)
==3459103==    by 0x24FE7D5B: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:388)
==3459103==    by 0x24FE7D5B: void Eigen::internal::generic_product_impl<Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> >, Eigen::Matrix<double, -1, 1, 0, -1, 1>, Eigen::DenseShape, Eigen::DenseShape, 7>::scaleAndAddTo<Eigen::Matrix<double, -1, 1, 0, -1, 1> >(Eigen::Matrix<double, -1, 1, 0, -1, 1>&, Eigen::Map<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 0, Eigen::Stride<0, 0> > const&, Eigen::Matrix<double, -1, 1, 0, -1, 1> const&, double const&) (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:376)
==3459103==    by 0x24FDD3BE: scaleAndAddTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:361)
==3459103==    by 0x24FDD3BE: evalTo<Eigen::Matrix<double, -1, 1> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:349)
==3459103==    by 0x24FDD3BE: run (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:148)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: run<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/ProductEvaluators.h:223)
==3459103==    by 0x24FDD3BE: call_assignment_no_alias<Eigen::Matrix<double, -1, 1>, Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> >, Eigen::internal::assign_op<double, double> > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/AssignEvaluator.h:890)
==3459103==    by 0x24FDD3BE: _set_noalias<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:797)
==3459103==    by 0x24FDD3BE: PlainObjectBase<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/PlainObjectBase.h:594)
==3459103==    by 0x24FDD3BE: Matrix<Eigen::CwiseBinaryOp<Eigen::internal::scalar_sum_op<double, double>, const Eigen::Product<Eigen::Map<Eigen::Matrix<double, -1, -1> >, Eigen::Matrix<double, -1, 1>, 0>, const Eigen::Product<Eigen::Transpose<const Eigen::Map<Eigen::SparseMatrix<double, 0, int> > >, Eigen::Matrix<double, -1, 1>, 0> > > (R-devel/site-library/RcppEigen/include/Eigen/src/Core/Matrix.h:423)
==3459103==    by 0x24FDD3BE: lme4::merPredD::linPred(double const&) const (packages/tests-vg/lme4/src/predModule.cpp:95)
==3459103==    by 0x24FB286F: internal_glmerWrkIter (packages/tests-vg/lme4/src/external.cpp:290)
==3459103==    by 0x24FB286F: pwrssUpdate (packages/tests-vg/lme4/src/external.cpp:330)
==3459103==    by 0x24FB8E6C: glmerLaplace (packages/tests-vg/lme4/src/external.cpp:386)
==3459103==    by 0x4A91E9: R_doDotCall (svn/R-devel/src/main/dotcode.c:770)
==3459103==    by 0x4E4283: bcEval_loop (svn/R-devel/src/main/eval.c:8700)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Address 0x1c44b2c0 is 96 bytes inside a block of size 8,240 free'd
==3459103==    at 0x4846B83: free (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:989)
==3459103==    by 0x5354C5: ReleaseLargeFreeVectors (svn/R-devel/src/main/memory.c:1167)
==3459103==    by 0x5354C5: RunGenCollect (svn/R-devel/src/main/memory.c:1951)
==3459103==    by 0x5354C5: R_gc_internal (svn/R-devel/src/main/memory.c:3237)
==3459103==    by 0x53830B: CONS_NR (svn/R-devel/src/main/memory.c:2519)
==3459103==    by 0x4F9870: Rf_evalList (svn/R-devel/src/main/eval.c:3685)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F9850: Rf_evalList (svn/R-devel/src/main/eval.c:3683)
==3459103==    by 0x4F4A95: Rf_eval (svn/R-devel/src/main/eval.c:1249)
==3459103==    by 0x4F80E4: do_begin (svn/R-devel/src/main/eval.c:3001)
==3459103==    by 0x4F49B2: Rf_eval (svn/R-devel/src/main/eval.c:1230)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F475A: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x4F475A: Rf_eval (svn/R-devel/src/main/eval.c:1278)
==3459103==  Block was alloc'd at
==3459103==    at 0x4843866: malloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:446)
==3459103==    by 0x539581: Rf_allocVector3 (svn/R-devel/src/main/memory.c:2894)
==3459103==    by 0x44084B: Rf_allocVector (svn/R-devel/src/include/Rinlinedfuns.h:609)
==3459103==    by 0x44084B: Rf_allocMatrix (svn/R-devel/src/main/array.c:236)
==3459103==    by 0x162E0DA3: modelmatrix (svn/R-devel/src/library/stats/src/model.c:676)
==3459103==    by 0x4A76CE: do_External (svn/R-devel/src/main/dotcode.c:573)
==3459103==    by 0x4E63CC: bcEval_loop (svn/R-devel/src/main/eval.c:8150)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103==    by 0x4F729F: applyClosure_core (svn/R-devel/src/main/eval.c:2314)
==3459103==    by 0x4F7D21: Rf_applyClosure (svn/R-devel/src/main/eval.c:2333)
==3459103==    by 0x53FD2B: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==3459103== 
Error: Downdated VtV is not positive definite
Execution halted
==3459103== 
==3459103== HEAP SUMMARY:
==3459103==     in use at exit: 209,898,119 bytes in 44,796 blocks
==3459103==   total heap usage: 240,672 allocs, 195,876 frees, 394,460,329 bytes allocated
==3459103== 
==3459103== 416 bytes in 1 blocks are possibly lost in loss record 232 of 2,589
==3459103==    at 0x484B133: calloc (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:1675)
==3459103==    by 0x4011F63: UnknownInlinedFun (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../include/rtld-malloc.h:44)
==3459103==    by 0x4011F63: allocate_dtv (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../elf/dl-tls.c:395)
==3459103==    by 0x4012A61: _dl_allocate_tls (/usr/src/debug/glibc-2.39-38.fc40.x86_64/elf/../elf/dl-tls.c:673)
==3459103==    by 0x557CC03: allocate_stack (/usr/src/debug/glibc-2.39-38.fc40.x86_64/nptl/allocatestack.c:431)
==3459103==    by 0x557CC03: pthread_create@@GLIBC_2.34 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/nptl/pthread_create.c:660)
==3459103==    by 0x20968DB6: cli__start_thread (/tmp/RtmpOgXr6R/R.INSTALL3b00f1391d945e/cli/src/thread.c:54)
==3459103==    by 0x20968DF7: clic_start_thread (/tmp/RtmpOgXr6R/R.INSTALL3b00f1391d945e/cli/src/thread.c:78)
==3459103==    by 0x4A924D: R_doDotCall (svn/R-devel/src/main/dotcode.c:760)
==3459103==    by 0x4A97A3: do_dotcall (svn/R-devel/src/main/dotcode.c:1437)
==3459103==    by 0x4EEA31: bcEval_loop (svn/R-devel/src/main/eval.c:8085)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7533)
==3459103==    by 0x4F4317: bcEval (svn/R-devel/src/main/eval.c:7518)
==3459103==    by 0x4F464A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==3459103==    by 0x4F65DD: R_execClosure (svn/R-devel/src/main/eval.c:2398)
==3459103== 
==3459103== LEAK SUMMARY:
==3459103==    definitely lost: 0 bytes in 0 blocks
==3459103==    indirectly lost: 0 bytes in 0 blocks
==3459103==      possibly lost: 416 bytes in 1 blocks
==3459103==    still reachable: 209,897,703 bytes in 44,795 blocks
==3459103==         suppressed: 0 bytes in 0 blocks
==3459103== Reachable blocks (those to which a pointer was found) are not shown.
==3459103== To see them, rerun with: --leak-check=full --show-leak-kinds=all
==3459103== 
==3459103== For lists of detected and suppressed errors, rerun with: -s
==3459103== ERROR SUMMARY: 1844068 errors from 6 contexts (suppressed: 0 from 0)