Outdated tutorial

[blegat]

This tutorial seems outdated.
This line:

prob_femheat_diffusionconstants = HeatProblem(u₀,f,D=[.01 .001])

throws

ERROR: MethodError: no method matching DiffEqPDEBase.HeatProblem{islinear,isstochastic,MeshType,F,F2,F3,F4,F5,F6,F7,DiffType}(::##3#4, ::##1#2; D=[0.01 0.001])
Closest candidates are:
  DiffEqPDEBase.HeatProblem{islinear,isstochastic,MeshType,F,F2,F3,F4,F5,F6,F7,DiffType}(::Any, ::Any, ::Any; noisetype, gD, gN, σ, numvars, D) at /home/blegat/.julia/v0.5/DiffEqPDEBase/src/problems.jl:41
  DiffEqPDEBase.HeatProblem{islinear,isstochastic,MeshType,F,F2,F3,F4,F5,F6,F7,DiffType}(::Any, ::Any, ::Any, ::Any; noisetype, gN, σ, numvars, D) at /home/blegat/.julia/v0.5/DiffEqPDEBase/src/problems.jl:16
  DiffEqPDEBase.HeatProblem{islinear,isstochastic,MeshType,F,F2,F3,F4,F5,F6,F7,DiffType}{T}(::Any) at sysimg.jl:53 got unsupported keyword argument "D"
 in eval_user_input(::Any, ::Base.REPL.REPLBackend) at ./REPL.jl:64
 in macro expansion at ./REPL.jl:95 [inlined]
 in (::Base.REPL.##3#4{Base.REPL.REPLBackend})() at ./event.jl:68

[ChrisRackauckas]

Yeah, the code here is quite old. The mesh needs to be passed to the HeatProblem, and solve needs to be in the more modern form using type dispatch.

using DifferentialEquations, Plots
f = (t,x,u)  -> [zeros(size(x,1))    zeros(size(x,1))]
u₀ = (x) -> [float((abs(x[:,1]-.5) .< 1e-6) & (abs(x[:,2]-.5) .< 1e-6)) float((abs(x[:,1]-.5) .< 1e-6) & (abs(x[:,2]-.5) .< 1e-6))]  # size (x,2), 2 meaning 2 variables

## Build the mesh
T = 30
Δx = 1//2^(4)
Δt = 1//2^(2)
fem_mesh = parabolic_squaremesh([0 1 0 1],Δx,Δt,T,:neumann)

## Build the problem type
prob_femheat_diffusionconstants = HeatProblem(u₀,f,fem_mesh,D=[.01 .001])

sol = solve(prob,alg=FEMDiffEqHeatImplicitEuler())
animate(sol,zlim=(0,.01),cbar=false,filename="diffconsts.gif")

I didn't give it a try, but that should be it. This tutorial I know is up-to-date:

http://docs.juliadiffeq.org/latest/tutorials/femheat_example.html

Could you run it again with this?

[blegat]

Your code code give the following error:

ERROR: LoadError: LoadError: UndefVarError: @compat not defined
 in macro expansion; at ./none:2 [inlined]
 in anonymous at ./<missing>:?
while loading /home/blegat/.julia/v0.5/StochasticDiffEq/src/algorithms.jl, in expression starting on line 5
while loading /home/blegat/.julia/v0.5/StochasticDiffEq/src/StochasticDiffEq.jl, in expression starting on line 21
ERROR: LoadError: LoadError: Failed to precompile StochasticDiffEq to /home/blegat/.julia/lib/v0.5/StochasticDiffEq.ji.
while loading /home/blegat/.julia/v0.5/DifferentialEquations/src/DifferentialEquations.jl, in expression starting on line 9
while loading /home/blegat/tmp.jl, in expression starting on line 1

The package StochasticDiffEq fails to compile :(

[ChrisRackauckas]

Yeah, this'll have to wait until the fix for that goes through. Sorry about the problem... hasn't been a good day for DiffEq I guess.

[ChrisRackauckas]

Should be able to run now.

[ChrisRackauckas]

updated.

[blegat]

Thanks !