Skip to content
Snippets Groups Projects
Commit f38a1eda authored by Johannes Pfeifer's avatar Johannes Pfeifer
Browse files

posterior_sampler.m: trap case where record.ProposalCovariance is not filled (slice sampler)

parent 765109ee
1 merge request!1869posterior_sampler.m: trap case where record.ProposalCovariance is not filled (slice sampler)
...@@ -156,10 +156,12 @@ else ...@@ -156,10 +156,12 @@ else
record.FunctionEvalPerIteration(offset+1:sum(nBlockPerCPU(1:j)))=fout(j).record.FunctionEvalPerIteration(offset+1:sum(nBlockPerCPU(1:j))); record.FunctionEvalPerIteration(offset+1:sum(nBlockPerCPU(1:j)))=fout(j).record.FunctionEvalPerIteration(offset+1:sum(nBlockPerCPU(1:j)));
record.LastSeeds(offset+1:sum(nBlockPerCPU(1:j)))=fout(j).record.LastSeeds(offset+1:sum(nBlockPerCPU(1:j))); record.LastSeeds(offset+1:sum(nBlockPerCPU(1:j)))=fout(j).record.LastSeeds(offset+1:sum(nBlockPerCPU(1:j)));
if j==1 if j==1
if isfield(fout(j).record,'ProposalCovariance') && isfield(fout(j).record,'ProposalScaleVec')
record.ProposalCovariance=fout(j).record.ProposalCovariance; record.ProposalCovariance=fout(j).record.ProposalCovariance;
record.ProposalScaleVec=fout(j).record.ProposalScaleVec; record.ProposalScaleVec=fout(j).record.ProposalScaleVec;
end end
end end
end
options_.parallel_info.parallel_recover = 0; options_.parallel_info.parallel_recover = 0;
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment