diff --git a/src/solutions/pde_solutions.jl b/src/solutions/pde_solutions.jl
index 7677d6f69..bfc9b2f30 100644
--- a/src/solutions/pde_solutions.jl
+++ b/src/solutions/pde_solutions.jl
@@ -29,7 +29,7 @@ Solution to a PDE, solved from an ODEProblem generated by a discretizer.
   exited due to an error. For more details, see
   [the return code documentation](https://docs.sciml.ai/SciMLBase/stable/interfaces/Solutions/#retcodes).
 """
-struct PDETimeSeriesSolution{T, N, uType, Disc, Sol, DType, tType, domType, ivType, dvType,
+struct PDETimeSeriesSolution{T, N, uType, Disc, isCpx, Sol, DType, tType, domType, ivType, dvType,
                              P, A,
                              IType} <: AbstractPDETimeSeriesSolution{T, N, uType, Disc}
     u::uType
@@ -77,7 +77,7 @@ Solution to a PDE, solved from an NonlinearProblem generated by a discretizer.
   callback (`sol.retcode === ReturnCode.Terminated`), or whether it exited due to an error. For more
   details, see the return code section of the ODEProblem.jl documentation.
 """
-struct PDENoTimeSolution{T, N, uType, Disc, Sol, domType, ivType, dvType, P, A,
+struct PDENoTimeSolution{T, N, uType, Disc, isCpx, Sol, domType, ivType, dvType, P, A,
                          IType} <: AbstractPDENoTimeSolution{T, N, uType, Disc}
     u::uType
     original_sol::Sol
@@ -115,6 +115,10 @@ function SciMLBase.wrap_sol(sol,
     end
 end
 
+iscomplex(sol::PDETimeSeriesSolution{T, N, U, D, C}) where {T, N, U, D, C} = C <: Val{true} ? true : false
+iscomplex(sol::PDENoTimeSolution{T, N, U, D, C}) where {T, N, U, D, C} = C <: Val{true} ? true : false
+
+
 function Base.show(io::IO, m::MIME"text/plain", A::PDETimeSeriesSolution)
     println(io, string("retcode: ", A.retcode))
     println(io, string("Interpolation: "), typeof(A.interp))