metal_small example failing at make_ast_inputs

[crobinett1]

In attempting to create an input AST list from the metal_small example, I am experiencing the following error which appears to be due to differing array dimensions in pick_models_toothpick_style.

---

ValueError Traceback (most recent call last) Cell In[4], line 1 ----> 1 make_ast_inputs.make_ast_inputs('beast_settings.txt')

File ~/miniconda3/lib/python3.10/site-packages/beast/tools/run/make_ast_inputs.py:79, in make_ast_inputs(beast_settings_info, pick_method) 76 min_N_per_flux = settings.ast_n_per_flux_bin 77 bins_outfile = "./{0}/{0}_ASTfluxbins.txt".format(settings.project) ---> 79 chosen_seds = pick_models_toothpick_style( 80 modelsedgrid_filename, 81 settings.filters, 82 N_fluxes, 83 min_N_per_flux, 84 outfile=outfile_seds, 85 outfile_params=outfile_params, 86 bins_outfile=bins_outfile, 87 ) 89 if pick_method == "random_pick": 90 91 # construct magnitude cuts 92 mag_cuts = settings.ast_maglimit

File ~/miniconda3/lib/python3.10/site-packages/beast/observationmodel/ast/make_ast_input_list.py:188, in pick_models_toothpick_style(sedgrid_fname, filters, N_fluxes, min_N_per_flux, outfile, outfile_params, bins_outfile, bright_cut) 182 add_these[r] = True 184 # If all these bins are full... 185 else: 186 # ... do not include this model again, since we will reject it 187 # anyway. --> 188 include_mask[idxs == rand_idx] = False 190 # Add the approved models 191 chosen_idxs.extend(rand_idx[add_these])

ValueError: operands could not be broadcast together with shapes (493713,) (100000,)

[karllark]

Hummm.... I've tested most of the rest of the metal_small example in the past. But have to admit I have not been testing the make_ast_inputs part. Can you confirm that you are using the as provided beast_settings.txt? If so, I'll work on testing this and seeing if I can find a fix.

Thanks for opening an issue - much appreciated.

[crobinett1]

The beast_settings.txt I am using is an exact copy of the one provided in the beast-examples/metal_small folder in the beast-examples repository.

Thanks for looking into this.