Daniel VandenHeuvel

Indexing and iterating in ODEProblem and modelingtoolkitize / error messages

3

I am working on some complicated 2D finite volume code and am trying to refactor it to allow for `modelingtoolkitize`. In doing so I am having some difficulties dealing with...

Possibly missing documentation for boundary conditions

2

At http://juliamath.github.io/Interpolations.jl/latest/control/, it is said > The following boundary conditions are implemented: Flat, Line (alternatively, Natural), Free, Periodic and Reflect; **their mathematical implications are described in detail in their docstrings**....

NOMADOpt: libnomadInterface not defined in doc example

3

The `NOMADOpt()` example in http://optimization.sciml.ai/dev/optimization_packages/nomad/ doesn't seem to work: ``` using Optimization using OptimizationNOMAD rosenbrock(x, p) = (p[1] - x[1])^2 + p[2] * (x[2] - x[1]^2)^2 x0 = zeros(2) p...

Support for scalar inputs

1

There is currently no (little? Haven't tried all the solvers) support for scalar inputs, e.g. ``` julia> f = (u, p) -> u^2; u0 = 0.0; prob = OptimizationProblem(f, u0);...

DisjointSets with Any type never stops

1

The following code seems to never stop: ``` julia> using DataStructures julia> DisjointSets{Any}() # good DisjointSets{Any}(Dict{Any, Int64}(), Any[], IntDisjointSets{Int64}(Int64[], Int64[], 0)) julia> DisjointSets{Any}([]) ERROR: InterruptException: Stacktrace: [1] DisjointSets{Any}(xs::Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Tuple{Vector{Any}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}) (repeats 204...

More types of infinity transformations

2

Currently the only type of transformation used for handling infinite bounds is $u \mapsto a + t/(1-t)$ and $u \mapsto t/(1-t)^2$. It could be nice to eventually support other transformations....

Issues with KLUFactorization() and sparse matrices with SplitODEFunction

1

There are some issues when using a `SplitODEFunction` with `KLUFactorization()`. If using a `MatrixOperator`: ```julia using OrdinaryDiffEq, LinearSolve, SparseArrays f2 = (du, u, p, t) -> u[2] A = MatrixOperator([1.0...

In https://github.com/JuliaDebug/Cthulhu.jl#viewing-the-internal-representation-of-julia-code, the section starts as > Anyone using Cthulhu to investigate the behavior of Julia's compiler will prefer to examine the While Cthulhu tries to place type-annotations... It looks...

Automatically change branch name if it already exists

4

I had made a PR to someones repo with this tool, and then had to follow up with another PR to include some other file types. Unfortunately, after going through...

GeneralDomain documentation is unclear

2

The documentation for `GeneralDomain` shows the docstring ``` GeneralDomain(g, u=nothing; nlsolve=NLSOLVEJL_SETUP(), save=true, abstol=nothing, scalefactor=nothing, autonomous=numargs(g)==2, nlopts=Dict(:ftol => 10*eps())) ``` which should now be ``` GeneralDomain(g, u = nothing; nlsolve =...