search

SciML / BoundaryValueDiffEq.jl

Boundary value problem (BVP) solvers for scientific machine learning (SciML)
Other
41 stars 31 forks source link

Shooting methods fail on flow in a channel tests #191

Closed ErikQQY closed 2 weeks ago

ErikQQY commented 2 weeks ago

Describe the bug 🐞

Flow in a channel failed

Expected behavior

Successfully solve and get a successful retcode

Minimal Reproducible Example 👇

Without MRE, we would only be able to help you to a limited extent, and attention to the issue would be limited. to know more about MRE refer to wikipedia and stackoverflow.

    ]activate --temp
    ]add BoundaryValueDiffEq
    using BoundaryValueDiffEq, LinearAlgebra, LinearSolve

    function flow_in_a_channel!(du, u, p, t)
        R, P = p
        A, f′′, f′, f, h′, h, θ′, θ = u
        du[1] = 0
        du[2] = R * (f′^2 - f * f′′) - R * A
        du[3] = f′′
        du[4] = f′
        du[5] = -R * f * h′ - 1
        du[6] = h′
        du[7] = -P * f * θ′
        du[8] = θ′
    end

    function bc_flow!(resid, sol, p, tspan)
        t₁, t₂ = extrema(tspan)
        solₜ₁ = sol(t₁)
        solₜ₂ = sol(t₂)
        resid[1] = solₜ₁[4]
        resid[2] = solₜ₁[3]
        resid[3] = solₜ₂[4] - 1
        resid[4] = solₜ₂[3]
        resid[5] = solₜ₁[6]
        resid[6] = solₜ₂[6]
        resid[7] = solₜ₁[8]
        resid[8] = solₜ₂[8] - 1
    end

    tspan = (0.0, 1.0)
    p = [10.0, 7.0]
    u0 = zeros(8)

    flow_bvp = BVProblem{true}(flow_in_a_channel!, bc_flow!, u0, tspan, p)
    solve(flow_bvp, Shooting(AutoTsit5(Rosenbrock23()),
            NewtonRaphson(; autodiff = AutoForwardDiff(; chunksize = 8))))

Error & Stacktrace ⚠️

