Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

WIP: Add Lobatto IIIa-b and Radau II tableaus 2.0 #115

Closed
wants to merge 108 commits into from
Closed
Changes from 1 commit
Commits
Show all changes
108 commits
Select commit Hold shift + click to select a range
2502733
First test that works
axla-io Sep 21, 2023
70f4aa0
Renamed MIRKCache to RKCache
axla-io Sep 21, 2023
f6366a6
Added RKTableau and Lobatto IIIb tableaus
axla-io Sep 21, 2023
c55b526
First test that works
axla-io Sep 21, 2023
411eba3
Renamed MIRKCache to RKCache
axla-io Sep 21, 2023
f098715
Added RKTableau and Lobatto IIIb tableaus
axla-io Sep 21, 2023
16a25bf
Merge branch 'al/lobattotableau' of https://github.com/axla-io/Bounda…
axla-io Oct 3, 2023
e47f05e
Fixing faulty rebase
axla-io Oct 3, 2023
7929300
Added LobattoIIIa
axla-io Oct 3, 2023
d6e1d3e
Radau IIa tables implemented
axla-io Oct 4, 2023
71dee6b
Full size jacobian
axla-io Oct 9, 2023
15b5ace
Full block jacobian
axla-io Oct 10, 2023
25ea9bc
revert to old banded comp
axla-io Oct 10, 2023
560d50b
Added nested nonlinear solve option
axla-io Oct 10, 2023
04b0e63
Added oop collocation
axla-io Oct 10, 2023
faa2b27
Added nested dispatch
axla-io Oct 10, 2023
29f81ef
Dispatch fix
axla-io Oct 10, 2023
cd1c23e
Nested solve working
axla-io Oct 10, 2023
9cca23b
Split up the 2Point BVP
avik-pal Oct 6, 2023
911df33
Update odeinterface_ex7.jl
ChrisRackauckas Oct 7, 2023
74d6ab7
Update odeinterface_ex7.jl
ChrisRackauckas Oct 7, 2023
936a955
Update Project.toml
ChrisRackauckas Oct 7, 2023
7eb87d0
Update Project.toml
ChrisRackauckas Oct 7, 2023
275da3c
Fast path for Two Point BVPs
avik-pal Sep 12, 2023
5d027ab
Allow non vector inputs for Single Shooting
avik-pal Sep 13, 2023
28dfba2
Working version of Multiple Shooting :tada:
avik-pal Sep 13, 2023
d06ae40
Tests
avik-pal Sep 23, 2023
5b90ab9
orbital multiple shooting
avik-pal Sep 23, 2023
c4265a2
Add type stability tests
avik-pal Sep 23, 2023
4de194d
Port shooting tests from NeuralBVP
avik-pal Sep 23, 2023
ad7ebae
Incorrect merge
avik-pal Sep 28, 2023
956e459
fix type stability??
avik-pal Sep 28, 2023
0e721e0
Update test/shooting/shooting_tests.jl
avik-pal Oct 1, 2023
0ac6dae
Fix formatting
avik-pal Oct 2, 2023
fa411ee
Mark Single Shooting test as borken
avik-pal Oct 2, 2023
898edb3
Fix alg
avik-pal Oct 2, 2023
cf18105
Use ensemble problem to speed up multiple shooting
avik-pal Oct 2, 2023
d27dfa2
Patch interpolation
avik-pal Oct 3, 2023
8baf3b0
Don't use custom jacobian
avik-pal Oct 5, 2023
ee8fb00
Use sparse Jacobian for Multiple Shooting
avik-pal Oct 5, 2023
9b7295a
xSetup colorvec construction
avik-pal Oct 5, 2023
a7e7683
Fix tests
avik-pal Oct 9, 2023
da513ff
Update the TwoPointBVP code
avik-pal Oct 9, 2023
e06a2fb
Initial Working version with colorring
avik-pal Oct 9, 2023
272b845
Generalize the Jacobian Computation for MultipleShooting
avik-pal Oct 10, 2023
c6381fc
Two Point Specialization for MultipleShooting
avik-pal Oct 10, 2023
7ea26d7
Cleanup sparse prototype creation
avik-pal Oct 10, 2023
cf85900
Move things around for MIRK
avik-pal Oct 11, 2023
cb0080d
Single Shooting Cleanup
avik-pal Oct 11, 2023
6bde443
Fix tests and make Multiple Shooting Type Stable
avik-pal Oct 11, 2023
bfa4a79
Split up the tests
avik-pal Oct 12, 2023
6caf472
Fix interpolant evaluation error
ErikQQY Oct 12, 2023
7b27f62
Added defect estimate functions
axla-io Oct 12, 2023
8c9aaf5
Added interpolation tables proper implementation
axla-io Oct 12, 2023
ab7d2cb
Cleanup after rebase
axla-io Oct 13, 2023
e8d78ad
Added poly_max
axla-io Oct 13, 2023
d3b03d0
updated adaptivity functions
axla-io Oct 13, 2023
11d1517
Started work on table coeffs
axla-io Oct 13, 2023
87189f3
Added coefficients
axla-io Oct 13, 2023
fc57806
Adding modification to table types for adaptivity.
axla-io Oct 16, 2023
4c0fa9b
Interpolation for RK methods
axla-io Oct 16, 2023
6a845b2
LobattoIIIa interpolation works.
axla-io Oct 19, 2023
d542397
tau_star for radau and lobatto
axla-io Oct 19, 2023
b8feb2c
Updated tables definition
axla-io Oct 19, 2023
e5eddaf
Fixed interpolation template
axla-io Oct 19, 2023
5735bc4
Adaptivity for two sides and doesn't fail
axla-io Oct 19, 2023
67e8701
Working adaptivity and interpolation
axla-io Oct 20, 2023
a0d226b
Debugged interpolation.
axla-io Oct 20, 2023
ab65390
Nested interpolation
axla-io Oct 20, 2023
5d9d380
Sketch implementation for FIRK sparse jacobian
axla-io Oct 24, 2023
02b9441
Sparse jacobian implementation for non-nested FIRK
axla-io Oct 25, 2023
d9598de
Sparse jac for M == 2
axla-io Oct 25, 2023
f323215
Generalized nested sparsity pattern
axla-io Oct 25, 2023
a2f141e
Simplest possible smart nested nlsolve
axla-io Dec 2, 2023
7ca14e7
Splitted MIRK and FIRK.
axla-io Dec 3, 2023
efd78c1
New default: nested
axla-io Dec 4, 2023
92563bc
Optimized allocations
axla-io Dec 4, 2023
b5139f3
Nicer initialization
axla-io Dec 4, 2023
525b7b7
removed adaptivity for RadauII1 and LobattoIIIb2
axla-io Dec 4, 2023
86d8459
Fixed uncommented line
axla-io Dec 4, 2023
848999e
User specified nested solve kwargs
axla-io Dec 4, 2023
e98dc0d
Added LobattoIIIC tableaus
axla-io Dec 4, 2023
ed9701c
Implemented interpolation polynomials for
axla-io Dec 4, 2023
2715132
Tested adaptivity for non nested
axla-io Dec 5, 2023
37cfd20
Removed unnecessary comments
axla-io Dec 5, 2023
a8879cf
Adaptivity for nested nonlinear solve now works
axla-io Dec 6, 2023
1fa33f9
Fixed merge conflicts for tests
axla-io Dec 9, 2023
a8a5f84
Nested nlsolve now works with AD
axla-io Dec 9, 2023
9950512
Nested nonlinear solve reduced allocs
axla-io Dec 10, 2023
5964bb6
Copy paste correct merge conflicts
axla-io Dec 10, 2023
f7a619c
More copy paste changes
axla-io Dec 10, 2023
dc6f7c4
It compiles
axla-io Dec 11, 2023
cdcc9d4
Everything except nested sparse jacobian works
axla-io Dec 11, 2023
d8af738
runs again but nested is slow
axla-io Dec 11, 2023
35c59d8
Adaptivity working.
axla-io Dec 11, 2023
afaf30f
Adaptivity and interpolation works again
axla-io Dec 11, 2023
66cdf8e
Adaptivity works for oop
axla-io Dec 11, 2023
24dadfd
Alloc reduce by storing transpose of a
axla-io Dec 12, 2023
d98d487
Started adding tests
axla-io Dec 12, 2023
8010305
Added test for non nested for radau
axla-io Dec 12, 2023
1184f83
Added lobatto convergence tests
axla-io Dec 12, 2023
09e68a5
Added docstrings for algorithms
axla-io Dec 12, 2023
b131096
permutedims instead of transpose
axla-io Dec 12, 2023
83a0547
Code fixed
axla-io Jan 2, 2024
1af08f9
Started cleaning up repo and reviewing tests
axla-io Feb 20, 2024
db1d8ee
Proper naming and order of methods
axla-io Mar 6, 2024
6962845
TP for expanded cache
axla-io Mar 7, 2024
d3fa4d3
Tests are working for SparseForwardDiff
axla-io Mar 7, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Update test/shooting/shooting_tests.jl
Co-authored-by: Qingyu Qu <52615090+ErikQQY@users.noreply.github.com>
  • Loading branch information
2 people authored and axla-io committed Oct 12, 2023
commit 0e721e0cfe93dcf0b6141e84db3d1383b29f94e5
12 changes: 12 additions & 0 deletions test/shooting/shooting_tests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -82,6 +82,18 @@ end

@testset "Shooting with Complex Values" begin
# Test for complex values
function f1!(du, u, p, t)
du[1] = u[2]
du[2] = -u[1]
return nothing
end

function bc1!(resid, sol, p, t)
t₀, t₁ = first(t), last(t)
resid[1] = sol(t₀)[1]
resid[2] = sol(t₁)[1] - 1
return nothing
end
u0 = [0.0, 1.0] .+ 1im
bvp = BVProblem(f1!, bc1!, u0, tspan)
resid_f = Array{ComplexF64}(undef, 2)
Expand Down