GSoC 5th week review
In this week, I continued the previous work on defect control adaptivity and added all of the BVP test problems and the others to the DiffEqProblemLibrary.jl, and closed the relating issue. All of the work in this week can be summarized as follows:
- Fixed the indexing on discrete stages
k_discrete
and interpolation coefficientsk_interp
. The original indexing method fork_discrete
andk_interp
is awkward, by using a new matrix structure, we can directly usek_discrete[ith_subinterval, ith_stage]
andk_interp[ith_subinterval, ith_stage]
to addressing discrete stages and interpolantion coefficients. - Fixed the new stages and solution addition in the routine
sum_stages
. After we scaled the new stages usingdt
, we need to add the new stages to the solution in the specific interval. - Add the nonlinear BVP test problems to DiffEqProblemLibrary.jl
- Add the regular nonlinear BVP test problems to DiffEqProblemLibrary.jl
Next week, I will continue to work on the problem that defect_norm
is always bigger than the threshold, by comparing with the MIRKDC routine and monitoring the computing flow, while the other routines are the same, I think the problem lies in the collocation part.