Please submit a bug report with steps to reproduce this fault, and any error messages that follow (in their entirety). Thanks.
Exception: EXCEPTION_ACCESS_VIOLATION at 0x7ff8367b728f -- gc_setmark_pool_ at C:/workdir/src\gc.c:874 [inlined]
gc_setmark_pool at C:/workdir/src\gc.c:893 [inlined]
gc_setmark at C:/workdir/src\gc.c:900 [inlined]
gc_setmark at C:/workdir/src\gc.c:896 [inlined]
gc_mark_outrefs at C:/workdir/src\gc.c:2618 [inlined]
gc_mark_loop_serial_ at C:/workdir/src\gc.c:2691
in expression starting at REPL[10]:1
gc_setmark_pool_ at C:/workdir/src\gc.c:874 [inlined]
gc_setmark_pool at C:/workdir/src\gc.c:893 [inlined]
gc_setmark at C:/workdir/src\gc.c:900 [inlined]
gc_setmark at C:/workdir/src\gc.c:896 [inlined]
gc_mark_outrefs at C:/workdir/src\gc.c:2618 [inlined]
gc_mark_loop_serial_ at C:/workdir/src\gc.c:2691
gc_mark_loop_serial at C:/workdir/src\gc.c:2714
gc_mark_loop at C:/workdir/src\gc.c:2827 [inlined]
_jl_gc_collect at C:/workdir/src\gc.c:3151
ijl_gc_collect at C:/workdir/src\gc.c:3448
maybe_collect at C:/workdir/src\gc.c:935 [inlined]
jl_gc_pool_alloc_inner at C:/workdir/src\gc.c:1291
jl_gc_pool_alloc_noinline at C:/workdir/src\gc.c:1348 [inlined]
jl_gc_alloc_ at C:/workdir/src\julia_internal.h:475 [inlined]
jl_gc_alloc at C:/workdir/src\gc.c:3500
_new_array_ at C:/workdir/src\array.c:134
_new_array at C:/workdir/src\array.c:198 [inlined]
ijl_alloc_array_1d at C:/workdir/src\array.c:436
Array at .\boot.jl:475 [inlined]
InstructionStream at .\compiler/ssair\ir.jl:196
NewNodeStream at .\compiler/ssair\ir.jl:278 [inlined]
NewNodeStream at .\compiler/ssair\ir.jl:278 [inlined]
IRCode at .\compiler/ssair\ir.jl:343 [inlined]
inflate_ir! at .\compiler/ssair\legacy.jl:45
inflate_ir! at .\compiler/ssair\legacy.jl:14
inflate_ir at .\compiler/ssair\legacy.jl:56 [inlined]
retrieve_ir_for_inlining at .\compiler/ssair\inlining.jl:988 [inlined]
#resolve_todo#457 at .\compiler/ssair\inlining.jl:908
resolve_todo at .\compiler/ssair\inlining.jl:871 [inlined]
#handle_const_prop_result!#465 at .\compiler/ssair\inlining.jl:1483
handle_const_prop_result! at .\compiler/ssair\inlining.jl:1473 [inlined]
#handle_any_const_result!#459 at .\compiler/ssair\inlining.jl:1326
handle_any_const_result! at .\compiler/ssair\inlining.jl:1312 [inlined]
compute_inlining_cases at .\compiler/ssair\inlining.jl:1404
handle_call! at .\compiler/ssair\inlining.jl:1450 [inlined]
assemble_inline_todo! at .\compiler/ssair\inlining.jl:1698
ssa_inlining_pass! at .\compiler/ssair\inlining.jl:78 [inlined]
run_passes at .\compiler\optimize.jl:489
run_passes at .\compiler\optimize.jl:504 [inlined]
optimize at .\compiler\optimize.jl:453 [inlined]
_typeinf at .\compiler\typeinfer.jl:272
typeinf at .\compiler\typeinfer.jl:216
typeinf_edge at .\compiler\typeinfer.jl:930
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618
abstract_eval_basic_statement at .\compiler\abstractinterpretation.jl:2883
typeinf_local at .\compiler\abstractinterpretation.jl:3092
typeinf_nocycle at .\compiler\abstractinterpretation.jl:3180
_typeinf at .\compiler\typeinfer.jl:247
typeinf at .\compiler\typeinfer.jl:216
typeinf_edge at .\compiler\typeinfer.jl:930
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618
abstract_eval_basic_statement at .\compiler\abstractinterpretation.jl:2907
typeinf_local at .\compiler\abstractinterpretation.jl:3092
typeinf_nocycle at .\compiler\abstractinterpretation.jl:3180
_typeinf at .\compiler\typeinfer.jl:247
typeinf at .\compiler\typeinfer.jl:216
typeinf_edge at .\compiler\typeinfer.jl:930
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618
abstract_eval_basic_statement at .\compiler\abstractinterpretation.jl:2907
typeinf_local at .\compiler\abstractinterpretation.jl:3092
typeinf_nocycle at .\compiler\abstractinterpretation.jl:3180
_typeinf at .\compiler\typeinfer.jl:247
typeinf at .\compiler\typeinfer.jl:216
typeinf_ext at .\compiler\typeinfer.jl:1051
typeinf_ext_toplevel at .\compiler\typeinfer.jl:1082
typeinf_ext_toplevel at .\compiler\typeinfer.jl:1078
jfptr_typeinf_ext_toplevel_38862.1 at D:\Julia\Julia-1.10.0-rc1\lib\julia\sys.dll (unknown line)
_jl_invoke at C:/workdir/src\gf.c:2892 [inlined]
ijl_apply_generic at C:/workdir/src\gf.c:3074 [inlined]
jl_apply at C:/workdir/src\julia.h:1976 [inlined]
jl_type_infer at C:/workdir/src\gf.c:394
jl_generate_fptr_impl at C:/workdir/src\jitlayers.cpp:502
jl_compile_method_internal at C:/workdir/src\gf.c:2478
jl_compile_method_internal at C:/workdir/src\gf.c:2370 [inlined]
_jl_invoke at C:/workdir/src\gf.c:2884 [inlined]
ijl_apply_generic at C:/workdir/src\gf.c:3074
GeneratedFunctionStub at .\boot.jl:600
jl_call_staged at C:/workdir/src\method.c:540
ijl_code_for_staged at C:/workdir/src\method.c:593
get_staged at .\compiler\utilities.jl:123
retrieve_code_info at .\compiler\utilities.jl:135 [inlined]
InferenceState at .\compiler\inferencestate.jl:430
typeinf_edge at .\compiler\typeinfer.jl:920
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618
abstract_eval_basic_statement at .\compiler\abstractinterpretation.jl:2907
typeinf_local at .\compiler\abstractinterpretation.jl:3092
typeinf_nocycle at .\compiler\abstractinterpretation.jl:3180
_typeinf at .\compiler\typeinfer.jl:247
typeinf at .\compiler\typeinfer.jl:216
typeinf_edge at .\compiler\typeinfer.jl:930
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618
abstract_eval_basic_statement at .\compiler\abstractinterpretation.jl:2883
typeinf_local at .\compiler\abstractinterpretation.jl:3092
typeinf_nocycle at .\compiler\abstractinterpretation.jl:3180
_typeinf at .\compiler\typeinfer.jl:247
typeinf at .\compiler\typeinfer.jl:216
typeinf_edge at .\compiler\typeinfer.jl:930
abstract_call_method at .\compiler\abstractinterpretation.jl:629
abstract_call_gf_by_type at .\compiler\abstractinterpretation.jl:95
abstract_call_known at .\compiler\abstractinterpretation.jl:2083
abstract_call at .\compiler\abstractinterpretation.jl:2165
abstract_call at .\compiler\abstractinterpretation.jl:2158
abstract_call at .\compiler\abstractinterpretation.jl:2348
abstract_eval_call at .\compiler\abstractinterpretation.jl:2364
abstract_eval_statement_expr at .\compiler\abstractinterpretation.jl:2374
abstract_eval_statement at .\compiler\abstractinterpretation.jl:2618

