model.eval() will use the running stats, which might be wrong or skewed depending on your dataset.
If you are concerned, if the batch norm layers give valid answers, could you please compare the same outputs with the same modes?
E.g. make sure that a fixed input creates the same output in model.eval().
Note that you will not get the same outputs for the fixed data in model.train(), as the running stats will be updated.