R Under development (unstable) (2025-04-18 r88159 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

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.

> library(unittest)
> 
> library(gadget3)
> 
> st_vbsimple <- g3_stock('st_vbsimple', seq(50, 100, by = 10)) |> g3s_age(1,5)
> st_multspec <- g3_stock('st_multspec', seq(50, 100, by = 10)) |> g3s_age(1,5)
> st_weightjones <- g3_stock('st_weightjones', seq(50, 100, by = 10)) |> g3s_age(1,5)
> temperature <- g3_timeareadata(
+     'temp',
+     data.frame(year = 2000, step=c(1,2), temp=c(10, 14)),
+     value_field = "temp" )
> 
> dummy_feeding_action <- function (stock, level_f = g3_parameterized('feedinglevel', value = 1, by_year = TRUE, by_step = TRUE)) {
+     action_name <- gadget3:::unique_action_name()
+ 
+     stock__feedinglevel <- g3_stock_instance(stock, desc = "Fraction of the available food that the predator is consuming")
+     out <- list()
+     out[[gadget3:::step_id(1, stock, action_name)]] <- gadget3:::g3_step(gadget3:::f_substitute(~stock_iterate(stock, {
+         stock_ss(stock__feedinglevel)[] <- 0.001 * stock__midlen + level_f
+     }), list(
+         end = NULL )))
+     return(out)
+ }
> 
> actions <- list(
+     g3a_time(2000, 2000, step_lengths = c(3, 9)),  # NB: Second step longer than first
+ 
+     g3a_otherfood(st_vbsimple, 1e6, 1e3),
+     g3a_growmature(st_vbsimple, g3a_grow_impl_bbinom(
+         g3a_grow_lengthvbsimple(),
+         g3a_grow_weightsimple(),
+         maxlengthgroupgrowth = 8 )),
+ 
+     g3a_otherfood(st_multspec, 1e6, 1e3),
+     dummy_feeding_action(st_multspec),
+     g3a_growmature(st_multspec, g3a_grow_impl_bbinom(
+         g3a_grow_length_multspec(temperature = temperature),
+         g3a_grow_weight_multspec(temperature = temperature),
+         maxlengthgroupgrowth = 8 )),
+     
+     g3a_otherfood(st_weightjones, 1e6, 1e3),
+     dummy_feeding_action(st_weightjones),
+     g3a_growmature(st_weightjones, g3a_grow_impl_bbinom(
+         g3a_grow_length_weightjones(temperature = temperature),
+         g3a_grow_weight_weightjones(temperature = temperature),
+         maxlengthgroupgrowth = 8 )),
+ 
+     # NB: Dummy parameter so model will compile in TMB
+     ~{nll <- nll + g3_param("x", value = 0)} )
> actions <- c(actions, list(
+     g3a_report_history(actions, var_re = "__num$|__wgt$", out_prefix = "late_"),
+     g3a_report_detail(actions) ))
> model_fn <- g3_to_r(actions)
> model_cpp <- g3_to_tmb(actions)
> 
> ok_group("Default params") ########
# Default params
> attr(model_fn, 'parameter_template') |>
+     g3_init_val("*.Linf", 180) |>
+     g3_init_val("*.K", 0.1) |>
+     g3_init_val("*.walpha", 0.01, optimise = FALSE) |>
+     g3_init_val("*.wbeta", 3, optimise = FALSE) |>
+     g3_init_val("st_weightjones.consumption.m0", 10) |>
+     g3_init_val("st_weightjones.consumption.m1", 0) |>
+     g3_init_val("st_weightjones.consumption.m2", 0) |>
+     g3_init_val("st_weightjones.consumption.m3", 1) |>
+     g3_init_val("st_weightjones.weightjones.q1", 0.001) |>
+     g3_init_val("st_weightjones.weightjones.q2", 0.001) |>
+     g3_init_val("st_weightjones.weightjones.q3", 1) |>
+     g3_init_val("st_weightjones.weightjones.q4", 1e-4) |>
+     g3_init_val("st_weightjones.weightjones.q5", -10) |>
+     identity() -> params
> 
> result <- model_fn(params)
> r <- attributes(result)
> 
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_vbsimple__num[,,time='2000-01']), '
+              age1      age2      age3      age4      age5
+ 50:60    961352.4  961352.4  961352.4  961352.4  961352.4
+ 60:70    964474.8  964474.8  964474.8  964474.8  964474.8
+ 70:80    967582.0  967582.0  967582.0  967582.0  967582.0
+ 80:90    970684.2  970684.2  970684.2  970684.2  970684.2
+ 90:100   973784.3  973784.3  973784.3  973784.3  973784.3
+ 100:Inf 1162122.3 1162122.3 1162122.3 1162122.3 1162122.3
+ ', tolerance = 1e-7), "late_st_vbsimple__num[,,time='2000-01']")
ok - late_st_vbsimple__num[,,time='2000-01']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_vbsimple__num[,,time='2000-02']), '
+              age1      age2      age3      age4      age5
+ 50:60    886896.4  886896.4  886896.4  886896.4  886896.4
+ 60:70    896034.1  896034.1  896034.1  896034.1  896034.1
+ 70:80    905127.4  905127.4  905127.4  905127.4  905127.4
+ 80:90    914206.1  914206.1  914206.1  914206.1  914206.1
+ 90:100   923278.8  923278.8  923278.8  923278.8  923278.8
+ 100:Inf 1474457.3 1474457.3 1474457.3 1474457.3 1474457.3
+ ', tolerance = 1e-7), "late_st_vbsimple__num[,,time='2000-02']")
ok - late_st_vbsimple__num[,,time='2000-02']
> 
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_vbsimple__num[,,time='2000-01']), '
+              age1      age2      age3      age4      age5
+ 50:60    961352.4  961352.4  961352.4  961352.4  961352.4
+ 60:70    964474.8  964474.8  964474.8  964474.8  964474.8
+ 70:80    967582.0  967582.0  967582.0  967582.0  967582.0
+ 80:90    970684.2  970684.2  970684.2  970684.2  970684.2
+ 90:100   973784.3  973784.3  973784.3  973784.3  973784.3
+ 100:Inf 1162122.3 1162122.3 1162122.3 1162122.3 1162122.3
+ ', tolerance = 1e-7), "late_st_vbsimple__num[,,time='2000-01']")
ok - late_st_vbsimple__num[,,time='2000-01']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_vbsimple__num[,,time='2000-02']), '
+              age1      age2      age3      age4      age5
+ 50:60    886896.4  886896.4  886896.4  886896.4  886896.4
+ 60:70    896034.1  896034.1  896034.1  896034.1  896034.1
+ 70:80    905127.4  905127.4  905127.4  905127.4  905127.4
+ 80:90    914206.1  914206.1  914206.1  914206.1  914206.1
+ 90:100   923278.8  923278.8  923278.8  923278.8  923278.8
+ 100:Inf 1474457.3 1474457.3 1474457.3 1474457.3 1474457.3
+ ', tolerance = 1e-7), "late_st_vbsimple__num[,,time='2000-02']")
ok - late_st_vbsimple__num[,,time='2000-02']
> 
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_multspec__num[,,time='2000-01']), '
+               age1       age2       age3       age4       age5
+ 50:60     816546.8   816546.8   816546.8   816546.8   816546.8
+ 60:70    1168617.6  1168617.6  1168617.6  1168617.6  1168617.6
+ 70:80    1526624.9  1526624.9  1526624.9  1526624.9  1526624.9
+ 80:90    1890891.9  1890891.9  1890891.9  1890891.9  1890891.9
+ 90:100   2261515.3  2261515.3  2261515.3  2261515.3  2261515.3
+ 100:Inf 18903686.1 18903686.1 18903686.1 18903686.1 18903686.1
+ ', tolerance = 1e-7), "late_st_multspec__num[,,time='2000-01']")
ok - late_st_multspec__num[,,time='2000-01']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_multspec__num[,,time='2000-02']), '
+             age1     age2     age3     age4     age5
+ 50:60    6629480  6629480  6629480  6629480  6629480
+ 60:70    8108169  8108169  8108169  8108169  8108169
+ 70:80    9611799  9611799  9611799  9611799  9611799
+ 80:90   11141721 11141721 11141721 11141721 11141721
+ 90:100  12698340 12698340 12698340 12698340 12698340
+ 100:Inf 82595455 82595455 82595455 82595455 82595455
+ ', tolerance = 1e-7), "late_st_multspec__num[,,time='2000-02']")
ok - late_st_multspec__num[,,time='2000-02']
> 
> ################### st_weightjones
> 
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_weightjones__num[,,time='2000-01']), '
+              age1      age2      age3      age4      age5
+ 50:60    548992.4  548992.4  548992.4  548992.4  548992.4
+ 60:70    462295.2  462295.2  462295.2  462295.2  462295.2
+ 70:80    373964.0  373964.0  373964.0  373964.0  373964.0
+ 80:90    284078.4  284078.4  284078.4  284078.4  284078.4
+ 90:100   192663.0  192663.0  192663.0  192663.0  192663.0
+ 100:Inf 4138007.0 4138007.0 4138007.0 4138007.0 4138007.0
+ ', tolerance = 1e-7), "late_st_weightjones__num[,,time='2000-01']")
ok - late_st_weightjones__num[,,time='2000-01']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_weightjones__num[,,time='2000-02']), '
+             age1     age2     age3     age4     age5
+ 50:60    3059086  3059086  3059086  3059086  3059086
+ 60:70    3845789  3845789  3845789  3845789  3845789
+ 70:80    4645760  4645760  4645760  4645760  4645760
+ 80:90    5459719  5459719  5459719  5459719  5459719
+ 90:100   6287882  6287882  6287882  6287882  6287882
+ 100:Inf 43474994 43474994 43474994 43474994 43474994
+ ', tolerance = 1e-7), "late_st_weightjones__num[,,time='2000-02']")
ok - late_st_weightjones__num[,,time='2000-02']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_weightjones__wgt[,,time='2000-01']), '
+             age1     age2     age3     age4     age5
+ 50:60   1036.016 1036.016 1036.016 1036.016 1036.016
+ 60:70   1042.962 1042.962 1042.962 1042.962 1042.962
+ 70:80   1050.028 1050.028 1050.028 1050.028 1050.028
+ 80:90   1057.206 1057.206 1057.206 1057.206 1057.206
+ 90:100  1064.480 1064.480 1064.480 1064.480 1064.480
+ 100:Inf 1056.993 1056.993 1056.993 1056.993 1056.993
+ ', tolerance = 1e-6), "late_st_weightjones__wgt[,,time='2000-01']")
ok - late_st_weightjones__wgt[,,time='2000-01']
> ok(gadget3:::ut_cmp_df(as.data.frame(r$late_st_weightjones__wgt[,,time='2000-02']), '
+             age1     age2     age3     age4     age5
+ 50:60   1324.144 1324.144 1324.144 1324.144 1324.144
+ 60:70   1386.658 1386.658 1386.658 1386.658 1386.658
+ 70:80   1450.290 1450.290 1450.290 1450.290 1450.290
+ 80:90   1515.036 1515.036 1515.036 1515.036 1515.036
+ 90:100  1580.892 1580.892 1580.892 1580.892 1580.892
+ 100:Inf 1532.325 1532.325 1532.325 1532.325 1532.325
+ ', tolerance = 1e-6), "late_st_weightjones__wgt[,,time='2000-02']")
ok - late_st_weightjones__wgt[,,time='2000-02']
> 
> gadget3:::ut_tmb_r_compare2(model_fn, model_cpp, params)
# skip: not running TMB tests
NULL
> ######## Default params
> 
> proc.time()
   user  system elapsed 
   3.68    0.12    3.79 
1..10
# Looks like you passed all 10 tests.