list_models() and get_model() differences

[mmundy42]

list_models() returns the number of reactions and compounds that it got from the workspace list() function which is not accurate when there are integrated gap fills. get_model() returns the data from the integrated model so the length of the lists in the model_data structure is different.

The number of genes as calculated by the workspace list() function is different than the length of the list returned in the model_data structure. For genome 226186.12 list_models() says there are 1,159 genes but get_model() returns a list of 9,871 genes.

Also, there is no interface to get both model stats and model data in one function call. If you call list_models() to get model stats you'll get stats for all of the user's models. And if you call get_model() to get model data you don't get information like the name, genome, etc.

[mmundy42]

Figured out that if I only count the genes that have a reaction list in the list of genes returned by get_model(), the count matches. So in the model for genome 226186.12 there are 1,159 genes with a reaction list.