My problem is this: Initialize some functions at every remake in PDEProblem
I’ve thought a possible solution would be to redefine
Temperature(t,x)
by defining a function:
function pde_solution(Tvec)
interp = LinearInterpolation(Tvec, [1.,2.,3.,4.,5.])
global Temperature(t, x) = interp(x)
solve(prob, TRBDF2())
end