fix: Set the simulation state count in each independent `modelLikelih…

…ood` class This is necessary so that, for example, the `modelLikelihoodGalaxyPopulation` class can store results to groups numbered by the state count.
galacticusorg · Dec 13, 2024 · 8557bc3 · 8557bc3
1 parent 07dff48
commit 8557bc3
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/source/models.likelihoods.independent_likelihoods.F90 b/source/models.likelihoods.independent_likelihoods.F90
@@ -307,6 +307,7 @@ Return the log-likelihood for the halo mass function likelihood function.
        end forall
        call modelLikelihood_%simulationState%update       (stateVectorMapped(1:size(modelLikelihood_%parameterMap)),logState=.false.,isConverged=.false.)
        call modelLikelihood_%simulationState%chainIndexSet(simulationState%chainIndex())
+       call modelLikelihood_%simulationState%countSet     (simulationState%count     ())
        ! Determine if the chain is already accepted - if it is we set the proposed prior to be impossible so that the model will not actually be evaluated.
        if (logLikelihoodCurrent > self%logLikelihoodAccept) then
           logPriorProposed_=logImpossible

diff --git a/source/models.likelihoods.independent_likelihoods.sequential.F90 b/source/models.likelihoods.independent_likelihoods.sequential.F90
@@ -262,6 +262,7 @@ Return the log-likelihood for the halo mass function likelihood function.
        end forall
        call modelLikelihood_%simulationState%update       (stateVectorMapped(1:size(modelLikelihood_%parameterMap)),logState=.false.,isConverged=.false.)
        call modelLikelihood_%simulationState%chainIndexSet(simulationState%chainIndex())
+       call modelLikelihood_%simulationState%countSet     (simulationState%count     ())
        ! Evaluate this likelihood
        timeEvaluate_=-1.0
        logLikelihood                                             =+modelLikelihood_%modelLikelihood_%evaluate(                                           &