Evaluate the rate constants in batch

[jiweiqi]

reactorch/reactorch/Solution.py

Line 240 in aa6d9d2

for i in range(self.n_reactions):

Currently, in order to avoid in-place operation, we iteratively evaluate the forward rates by looping the reactions. It seems that we can evaluate it in matrix form by some careful coding without in-place operation. Initial thoughts are having a mask to evaluate the normal reactions in batch and the pressure-dependent reaction iteratively. Since the pressure-dependent reactions are usually only a small portion, and related to C0-C1 core mechanism.