fix issue with ReverseDiff.TrackedArray treatment

the original function was actually effectively removing
the tracking with kfs = zeros(Q.parameters[1],length(ph.reactions))
resulting in a Array{Float64,1}

This behavior was originally hidden by a bug in ReverseDiff.jl that
was subsequently fixed making us incompatible with recent versions
this commit should make us compatible with newer versions of ReverseDiff

Author: mjohnson541

src/PhaseState.jl

@@ -2,6 +2,7 @@ using Parameters
 using SpecialFunctions
 using LinearAlgebra
 using Tracker
+using ReverseDiff
 using RecursiveArrayTools
 
 @inline function calcgibbs(ph::U,T::W) where {U<:IdealPhase,W<:Real}
@@ -56,7 +57,7 @@ end
 export getkf
 
 @inline function getkfs(ph::U,T::W1,P::W2,C::W3,ns::Q,V::W4,phi) where {U,W1,W2,W3,W4<:Real,Q<:AbstractArray}
-    kfs = zeros(Q.parameters[1],length(ph.reactions))
+    kfs = similar(ns,length(ph.reactions))
     i = 1
     oldind = 1
     ind = 0
@@ -71,6 +72,7 @@ export getkf
     end
     return kfs
 end
+
 export getkfs
 
 """