Package evaluation to test QuantumCollocation on Julia 1.14.0-DEV.3141 (fbd24f51d0*) started at 2026-09-11T03:38:01.449 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 15.4s ################################################################################ # Installation # Installing QuantumCollocation... Resolving package versions... Updating `~/.julia/environments/v1.14/Project.toml` [0dc23a59] + QuantumCollocation v0.10.2 Updating `~/.julia/environments/v1.14/Manifest.toml` [47edcb42] + ADTypes v1.24.0 ⌃ [c29ec348] + AbstractDifferentiation v0.4.4 [7d9f7c33] + Accessors v0.1.45 [79e6a3ab] + Adapt v4.7.0 [66dad0bd] + AliasTables v1.1.3 [4fba245c] + ArrayInterface v7.30.1 [4c555306] + ArrayLayouts v1.12.2 [62783981] + BitTwiddlingConvenienceFunctions v0.1.6 ⌃ [70df07ce] + BracketingNonlinearSolve v1.12.1 [2a0fbf3d] + CPUSummary v0.2.7 [d360d2e6] + ChainRulesCore v1.26.1 [0b6fb165] + ChunkCodecCore v1.0.2 [4c0bbee4] + ChunkCodecLibZlib v1.1.0 [55437552] + ChunkCodecLibZstd v1.0.0 [fb6a15b2] + CloseOpenIntervals v0.1.13 [523fee87] + CodecBzip2 v0.8.5 [944b1d66] + CodecZlib v0.7.9 [38540f10] + CommonSolve v0.2.14 [bbf7d656] + CommonSubexpressions v0.3.1 [f70d9fcc] + CommonWorldInvalidations v1.2.2 [34da2185] + Compat v4.18.1 [a33af91c] + CompositionsBase v0.1.2 [95dc2771] + ComputePipeline v0.1.8 [2569d6c7] + ConcreteStructs v0.2.8 [187b0558] + ConstructionBase v1.6.0 [adafc99b] + CpuId v0.3.1 [a8cc5b0e] + Crayons v4.2.0 [9a962f9c] + DataAPI v1.16.0 ⌅ [82cc6244] + DataInterpolations v8.10.0 [864edb3b] + DataStructures v0.19.6 [e2d170a0] + DataValueInterfaces v1.0.0 ⌅ [2b5f629d] + DiffEqBase v6.218.0 [163ba53b] + DiffResults v1.1.0 [b552c78f] + DiffRules v1.16.0 [a0c0ee7d] + DifferentiationInterface v0.7.21 ⌅ [c823fa1f] + DirectTrajOpt v0.8.11 [31c24e10] + Distributions v0.25.131 [ffbed154] + DocStringExtensions v0.9.5 [4e289a0a] + EnumX v1.0.7 [f151be2c] + EnzymeCore v0.8.21 [e24c0720] + ExponentialAction v0.2.11 ⌃ [d4d017d3] + ExponentialUtilities v1.31.0 [e2ba6199] + ExprTools v0.1.11 [55351af7] + ExproniconLite v0.10.14 [7034ab61] + FastBroadcast v1.4.0 [9aa1b823] + FastClosures v0.3.2 [a4df4552] + FastPower v1.5.0 [5789e2e9] + FileIO v1.20.0 [1a297f60] + FillArrays v1.17.0 ⌅ [64ca27bc] + FindFirstFunctions v1.8.0 [6a86dc24] + FiniteDiff v2.33.0 [f6369f11] + ForwardDiff v1.4.5 [069b7b12] + FunctionWrappers v1.1.3 [77dc65aa] + FunctionWrappersWrappers v1.13.0 [46192b85] + GPUArraysCore v0.2.0 [a0844989] + Gamma v1.2.0 [c145ed77] + GenericSchur v0.5.8 [076d061b] + HashArrayMappedTries v0.2.0 [34004b35] + HypergeometricFunctions v0.3.30 [615f187c] + IfElse v0.1.1 [3587e190] + InverseFunctions v0.1.17 [b6b21f68] + Ipopt v1.16.0 [92d709cd] + IrrationalConstants v0.2.6 [82899510] + IteratorInterfaceExtensions v1.0.0 [033835bb] + JLD2 v0.6.6 [692b3bcd] + JLLWrappers v1.8.0 [682c06a0] + JSON v1.8.0 [ae98c720] + Jieko v0.2.1 [b964fa9f] + LaTeXStrings v1.4.1 [10f19ff3] + LayoutPointers v0.1.17 [5078a376] + LazyArrays v2.13.0 [2ab3a3ac] + LogExpFunctions v1.0.1 [e6f89c97] + LoggingExtras v1.2.0 [1914dd2f] + MacroTools v0.5.16 [d125e4d3] + ManualMemory v0.1.8 [b8f27783] + MathOptInterface v1.53.0 [a3b82374] + MatrixFactorizations v3.1.3 [bb5d69b7] + MaybeInplace v0.1.8 [e1d29d7a] + Missings v1.2.0 [2e0e35c7] + Moshi v0.3.12 [46d2c3a1] + MuladdMacro v0.2.7 [d8a4904e] + MutableArithmetics v1.8.0 [77ba4419] + NaNMath v1.1.4 ⌅ [538bc3a1] + NamedTrajectories v0.8.4 ⌃ [be0214bd] + NonlinearSolveBase v2.30.3 [510215fc] + Observables v0.5.5 ⌅ [bac558e1] + OrderedCollections v1.8.2 ⌅ [bbf590c4] + OrdinaryDiffEqCore v3.33.1 ⌅ [521117fe] + OrdinaryDiffEqLinear v1.12.0 ⌅ [b1df2697] + OrdinaryDiffEqTsit5 v1.12.0 [90014a1f] + PDMats v0.11.41 [69de0a69] + Parsers v3.0.0 [5a402ddf] + PiccoloQuantumObjects v0.10.2 [f517fe37] + Polyester v0.7.19 [1d0040c9] + PolyesterWeave v0.2.2 [d236fae5] + PreallocationTools v1.7.1 [aea7be01] + PrecompileTools v1.3.4 [21216c6a] + Preferences v1.5.2 [08abe8d2] + PrettyTables v3.4.8 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.4.0 [1fd47b50] + QuadGK v2.11.3 [0dc23a59] + QuantumCollocation v0.10.2 [3cdcf5f2] + RecipesBase v1.3.4 ⌅ [731186ca] + RecursiveArrayTools v3.54.0 [189a3867] + Reexport v1.2.2 [ae029012] + Requires v1.3.1 [79098fc4] + Rmath v0.9.0 [f2b01f46] + Roots v3.0.8 [7e49a35a] + RuntimeGeneratedFunctions v0.5.26 [94e857df] + SIMDTypes v0.1.0 ⌅ [0bca4576] + SciMLBase v2.155.2 ⌃ [19f34311] + SciMLJacobianOperators v0.1.17 ⌅ [a6db7da4] + SciMLLogging v1.10.1 [c0aeaf25] + SciMLOperators v1.30.0 [431bcebd] + SciMLPublic v1.3.0 [53ae85a6] + SciMLStructures v1.10.5 [7e506255] + ScopedValues v1.6.2 [efcf1570] + Setfield v1.1.2 [a2af1166] + SortingAlgorithms v1.2.3 [276daf66] + SpecialFunctions v2.9.0 [aedffcd0] + Static v1.4.6 [0d7ed370] + StaticArrayInterface v1.10.0 [90137ffa] + StaticArrays v1.9.20 [1e83bf80] + StaticArraysCore v1.4.4 [10745b16] + Statistics v1.11.5 [82ae8749] + StatsAPI v1.8.0 [2913bbd2] + StatsBase v0.34.13 [4c63d2b9] + StatsFuns v2.2.1 [7792a7ef] + StrideArraysCore v0.5.9 ⌅ [892a3eda] + StringManipulation v0.5.0 [ec057cc2] + StructUtils v2.8.5 [2efcf032] + SymbolicIndexingInterface v0.3.55 [3783bdb8] + TableTraits v1.0.1 [bd369af6] + Tables v1.14.0 [f8b46487] + TestItemRunner v1.3.2 [1c621080] + TestItems v1.1.0 [8290d209] + ThreadingUtilities v0.5.6 ⌅ [a759f4b9] + TimerOutputs v0.5.29 [6dad8b7f] + TrajectoryIndexingUtils v0.1.0 [3bb67fe8] + TranscodingStreams v0.11.3 [781d530d] + TruncatedStacktraces v1.4.0 [ae81ac8f] + ASL_jll v0.1.5+0 [6e34b625] + Bzip2_jll v1.0.9+0 [e33a78d0] + Hwloc_jll v2.14.0+0 ⌅ [9cc047cb] + Ipopt_jll v300.1400.1902+0 [94ce4f54] + Libiconv_jll v1.18.0+0 [d00139f3] + METIS_jll v5.1.3+0 [d7ed1dd3] + MUMPS_seq_jll v500.900.100+0 [656ef2d0] + OpenBLAS32_jll v0.3.34+0 [efe28fd5] + OpenSpecFun_jll v0.5.6+0 [f50d1b31] + Rmath_jll v0.5.2+0 [319450e9] + SPRAL_jll v2025.9.18+1 ⌅ [02c8fc9c] + XML2_jll v2.13.9+0 [a65dc6b1] + Xorg_libpciaccess_jll v0.19.0+0 [0dad84c5] + ArgTools v1.2.0 [56f22d72] + Artifacts v1.11.0 [2a0f44e3] + Base64 v1.11.0 [ade2ca70] + Dates v1.11.0 [8ba89e20] + Distributed v1.12.0 [f43a241f] + Downloads v1.7.0 [7b1f6079] + FileWatching v1.11.0 [9fa8497b] + Future v1.11.0 [b77e0a4c] + InteractiveUtils v1.11.0 [ac6e5ff7] + JuliaSyntaxHighlighting v1.13.0 [b27032c2] + LibCURL v1.0.0 [76f85450] + LibGit2 v1.11.0 [8f399da3] + Libdl v1.11.0 [37e2e46d] + LinearAlgebra v1.14.0 [56ddb016] + Logging v1.11.0 [d6f4376e] + Markdown v1.11.0 [a63ad114] + Mmap v1.11.0 [ca575930] + NetworkOptions v1.3.0 [44cfe95a] + Pkg v1.14.0 [de0858da] + Printf v1.11.0 [3fa0cd96] + REPL v1.11.0 [9a3f8284] + Random v1.11.0 [ea8e919c] + SHA v1.13.0 [9e88b42a] + Serialization v1.11.0 [6462fe0b] + Sockets v1.11.0 [2f01184e] + SparseArrays v1.13.0 [f489334b] + StyledStrings v1.13.0 [4607b0f0] + SuiteSparse [fa267f1f] + TOML v1.0.3 [a4e569a6] + Tar v1.10.0 [8dfed614] + Test v1.11.0 [cf7118a7] + UUIDs v1.11.0 [4ec0a83e] + Unicode v1.11.0 [e66e0078] + CompilerSupportLibraries_jll v1.5.7+0 [deac9b47] + LibCURL_jll v8.22.0+0 [e37daf67] + LibGit2_jll v1.9.7+0 [29816b5a] + LibSSH2_jll v1.11.104+0 [14a3606d] + MozillaCACerts_jll v2026.8.13 [4536629a] + OpenBLAS_jll v0.3.34+0 [05823500] + OpenLibm_jll v0.8.8+0 [458c3c95] + OpenSSL_jll v3.5.8+0 [efcefdf7] + PCRE2_jll v10.48.0+0 [bea87d4a] + SuiteSparse_jll v7.10.1+0 [83775a58] + Zlib_jll v1.3.2+0 [3161d3a3] + Zstd_jll v1.5.7+1 [8e850b90] + libblastrampoline_jll v5.15.0+0 [8e850ede] + nghttp2_jll v1.70.0+0 [3f19e933] + p7zip_jll v17.8.2+0 Info Packages marked with ⌃ and ⌅ have new versions available. Those with ⌃ may be upgradable, but those with ⌅ are restricted by compatibility constraints from upgrading. To see why use `status --outdated -m` Installation completed after 4.84s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling project... 10.2 s ✓ NamedTrajectories → InterpolationsExt 25.9 s ✓ DirectTrajOpt 27.6 s ✓ PiccoloQuantumObjects 36.8 s ✓ QuantumCollocation 4 dependencies successfully precompiled in 103 seconds. 274 already precompiled. Precompilation completed after 117.04s ################################################################################ # Testing # Testing QuantumCollocation Status `/tmp/jl_uK67F4/Project.toml` ⌅ [c823fa1f] DirectTrajOpt v0.8.11 [31c24e10] Distributions v0.25.131 [e24c0720] ExponentialAction v0.2.11 ⌅ [538bc3a1] NamedTrajectories v0.8.4 [5a402ddf] PiccoloQuantumObjects v0.10.2 [0dc23a59] QuantumCollocation v0.10.2 [189a3867] Reexport v1.2.2 [f8b46487] TestItemRunner v1.3.2 [1c621080] TestItems v1.1.0 [6dad8b7f] TrajectoryIndexingUtils v0.1.0 [8f399da3] Libdl v1.11.0 [37e2e46d] LinearAlgebra v1.14.0 [9a3f8284] Random v1.11.0 [2f01184e] SparseArrays v1.13.0 [8dfed614] Test v1.11.0 Status `/tmp/jl_uK67F4/Manifest.toml` [47edcb42] ADTypes v1.24.0 ⌃ [c29ec348] AbstractDifferentiation v0.4.4 [7d9f7c33] Accessors v0.1.45 [79e6a3ab] Adapt v4.7.0 [66dad0bd] AliasTables v1.1.3 [4fba245c] ArrayInterface v7.30.1 [4c555306] ArrayLayouts v1.12.2 [62783981] BitTwiddlingConvenienceFunctions v0.1.6 ⌃ [70df07ce] BracketingNonlinearSolve v1.12.1 [2a0fbf3d] CPUSummary v0.2.7 [d360d2e6] ChainRulesCore v1.26.1 [0b6fb165] ChunkCodecCore v1.0.2 [4c0bbee4] ChunkCodecLibZlib v1.1.0 [55437552] ChunkCodecLibZstd v1.0.0 [fb6a15b2] CloseOpenIntervals v0.1.13 [523fee87] CodecBzip2 v0.8.5 [944b1d66] CodecZlib v0.7.9 [38540f10] CommonSolve v0.2.14 [bbf7d656] CommonSubexpressions v0.3.1 [f70d9fcc] CommonWorldInvalidations v1.2.2 [34da2185] Compat v4.18.1 [a33af91c] CompositionsBase v0.1.2 [95dc2771] ComputePipeline v0.1.8 [2569d6c7] ConcreteStructs v0.2.8 [187b0558] ConstructionBase v1.6.0 [adafc99b] CpuId v0.3.1 [a8cc5b0e] Crayons v4.2.0 [9a962f9c] DataAPI v1.16.0 ⌅ [82cc6244] DataInterpolations v8.10.0 [864edb3b] DataStructures v0.19.6 [e2d170a0] DataValueInterfaces v1.0.0 ⌅ [2b5f629d] DiffEqBase v6.218.0 [163ba53b] DiffResults v1.1.0 [b552c78f] DiffRules v1.16.0 [a0c0ee7d] DifferentiationInterface v0.7.21 ⌅ [c823fa1f] DirectTrajOpt v0.8.11 [31c24e10] Distributions v0.25.131 [ffbed154] DocStringExtensions v0.9.5 [4e289a0a] EnumX v1.0.7 [f151be2c] EnzymeCore v0.8.21 [e24c0720] ExponentialAction v0.2.11 ⌃ [d4d017d3] ExponentialUtilities v1.31.0 [e2ba6199] ExprTools v0.1.11 [55351af7] ExproniconLite v0.10.14 [7034ab61] FastBroadcast v1.4.0 [9aa1b823] FastClosures v0.3.2 [a4df4552] FastPower v1.5.0 [5789e2e9] FileIO v1.20.0 [1a297f60] FillArrays v1.17.0 ⌅ [64ca27bc] FindFirstFunctions v1.8.0 [6a86dc24] FiniteDiff v2.33.0 [f6369f11] ForwardDiff v1.4.5 [069b7b12] FunctionWrappers v1.1.3 [77dc65aa] FunctionWrappersWrappers v1.13.0 [46192b85] GPUArraysCore v0.2.0 [a0844989] Gamma v1.2.0 [c145ed77] GenericSchur v0.5.8 [076d061b] HashArrayMappedTries v0.2.0 [34004b35] HypergeometricFunctions v0.3.30 [615f187c] IfElse v0.1.1 [3587e190] InverseFunctions v0.1.17 [b6b21f68] Ipopt v1.16.0 [92d709cd] IrrationalConstants v0.2.6 [82899510] IteratorInterfaceExtensions v1.0.0 [033835bb] JLD2 v0.6.6 [692b3bcd] JLLWrappers v1.8.0 [682c06a0] JSON v1.8.0 [ae98c720] Jieko v0.2.1 [b964fa9f] LaTeXStrings v1.4.1 [10f19ff3] LayoutPointers v0.1.17 [5078a376] LazyArrays v2.13.0 [2ab3a3ac] LogExpFunctions v1.0.1 [e6f89c97] LoggingExtras v1.2.0 [1914dd2f] MacroTools v0.5.16 [d125e4d3] ManualMemory v0.1.8 [b8f27783] MathOptInterface v1.53.0 [a3b82374] MatrixFactorizations v3.1.3 [bb5d69b7] MaybeInplace v0.1.8 [e1d29d7a] Missings v1.2.0 [2e0e35c7] Moshi v0.3.12 [46d2c3a1] MuladdMacro v0.2.7 [d8a4904e] MutableArithmetics v1.8.0 [77ba4419] NaNMath v1.1.4 ⌅ [538bc3a1] NamedTrajectories v0.8.4 ⌃ [be0214bd] NonlinearSolveBase v2.30.3 [510215fc] Observables v0.5.5 ⌅ [bac558e1] OrderedCollections v1.8.2 ⌅ [bbf590c4] OrdinaryDiffEqCore v3.33.1 ⌅ [521117fe] OrdinaryDiffEqLinear v1.12.0 ⌅ [b1df2697] OrdinaryDiffEqTsit5 v1.12.0 [90014a1f] PDMats v0.11.41 [69de0a69] Parsers v3.0.0 [5a402ddf] PiccoloQuantumObjects v0.10.2 [f517fe37] Polyester v0.7.19 [1d0040c9] PolyesterWeave v0.2.2 [d236fae5] PreallocationTools v1.7.1 [aea7be01] PrecompileTools v1.3.4 [21216c6a] Preferences v1.5.2 [08abe8d2] PrettyTables v3.4.8 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.4.0 [1fd47b50] QuadGK v2.11.3 [0dc23a59] QuantumCollocation v0.10.2 [3cdcf5f2] RecipesBase v1.3.4 ⌅ [731186ca] RecursiveArrayTools v3.54.0 [189a3867] Reexport v1.2.2 [ae029012] Requires v1.3.1 [79098fc4] Rmath v0.9.0 [f2b01f46] Roots v3.0.8 [7e49a35a] RuntimeGeneratedFunctions v0.5.26 [94e857df] SIMDTypes v0.1.0 ⌅ [0bca4576] SciMLBase v2.155.2 ⌃ [19f34311] SciMLJacobianOperators v0.1.17 ⌅ [a6db7da4] SciMLLogging v1.10.1 [c0aeaf25] SciMLOperators v1.30.0 [431bcebd] SciMLPublic v1.3.0 [53ae85a6] SciMLStructures v1.10.5 [7e506255] ScopedValues v1.6.2 [efcf1570] Setfield v1.1.2 [a2af1166] SortingAlgorithms v1.2.3 [276daf66] SpecialFunctions v2.9.0 [aedffcd0] Static v1.4.6 [0d7ed370] StaticArrayInterface v1.10.0 [90137ffa] StaticArrays v1.9.20 [1e83bf80] StaticArraysCore v1.4.4 [10745b16] Statistics v1.11.5 [82ae8749] StatsAPI v1.8.0 [2913bbd2] StatsBase v0.34.13 [4c63d2b9] StatsFuns v2.2.1 [7792a7ef] StrideArraysCore v0.5.9 ⌅ [892a3eda] StringManipulation v0.5.0 [ec057cc2] StructUtils v2.8.5 [2efcf032] SymbolicIndexingInterface v0.3.55 [3783bdb8] TableTraits v1.0.1 [bd369af6] Tables v1.14.0 [f8b46487] TestItemRunner v1.3.2 [1c621080] TestItems v1.1.0 [8290d209] ThreadingUtilities v0.5.6 ⌅ [a759f4b9] TimerOutputs v0.5.29 [6dad8b7f] TrajectoryIndexingUtils v0.1.0 [3bb67fe8] TranscodingStreams v0.11.3 [781d530d] TruncatedStacktraces v1.4.0 [ae81ac8f] ASL_jll v0.1.5+0 [6e34b625] Bzip2_jll v1.0.9+0 [e33a78d0] Hwloc_jll v2.14.0+0 ⌅ [9cc047cb] Ipopt_jll v300.1400.1902+0 [94ce4f54] Libiconv_jll v1.18.0+0 [d00139f3] METIS_jll v5.1.3+0 [d7ed1dd3] MUMPS_seq_jll v500.900.100+0 [656ef2d0] OpenBLAS32_jll v0.3.34+0 [efe28fd5] OpenSpecFun_jll v0.5.6+0 [f50d1b31] Rmath_jll v0.5.2+0 [319450e9] SPRAL_jll v2025.9.18+1 ⌅ [02c8fc9c] XML2_jll v2.13.9+0 [a65dc6b1] Xorg_libpciaccess_jll v0.19.0+0 [0dad84c5] ArgTools v1.2.0 [56f22d72] Artifacts v1.11.0 [2a0f44e3] Base64 v1.11.0 [ade2ca70] Dates v1.11.0 [8ba89e20] Distributed v1.12.0 [f43a241f] Downloads v1.7.0 [7b1f6079] FileWatching v1.11.0 [9fa8497b] Future v1.11.0 [b77e0a4c] InteractiveUtils v1.11.0 [ac6e5ff7] JuliaSyntaxHighlighting v1.13.0 [b27032c2] LibCURL v1.0.0 [76f85450] LibGit2 v1.11.0 [8f399da3] Libdl v1.11.0 [37e2e46d] LinearAlgebra v1.14.0 [56ddb016] Logging v1.11.0 [d6f4376e] Markdown v1.11.0 [a63ad114] Mmap v1.11.0 [ca575930] NetworkOptions v1.3.0 [44cfe95a] Pkg v1.14.0 [de0858da] Printf v1.11.0 [3fa0cd96] REPL v1.11.0 [9a3f8284] Random v1.11.0 [ea8e919c] SHA v1.13.0 [9e88b42a] Serialization v1.11.0 [6462fe0b] Sockets v1.11.0 [2f01184e] SparseArrays v1.13.0 [f489334b] StyledStrings v1.13.0 [4607b0f0] SuiteSparse [fa267f1f] TOML v1.0.3 [a4e569a6] Tar v1.10.0 [8dfed614] Test v1.11.0 [cf7118a7] UUIDs v1.11.0 [4ec0a83e] Unicode v1.11.0 [e66e0078] CompilerSupportLibraries_jll v1.5.7+0 [deac9b47] LibCURL_jll v8.22.0+0 [e37daf67] LibGit2_jll v1.9.7+0 [29816b5a] LibSSH2_jll v1.11.104+0 [14a3606d] MozillaCACerts_jll v2026.8.13 [4536629a] OpenBLAS_jll v0.3.34+0 [05823500] OpenLibm_jll v0.8.8+0 [458c3c95] OpenSSL_jll v3.5.8+0 [efcefdf7] PCRE2_jll v10.48.0+0 [bea87d4a] SuiteSparse_jll v7.10.1+0 [83775a58] Zlib_jll v1.3.2+0 [3161d3a3] Zstd_jll v1.5.7+1 [8e850b90] libblastrampoline_jll v5.15.0+0 [8e850ede] nghttp2_jll v1.70.0+0 [3f19e933] p7zip_jll v17.8.2+0 Info Packages marked with ⌃ and ⌅ have new versions available. Those with ⌃ may be upgradable, but those with ⌅ are restricted by compatibility constraints from upgrading. Testing Running tests... added GlobalBoundsConstraint for :δ with bounds 0.5 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt ┌ Warning: `disable_sigint` no longer defers Ctrl-C: SIGINT is delivered as a cancellation of the current ^C scope and observed at cancellation points regardless of the sigatomic region this establishes. Shield a region from cancellation by scoping `Base.CANCEL_TOKEN => nothing` over it instead. │ caller = IpoptSolve(prob::Ipopt.IpoptProblem) at C_wrapper.jl:407 [inlined] └ @ Ipopt ~/.julia/packages/Ipopt/JN6dO/src/C_wrapper.jl:407 ****************************************************************************** This program contains Ipopt, a library for large-scale nonlinear optimization. Ipopt is released as open source code under the Eclipse Public License (EPL). For more information visit https://github.com/coin-or/Ipopt ****************************************************************************** constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}... final fidelity constraint: 0.95 minimum-time weight D: 100.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.KetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}... final fidelity constraint: 0.9 minimum-time weight D: 50.0 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}... final fidelity constraint: 0.95 minimum-time weight D: 100.0 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt ┌ Warning: Trajectory has timestep variable :Δt but no bounds on it. │ Adding default lower bound of 0 to prevent negative timesteps. │ │ Recommended: Add explicit bounds when creating the trajectory: │ NamedTrajectory(...; Δt_bounds=(min, max)) │ Example: │ NamedTrajectory(qtraj, N; Δt_bounds=(1e-3, 0.5)) │ │ Or use timesteps_all_equal=true in problem options to fix timesteps. └ @ DirectTrajOpt.Problems ~/.julia/packages/DirectTrajOpt/OwEJV/src/problems.jl:66 constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}... final fidelity constraint: 0.99 minimum-time weight D: 100.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.KetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}... final fidelity constraint: 0.99 minimum-time weight D: 100.0 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.SamplingTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}}... final fidelity constraint: 0.9 minimum-time weight D: 50.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.SamplingTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, PiccoloQuantumObjects.QuantumTrajectories.KetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}... final fidelity constraint: 0.85 minimum-time weight D: 30.0 constructing SmoothPulseProblem for MultiKetTrajectory (2 states)... ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 applying timesteps_all_equal constraint: Δt ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}}... final fidelity constraint: 0.9 minimum-time weight D: 50.0 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##155".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##155".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##155".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##155".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}... final fidelity constraint: 0.85 minimum-time weight D: 50.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.KetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##156".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##156".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##156".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##156".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}... final fidelity constraint: 0.85 minimum-time weight D: 50.0 constructing SmoothPulseProblem for MultiKetTrajectory (2 states)... ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 applying timesteps_all_equal constraint: Δt ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##157".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##157".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}}... final fidelity constraint: 0.8 minimum-time weight D: 50.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.SamplingTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, PiccoloQuantumObjects.QuantumTrajectories.KetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##158".H1)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##158".H1)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##158".H1)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##158".H1)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}... final fidelity constraint: 0.6 minimum-time weight D: 50.0 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 29666 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 26679 Total number of variables............................: 933 variables with only lower bounds: 50 variables with lower and upper bounds: 834 variables with only upper bounds: 0 Total number of equality constraints.................: 833 Total number of inequality constraints...............: 0 inequality constraints with only lower bounds: 0 inequality constraints with lower and upper bounds: 0 inequality constraints with only upper bounds: 0 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 0 1.9999965e+02 2.05e-02 1.00e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 1.9996912e+02 1.50e-02 3.92e+01 -0.3 5.01e-01 2.0 1.00e+00 2.67e-01f 1 2 1.9984571e+02 1.36e-02 9.68e+01 -4.0 5.30e-01 1.5 2.92e-01 9.08e-02h 1 3 1.9842898e+02 6.59e-02 4.25e+02 0.0 5.68e-01 1.9 1.00e+00 7.19e-01f 1 4 1.9654415e+02 3.99e-02 2.90e+03 0.8 5.57e-01 2.4 9.11e-01 3.94e-01f 1 5 1.9561665e+02 5.67e-03 5.48e+02 0.5 1.20e-01 2.8 9.99e-01 1.00e+00h 1 6 1.9273236e+02 1.06e-03 2.30e+02 0.3 6.19e-02 2.3 9.93e-01 1.00e+00f 1 7 1.7733543e+02 8.89e-03 1.46e+02 0.0 1.96e-01 1.8 9.97e-01 1.00e+00f 1 8 9.4246162e+01 3.23e-01 1.13e+02 -0.3 3.85e+00 1.4 1.31e-01 2.79e-01f 1 9 7.9813472e+01 3.62e-02 1.14e+02 0.3 6.26e-01 1.8 8.70e-01 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 7.4729921e+01 4.36e-03 1.84e+01 -0.2 1.16e-01 2.2 9.95e-01 1.00e+00f 1 Number of Iterations....: 10 (scaled) (unscaled) Objective...............: 7.4729921182783940e+01 7.4729921182783940e+01 Dual infeasibility......: 1.8437859625225624e+01 1.8437859625225624e+01 Constraint violation....: 4.3605052884085893e-03 4.3605052884085893e-03 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 7.1725292363239452e+00 7.1725292363239452e+00 Overall NLP error.......: 1.8437859625225624e+01 1.8437859625225624e+01 Number of objective function evaluations = 11 Number of objective gradient evaluations = 11 Number of equality constraint evaluations = 11 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 11 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 10 Total seconds in IPOPT = 4.335 EXIT: Maximum Number of Iterations Exceeded. constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 3 systems constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing MinimumTimeProblem from QuantumControlProblem{PiccoloQuantumObjects.QuantumTrajectories.SamplingTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, PiccoloQuantumObjects.QuantumTrajectories.UnitaryTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.ODESolution{ComplexF64, 3, Vector{Matrix{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Matrix{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Matrix{ComplexF64}}, Vector{Float64}, Vector{Vector{Matrix{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Matrix{ComplexF64}, Matrix{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}, Matrix{ComplexF64}}}}... final fidelity constraint: 0.9 minimum-time weight D: 50.0 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 9580 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 10762 Total number of variables............................: 591 variables with only lower bounds: 50 variables with lower and upper bounds: 392 variables with only upper bounds: 0 Total number of equality constraints.................: 490 Total number of inequality constraints...............: 0 inequality constraints with only lower bounds: 0 inequality constraints with lower and upper bounds: 0 inequality constraints with only upper bounds: 0 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 0 1.6167171e+01 2.96e+01 1.73e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 2.2731880e+01 8.91e+00 4.27e+02 -0.0 1.83e+01 - 1.09e-01 7.00e-01f 1 2 3.3654289e+01 2.83e+00 5.64e+02 0.2 5.34e+00 2.0 1.00e+00 6.83e-01f 1 3 3.4033818e+01 6.50e-03 3.11e+02 0.7 1.72e+00 1.5 9.49e-01 1.00e+00f 1 4 2.5997088e+01 1.39e-03 2.00e+01 0.1 1.21e-01 1.0 9.96e-01 1.00e+00f 1 5 5.5184170e+00 5.87e-03 7.38e+00 -0.3 3.74e-01 0.6 9.66e-01 1.00e+00f 1 6 1.3895263e+00 3.79e-03 1.88e+02 -0.5 2.79e-01 0.1 8.62e-01 1.00e+00f 1 7 3.9068827e-01 2.16e-03 4.73e+01 -0.8 9.01e-01 - 6.43e-01 2.50e-01f 3 8 1.6104250e+00 2.31e-03 1.94e+02 -1.1 1.83e-01 0.5 9.96e-01 1.00e+00h 1 9 7.8832818e-01 2.70e-03 1.90e+02 -1.7 4.08e-01 - 8.46e-01 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 7.3642736e-01 3.72e-05 2.10e+00 -1.7 2.18e-02 1.8 1.00e+00 1.00e+00h 1 11 5.0891791e-01 2.24e-05 3.88e-01 -3.1 1.66e-02 1.4 9.99e-01 1.00e+00f 1 12 1.4817906e-01 8.20e-05 2.45e-01 -4.0 3.15e-02 0.9 1.00e+00 1.00e+00f 1 13 5.6027303e-02 1.00e-04 1.93e+02 -4.0 3.59e-02 0.4 1.00e+00 1.00e+00f 1 14 1.4358570e-02 5.55e-05 9.65e+01 -4.0 1.90e-02 -0.1 1.00e+00 5.00e-01f 2 15 3.4099223e-02 4.34e-06 1.70e-01 -4.0 7.18e-03 1.3 1.00e+00 1.00e+00h 1 16 7.7144846e-03 7.53e-06 6.19e-02 -4.0 1.00e-02 0.8 1.00e+00 1.00e+00f 1 17 5.1015376e-03 7.55e-06 1.93e+02 -4.0 1.02e-02 0.3 1.00e+00 1.00e+00h 1 18 4.7231790e-03 1.32e-06 1.93e+02 -4.0 4.58e-03 -0.2 1.00e+00 1.00e+00h 1 19 1.4564086e-03 1.44e-06 2.85e-01 -4.0 3.86e-03 1.2 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 20 8.4046867e-04 5.34e-07 1.31e-02 -4.0 2.69e-03 0.7 1.00e+00 1.00e+00h 1 21 8.7085916e-04 4.20e-07 1.94e+02 -4.0 2.37e-03 0.2 1.00e+00 1.00e+00h 1 22 7.8591847e-04 5.90e-08 1.94e+02 -4.1 9.27e-04 -0.3 1.00e+00 1.00e+00h 1 23 6.8297371e-04 1.24e-07 1.94e+02 -4.0 1.07e-03 1.1 1.00e+00 1.00e+00h 1 24 7.2259234e-04 1.84e-08 1.94e+02 -4.0 4.64e-04 0.6 1.00e+00 1.00e+00h 1 25 6.8961278e-04 9.69e-08 1.94e+02 -4.1 8.97e-04 1.0 1.00e+00 1.00e+00h 1 26 6.6919971e-04 8.33e-09 1.94e+02 -4.1 2.92e-04 0.5 1.00e+00 1.00e+00h 1 27 6.6964227e-04 2.33e-10 3.91e-03 -4.0 4.30e-05 1.9 1.00e+00 1.00e+00h 1 28 6.6139848e-04 1.20e-09 2.41e-03 -4.0 9.91e-05 1.4 1.00e+00 1.00e+00h 1 29 6.4756129e-04 4.77e-09 1.62e-03 -4.0 2.00e-04 0.9 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 30 6.4074352e-04 6.87e-09 1.94e+02 -4.0 2.43e-04 0.4 1.00e+00 1.00e+00h 1 31 6.4508904e-04 2.75e-09 1.94e+02 -4.1 2.31e-04 -0.0 1.00e+00 1.00e+00h 1 32 6.4332775e-04 7.68e-10 6.43e-03 -4.1 9.00e-05 1.3 1.00e+00 1.00e+00h 1 33 6.4082539e-04 9.28e-10 6.32e-04 -4.0 9.87e-05 0.8 1.00e+00 1.00e+00h 1 34 6.3948673e-04 1.09e-09 1.94e+02 -4.0 1.12e-04 0.3 1.00e+00 1.00e+00h 1 35 6.4079950e-04 9.21e-10 1.94e+02 -4.1 2.75e-04 -0.1 1.00e+00 1.00e+00h 1 36 6.4013223e-04 2.99e-10 2.68e-03 -4.1 6.17e-05 1.2 1.00e+00 1.00e+00h 1 37 6.3972786e-04 1.56e-10 2.28e-04 -4.0 4.51e-05 0.7 1.00e+00 1.00e+00h 1 38 6.3941702e-04 2.74e-10 1.90e-04 -4.0 1.13e-04 0.2 1.00e+00 1.00e+00h 1 39 6.3937292e-04 1.19e-09 1.94e+02 -4.0 3.63e-04 -0.2 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 40 6.5459838e-04 2.17e-03 1.73e+02 -4.1 2.73e+00 - 3.35e-01 1.07e-01f 3 41 3.4234182e-01 2.33e-03 1.07e+02 -4.1 3.46e-01 -0.7 8.66e-01 5.00e-01h 2 42 9.4360098e-03 1.63e-05 6.99e-01 -4.1 1.49e-02 1.5 1.00e+00 1.00e+00h 1 43 3.1513436e-03 7.37e-07 3.39e-02 -4.1 3.18e-03 1.0 1.00e+00 1.00e+00h 1 44 1.1055245e-03 1.24e-06 1.94e+02 -4.1 4.23e-03 0.6 1.00e+00 1.00e+00h 1 45 2.1551742e-03 3.87e-07 1.94e+02 -4.1 2.59e-03 0.1 1.00e+00 1.00e+00h 1 46 1.6572421e-03 1.37e-07 6.49e-02 -4.1 1.31e-03 1.4 1.00e+00 1.00e+00h 1 47 9.8082934e-04 1.71e-07 1.35e-02 -4.1 1.60e-03 0.9 1.00e+00 1.00e+00h 1 48 7.7496919e-04 2.26e-07 1.94e+02 -4.1 1.86e-03 0.4 1.00e+00 1.00e+00h 1 49 8.3045228e-04 5.24e-08 1.94e+02 -4.1 9.15e-04 -0.0 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 50 7.4022483e-04 2.95e-08 3.42e-02 -4.1 6.16e-04 1.3 1.00e+00 1.00e+00h 1 51 6.7586917e-04 2.24e-08 3.91e-03 -4.1 5.88e-04 0.8 1.00e+00 1.00e+00h 1 52 6.7037405e-04 2.31e-08 1.93e+02 -4.1 5.96e-04 0.3 1.00e+00 1.00e+00h 1 53 6.6514682e-04 4.00e-09 1.93e+02 -4.1 2.42e-04 -0.1 1.00e+00 1.00e+00h 1 54 6.5738650e-04 3.28e-09 1.35e-02 -4.1 2.00e-04 1.2 1.00e+00 1.00e+00h 1 55 6.5572279e-04 1.33e-09 7.22e-04 -4.1 1.37e-04 0.7 1.00e+00 1.00e+00h 1 56 6.5607549e-04 9.46e-10 1.93e+02 -6.1 1.03e-04 0.2 1.00e+00 1.00e+00h 1 57 6.5577861e-04 1.07e-10 1.93e+02 -6.1 4.28e-05 -0.2 1.00e+00 1.00e+00h 1 58 6.5559404e-04 1.56e-10 1.93e+02 -6.1 3.87e-05 1.1 1.00e+00 1.00e+00h 1 59 6.5565581e-04 2.57e-11 1.93e+02 -6.1 1.50e-05 0.6 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 60 6.5558724e-04 1.10e-10 1.93e+02 -6.1 3.18e-05 1.0 1.00e+00 1.00e+00h 1 61 6.5556428e-04 1.39e-11 1.93e+02 -6.1 1.18e-05 0.6 1.00e+00 1.00e+00h 1 62 6.5555321e-04 5.04e-11 1.93e+02 -6.1 2.16e-05 1.0 1.00e+00 1.00e+00h 1 63 6.5552543e-04 6.42e-12 1.93e+02 -6.1 9.33e-06 0.5 1.00e+00 1.00e+00h 1 64 6.5552733e-04 3.41e-14 5.03e-05 -6.1 7.06e-07 1.8 1.00e+00 1.00e+00h 1 65 6.5552715e-04 1.53e-12 1.73e-04 -4.0 7.39e-06 1.4 1.00e+00 1.00e+00h 1 66 6.6806946e-04 9.12e-06 1.93e+02 -4.1 1.04e+00 - 8.23e-01 1.46e-02f 7 67 6.7689603e-04 3.96e-05 5.08e-01 -4.1 1.42e+01 - 5.15e-02 2.62e-03h 6 68 6.5666298e-04 3.12e-09 1.56e-03 -4.1 2.00e-04 0.9 1.00e+00 1.00e+00h 1 69 6.5583248e-04 2.17e-10 1.93e+02 -4.1 6.56e-05 0.4 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 70 6.7687768e-04 1.07e-02 1.10e+02 -4.1 1.07e+00 - 1.00e+00 4.33e-01f 2 71 7.2872569e-04 1.07e-02 8.38e+01 -4.1 3.84e-02 0.8 1.00e+00 2.44e-04h 13 72 7.6404587e-04 1.07e-02 8.38e+01 -4.1 3.57e-02 1.3 1.00e+00 4.88e-04h 12 73 7.7201857e-04 1.06e-02 8.38e+01 -4.1 4.10e-02 0.8 1.00e+00 2.44e-04h 13 74 7.7306970e-04 1.06e-02 8.38e+01 -4.1 3.58e-02 1.2 1.00e+00 1.53e-05h 17 75 7.0751029e-02 1.03e-04 1.96e+02 -4.1 4.43e-02 0.7 1.00e+00 1.00e+00s 22 76 8.7613825e-03 2.39e-06 1.96e+02 -4.1 6.19e-03 0.3 1.00e+00 1.00e+00s 22 77 7.5483173e-03 4.05e-07 1.38e-01 -4.1 1.85e-03 1.6 1.00e+00 1.00e+00s 22 78 3.9528240e-03 1.07e-06 3.95e-02 -4.1 3.05e-03 1.1 1.00e+00 1.00e+00s 22 79 8.0982483e-04 2.17e-06 1.96e+02 -4.1 4.35e-03 0.6 1.00e+00 0.00e+00S 22 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 80 2.8845798e-03 1.05e-06 1.96e+02 -4.1 3.11e-03 0.2 1.00e+00 1.00e+00h 1 81 2.4182427e-03 1.88e-07 9.49e-02 -4.1 1.27e-03 1.5 1.00e+00 1.00e+00h 1 82 1.3705872e-03 3.72e-07 1.86e-02 -4.1 1.82e-03 1.0 1.00e+00 1.00e+00h 1 83 7.7930815e-04 6.01e-07 1.96e+02 -4.1 2.35e-03 0.5 1.00e+00 1.00e+00h 1 84 1.0931263e-03 2.15e-07 1.96e+02 -4.1 1.44e-03 0.1 1.00e+00 1.00e+00h 1 85 9.5796762e-04 5.84e-08 5.26e-02 -4.1 7.27e-04 1.4 1.00e+00 1.00e+00h 1 86 7.5938690e-04 8.28e-08 7.16e-03 -4.1 8.85e-04 0.9 1.00e+00 1.00e+00h 1 87 6.9043327e-04 1.08e-07 1.96e+02 -4.1 1.01e-03 0.4 1.00e+00 1.00e+00h 1 88 7.1699025e-04 2.87e-08 1.96e+02 -4.1 5.44e-04 -0.0 1.00e+00 1.00e+00h 1 89 6.9217543e-04 1.17e-08 2.30e-02 -4.1 3.34e-04 1.3 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 90 6.6871990e-04 1.13e-08 2.15e-03 -4.1 3.36e-04 0.8 1.00e+00 1.00e+00h 1 91 6.6391048e-04 1.16e-08 1.96e+02 -4.1 3.44e-04 0.3 1.00e+00 1.00e+00h 1 92 6.6543763e-04 3.03e-09 1.96e+02 -4.1 1.74e-04 -0.1 1.00e+00 1.00e+00h 1 93 6.6224857e-04 1.51e-09 6.77e-03 -4.1 1.25e-04 1.2 1.00e+00 1.00e+00h 1 94 6.6057362e-04 8.68e-10 5.04e-04 -4.1 9.98e-05 0.7 1.00e+00 1.00e+00h 1 95 6.6036729e-04 7.30e-10 1.96e+02 -6.1 9.49e-05 0.2 1.00e+00 1.00e+00h 1 96 6.6025496e-04 3.00e-10 1.96e+02 -6.1 5.60e-05 -0.3 1.00e+00 1.00e+00h 1 97 6.6005404e-04 3.20e-10 1.96e+02 -6.1 5.38e-05 1.1 1.00e+00 1.00e+00h 1 98 6.6013425e-04 5.40e-11 1.96e+02 -6.1 2.13e-05 0.6 1.00e+00 1.00e+00h 1 99 6.6001270e-04 1.64e-12 2.62e-02 -6.1 1.34e-04 1.0 1.00e+00 1.00e+00H 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 100 6.5996695e-04 5.07e-11 1.96e+02 -6.1 2.11e-05 0.6 1.00e+00 1.00e+00h 1 Number of Iterations....: 100 (scaled) (unscaled) Objective...............: 6.5996694867133062e-04 6.5996694867133062e-04 Dual infeasibility......: 1.9554491748522474e+02 1.9554491748522474e+02 Constraint violation....: 5.0684123564792571e-11 5.0684123564792571e-11 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 7.1554219627044421e-07 7.1554219627044421e-07 Overall NLP error.......: 1.9554491748522474e+02 1.9554491748522474e+02 Number of objective function evaluations = 203 Number of objective gradient evaluations = 101 Number of equality constraint evaluations = 203 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 101 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 100 Total seconds in IPOPT = 41.554 EXIT: Maximum Number of Iterations Exceeded. ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 18904 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 18984 Total number of variables............................: 787 variables with only lower bounds: 50 variables with lower and upper bounds: 588 variables with only upper bounds: 0 Total number of equality constraints.................: 686 Total number of inequality constraints...............: 0 inequality constraints with only lower bounds: 0 inequality constraints with lower and upper bounds: 0 inequality constraints with only upper bounds: 0 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 0 5.5146996e+01 6.86e+00 1.66e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 4.5126323e+01 2.52e+00 3.13e+01 -0.6 4.22e+00 - 3.10e-01 6.32e-01f 1 2 1.0902812e+01 2.23e+00 5.25e+01 0.2 1.54e+01 0.0 1.42e-01 1.17e-01f 1 3 1.6651982e+01 8.82e-01 3.47e+01 -0.1 2.29e+00 -0.5 5.08e-01 6.04e-01f 1 4 1.7448499e+01 1.01e-03 3.88e+01 -0.3 5.13e-01 0.9 6.00e-01 1.00e+00f 1 5 9.1488798e+00 6.54e-03 1.57e+01 -0.2 1.56e+00 - 6.23e-01 6.07e-01f 1 6 4.2552813e+00 1.49e-03 8.33e+00 -0.7 5.14e-01 - 9.91e-01 1.00e+00f 1 7 2.3802936e+00 4.55e-04 1.97e+00 -1.2 3.35e-01 - 9.92e-01 1.00e+00f 1 8 1.2255901e+00 2.56e-04 3.13e-01 -1.8 3.07e-01 - 9.97e-01 1.00e+00f 1 9 6.4149402e-01 1.48e-04 1.91e-01 -2.7 2.93e-01 - 9.99e-01 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 4.7746310e-01 2.82e-05 4.69e-02 -3.8 1.23e-01 - 1.00e+00 1.00e+00f 1 11 4.4806848e-01 2.06e-06 3.06e-03 -4.0 3.70e-02 - 1.00e+00 1.00e+00h 1 12 4.4647094e-01 5.58e-08 1.40e-05 -4.0 1.66e-02 - 1.00e+00 1.00e+00h 1 13 4.4646514e-01 8.22e-14 1.30e-08 -4.0 6.67e-06 - 1.00e+00 1.00e+00h 1 14 4.4646514e-01 5.83e-16 1.00e-09 -4.0 4.25e-09 - 1.00e+00 1.00e+00h 1 Number of Iterations....: 14 (scaled) (unscaled) Objective...............: 4.4646514482651906e-01 4.4646514482651906e-01 Dual infeasibility......: 1.0001881193925921e-09 1.0001881193925921e-09 Constraint violation....: 5.8286708792820718e-16 5.8286708792820718e-16 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 1.0000000000001254e-04 1.0000000000001254e-04 Overall NLP error.......: 1.0001881193925921e-09 1.0001881193925921e-09 Number of objective function evaluations = 15 Number of objective gradient evaluations = 15 Number of equality constraint evaluations = 15 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 15 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 14 Total seconds in IPOPT = 21.964 EXIT: Optimal Solution Found. constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 29666 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 26679 Total number of variables............................: 933 variables with only lower bounds: 50 variables with lower and upper bounds: 834 variables with only upper bounds: 0 Total number of equality constraints.................: 833 Total number of inequality constraints...............: 0 inequality constraints with only lower bounds: 0 inequality constraints with lower and upper bounds: 0 inequality constraints with only upper bounds: 0 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 0 1.9870369e+02 6.23e-02 2.64e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 1.9248081e+02 1.36e-02 2.05e+02 -0.5 4.37e-01 2.0 8.29e-01 9.88e-01f 1 2 1.9177184e+02 9.39e-04 9.76e+01 -0.4 8.57e-02 2.4 1.00e+00 1.00e+00f 1 3 1.8897677e+02 2.55e-03 7.36e+00 -0.9 7.57e-02 1.9 9.76e-01 1.00e+00f 1 4 1.7651478e+02 2.48e-02 2.66e+01 -0.0 5.40e+00 1.5 4.46e-01 3.76e-02f 1 5 1.6591383e+02 1.07e-02 4.14e+01 -0.0 1.03e-01 1.9 7.27e-01 1.00e+00f 1 6 1.3022950e+02 3.52e-02 1.66e+01 -0.1 4.07e-01 1.4 8.77e-01 5.22e-01f 1 7 1.2064583e+02 5.23e-03 1.07e+01 -0.3 8.91e-02 1.8 7.54e-01 1.00e+00f 1 8 8.9121785e+01 1.67e-02 8.59e+00 -0.7 2.20e-01 1.4 9.72e-01 6.89e-01f 1 9 7.8744481e+01 3.41e-03 4.74e+00 -1.1 8.10e-02 1.8 9.94e-01 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 7.9243026e+01 4.75e-02 8.27e+00 0.6 4.65e+00 1.3 9.94e-02 7.68e-02f 1 11 2.5984640e+01 3.74e-02 1.39e+01 -0.5 4.52e-01 0.8 7.92e-01 7.60e-01f 1 12 6.4331497e+00 1.45e-01 3.68e+01 -0.5 4.62e-01 0.4 4.22e-01 1.00e+00f 1 13 1.0658976e+01 1.97e-02 5.70e+00 -0.6 2.64e-01 0.8 9.92e-01 1.00e+00f 1 14 1.0829314e+01 7.13e-04 2.78e+00 -1.6 3.96e-02 1.2 9.97e-01 1.00e+00h 1 15 9.5279805e+00 9.03e-04 6.36e-01 -1.5 9.03e-02 0.7 8.64e-01 1.00e+00f 1 16 6.4901301e+00 4.22e-03 2.53e+00 -2.1 1.85e-01 0.3 9.93e-01 7.67e-01f 1 17 5.5894377e+00 6.45e-04 4.11e-01 -2.3 6.70e-02 0.7 9.97e-01 1.00e+00f 1 18 4.6053092e+00 1.82e-03 1.22e+00 -2.9 1.82e-01 0.2 9.98e-01 4.03e-01f 1 19 4.1414306e+00 8.23e-04 5.03e-01 -3.6 6.55e-02 0.6 1.00e+00 5.58e-01f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 20 3.7116517e+00 1.08e-03 3.64e-01 -4.0 1.20e-01 0.2 1.00e+00 2.61e-01f 1 21 3.4546084e+00 7.84e-04 5.57e-01 -4.0 5.20e-02 0.6 1.00e+00 4.45e-01f 1 22 3.0372656e+00 1.04e-03 4.42e-01 -4.0 1.22e-01 0.1 1.00e+00 3.65e-01f 1 23 2.6828209e+00 7.17e-04 1.75e-01 -4.0 4.92e-02 0.5 1.00e+00 8.97e-01f 1 24 2.6643767e+00 6.83e-04 6.74e-01 -3.7 4.67e-01 0.1 1.00e+00 2.02e-02h 1 25 2.5533378e+00 5.25e-04 6.42e-01 -4.0 4.82e-02 0.5 1.00e+00 3.59e-01h 1 26 1.9874806e+00 3.59e-03 5.41e-01 -4.0 1.33e-01 0.0 1.00e+00 7.40e-01f 1 27 1.9310468e+00 5.95e-02 2.72e+01 -4.0 2.08e+01 -0.5 8.54e-03 1.64e-02f 1 28 1.9005720e+00 5.91e-02 2.70e+01 -2.7 5.30e-01 0.9 1.00e+00 6.44e-03h 1 29 1.5056872e+00 1.25e-03 1.06e+00 -3.0 3.21e-02 1.3 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 30 1.8535069e+00 1.00e-04 1.29e-01 -3.5 1.98e-02 0.8 1.00e+00 1.00e+00h 1 31 1.7179410e+00 2.38e-04 1.24e-01 -4.0 5.59e-02 0.3 1.00e+00 4.80e-01f 1 32 1.3738899e+00 1.66e-03 3.52e-01 -4.0 1.72e-01 -0.1 1.00e+00 5.72e-01f 1 33 1.1394349e+00 1.12e-02 4.80e+00 -2.4 9.04e-01 -0.6 4.37e-01 2.77e-01f 1 34 1.0170549e+00 1.26e-04 1.25e-01 -3.6 2.29e-02 0.7 1.00e+00 1.00e+00h 1 35 8.9706439e-01 3.74e-04 1.03e-01 -4.0 4.83e-02 0.2 1.00e+00 7.41e-01f 1 36 5.3893116e-01 3.84e-03 5.26e-01 -4.0 1.10e-01 -0.2 1.00e+00 1.00e+00f 1 37 4.8376229e-01 3.23e-04 9.27e-02 -4.0 3.78e-02 0.2 1.00e+00 1.00e+00h 1 38 3.7337670e-01 6.75e-04 1.52e-01 -4.0 8.85e-02 -0.3 1.00e+00 5.09e-01f 1 39 1.7658505e-01 2.33e-03 5.82e-01 -4.0 2.18e-01 -0.8 1.00e+00 5.54e-01f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 40 1.1700070e-01 6.23e-04 1.55e-01 -4.0 4.70e-02 -0.4 1.00e+00 8.84e-01h 1 41 1.5329518e-02 1.39e-03 1.68e-01 -4.0 6.17e-02 -0.8 1.00e+00 9.06e-01h 1 42 4.5771245e-02 1.66e-04 8.15e-01 -4.0 2.17e-02 -0.4 1.00e+00 1.00e+00h 1 43 2.6339141e-02 1.69e-05 8.35e-02 -4.0 5.83e-03 0.0 1.00e+00 1.00e+00h 1 44 1.9572110e-02 9.12e-05 3.72e-02 -4.0 1.30e-02 -0.5 1.00e+00 1.00e+00h 1 45 1.4048365e-02 2.42e-04 4.59e-01 -4.0 2.06e-02 -0.9 1.00e+00 1.00e+00h 1 46 1.5467270e-02 4.69e-05 7.17e-01 -4.0 8.73e-03 -0.5 1.00e+00 1.00e+00h 1 47 1.4794845e-02 3.66e-06 3.69e-02 -4.0 2.45e-03 -0.1 1.00e+00 1.00e+00h 1 48 1.4271451e-02 2.60e-05 1.38e-01 -4.0 6.16e-03 -0.6 1.00e+00 1.00e+00h 1 49 1.3606551e-02 9.16e-05 7.15e-01 -4.0 1.40e-02 -1.0 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 50 1.3687554e-02 2.93e-07 3.51e-03 -4.0 6.36e-04 0.3 1.00e+00 1.00e+00h 1 51 1.3611985e-02 2.47e-06 1.35e-02 -4.0 1.96e-03 -0.2 1.00e+00 1.00e+00h 1 52 1.3442763e-02 1.82e-05 6.18e-02 -4.0 5.28e-03 -0.7 1.00e+00 1.00e+00h 1 53 1.3027576e-02 8.98e-05 3.60e-01 -4.0 1.48e-02 -1.1 1.00e+00 1.00e+00h 1 54 1.2924067e-02 5.70e-05 2.26e-02 -4.0 9.63e-03 -0.7 1.00e+00 1.00e+00h 1 55 1.2709472e-02 1.13e-04 2.35e-01 -4.0 1.44e-02 -1.2 1.00e+00 1.00e+00h 1 56 1.2434071e-02 2.43e-05 3.04e-01 -4.0 6.35e-03 -0.8 1.00e+00 1.00e+00h 1 57 1.2374840e-02 1.94e-06 3.17e-02 -4.0 2.09e-03 -0.3 1.00e+00 1.00e+00h 1 58 1.2294398e-02 1.59e-05 9.20e-02 -4.0 5.60e-03 -0.8 1.00e+00 1.00e+00h 1 59 1.1872109e-02 6.83e-04 3.19e-01 -4.0 3.67e-02 -1.3 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 60 1.1808445e-02 3.90e-04 6.49e-02 -4.0 2.87e-02 -0.9 1.00e+00 1.00e+00h 1 61 1.1617932e-02 1.11e-04 1.41e-01 -4.0 1.71e-02 -1.3 1.00e+00 1.00e+00h 1 62 1.1459548e-02 2.06e-05 8.46e-02 -4.0 6.08e-03 -0.9 1.00e+00 1.00e+00h 1 63 1.0994783e-02 2.00e-03 1.10e+00 -4.0 6.13e-02 -1.4 1.00e+00 1.00e+00h 1 64 1.0998497e-02 1.20e-06 7.54e-03 -4.0 1.48e-03 -0.1 1.00e+00 1.00e+00h 1 65 1.1035521e-02 2.71e-06 7.41e-02 -4.0 2.51e-03 -0.5 1.00e+00 1.00e+00h 1 66 1.0800728e-02 2.56e-05 5.75e-02 -4.1 5.99e-03 -1.0 1.00e+00 1.00e+00h 1 67 1.0477321e-02 2.90e-04 7.43e-01 -4.0 2.61e-02 -1.5 1.00e+00 1.00e+00h 1 68 1.0515141e-02 2.86e-06 4.00e-02 -4.0 2.23e-03 -0.2 1.00e+00 1.00e+00h 1 69 1.0501138e-02 3.42e-06 1.90e-02 -4.0 2.38e-03 -0.6 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 70 1.0337437e-02 3.10e-05 8.99e-02 -4.0 7.12e-03 -1.1 1.00e+00 1.00e+00h 1 71 1.0275067e-02 3.70e-06 1.61e-02 -4.0 2.91e-03 -0.7 1.00e+00 1.00e+00h 1 72 1.0136405e-02 3.40e-05 2.23e-02 -4.0 8.08e-03 -1.2 1.00e+00 1.00e+00h 1 73 9.6950463e-03 3.07e-04 2.28e-01 -4.0 2.42e-02 -1.6 1.00e+00 1.00e+00h 1 74 9.6662516e-03 5.92e-07 6.19e-03 -4.0 1.33e-03 -0.3 1.00e+00 1.00e+00h 1 75 9.6469701e-03 4.93e-06 4.71e-02 -4.0 3.83e-03 -0.8 1.00e+00 1.00e+00h 1 76 9.4749637e-03 5.17e-05 5.61e-02 -4.0 1.07e-02 -1.3 1.00e+00 1.00e+00h 1 77 9.3960251e-03 7.45e-06 4.45e-02 -4.0 4.42e-03 -0.8 1.00e+00 1.00e+00h 1 78 9.2451747e-03 7.95e-05 9.68e-02 -4.0 1.56e-02 -1.3 1.00e+00 1.00e+00h 1 79 9.1849284e-03 7.61e-06 2.19e-03 -4.0 4.61e-03 -0.9 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 80 8.9936415e-03 8.01e-05 7.09e-02 -4.0 1.55e-02 -1.4 1.00e+00 1.00e+00h 1 81 8.9272783e-03 1.25e-05 8.76e-04 -4.0 6.17e-03 -0.9 1.00e+00 1.00e+00h 1 82 8.7454641e-03 1.03e-04 3.73e-02 -4.0 1.80e-02 -1.4 1.00e+00 1.00e+00h 1 83 8.2001682e-03 1.13e-03 1.80e-01 -4.0 4.44e-02 -1.9 1.00e+00 1.00e+00h 1 84 8.1683081e-03 2.31e-06 6.24e-03 -4.0 2.40e-03 -0.6 1.00e+00 1.00e+00h 1 85 8.1006864e-03 1.71e-05 1.61e-03 -4.0 7.21e-03 -1.1 1.00e+00 1.00e+00h 1 86 7.9106718e-03 1.51e-04 1.65e-02 -4.0 2.13e-02 -1.5 1.00e+00 1.00e+00h 1 87 7.4439437e-03 5.14e-04 5.22e-02 -4.0 4.32e-02 -2.0 1.00e+00 1.00e+00h 1 88 7.2830388e-03 6.26e-05 1.26e-02 -4.0 1.90e-02 -1.6 1.00e+00 1.00e+00h 1 89 6.8179277e-03 1.00e-03 2.94e-01 -4.0 5.08e-02 -2.1 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 90 6.7859774e-03 2.99e-06 8.69e-04 -4.0 3.99e-03 -0.7 1.00e+00 1.00e+00h 1 91 6.7093442e-03 1.44e-05 6.57e-02 -4.0 7.50e-03 -1.2 1.00e+00 1.00e+00h 1 92 6.5387772e-03 2.74e-04 2.45e-02 -4.0 2.68e-02 -1.7 1.00e+00 1.00e+00h 1 93 6.2178574e-03 1.30e-03 3.44e-01 -4.0 6.59e-02 -2.2 1.00e+00 1.00e+00h 1 94 6.1945015e-03 8.80e-06 3.85e-03 -4.0 4.69e-03 -0.8 1.00e+00 1.00e+00h 1 95 6.1355563e-03 1.61e-05 2.69e-02 -4.0 9.19e-03 -1.3 1.00e+00 1.00e+00h 1 96 5.9818834e-03 2.77e-04 2.56e-02 -4.0 2.81e-02 -1.8 1.00e+00 1.00e+00h 1 97 5.6450658e-03 1.09e-02 2.59e-01 -4.0 1.40e-01 -2.3 1.00e+00 1.00e+00h 1 98 5.6416046e-03 6.43e-05 2.74e-02 -4.0 1.29e-02 -0.0 1.00e+00 1.00e+00h 1 99 5.6248013e-03 5.13e-07 7.51e-03 -4.0 1.35e-03 -0.5 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 100 5.6156659e-03 3.49e-06 1.94e-02 -4.0 3.72e-03 -1.0 1.00e+00 1.00e+00h 1 101 5.5616981e-03 3.63e-05 3.97e-03 -4.0 1.14e-02 -1.5 1.00e+00 1.00e+00h 1 102 5.3969903e-03 2.13e-04 2.50e-01 -4.0 2.91e-02 -1.9 1.00e+00 1.00e+00h 1 103 5.4048973e-03 1.52e-06 1.83e-02 -4.0 2.63e-03 -0.6 1.00e+00 1.00e+00h 1 104 5.4036818e-03 4.29e-06 2.39e-02 -4.0 4.10e-03 -1.1 1.00e+00 1.00e+00h 1 105 5.3387104e-03 5.20e-05 7.85e-02 -4.0 1.42e-02 -1.6 1.00e+00 1.00e+00h 1 106 5.3116005e-03 5.01e-06 1.49e-02 -4.0 5.22e-03 -1.1 1.00e+00 1.00e+00h 1 107 5.2618862e-03 4.61e-05 4.47e-02 -4.0 1.53e-02 -1.6 1.00e+00 1.00e+00h 1 108 5.2376462e-03 7.34e-06 9.06e-03 -4.0 5.94e-03 -1.2 1.00e+00 1.00e+00h 1 109 5.1771411e-03 5.25e-05 3.81e-03 -4.0 1.74e-02 -1.7 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 110 5.0192558e-03 3.47e-04 4.10e-02 -4.0 4.90e-02 -2.1 1.00e+00 1.00e+00h 1 111 4.9588987e-03 8.28e-05 2.39e-03 -4.0 2.43e-02 -1.7 1.00e+00 1.00e+00h 1 112 4.8711307e-03 3.49e-04 4.86e-01 -4.0 4.07e-02 -2.2 1.00e+00 1.00e+00h 1 113 4.8578362e-03 2.99e-06 5.76e-03 -4.0 2.96e-03 -0.9 1.00e+00 1.00e+00h 1 114 4.7630540e-03 8.00e-06 1.57e-01 -4.0 8.75e-03 -1.3 1.00e+00 1.00e+00h 1 115 4.8102608e-03 8.02e-08 1.26e-02 -4.0 5.37e-04 -0.0 1.00e+00 1.00e+00h 1 116 4.8061927e-03 1.54e-07 3.01e-03 -4.0 1.06e-03 -0.5 1.00e+00 1.00e+00h 1 117 4.7993788e-03 1.32e-06 3.26e-04 -4.0 3.17e-03 -1.0 1.00e+00 1.00e+00h 1 118 4.7763464e-03 1.17e-05 4.56e-03 -4.0 9.50e-03 -1.4 1.00e+00 1.00e+00h 1 119 4.6105601e-03 7.03e-05 1.04e-02 -4.1 2.85e-02 -1.9 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 120 4.5646552e-03 1.33e-03 8.54e-02 -4.0 8.26e-02 -2.4 1.00e+00 1.00e+00h 1 121 4.5543593e-03 2.92e-06 2.02e-03 -4.0 3.36e-03 -1.1 1.00e+00 1.00e+00h 1 122 4.5140119e-03 1.30e-05 5.53e-02 -4.0 1.10e-02 -1.5 1.00e+00 1.00e+00h 1 123 4.5157385e-03 2.18e-06 1.39e-02 -4.0 3.98e-03 -1.1 1.00e+00 1.00e+00h 1 124 4.4936431e-03 1.25e-05 4.37e-03 -4.0 1.15e-02 -1.6 1.00e+00 1.00e+00h 1 125 4.4262577e-03 1.23e-04 2.72e-02 -4.0 3.30e-02 -2.1 1.00e+00 1.00e+00h 1 126 4.4062316e-03 2.37e-05 1.22e-02 -4.0 1.25e-02 -1.6 1.00e+00 1.00e+00h 1 127 4.3557474e-03 2.00e-04 5.00e-02 -4.0 3.51e-02 -2.1 1.00e+00 1.00e+00h 1 128 4.3295634e-03 1.65e-05 1.81e-02 -4.0 9.82e-03 -1.7 1.00e+00 1.00e+00h 1 129 4.3003279e-03 3.39e-06 2.75e-02 -4.0 5.00e-03 -1.3 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 130 4.2883146e-03 4.32e-05 6.84e-02 -4.0 1.66e-02 -1.8 1.00e+00 1.00e+00h 1 131 4.2889477e-03 7.46e-06 1.32e-02 -4.0 5.38e-03 -1.3 1.00e+00 1.00e+00h 1 132 4.2487196e-03 3.37e-05 6.34e-02 -4.0 1.52e-02 -1.8 1.00e+00 1.00e+00h 1 133 4.2422585e-03 9.24e-06 1.03e-02 -4.0 7.34e-03 -1.4 1.00e+00 1.00e+00h 1 134 4.2236244e-03 4.08e-05 2.22e-02 -4.0 1.77e-02 -1.9 1.00e+00 1.00e+00h 1 135 4.2035477e-03 4.77e-06 1.72e-02 -4.0 5.85e-03 -1.4 1.00e+00 1.00e+00h 1 136 4.1741054e-03 1.12e-04 3.52e-02 -4.0 2.73e-02 -1.9 1.00e+00 1.00e+00h 1 137 4.1748135e-03 8.39e-06 1.65e-02 -4.0 8.24e-03 -1.5 1.00e+00 1.00e+00h 1 138 4.1264254e-03 8.56e-05 1.09e-01 -4.0 1.96e-02 -2.0 1.00e+00 1.00e+00h 1 139 4.1338726e-03 4.57e-07 6.18e-03 -4.0 1.46e-03 -0.6 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 140 4.1328271e-03 1.10e-06 1.31e-03 -4.0 3.19e-03 -1.1 1.00e+00 1.00e+00h 1 141 4.1172152e-03 1.01e-05 9.68e-03 -4.0 9.49e-03 -1.6 1.00e+00 1.00e+00h 1 142 4.0767978e-03 2.59e-04 2.07e-02 -4.0 4.03e-02 -2.1 1.00e+00 1.00e+00h 1 143 4.0738885e-03 1.31e-05 1.13e-02 -4.0 1.15e-02 -1.6 1.00e+00 1.00e+00h 1 144 4.0147501e-03 1.44e-04 1.46e-01 -4.0 2.76e-02 -2.1 1.00e+00 1.00e+00h 1 145 4.0195730e-03 9.02e-07 6.71e-03 -4.0 2.23e-03 -0.8 1.00e+00 1.00e+00h 1 146 4.0230333e-03 2.69e-06 8.77e-03 -4.0 4.58e-03 -1.3 1.00e+00 1.00e+00h 1 147 3.9142958e-03 2.19e-05 4.06e-03 -4.1 1.27e-02 -1.7 1.00e+00 1.00e+00h 1 148 3.9600838e-03 4.19e-04 6.28e-02 -4.0 4.85e-02 -2.2 1.00e+00 1.00e+00h 1 149 3.8716535e-03 5.85e-07 5.01e-04 -4.1 2.19e-03 -0.9 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 150 3.9260332e-03 5.46e-06 5.50e-02 -4.0 5.48e-03 -1.4 1.00e+00 1.00e+00h 1 Number of Iterations....: 150 (scaled) (unscaled) Objective...............: 3.9260331977887470e-03 3.9260331977887470e-03 Dual infeasibility......: 5.5015966222676085e-02 5.5015966222676085e-02 Constraint violation....: 5.4557546105837851e-06 5.4557546105837851e-06 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 1.0092118189813168e-04 1.0092118189813168e-04 Overall NLP error.......: 5.5015966222676085e-02 5.5015966222676085e-02 Number of objective function evaluations = 151 Number of objective gradient evaluations = 151 Number of equality constraint evaluations = 151 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 151 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 150 Total seconds in IPOPT = 65.406 EXIT: Maximum Number of Iterations Exceeded. constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing SmoothPulseProblem for MultiKetTrajectory (2 states)... ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##175".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##175".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##175".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##175".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 applying timesteps_all_equal constraint: Δt ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.ZeroOrderPulse{DataInterpolations.ConstantInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##175".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##175".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##175".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##175".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SamplingProblem... using 2 systems constructing SmoothPulseProblem for UnitaryTrajectory... constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt added DerivativeIntegrator for LinearSplinePulse constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.CubicSplinePulse{DataInterpolations.CubicHermiteSpline{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Matrix{Float64}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.CubicSplinePulse{DataInterpolations.CubicHermiteSpline{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Matrix{Float64}, Vector{Vector{Float64}}, Float64}}... set du bounds to ±5.0 for CubicSplinePulse applying timesteps_all_equal constraint: Δt constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.CubicSplinePulse{DataInterpolations.CubicHermiteSpline{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Matrix{Float64}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt added DerivativeIntegrator for LinearSplinePulse constructing SplinePulseProblem for MultiKetTrajectory with PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}... with 2 state transfers ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#26#27"{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}, Vector{SparseArrays.SparseMatrixCSC{ComplexF64, Int64}}, Int64}, PiccoloQuantumObjects.QuantumSystems.var"#28#29"{Vector{SparseArrays.SparseMatrixCSC{Float64, Int64}}, Int64, SparseArrays.SparseMatrixCSC{Float64, Int64}}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 applying timesteps_all_equal constraint: Δt added DerivativeIntegrator for LinearSplinePulse constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt added DerivativeIntegrator for LinearSplinePulse constructing SamplingProblem... using 2 systems constructing SplinePulseProblem with PiccoloQuantumObjects.Pulses.CubicSplinePulse{DataInterpolations.CubicHermiteSpline{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Matrix{Float64}, Vector{Vector{Float64}}, Float64}}... applying timesteps_all_equal constraint: Δt initializing optimizer... building evaluator: 1 integrators, 0 nonlinear constraints dynamics constraints: 40, nonlinear constraints: 0 integrator 1 jacobian structure: 0.082s jacobian structure: 560 nonzeros (0.082s) integrator 1 hessian structure: 0.08s computing objective hessian structure (QuadraticRegularizer)... objective hessian structure: 0.105s hessian structure: 798 nonzeros (0.186s) linear index maps built (0.0s) evaluator ready (total: 0.269s) evaluator created (0.269s) NL constraint bounds extracted (0.0s) NLP block data built (0.0s) Ipopt optimizer configured (0.0s) variables set (0.185s) applying constraint: initial value of ψ̃ applying constraint: initial value of u applying constraint: final value of u applying constraint: bounds on ψ̃ applying constraint: bounds on u applying constraint: bounds on Δt applying constraint: time consistency constraint (t_{k+1} = t_k + Δt_k) applying constraint: initial time t₁ = 0 linear constraints added: 8 (1.99s) optimizer initialization complete (total: 2.445s) This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 561 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 715 Number of Iterations....: 0 Number of objective function evaluations = 0 Number of objective gradient evaluations = 1 Number of equality constraint evaluations = 0 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 1 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 0 Total seconds in IPOPT = 9.683 EXIT: Invalid number in NLP function or derivative detected. initializing optimizer... building evaluator: 1 integrators, 0 nonlinear constraints dynamics constraints: 40, nonlinear constraints: 0 integrator 1 jacobian structure: 0.0s jacobian structure: 560 nonzeros (0.0s) integrator 1 hessian structure: 0.0s computing objective hessian structure (QuadraticRegularizer)... objective hessian structure: 0.0s hessian structure: 798 nonzeros (0.001s) linear index maps built (0.0s) evaluator ready (total: 0.002s) evaluator created (0.002s) NL constraint bounds extracted (0.0s) NLP block data built (0.0s) Ipopt optimizer configured (0.0s) variables set (0.0s) applying constraint: initial value of ψ̃ applying constraint: initial value of u applying constraint: final value of u applying constraint: bounds on ψ̃ applying constraint: bounds on u applying constraint: bounds on Δt applying constraint: time consistency constraint (t_{k+1} = t_k + Δt_k) applying constraint: initial time t₁ = 0 linear constraints added: 8 (0.0s) optimizer initialization complete (total: 0.003s) This is Ipopt version 3.14.19, running with linear solver MUMPS 5.9.0. Number of nonzeros in equality constraint Jacobian...: 561 Number of nonzeros in inequality constraint Jacobian.: 0 Number of nonzeros in Lagrangian Hessian.............: 715 Total number of variables............................: 70 variables with only lower bounds: 11 variables with lower and upper bounds: 49 variables with only upper bounds: 0 Total number of equality constraints.................: 50 Total number of inequality constraints...............: 0 inequality constraints with only lower bounds: 0 inequality constraints with lower and upper bounds: 0 inequality constraints with only upper bounds: 0 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 0 0.0000000e+00 1.00e-02 1.00e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 Number of Iterations....: 0 (scaled) (unscaled) Objective...............: 0.0000000000000000e+00 0.0000000000000000e+00 Dual infeasibility......: 1.0000000000000000e+00 1.0000000000000000e+00 Constraint violation....: 9.9999901000000779e-03 9.9999901000000779e-03 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 2.0000000199999999e+00 2.0000000199999999e+00 Overall NLP error.......: 1.9999000199999999e+00 1.9999000199999999e+00 Number of objective function evaluations = 1 Number of objective gradient evaluations = 1 Number of equality constraint evaluations = 1 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 1 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 0 Total seconds in IPOPT = 0.215 EXIT: Maximum Number of Iterations Exceeded. [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=88, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=220, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=220, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=77, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##198".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##198".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##198".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##198".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 ┌ Warning: `Base.getindex(VA::AbstractVectorOfArray{T, N, A}, I::Int) where {T, N, A <: Union{AbstractArray, AbstractVectorOfArray}}` is deprecated, use `VA.u[I]` instead. │ caller = getindex(traj::PiccoloQuantumObjects.QuantumTrajectories.MultiKetTrajectory{PiccoloQuantumObjects.Pulses.LinearSplinePulse{DataInterpolations.LinearInterpolation{Matrix{Float64}, Vector{Float64}, Vector{Union{}}, Vector{Vector{Float64}}, Float64}}, SciMLBase.EnsembleSolution{ComplexF64, 3, Vector{SciMLBase.ODESolution{ComplexF64, 2, Vector{Vector{ComplexF64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, SciMLBase.ODEProblem{Vector{ComplexF64}, Tuple{Float64, Float64}, true, SciMLBase.NullParameters, SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##198".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##198".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Base.Pairs{Symbol, Vector{Float64}, Nothing, @NamedTuple{tstops::Vector{Float64}, saveat::Vector{Float64}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEqLinear.MagnusGL4, OrdinaryDiffEqCore.InterpolationData{SciMLBase.ODEFunction{true, SciMLBase.FullSpecialize, SciMLOperators.MatrixOperator{ComplexF64, Matrix{ComplexF64}, SciMLOperators.FilterKwargs{Nothing, Val{()}}, SciMLOperators.FilterKwargs{PiccoloQuantumObjects.Rollouts.var"#update!#_construct_operator##2"{PiccoloQuantumObjects.QuantumSystems.QuantumSystem{PiccoloQuantumObjects.QuantumSystems.var"#10#11"{typeof(Main.var"##198".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##198".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing, PiccoloQuantumObjects.Rollouts.PiccoloRolloutSystem{Union{Int64, AbstractVector{Int64}, CartesianIndex, CartesianIndices}}, Nothing, Nothing}, Vector{Vector{ComplexF64}}, Vector{Float64}, Vector{Vector{Vector{ComplexF64}}}, Nothing, OrdinaryDiffEqLinear.MagnusGL4Cache{Vector{ComplexF64}, Vector{ComplexF64}, Matrix{ComplexF64}, Nothing}, Nothing}, SciMLBase.DEStats, Nothing, Nothing, Nothing, Nothing}}}}, i::Int64) at ensemble_trajectory.jl:109 [inlined] └ @ PiccoloQuantumObjects.QuantumTrajectories ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=40, datavec_len=132, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=209, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! [test_integrator] Starting... [test_integrator] constraint_dim=80, datavec_len=209, global_len=0 [test_integrator] Testing evaluate!... [test_integrator] Computing analytic Jacobian... [test_integrator] Computing finite diff Jacobian... [test_integrator] Testing Jacobian... [test_integrator] Computing analytic Hessian... [test_integrator] Computing finite diff Hessian... [test_integrator] Testing Hessian... [test_integrator] Done! Test Summary: | Pass Broken Total Time QuantumCollocation | 221 4 225 27m26.3s Testing QuantumCollocation tests passed Testing completed after 1688.39s PkgEval succeeded after 1863.3s