The best way is to use the .eval() method. This allows you to eagerly compute a single variable in your model. You should check that all of these different components have the correct shapes with e.g. (b_time * time).shape.eval(), (b_hematocrit*hematocrit).shape.eval(), and so on.