when wkr_size_sd is small, NaN introduced in model matrix and lambda() fails about bumbl HOT 2 CLOSED

foraging_index <- match(FALSE, cum_prop < prop_nforaging)
if (is.na(foraging_index)) { #i.e., no foraging workers
    wkr_larv[] <- 0
  } else if (foraging_index > 1) {
    wkr_larv[1:foraging_index - 1] <- 0
    border_case_correction <-
      1 - (prop_nforaging - max(cum_prop[1:foraging_index - 1])) / prop_wkr_size[foraging_index]
  } else if (foraging_index == 1) { #i.e only part of the smallest size class doesn't forage
    border_case_correction <-
      1 - (prop_nforaging / prop_wkr_size[1]) #if prop_foraging = 1 this equals 1
  }

  if (!is.na(foraging_index)) {
    wkr_larv[foraging_index] <-
      wkr_larv[foraging_index] * border_case_correction
  }