if we would use this for DIRECT source terms, we would have ```_mesh->getNNodes() > prop_vals.size()```. I cannot think of a case where we would need *more* prop_vals than nodes but having fewer seems plausible. 
If you save vectors/matrices as linear double values, you might also have more prop_vals than nodes.