diff --git a/make/_make.R b/make/_make.R
index d07615a..efe990c 100755
--- a/make/_make.R
+++ b/make/_make.R
@@ -97,3 +97,6 @@ run_commit_push("docs_bayescomp-poster")
 run_commit_push("docs_bioinference-poster")
 run_commit_push("plot-tikz_algorithm-flowchart")
 run_commit_push("plot-tikz_simplified-naomi")
+
+#' Checks
+run_commit_push("check_sd-estimation")
diff --git a/src/check_sd-estimation/orderly.yml b/src/check_sd-estimation/orderly.yml
index f9786e8..520b39f 100755
--- a/src/check_sd-estimation/orderly.yml
+++ b/src/check_sd-estimation/orderly.yml
@@ -8,6 +8,9 @@ artefacts:
         - mode-locations.pdf
         - reduced-hessian-draws.pdf
 
+resources:
+  - naomi_simple.cpp
+
 packages:
   - dplyr
   - forcats
diff --git a/src/check_sd-estimation/script.R b/src/check_sd-estimation/script.R
index 4c2c92a..264a18e 100755
--- a/src/check_sd-estimation/script.R
+++ b/src/check_sd-estimation/script.R
@@ -80,6 +80,7 @@ aghq_cov <- solve(aghq_hess)
 tmb_mode <- tmb$fit$par.full
 j2 <- which(names(tmb_mode) == par)[i]
 
+TMB::compile("naomi_simple.cpp")
 dyn.load(TMB::dynlib("naomi_simple"))
 tmb_hess <- naomi:::sdreport_joint_precision(tmb$fit$obj, tmb$fit$par.fixed)
 dim(tmb_hess) #' 491 x 491