diff --git a/Project.toml b/Project.toml
index 6da65e226..17debf801 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "VoronoiFVM"
 uuid = "82b139dc-5afc-11e9-35da-9b9bdfd336f3"
 authors = ["Jürgen Fuhrmann <juergen.fuhrmann@wias-berlin.de>", "Dilara Abdel", "Jan Weidner", "Alexander Seiler", "Patricio Farrell", "Matthias Liero"]
-version = "1.16.1"
+version = "1.17.0"
 
 [deps]
 BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
@@ -42,7 +42,7 @@ LinearAlgebra = "1.6"
 LinearSolve = "2.2"
 Printf = "1.6"
 Random = "1.6"
-RecursiveArrayTools = "2.14.2 - 3.4.1"
+RecursiveArrayTools = "2,3"
 RecursiveFactorization = "0.2.18"
 SparseArrays = "1.6"
 SparseDiffTools = "^1.19, 2"
diff --git a/src/vfvm_transientsolution.jl b/src/vfvm_transientsolution.jl
index 1947649df..d65d7b066 100644
--- a/src/vfvm_transientsolution.jl
+++ b/src/vfvm_transientsolution.jl
@@ -46,7 +46,7 @@ mutable struct TransientSolution{T, N, A, B} <: AbstractTransientSolution{T, N,
 end
 
 function TransientSolution(vec::AbstractVector{T}, ts, ::NTuple{N}) where {T, N}
-    TransientSolution{T, N, typeof(vec), typeof(ts)}(vec, ts, TransientSolverHistory())
+    TransientSolution{eltype(T), N, typeof(vec), typeof(ts)}(vec, ts, TransientSolverHistory())
 end
 
 TransientSolution(vec::AbstractVector, ts::AbstractVector) = TransientSolution(vec, ts, (size(vec[1])..., length(vec)))