Environment (please complete the following information):

[764a87c0] BoundaryValueDiffEq v5.8.0
  [47edcb42] ADTypes v1.3.0
  [7d9f7c33] Accessors v0.1.36
  [79e6a3ab] Adapt v4.0.4
  [ec485272] ArnoldiMethod v0.4.0
  [4fba245c] ArrayInterface v7.11.0
  [4c555306] ArrayLayouts v1.9.4
  [aae01518] BandedMatrices v1.7.2
  [62783981] BitTwiddlingConvenienceFunctions v0.1.5
  [764a87c0] BoundaryValueDiffEq v5.8.0
  [2a0fbf3d] CPUSummary v0.2.5
  [d360d2e6] ChainRulesCore v1.24.0
  [fb6a15b2] CloseOpenIntervals v0.1.12
  [38540f10] CommonSolve v0.2.4
  [bbf7d656] CommonSubexpressions v0.3.0
  [34da2185] Compat v4.15.0
  [a33af91c] CompositionsBase v0.1.2
  [2569d6c7] ConcreteStructs v0.2.3
  [187b0558] ConstructionBase v1.5.5
  [adafc99b] CpuId v0.3.1
  [9a962f9c] DataAPI v1.16.0
  [864edb3b] DataStructures v0.18.20
  [e2d170a0] DataValueInterfaces v1.0.0
  [2b5f629d] DiffEqBase v6.151.2
  [163ba53b] DiffResults v1.1.0
  [b552c78f] DiffRules v1.15.1
  [a0c0ee7d] DifferentiationInterface v0.5.5
  [ffbed154] DocStringExtensions v0.9.3
  [4e289a0a] EnumX v1.0.4
  [f151be2c] EnzymeCore v0.7.4
  [d4d017d3] ExponentialUtilities v1.26.1
  [e2ba6199] ExprTools v0.1.10
  [9d29842c] FastAlmostBandedMatrices v0.1.2
  [7034ab61] FastBroadcast v0.3.3
  [9aa1b823] FastClosures v0.3.2
  [29a986be] FastLapackInterface v2.0.4
  [1a297f60] FillArrays v1.11.0
  [6a86dc24] FiniteDiff v2.23.1
  [f6369f11] ForwardDiff v0.10.36
  [069b7b12] FunctionWrappers v1.1.3
  [77dc65aa] FunctionWrappersWrappers v0.1.3
  [46192b85] GPUArraysCore v0.1.6
  [c145ed77] GenericSchur v0.5.4
  [86223c79] Graphs v1.11.1
  [3e5b6fbb] HostCPUFeatures v0.1.16
  [615f187c] IfElse v0.1.1
  [d25df0c9] Inflate v0.1.5
  [3587e190] InverseFunctions v0.1.14
  [92d709cd] IrrationalConstants v0.2.2
  [82899510] IteratorInterfaceExtensions v1.0.0
  [692b3bcd] JLLWrappers v1.5.0
  [ef3ab10e] KLU v0.6.0
  [ba0b0d4f] Krylov v0.9.6
  [10f19ff3] LayoutPointers v0.1.15
  [5078a376] LazyArrays v2.0.4
  [d3d80556] LineSearches v7.2.0
  [7ed4a6bd] LinearSolve v2.30.1
  [2ab3a3ac] LogExpFunctions v0.3.28
  [bdcacae8] LoopVectorization v0.12.170
  [1914dd2f] MacroTools v0.5.13
  [d125e4d3] ManualMemory v0.1.8
  [a3b82374] MatrixFactorizations v3.0.0
  [bb5d69b7] MaybeInplace v0.1.3
  [46d2c3a1] MuladdMacro v0.2.4
  [d41bc354] NLSolversBase v7.8.3
  [77ba4419] NaNMath v1.0.2
  [8913a72c] NonlinearSolve v3.12.4
  [6fe1bfb0] OffsetArrays v1.14.0
  [bac558e1] OrderedCollections v1.6.3
  [1dea7af3] OrdinaryDiffEq v6.82.0
  [65ce6f38] PackageExtensionCompat v1.0.2
  [d96e819e] Parameters v0.12.3
  [f517fe37] Polyester v0.7.14
  [1d0040c9] PolyesterWeave v0.2.1
  [d236fae5] PreallocationTools v0.4.22
  [aea7be01] PrecompileTools v1.2.1
  [21216c6a] Preferences v1.4.3
  [3cdcf5f2] RecipesBase v1.3.4
  [731186ca] RecursiveArrayTools v3.23.1
  [f2c3362d] RecursiveFactorization v0.2.23
  [189a3867] Reexport v1.2.2
  [ae029012] Requires v1.3.0
  [7e49a35a] RuntimeGeneratedFunctions v0.5.13
  [94e857df] SIMDTypes v0.1.0
  [476501e8] SLEEFPirates v0.6.42
  [0bca4576] SciMLBase v2.41.1
  [c0aeaf25] SciMLOperators v0.3.8
  [53ae85a6] SciMLStructures v1.3.0
  [efcf1570] Setfield v1.1.1
  [727e6d20] SimpleNonlinearSolve v1.10.0
  [699a6c99] SimpleTraits v0.9.4
  [ce78b400] SimpleUnPack v1.1.0
  [47a9eef4] SparseDiffTools v2.19.0
  [0a514795] SparseMatrixColorings v0.3.3
  [e56a9233] Sparspak v0.3.9
  [276daf66] SpecialFunctions v2.4.0
  [aedffcd0] Static v0.8.10
  [0d7ed370] StaticArrayInterface v1.5.0
  [90137ffa] StaticArrays v1.9.5
  [1e83bf80] StaticArraysCore v1.4.3
  [7792a7ef] StrideArraysCore v0.5.6
  [2efcf032] SymbolicIndexingInterface v0.3.22
  [3783bdb8] TableTraits v1.0.1
  [bd369af6] Tables v1.11.1
  [8290d209] ThreadingUtilities v0.5.2
  [a759f4b9] TimerOutputs v0.5.24
  [d5829a12] TriangularSolve v0.2.0
  [410a4b4d] Tricks v0.1.8
  [781d530d] TruncatedStacktraces v1.4.0
  [3a884ed6] UnPack v1.0.2
  [3d5dd08c] VectorizationBase v0.21.68
  [19fa3120] VertexSafeGraphs v0.2.0
  [1d5cc7b8] IntelOpenMP_jll v2024.1.0+0
  [856f044c] MKL_jll v2024.1.0+0
  [efe28fd5] OpenSpecFun_jll v0.5.5+0
  [1317d2d5] oneTBB_jll v2021.12.0+0
  [0dad84c5] ArgTools v1.1.1
  [56f22d72] Artifacts
  [2a0f44e3] Base64
  [ade2ca70] Dates
  [8ba89e20] Distributed
  [f43a241f] Downloads v1.6.0
  [7b1f6079] FileWatching
  [9fa8497b] Future
  [b77e0a4c] InteractiveUtils
  [4af54fe1] LazyArtifacts
  [b27032c2] LibCURL v0.6.4
  [76f85450] LibGit2
  [8f399da3] Libdl
  [37e2e46d] LinearAlgebra
  [56ddb016] Logging
  [d6f4376e] Markdown
  [a63ad114] Mmap
  [ca575930] NetworkOptions v1.2.0
  [44cfe95a] Pkg v1.10.0
  [de0858da] Printf
  [3fa0cd96] REPL
  [9a3f8284] Random
  [ea8e919c] SHA v0.7.0
  [9e88b42a] Serialization
  [1a1011a3] SharedArrays
  [6462fe0b] Sockets
  [2f01184e] SparseArrays v1.10.0
  [10745b16] Statistics v1.10.0
  [4607b0f0] SuiteSparse
  [fa267f1f] TOML v1.0.3
  [a4e569a6] Tar v1.10.0
  [8dfed614] Test
  [cf7118a7] UUIDs
  [4ec0a83e] Unicode
  [e66e0078] CompilerSupportLibraries_jll v1.0.5+1
  [deac9b47] LibCURL_jll v8.4.0+0
  [e37daf67] LibGit2_jll v1.6.4+0
  [29816b5a] LibSSH2_jll v1.11.0+1
  [c8ffd9c3] MbedTLS_jll v2.28.2+1
  [14a3606d] MozillaCACerts_jll v2023.1.10
  [4536629a] OpenBLAS_jll v0.3.23+2
  [05823500] OpenLibm_jll v0.8.1+2
  [bea87d4a] SuiteSparse_jll v7.2.1+1
  [83775a58] Zlib_jll v1.2.13+1
  [8e850b90] libblastrampoline_jll v5.8.0+1
  [8e850ede] nghttp2_jll v1.52.0+1
  [3f19e933] p7zip_jll v17.4.0+2
Julia Version 1.10.0-rc1
Commit 5aaa948543 (2023-11-03 07:44 UTC)
Build Info:
  Official https://julialang.org/ release
Platform Info:
  OS: Windows (x86_64-w64-mingw32)
  CPU: 8 × Intel(R) Core(TM) i5-10210U CPU @ 1.60GHz
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-15.0.7 (ORCJIT, skylake)
  Threads: 1 on 8 virtual cores
Environment:
  JULIA_DEPOT_PATH = D:\Julia\.julia
  JULIA_PKGDIR = D:\Julia\JuliaPKG
  JULIA_PKG_SERVER = us-west.pkg.julialang.org
ChrisRackauckas commented 2 weeks ago

Try Julia v1.10.4?

ChrisRackauckas commented 2 weeks ago

This seems to be an issue with the RC. v1.10.4 is fine.