@@ -9,24 +9,30 @@ xt::xarray<double> forward()
9
9
10
10
pyprob_cpp::setDefaultControl (true );
11
11
pyprob_cpp::setDefaultReplace (false );
12
-
13
12
auto normal1 = pyprob_cpp::distributions::Normal (prior_mean, prior_stddev);
14
- auto mu1 = pyprob_cpp::sample (normal1, " mu1" );
15
- mu1 = pyprob_cpp::sample (normal1, " mu1" );
13
+ xt::xarray<double > mu1;
14
+ for (int i = 0 ; i < 2 ; i++)
15
+ {
16
+ mu1 = pyprob_cpp::sample (normal1, " mu1" );
17
+ }
16
18
17
19
pyprob_cpp::setDefaultControl (true );
18
20
pyprob_cpp::setDefaultReplace (true );
19
-
20
21
auto normal2 = pyprob_cpp::distributions::Normal (mu1, prior_stddev);
21
- auto mu2 = pyprob_cpp::sample (normal2, " mu2" );
22
- mu2 = pyprob_cpp::sample (normal2, " mu2" );
22
+ xt::xarray<double > mu2;
23
+ for (int i = 0 ; i < 2 ; i++)
24
+ {
25
+ mu2 = pyprob_cpp::sample (normal2, " mu2" );
26
+ }
23
27
24
28
pyprob_cpp::setDefaultControl (false );
25
29
pyprob_cpp::setDefaultReplace (false );
26
-
27
30
auto normal3 = pyprob_cpp::distributions::Normal (mu2, prior_stddev);
28
- auto mu3 = pyprob_cpp::sample (normal3, " mu3" );
29
- mu3 = pyprob_cpp::sample (normal3, " mu3" );
31
+ xt::xarray<double > mu3;
32
+ for (int i = 0 ; i < 2 ; i++)
33
+ {
34
+ mu3 = pyprob_cpp::sample (normal3, " mu3" );
35
+ }
30
36
31
37
auto likelihood = pyprob_cpp::distributions::Normal (mu3, likelihood_stddev);
32
38
pyprob_cpp::observe (likelihood, " obs" );
0 commit comments