Skip to content

Commit a658a24

Browse files
update since FEM is gone
1 parent e86426f commit a658a24

File tree

3 files changed

+1
-20
lines changed

3 files changed

+1
-20
lines changed

REQUIRE

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -2,7 +2,7 @@ julia 0.6
22
RecipesBase 0.1.0
33
DiffEqBase 3.0.0
44
RecursiveArrayTools 0.4.0
5-
DiffEqPDEBase
5+
DiffEqPDEBase 0.4.0
66
DiffEqMonteCarlo
77
DiffEqNoiseProcess
88
Juno

src/convergence.jl

Lines changed: 0 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -107,19 +107,6 @@ function analyticless_test_convergence(dts::AbstractArray,prob::AbstractODEProbl
107107
ConvergenceSimulation(solutions,dts,auxdata=auxdata)
108108
end
109109

110-
#=
111-
function test_convergence(dts::AbstractArray,dxs::AbstractArray,prob::AbstractHeatProblem,convergence_axis;T=1,alg=:Euler)
112-
if length(dts)!=length(dxs) error("Lengths of dts!=dxs. Invalid convergence simulation") end
113-
solutions = [solve(parabolic_squaremesh([0 1 0 1],dxs[i],dts[i],T,:dirichlet),prob,alg=alg) for i in eachindex(dts)]
114-
auxdata = Dict(
115-
:dts => [sol.fem_mesh.dt for sol in solutions],
116-
:dxs => [sol.fem_mesh.dx for sol in solutions],
117-
:Δμs => [sol.fem_mesh.μ for sol in solutions],
118-
:Δνs => [sol.fem_mesh.ν for sol in solutions])
119-
return(ConvergenceSimulation(solutions,convergence_axis,auxdata=auxdata))
120-
end
121-
=#
122-
123110
function test_convergence(probs,convergence_axis,alg;kwargs...)
124111
ConvergenceSimulation([solve(prob,alg;kwargs...) for prob in probs],convergence_axis)
125112
end

src/test_solution.jl

Lines changed: 0 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -57,12 +57,6 @@ function appxtrue(sol::AbstractODESolution,sol2::TestSolution)
5757
build_solution(sol,timeseries_analytic,errors)
5858
end
5959

60-
function appxtrue(sol::AbstractFEMSolution,sol2::AbstractFEMSolution)
61-
u_analytic = sol2[end]
62-
errors = Dict(:l∞=>maximum(abs.(sol[end]-u_analytic)),:l2=>norm(sol[end]-u_analytic,2))
63-
FEMSolution(sol,u_analytic,errors)
64-
end
65-
6660
"""
6761
`appxtrue(sol::AbstractODESolution,sol2::AbstractODESolution)`
6862

0 commit comments

Comments
 (0)