Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Batching Nodally defined model parameters #169

Open
odunbar opened this issue Jul 22, 2020 · 0 comments
Open

Batching Nodally defined model parameters #169

odunbar opened this issue Jul 22, 2020 · 0 comments

Comments

@odunbar
Copy link
Collaborator

odunbar commented Jul 22, 2020

We have a batching procedure for observed states in the DataAssimilator. This only batches the states and not model parameters. It is inefficient when dealing with the TransitionRates type model parameters which can be nodally defined.

Consider observing 100 obseved statess (and only 1 state per node) in 4 batches of 25. Assume the TransitionRates clinical parameters are made up of 3 constant parameters and one nodally defined parameter, in the order [const, const, nodal, const], the total size of the parameters is 1+1+100+1 = 103. In the current setup we will do 4 batches of 25 states and 103 parameters for each EAKF batch. If we found the correct nodes that relate to the batched states we could do 4 batches of 25 states and 1+1+25+1=28 parameters.

The problem is to match which model parameters will go into which batch as the clinical parameters, when loaded into the data assimilator are all concatenated together into an np.array.

Coauthored with @jinlong83

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant