Running 'doRUnit.R' [56s/61s]
Running the tests in 'tests/doRUnit.R' failed.
Complete output:
  > #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  > #+++ call testthat unit tests in directory inst\tests +++
  > #++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  > 
  > if (require("testthat", quietly = TRUE)) {
  + 	pkg <- "REddyProc"
  + 	#library(pkg, character.only = TRUE)
  + 	#test_package(pkg)
  + 	test_check(pkg)
  + } else {
  + 	warning("cannot run unit tests -- package testthat is not available")
  + }
  Loading required package: REddyProc
  -- 1. Error: partGLFitLRCWindows outputs are in accepted range (@test_partGL.R#4
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds, nRecInDay = 48L, controlGLPart = partGLControl()) at testthat/test_partGL.R:460
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 2. Error: partGLFitLRCWindows error on low uncertainty (@test_partGL.R#506)  
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds, nRecInDay = 48L, controlGLPart = partGLControl()) at testthat/test_partGL.R:506
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 3. Error: partitionNEEGL (@test_partGL.R#591)  ------------------------------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: partitionNEEGL(dsNEE1, RadVar = "Rg_f") at testthat/test_partGL.R:591
  2: partGLFitNightTimeTRespSens(dsR, nRecInDay = nRecInDay, controlGLPart = controlGLPart)
  3: partGLSmoothTempSens(resNight)
  4: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  5: FUN(X[[i]], ...)
  6: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  7: evalVis(expr)
  8: withVisible(eval(expr, pf))
  9: eval(expr, pf)
  10: eval(expr, pf)
  11: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  12: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  13: FUN(X[[i]], ...)
  
  -- 4. Error: partitionNEEGL sparse data (@test_partGL.R#706)  ------------------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds, nRecInDay = 48L, controlGLPart = partGLControl()) at testthat/test_partGL.R:706
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 5. Error: partitionNEEGL isNeglectVPDEffect (@test_partGL.R#748)  -----------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds, nRecInDay = 48L, controlGLPart = ctrl) at testthat/test_partGL.R:748
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 6. Error: partGLPartitionFluxes missing night time data (@test_partGL.R#804) 
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds, nRecInDay = 48L, controlGLPart = partGLControl(), 
         winSizeNight = 4L, winExtendSizes = c()) at testthat/test_partGL.R:804
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 7. Error: partGLPartitionFluxes filter Meteo flag not enough without VPD (@te
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: partitionNEEGL(ds, controlGLPart = partGLControl(isFilterMeteoQualityFlag = TRUE)) at testthat/test_partGL.R:834
  2: partGLFitNightTimeTRespSens(dsR, nRecInDay = nRecInDay, controlGLPart = controlGLPart)
  3: partGLSmoothTempSens(resNight)
  4: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  5: FUN(X[[i]], ...)
  6: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  7: evalVis(expr)
  8: withVisible(eval(expr, pf))
  9: eval(expr, pf)
  10: eval(expr, pf)
  11: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  12: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  13: FUN(X[[i]], ...)
  
  -- 8. Error: partGLPartitionFluxes missing prediction VPD (@test_partGL.R#874)  
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: expect_warning(tmp <- partitionNEEGL(ds, controlGLPart = partGLControl(isFilterMeteoQualityFlag = TRUE, 
         isRefitMissingVPDWithNeglectVPDEffect = FALSE))) at testthat/test_partGL.R:874
  2: quasi_capture(enquo(object), label, capture_warnings)
  3: .capture(act$val <- eval_bare(get_expr(.quo), get_env(.quo)), ...)
  4: withCallingHandlers(code, warning = function(condition) {
         out$push(condition)
         invokeRestart("muffleWarning")
     })
  5: eval_bare(get_expr(.quo), get_env(.quo))
  6: partitionNEEGL(ds, controlGLPart = partGLControl(isFilterMeteoQualityFlag = TRUE, 
         isRefitMissingVPDWithNeglectVPDEffect = FALSE))
  7: partGLFitNightTimeTRespSens(dsR, nRecInDay = nRecInDay, controlGLPart = controlGLPart)
  8: partGLSmoothTempSens(resNight)
  9: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  10: FUN(X[[i]], ...)
  11: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  12: evalVis(expr)
  13: withVisible(eval(expr, pf))
  14: eval(expr, pf)
  15: eval(expr, pf)
  16: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  17: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  18: FUN(X[[i]], ...)
  
  -- 9. Error: partitionNEEGL fixed tempSens (@test_partGL.R#912)  ---------------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: REddyProc:::partGLFitNightTimeTRespSens(ds) at testthat/test_partGL.R:912
  2: partGLSmoothTempSens(resNight)
  3: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  4: FUN(X[[i]], ...)
  5: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  6: evalVis(expr)
  7: withVisible(eval(expr, pf))
  8: eval(expr, pf)
  9: eval(expr, pf)
  10: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  11: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  12: FUN(X[[i]], ...)
  
  -- 10. Error: partitionNEEGL: missing sdNEE (@test_partGL.R#940)  --------------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: partitionNEEGL(dsNEE1, RadVar = "Rg_f", controlGLPart = partGLControl()) at testthat/test_partGL.R:940
  2: partGLFitNightTimeTRespSens(dsR, nRecInDay = nRecInDay, controlGLPart = controlGLPart)
  3: partGLSmoothTempSens(resNight)
  4: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  5: FUN(X[[i]], ...)
  6: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  7: evalVis(expr)
  8: withVisible(eval(expr, pf))
  9: eval(expr, pf)
  10: eval(expr, pf)
  11: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  12: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  13: FUN(X[[i]], ...)
  
  -- 11. Error: partitionNEETK (@test_partGL.R#974)  -----------------------------
  BLAS/LAPACK routine 'DPPTRS' gave error code -1
  1: partitionNEEGL(dsNEE1, RadVar = "Rg_f", controlGLPart = partGLControl(useNightimeBasalRespiration = TRUE)) at testthat/test_partGL.R:974
  2: partGLFitNightTimeTRespSens(dsR, nRecInDay = nRecInDay, controlGLPart = controlGLPart)
  3: partGLSmoothTempSens(resNight)
  4: lapply(unique(E0Win$year), function(yr) {
         E0WinYr <- E0Win[E0Win$year == yr, , drop = FALSE]
         isFiniteE0 <- is.finite(E0WinYr$E0)
         E0WinFinite <- E0WinYr[isFiniteE0, ]
         if (!nrow(E0WinFinite)) {
             warning("No respiration-temperature relationship for any period of the ", 
                 yr, "th year. ", "Using mean temperature sensitivity of other years for this year.")
         }
         else {
             output <- capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, 
                 nugget = E0WinFinite$sdE0^2))
             pred1 <- predict(gpFit, matrix(E0WinYr$iCentralRec, ncol = 1), se.fit = TRUE)
             nuggetNewObs <- quantile(gpFit$nugget, 0.9)
             nugget <- rep(nuggetNewObs, nrow(E0WinYr))
             nugget[isFiniteE0] <- gpFit$nugget
             E0WinYr$E0 <- as.vector(pred1$fit)
             E0WinYr$sdE0 <- as.vector(pred1$se.fit) + unname(sqrt(nugget))
         }
         E0WinYr
     })
  5: FUN(X[[i]], ...)
  6: capture.output(gpFit <- mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2))
  7: evalVis(expr)
  8: withVisible(eval(expr, pf))
  9: eval(expr, pf)
  10: eval(expr, pf)
  11: mlegp(X = E0WinFinite$iCentralRec, Z = E0WinFinite$E0, nugget = E0WinFinite$sdE0^2)
  12: lapply(Zlist, mlegp2, XX = X, orig.XX = orig.X, nugget = nugget, nugget.known = nugget.known, 
         reps = reps, un.sum = un.sum, numEstimates = numEstimates, constantMean = constantMean, 
         simplex.ntries = simplex.ntries, simplex.maxiter = simplex.maxiter, simplex.reltol = simplex.reltol, 
         BFGS.maxiter = BFGS.maxiter, BFGS.tol = BFGS.tol, BFGS.h = BFGS.h, min.nugget = min.nugget, 
         verbose = verbose, param.names = param.names)
  13: FUN(X[[i]], ...)
  
  == testthat results  ===========================================================
  OK: 262 SKIPPED: 16 WARNINGS: 2 FAILED: 11
  1. Error: partGLFitLRCWindows outputs are in accepted range (@test_partGL.R#460) 
  2. Error: partGLFitLRCWindows error on low uncertainty (@test_partGL.R#506) 
  3. Error: partitionNEEGL (@test_partGL.R#591) 
  4. Error: partitionNEEGL sparse data (@test_partGL.R#706) 
  5. Error: partitionNEEGL isNeglectVPDEffect (@test_partGL.R#748) 
  6. Error: partGLPartitionFluxes missing night time data (@test_partGL.R#804) 
  7. Error: partGLPartitionFluxes filter Meteo flag not enough without VPD (@test_partGL.R#834) 
  8. Error: partGLPartitionFluxes missing prediction VPD (@test_partGL.R#874) 
  9. Error: partitionNEEGL fixed tempSens (@test_partGL.R#912) 
  1. ...
  
  Error: testthat unit tests failed
  Execution halted