fix: Correct the size of the sampled state array

Was previously multiplying the grid size by the number of parameters, instead of raising to that power.
galacticusorg · Jan 13, 2025 · bfb8ed3 · bfb8ed3
1 parent 38044b1
commit bfb8ed3
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/source/posterior_sampling.state.samples.prior_grid.F90 b/source/posterior_sampling.state.samples.prior_grid.F90
@@ -99,8 +99,8 @@ subroutine priorGridSamples(self,simulationStates,modelParameters_)
     double precision                                           , dimension(size(modelParameters_))              :: stateVector
     type            (multiCounter              )                                                                :: counter
     integer         (c_size_t                  )                                                                :: i               , j
-        
-    allocate(simulationStates(self%countGrid*size(modelParameters_)))
+
+    allocate(simulationStates(self%countGrid**size(modelParameters_)))
     counter=multiCounter(spread(self%countGrid,1,size(modelParameters_)))
     i      =0
     do while (counter%increment())