Package evaluation to test QuantumCollocation on Julia 1.14.0-DEV.1621 (4d04bb6b3b*) started at 2026-01-28T07:37:21.984 ################################################################################ # Set-up # Installing PkgEval dependencies (TestEnv)... Activating project at `~/.julia/environments/v1.14` Set-up completed after 10.66s ################################################################################ # 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.21.0 ⌃ [c29ec348] + AbstractDifferentiation v0.4.4 [7d9f7c33] + Accessors v0.1.43 [79e6a3ab] + Adapt v4.4.0 [66dad0bd] + AliasTables v1.1.3 [4fba245c] + ArrayInterface v7.22.0 [4c555306] + ArrayLayouts v1.12.2 [6e4b80f9] + BenchmarkTools v1.6.3 [62783981] + BitTwiddlingConvenienceFunctions v0.1.6 [70df07ce] + BracketingNonlinearSolve v1.6.2 [2a0fbf3d] + CPUSummary v0.2.7 [d360d2e6] + ChainRulesCore v1.26.0 [0b6fb165] + ChunkCodecCore v1.0.1 [4c0bbee4] + ChunkCodecLibZlib v1.0.0 [55437552] + ChunkCodecLibZstd v1.0.0 [fb6a15b2] + CloseOpenIntervals v0.1.13 [523fee87] + CodecBzip2 v0.8.5 [944b1d66] + CodecZlib v0.7.8 [38540f10] + CommonSolve v0.2.6 [bbf7d656] + CommonSubexpressions v0.3.1 [f70d9fcc] + CommonWorldInvalidations v1.0.0 [34da2185] + Compat v4.18.1 [a33af91c] + CompositionsBase v0.1.2 [95dc2771] + ComputePipeline v0.1.6 [2569d6c7] + ConcreteStructs v0.2.3 [187b0558] + ConstructionBase v1.6.0 [adafc99b] + CpuId v0.3.1 [a8cc5b0e] + Crayons v4.1.1 [9a962f9c] + DataAPI v1.16.0 [82cc6244] + DataInterpolations v8.9.0 [864edb3b] + DataStructures v0.19.3 [e2d170a0] + DataValueInterfaces v1.0.0 [2b5f629d] + DiffEqBase v6.199.0 [163ba53b] + DiffResults v1.1.0 [b552c78f] + DiffRules v1.15.1 [a0c0ee7d] + DifferentiationInterface v0.7.15 [c823fa1f] + DirectTrajOpt v0.8.1 [31c24e10] + Distributions v0.25.123 [ffbed154] + DocStringExtensions v0.9.5 [4e289a0a] + EnumX v1.0.6 [f151be2c] + EnzymeCore v0.8.18 [e24c0720] + ExponentialAction v0.2.10 [d4d017d3] + ExponentialUtilities v1.30.0 [e2ba6199] + ExprTools v0.1.10 [55351af7] + ExproniconLite v0.10.14 [7034ab61] + FastBroadcast v0.3.5 [9aa1b823] + FastClosures v0.3.2 [a4df4552] + FastPower v1.3.0 [5789e2e9] + FileIO v1.17.1 [1a297f60] + FillArrays v1.16.0 [64ca27bc] + FindFirstFunctions v1.8.0 [6a86dc24] + FiniteDiff v2.29.0 [f6369f11] + ForwardDiff v1.3.1 [069b7b12] + FunctionWrappers v1.1.3 [77dc65aa] + FunctionWrappersWrappers v0.1.3 [46192b85] + GPUArraysCore v0.2.0 [c145ed77] + GenericSchur v0.5.6 [076d061b] + HashArrayMappedTries v0.2.0 [34004b35] + HypergeometricFunctions v0.3.28 [615f187c] + IfElse v0.1.1 [3587e190] + InverseFunctions v0.1.17 [b6b21f68] + Ipopt v1.14.0 [92d709cd] + IrrationalConstants v0.2.6 [82899510] + IteratorInterfaceExtensions v1.0.0 [033835bb] + JLD2 v0.6.3 [692b3bcd] + JLLWrappers v1.7.1 [682c06a0] + JSON v1.4.0 [0f8b85d8] + JSON3 v1.14.3 [ae98c720] + Jieko v0.2.1 [b964fa9f] + LaTeXStrings v1.4.0 [10f19ff3] + LayoutPointers v0.1.17 [5078a376] + LazyArrays v2.9.5 [2ab3a3ac] + LogExpFunctions v0.3.29 [e6f89c97] + LoggingExtras v1.2.0 [1914dd2f] + MacroTools v0.5.16 [d125e4d3] + ManualMemory v0.1.8 [b8f27783] + MathOptInterface v1.48.0 [bb5d69b7] + MaybeInplace v0.1.4 [e1d29d7a] + Missings v1.2.0 [2e0e35c7] + Moshi v0.3.7 [46d2c3a1] + MuladdMacro v0.2.4 [d8a4904e] + MutableArithmetics v1.6.7 [77ba4419] + NaNMath v1.1.3 [538bc3a1] + NamedTrajectories v0.8.2 [be0214bd] + NonlinearSolveBase v2.11.1 [510215fc] + Observables v0.5.5 [bac558e1] + OrderedCollections v1.8.1 [bbf590c4] + OrdinaryDiffEqCore v3.2.0 [521117fe] + OrdinaryDiffEqLinear v1.10.0 [b1df2697] + OrdinaryDiffEqTsit5 v1.9.0 [90014a1f] + PDMats v0.11.37 [69de0a69] + Parsers v2.8.3 [5a402ddf] + PiccoloQuantumObjects v0.10.2 [f517fe37] + Polyester v0.7.18 [1d0040c9] + PolyesterWeave v0.2.2 [d236fae5] + PreallocationTools v1.1.0 [aea7be01] + PrecompileTools v1.3.3 [21216c6a] + Preferences v1.5.1 [08abe8d2] + PrettyTables v3.1.2 [92933f4c] + ProgressMeter v1.11.0 [43287f4e] + PtrArrays v1.3.0 [1fd47b50] + QuadGK v2.11.2 [0dc23a59] + QuantumCollocation v0.10.2 [3cdcf5f2] + RecipesBase v1.3.4 [731186ca] + RecursiveArrayTools v3.46.0 [189a3867] + Reexport v1.2.2 [ae029012] + Requires v1.3.1 [79098fc4] + Rmath v0.9.0 [7e49a35a] + RuntimeGeneratedFunctions v0.5.16 [94e857df] + SIMDTypes v0.1.0 [0bca4576] + SciMLBase v2.135.0 [19f34311] + SciMLJacobianOperators v0.1.12 [a6db7da4] + SciMLLogging v1.8.0 [c0aeaf25] + SciMLOperators v1.14.1 [431bcebd] + SciMLPublic v1.0.1 [53ae85a6] + SciMLStructures v1.10.0 [7e506255] + ScopedValues v1.5.0 [efcf1570] + Setfield v1.1.2 [a2af1166] + SortingAlgorithms v1.2.2 [276daf66] + SpecialFunctions v2.6.1 [aedffcd0] + Static v1.3.1 [0d7ed370] + StaticArrayInterface v1.8.0 [90137ffa] + StaticArrays v1.9.16 [1e83bf80] + StaticArraysCore v1.4.4 [10745b16] + Statistics v1.11.1 [82ae8749] + StatsAPI v1.8.0 [2913bbd2] + StatsBase v0.34.10 [4c63d2b9] + StatsFuns v1.5.2 [7792a7ef] + StrideArraysCore v0.5.8 [892a3eda] + StringManipulation v0.4.2 [856f2bd8] + StructTypes v1.11.0 [ec057cc2] + StructUtils v2.6.2 [2efcf032] + SymbolicIndexingInterface v0.3.46 [3783bdb8] + TableTraits v1.0.1 [bd369af6] + Tables v1.12.1 [f8b46487] + TestItemRunner v1.1.4 [1c621080] + TestItems v1.0.0 [8290d209] + ThreadingUtilities v0.5.5 [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.3+0 [6e34b625] + Bzip2_jll v1.0.9+0 [e33a78d0] + Hwloc_jll v2.12.2+0 [9cc047cb] + Ipopt_jll v300.1400.1901+0 [94ce4f54] + Libiconv_jll v1.18.0+0 [d00139f3] + METIS_jll v5.1.3+0 [d7ed1dd3] + MUMPS_seq_jll v500.800.200+0 [656ef2d0] + OpenBLAS32_jll v0.3.30+0 [efe28fd5] + OpenSpecFun_jll v0.5.6+0 [f50d1b31] + Rmath_jll v0.5.1+0 [319450e9] + SPRAL_jll v2025.9.18+0 ⌅ [02c8fc9c] + XML2_jll v2.13.9+0 [a65dc6b1] + Xorg_libpciaccess_jll v0.18.1+0 [0dad84c5] + ArgTools v1.1.2 [56f22d72] + Artifacts v1.11.0 [2a0f44e3] + Base64 v1.11.0 [ade2ca70] + Dates v1.11.0 [8ba89e20] + Distributed v1.11.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.13.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 [9abbd945] + Profile v1.11.0 [3fa0cd96] + REPL v1.11.0 [9a3f8284] + Random v1.11.0 [ea8e919c] + SHA v1.0.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.3.0+1 [deac9b47] + LibCURL_jll v8.18.0+0 [e37daf67] + LibGit2_jll v1.9.2+0 [29816b5a] + LibSSH2_jll v1.11.3+1 [14a3606d] + MozillaCACerts_jll v2025.12.2 [4536629a] + OpenBLAS_jll v0.3.30+0 [05823500] + OpenLibm_jll v0.8.7+0 [458c3c95] + OpenSSL_jll v3.5.4+0 [efcefdf7] + PCRE2_jll v10.47.0+0 [bea87d4a] + SuiteSparse_jll v7.10.1+0 [83775a58] + Zlib_jll v1.3.1+2 [3161d3a3] + Zstd_jll v1.5.7+1 [8e850b90] + libblastrampoline_jll v5.15.0+0 [8e850ede] + nghttp2_jll v1.68.0+1 [3f19e933] + p7zip_jll v17.7.0+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 5.98s ################################################################################ # Precompilation # Precompiling PkgEval dependencies... Precompiling package dependencies... Precompiling packages... 26703.0 ms ✓ DirectTrajOpt 36154.7 ms ✓ QuantumCollocation 2 dependencies successfully precompiled in 74 seconds. 270 already precompiled. Precompilation completed after 87.83s ################################################################################ # Testing # Testing QuantumCollocation Status `/tmp/jl_w5vnFF/Project.toml` [c823fa1f] DirectTrajOpt v0.8.1 [31c24e10] Distributions v0.25.123 [e24c0720] ExponentialAction v0.2.10 [538bc3a1] NamedTrajectories v0.8.2 [5a402ddf] PiccoloQuantumObjects v0.10.2 [0dc23a59] QuantumCollocation v0.10.2 [189a3867] Reexport v1.2.2 [f8b46487] TestItemRunner v1.1.4 [1c621080] TestItems v1.0.0 [6dad8b7f] TrajectoryIndexingUtils v0.1.0 [8f399da3] Libdl v1.11.0 [37e2e46d] LinearAlgebra v1.13.0 [9a3f8284] Random v1.11.0 [2f01184e] SparseArrays v1.13.0 [8dfed614] Test v1.11.0 Status `/tmp/jl_w5vnFF/Manifest.toml` [47edcb42] ADTypes v1.21.0 ⌃ [c29ec348] AbstractDifferentiation v0.4.4 [7d9f7c33] Accessors v0.1.43 [79e6a3ab] Adapt v4.4.0 [66dad0bd] AliasTables v1.1.3 [4fba245c] ArrayInterface v7.22.0 [4c555306] ArrayLayouts v1.12.2 [6e4b80f9] BenchmarkTools v1.6.3 [62783981] BitTwiddlingConvenienceFunctions v0.1.6 [70df07ce] BracketingNonlinearSolve v1.6.2 [2a0fbf3d] CPUSummary v0.2.7 [d360d2e6] ChainRulesCore v1.26.0 [0b6fb165] ChunkCodecCore v1.0.1 [4c0bbee4] ChunkCodecLibZlib v1.0.0 [55437552] ChunkCodecLibZstd v1.0.0 [fb6a15b2] CloseOpenIntervals v0.1.13 [523fee87] CodecBzip2 v0.8.5 [944b1d66] CodecZlib v0.7.8 [38540f10] CommonSolve v0.2.6 [bbf7d656] CommonSubexpressions v0.3.1 [f70d9fcc] CommonWorldInvalidations v1.0.0 [34da2185] Compat v4.18.1 [a33af91c] CompositionsBase v0.1.2 [95dc2771] ComputePipeline v0.1.6 [2569d6c7] ConcreteStructs v0.2.3 [187b0558] ConstructionBase v1.6.0 [adafc99b] CpuId v0.3.1 [a8cc5b0e] Crayons v4.1.1 [9a962f9c] DataAPI v1.16.0 [82cc6244] DataInterpolations v8.9.0 [864edb3b] DataStructures v0.19.3 [e2d170a0] DataValueInterfaces v1.0.0 [2b5f629d] DiffEqBase v6.199.0 [163ba53b] DiffResults v1.1.0 [b552c78f] DiffRules v1.15.1 [a0c0ee7d] DifferentiationInterface v0.7.15 [c823fa1f] DirectTrajOpt v0.8.1 [31c24e10] Distributions v0.25.123 [ffbed154] DocStringExtensions v0.9.5 [4e289a0a] EnumX v1.0.6 [f151be2c] EnzymeCore v0.8.18 [e24c0720] ExponentialAction v0.2.10 [d4d017d3] ExponentialUtilities v1.30.0 [e2ba6199] ExprTools v0.1.10 [55351af7] ExproniconLite v0.10.14 [7034ab61] FastBroadcast v0.3.5 [9aa1b823] FastClosures v0.3.2 [a4df4552] FastPower v1.3.0 [5789e2e9] FileIO v1.17.1 [1a297f60] FillArrays v1.16.0 [64ca27bc] FindFirstFunctions v1.8.0 [6a86dc24] FiniteDiff v2.29.0 [f6369f11] ForwardDiff v1.3.1 [069b7b12] FunctionWrappers v1.1.3 [77dc65aa] FunctionWrappersWrappers v0.1.3 [46192b85] GPUArraysCore v0.2.0 [c145ed77] GenericSchur v0.5.6 [076d061b] HashArrayMappedTries v0.2.0 [34004b35] HypergeometricFunctions v0.3.28 [615f187c] IfElse v0.1.1 [3587e190] InverseFunctions v0.1.17 [b6b21f68] Ipopt v1.14.0 [92d709cd] IrrationalConstants v0.2.6 [82899510] IteratorInterfaceExtensions v1.0.0 [033835bb] JLD2 v0.6.3 [692b3bcd] JLLWrappers v1.7.1 [682c06a0] JSON v1.4.0 [0f8b85d8] JSON3 v1.14.3 [ae98c720] Jieko v0.2.1 [b964fa9f] LaTeXStrings v1.4.0 [10f19ff3] LayoutPointers v0.1.17 [5078a376] LazyArrays v2.9.5 [2ab3a3ac] LogExpFunctions v0.3.29 [e6f89c97] LoggingExtras v1.2.0 [1914dd2f] MacroTools v0.5.16 [d125e4d3] ManualMemory v0.1.8 [b8f27783] MathOptInterface v1.48.0 [bb5d69b7] MaybeInplace v0.1.4 [e1d29d7a] Missings v1.2.0 [2e0e35c7] Moshi v0.3.7 [46d2c3a1] MuladdMacro v0.2.4 [d8a4904e] MutableArithmetics v1.6.7 [77ba4419] NaNMath v1.1.3 [538bc3a1] NamedTrajectories v0.8.2 [be0214bd] NonlinearSolveBase v2.11.1 [510215fc] Observables v0.5.5 [bac558e1] OrderedCollections v1.8.1 [bbf590c4] OrdinaryDiffEqCore v3.2.0 [521117fe] OrdinaryDiffEqLinear v1.10.0 [b1df2697] OrdinaryDiffEqTsit5 v1.9.0 [90014a1f] PDMats v0.11.37 [69de0a69] Parsers v2.8.3 [5a402ddf] PiccoloQuantumObjects v0.10.2 [f517fe37] Polyester v0.7.18 [1d0040c9] PolyesterWeave v0.2.2 [d236fae5] PreallocationTools v1.1.0 [aea7be01] PrecompileTools v1.3.3 [21216c6a] Preferences v1.5.1 [08abe8d2] PrettyTables v3.1.2 [92933f4c] ProgressMeter v1.11.0 [43287f4e] PtrArrays v1.3.0 [1fd47b50] QuadGK v2.11.2 [0dc23a59] QuantumCollocation v0.10.2 [3cdcf5f2] RecipesBase v1.3.4 [731186ca] RecursiveArrayTools v3.46.0 [189a3867] Reexport v1.2.2 [ae029012] Requires v1.3.1 [79098fc4] Rmath v0.9.0 [7e49a35a] RuntimeGeneratedFunctions v0.5.16 [94e857df] SIMDTypes v0.1.0 [0bca4576] SciMLBase v2.135.0 [19f34311] SciMLJacobianOperators v0.1.12 [a6db7da4] SciMLLogging v1.8.0 [c0aeaf25] SciMLOperators v1.14.1 [431bcebd] SciMLPublic v1.0.1 [53ae85a6] SciMLStructures v1.10.0 [7e506255] ScopedValues v1.5.0 [efcf1570] Setfield v1.1.2 [a2af1166] SortingAlgorithms v1.2.2 [276daf66] SpecialFunctions v2.6.1 [aedffcd0] Static v1.3.1 [0d7ed370] StaticArrayInterface v1.8.0 [90137ffa] StaticArrays v1.9.16 [1e83bf80] StaticArraysCore v1.4.4 [10745b16] Statistics v1.11.1 [82ae8749] StatsAPI v1.8.0 [2913bbd2] StatsBase v0.34.10 [4c63d2b9] StatsFuns v1.5.2 [7792a7ef] StrideArraysCore v0.5.8 [892a3eda] StringManipulation v0.4.2 [856f2bd8] StructTypes v1.11.0 [ec057cc2] StructUtils v2.6.2 [2efcf032] SymbolicIndexingInterface v0.3.46 [3783bdb8] TableTraits v1.0.1 [bd369af6] Tables v1.12.1 [f8b46487] TestItemRunner v1.1.4 [1c621080] TestItems v1.0.0 [8290d209] ThreadingUtilities v0.5.5 [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.3+0 [6e34b625] Bzip2_jll v1.0.9+0 [e33a78d0] Hwloc_jll v2.12.2+0 [9cc047cb] Ipopt_jll v300.1400.1901+0 [94ce4f54] Libiconv_jll v1.18.0+0 [d00139f3] METIS_jll v5.1.3+0 [d7ed1dd3] MUMPS_seq_jll v500.800.200+0 [656ef2d0] OpenBLAS32_jll v0.3.30+0 [efe28fd5] OpenSpecFun_jll v0.5.6+0 [f50d1b31] Rmath_jll v0.5.1+0 [319450e9] SPRAL_jll v2025.9.18+0 ⌅ [02c8fc9c] XML2_jll v2.13.9+0 [a65dc6b1] Xorg_libpciaccess_jll v0.18.1+0 [0dad84c5] ArgTools v1.1.2 [56f22d72] Artifacts v1.11.0 [2a0f44e3] Base64 v1.11.0 [ade2ca70] Dates v1.11.0 [8ba89e20] Distributed v1.11.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.13.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 [9abbd945] Profile v1.11.0 [3fa0cd96] REPL v1.11.0 [9a3f8284] Random v1.11.0 [ea8e919c] SHA v1.0.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.3.0+1 [deac9b47] LibCURL_jll v8.18.0+0 [e37daf67] LibGit2_jll v1.9.2+0 [29816b5a] LibSSH2_jll v1.11.3+1 [14a3606d] MozillaCACerts_jll v2025.12.2 [4536629a] OpenBLAS_jll v0.3.30+0 [05823500] OpenLibm_jll v0.8.7+0 [458c3c95] OpenSSL_jll v3.5.4+0 [efcefdf7] PCRE2_jll v10.47.0+0 [bea87d4a] SuiteSparse_jll v7.10.1+0 [83775a58] Zlib_jll v1.3.1+2 [3161d3a3] Zstd_jll v1.5.7+1 [8e850b90] libblastrampoline_jll v5.15.0+0 [8e850ede] nghttp2_jll v1.68.0+1 [3f19e933] p7zip_jll v17.7.0+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... [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! ┌ 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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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=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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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! constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt ****************************************************************************** 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, 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, 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, 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, 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, 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, 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/4TeZc/src/problems.jl:65 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, 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, 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, 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, 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, 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, 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, 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, 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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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, 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, 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"##165".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##165".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##165".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##165".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##166".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##166".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##166".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##166".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.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"##167".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##167".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##167".H)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##167".H)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##168".H1)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##168".H1)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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"##168".H1)}, PiccoloQuantumObjects.QuantumSystems.var"#12#13"{typeof(Main.var"##168".H1)}, @NamedTuple{}}}, Val{()}}}, LinearAlgebra.UniformScaling{Bool}, 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.8.2. 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.9924628e+02 2.07e-02 3.17e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 1.9910906e+02 1.05e-02 6.91e+01 -0.3 4.96e-01 2.0 1.00e+00 5.00e-01f 1 2 1.4041109e+02 5.36e-02 2.10e+02 0.1 6.98e-01 1.5 3.86e-01 1.00e+00f 1 3 1.5793327e+02 1.74e-03 1.63e+02 0.6 1.54e-01 2.9 9.45e-01 1.00e+00f 1 4 1.5276644e+02 9.20e-04 6.53e+01 -0.0 3.49e-02 2.4 9.81e-01 1.00e+00f 1 5 1.2837982e+02 2.39e-02 1.45e+01 -0.5 2.05e-01 1.9 9.98e-01 1.00e+00f 1 6 1.2416234e+02 1.05e-02 1.43e+01 -1.2 6.25e-02 2.3 9.97e-01 5.75e-01f 1 7 1.0291749e+02 1.10e-02 8.95e+00 -1.8 1.35e-01 1.8 9.99e-01 1.00e+00f 1 8 9.8533531e+01 1.04e-03 9.08e+00 -2.3 4.94e-02 2.3 9.99e-01 1.00e+00f 1 9 9.7504158e+01 9.94e-04 7.35e+00 -3.1 9.08e-02 1.8 1.00e+00 6.91e-02f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 9.1774919e+01 9.02e-04 6.02e+00 -3.8 3.62e-02 2.2 1.00e+00 1.00e+00f 1 Number of Iterations....: 10 (scaled) (unscaled) Objective...............: 9.1774919330416353e+01 9.1774919330416353e+01 Dual infeasibility......: 6.0244164151651320e+00 6.0244164151651320e+00 Constraint violation....: 9.0205964343931244e-04 9.0205964343931244e-04 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 1.6384965913018711e-03 1.6384965913018711e-03 Overall NLP error.......: 6.0244164151651320e+00 6.0244164151651320e+00 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 = 5.530 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, 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, 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 added GlobalBoundsConstraint for :δ with bounds 0.5 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 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 constructing SmoothPulseProblem for KetTrajectory... applying timesteps_all_equal constraint: Δt This is Ipopt version 3.14.19, running with linear solver MUMPS 5.8.2. 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 3.4447362e+01 4.26e+01 1.72e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 3.7700980e+01 2.66e+01 1.40e+02 -0.3 2.28e+01 - 3.20e-02 3.77e-01f 1 2 3.7521900e+01 2.36e+01 4.33e+02 0.2 1.40e+01 2.0 9.00e-01 1.10e-01h 1 3 3.1342024e+01 1.88e+01 3.50e+02 -4.0 1.26e+01 - 9.50e-02 2.03e-01h 1 4 3.8717397e+01 7.06e+00 3.60e+02 0.5 1.02e+01 1.5 6.84e-01 6.26e-01f 1 5 4.3731469e+01 2.31e-02 4.11e+02 -4.0 3.91e+00 - 3.74e-01 1.00e+00h 1 6 3.4002363e+01 1.02e-02 1.06e+02 0.5 5.61e-01 1.0 1.00e+00 1.00e+00f 1 7 3.2684753e+01 8.87e-04 6.04e+01 0.2 1.74e-01 1.5 8.56e-01 1.00e+00f 1 8 2.5253281e+01 1.26e-03 2.90e+00 -0.3 1.45e-01 1.0 9.46e-01 1.00e+00f 1 9 3.1962339e-02 9.29e-03 3.02e+00 -0.7 4.96e-01 0.5 9.91e-01 9.51e-01f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 3.5540528e-01 5.29e-03 2.05e+02 -1.3 2.42e-01 0.0 9.95e-01 4.49e-01h 1 11 1.8928230e+00 1.31e-02 1.35e+01 0.7 1.79e+01 -0.4 1.60e-01 6.96e-02f 1 12 2.8776425e+00 2.29e-02 1.65e+02 -0.2 9.43e-01 - 8.76e-01 1.00e+00f 1 13 5.7985166e+00 1.41e-03 1.46e+02 -0.2 3.33e-01 0.9 3.63e-01 1.00e+00f 1 14 3.4595365e+00 2.38e-04 2.99e+00 -1.1 6.22e-02 1.3 1.00e+00 1.00e+00f 1 15 8.6588791e-01 5.60e-04 6.27e-01 -2.3 9.31e-02 0.8 1.00e+00 1.00e+00f 1 16 4.6865409e-01 6.62e-04 1.64e+02 -3.2 9.79e-02 0.4 1.00e+00 1.00e+00f 1 17 4.7491372e-01 1.62e-04 1.63e+02 -2.1 5.17e-02 -0.1 1.00e+00 1.00e+00f 1 18 1.5211365e-01 1.10e-04 2.63e+00 -3.4 3.85e-02 1.2 1.00e+00 1.00e+00f 1 19 4.0523045e-02 5.41e-05 1.46e-01 -4.0 2.67e-02 0.7 1.00e+00 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 20 3.1620410e-02 5.05e-05 1.60e+02 -4.0 2.64e-02 0.3 1.00e+00 1.00e+00f 1 21 2.5278562e-02 8.44e-06 1.60e+02 -4.0 1.13e-02 -0.2 1.00e+00 1.00e+00f 1 22 1.5735019e-03 1.09e-05 1.10e+00 -4.0 1.09e-02 1.1 1.00e+00 1.00e+00f 1 23 1.5599616e-03 3.12e-06 2.75e-02 -4.0 6.35e-03 0.6 1.00e+00 1.00e+00h 1 24 1.8940242e-03 2.34e-06 1.59e+02 -4.0 5.48e-03 0.2 1.00e+00 1.00e+00h 1 25 1.3311891e-03 3.18e-07 1.59e+02 -4.0 2.07e-03 -0.3 1.00e+00 1.00e+00h 1 26 7.8475518e-04 1.05e-08 7.04e-01 -4.0 3.73e-03 1.0 1.00e+00 1.00e+00H 1 27 6.3047713e-04 8.07e-08 1.59e+02 -4.0 9.35e-04 0.5 1.00e+00 1.00e+00h 1 28 7.1598583e-04 3.93e-08 1.59e+02 -4.0 6.28e-04 0.1 1.00e+00 1.00e+00h 1 29 6.9433011e-04 7.02e-09 2.29e-02 -4.1 2.75e-04 1.4 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.4889587e-04 1.33e-08 3.11e-03 -4.0 3.83e-04 0.9 1.00e+00 1.00e+00h 1 31 6.2420391e-04 2.07e-08 1.59e+02 -4.0 4.89e-04 0.4 1.00e+00 1.00e+00h 1 32 6.3818771e-04 9.17e-09 1.59e+02 -4.1 4.97e-04 -0.0 1.00e+00 1.00e+00h 1 33 6.3242642e-04 2.10e-09 1.07e-02 -4.1 1.51e-04 1.3 1.00e+00 1.00e+00h 1 34 6.2371307e-04 2.88e-09 1.16e-03 -4.0 1.81e-04 0.8 1.00e+00 1.00e+00h 1 35 6.2014789e-04 3.40e-09 1.59e+02 -4.0 2.12e-04 0.3 1.00e+00 1.00e+00h 1 36 6.2256665e-04 2.48e-09 1.59e+02 -4.1 4.86e-04 -0.1 1.00e+00 1.00e+00h 1 37 6.2098431e-04 6.34e-10 5.02e-03 -4.1 8.99e-05 1.2 1.00e+00 1.00e+00h 1 38 6.1984785e-04 3.95e-10 3.83e-04 -4.0 7.56e-05 0.7 1.00e+00 1.00e+00h 1 39 6.1913572e-04 3.02e-10 2.28e-04 -4.0 1.35e-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.1916362e-04 6.45e-10 3.04e-04 -4.0 5.41e-04 -0.2 1.00e+00 1.00e+00h 1 41 6.2034252e-04 1.05e-07 1.59e+02 -4.1 8.42e+00 - 9.86e-02 2.46e-04h 10 42 6.1956114e-04 3.50e-09 1.59e+02 -4.1 1.44e-03 -0.7 1.00e+00 1.00e+00h 1 43 6.1973243e-04 1.53e-11 9.90e-04 -4.1 1.76e-05 1.5 1.00e+00 1.00e+00h 1 44 6.1946426e-04 2.29e-11 2.46e-04 -4.0 2.31e-05 1.0 1.00e+00 1.00e+00h 1 45 6.1918148e-04 5.33e-11 2.30e-04 -4.0 6.46e-05 0.6 1.00e+00 1.00e+00h 1 46 6.1907069e-04 1.95e-10 2.92e-04 -4.0 2.46e-04 0.1 1.00e+00 1.00e+00h 1 47 6.1934313e-04 2.54e-08 1.59e+02 -4.0 1.83e+01 - 4.64e-02 5.59e-05f 11 48 6.1949374e-04 1.44e-09 1.59e+02 -4.1 6.94e-04 -0.4 1.00e+00 1.00e+00h 1 49 6.1979997e-04 2.75e-12 4.96e-04 -4.1 7.79e-06 1.8 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 50 6.1950241e-04 2.33e-11 4.92e-04 -6.2 2.19e-05 1.4 1.00e+00 1.00e+00h 1 51 6.1906794e-04 7.56e-11 2.97e-04 -6.2 3.96e-05 0.9 1.00e+00 1.00e+00h 1 52 6.1903645e-04 5.04e-11 2.75e-04 -4.0 1.10e-04 0.4 1.00e+00 1.00e+00h 1 53 6.1896941e-04 2.73e-10 3.00e-04 -4.0 3.60e-04 -0.1 1.00e+00 1.00e+00h 1 54 6.2056736e-04 1.08e-07 1.59e+02 -4.0 7.00e+00 - 1.22e-01 2.95e-04f 10 55 6.1930878e-04 1.59e-09 1.59e+02 -4.1 9.55e-04 -0.6 1.00e+00 1.00e+00h 1 56 6.1955418e-04 5.57e-12 5.77e-04 -4.1 1.10e-05 1.7 1.00e+00 1.00e+00h 1 57 6.1921218e-04 3.54e-11 4.28e-04 -6.2 2.71e-05 1.2 1.00e+00 1.00e+00h 1 58 6.1882429e-04 8.22e-11 2.21e-04 -6.2 4.20e-05 0.7 1.00e+00 1.00e+00h 1 59 6.1887990e-04 6.35e-11 2.92e-04 -4.0 1.66e-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 60 6.2014266e-04 1.06e-07 1.59e+02 -4.1 7.27e+00 - 1.14e-01 2.85e-04f 10 61 6.1922129e-04 1.07e-09 1.59e+02 -4.1 4.83e-04 -0.2 1.00e+00 1.00e+00h 1 62 6.1872346e-04 3.38e-10 1.59e+02 -4.1 6.85e-05 1.1 1.00e+00 1.00e+00h 1 63 6.1931801e-04 2.38e-10 1.59e+02 -4.1 1.00e-04 0.6 1.00e+00 1.00e+00h 1 64 6.1883864e-04 7.13e-14 1.51e-02 -4.1 8.46e-05 1.0 1.00e+00 1.00e+00H 1 65 6.1879398e-04 1.57e-11 2.18e-04 -4.1 5.91e-05 0.6 1.00e+00 1.00e+00h 1 66 6.4573290e-04 1.83e-03 1.59e+02 -4.0 1.79e+02 - 4.87e-03 1.61e-03f 3 67 1.1684422e-01 1.38e-04 1.58e+02 -4.2 5.25e-02 0.1 1.00e+00 1.00e+00h 1 68 6.2282787e-03 6.62e-07 1.89e-01 -4.2 3.25e-03 1.4 1.00e+00 1.00e+00f 1 69 2.4207717e-03 8.09e-07 3.12e-02 -4.2 3.56e-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 70 1.2765125e-03 1.05e-06 1.59e+02 -4.2 3.88e-03 0.5 1.00e+00 1.00e+00h 1 71 1.6214187e-03 2.19e-07 1.59e+02 -4.2 2.04e-03 -0.0 1.00e+00 1.00e+00h 1 72 1.1022815e-03 1.41e-07 9.16e-02 -4.2 1.36e-03 1.3 1.00e+00 1.00e+00h 1 73 7.4115577e-04 1.06e-07 8.80e-03 -4.2 1.27e-03 0.8 1.00e+00 1.00e+00h 1 74 6.9315406e-04 1.14e-07 1.59e+02 -4.2 1.34e-03 0.4 1.00e+00 1.00e+00h 1 75 6.8440944e-04 1.73e-08 1.59e+02 -4.2 5.72e-04 -0.1 1.00e+00 1.00e+00h 1 76 6.3372703e-04 1.81e-08 3.67e-02 -4.2 4.79e-04 1.2 1.00e+00 1.00e+00h 1 77 6.2027607e-04 8.26e-09 1.93e-03 -4.2 3.53e-04 0.7 1.00e+00 1.00e+00h 1 78 6.2088467e-04 7.68e-09 1.59e+02 -4.2 3.29e-04 0.3 1.00e+00 1.00e+00h 1 79 6.1825584e-04 7.21e-10 1.59e+02 -4.2 1.33e-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 80 6.1617354e-04 1.06e-09 1.59e+02 -4.2 1.11e-04 1.1 1.00e+00 1.00e+00h 1 81 6.1753392e-04 2.70e-10 1.59e+02 -4.2 5.81e-05 0.6 1.00e+00 1.00e+00h 1 82 6.1625176e-04 7.75e-10 1.59e+02 -4.2 9.15e-05 1.1 1.00e+00 1.00e+00h 1 83 6.1686273e-04 2.05e-10 1.59e+02 -4.2 5.04e-05 0.6 1.00e+00 1.00e+00h 1 84 6.1634623e-04 6.06e-10 1.59e+02 -4.2 8.23e-05 1.0 1.00e+00 1.00e+00h 1 85 6.1631331e-04 1.02e-10 1.59e+02 -4.2 4.24e-05 0.5 1.00e+00 1.00e+00h 1 86 6.1634739e-04 1.50e-12 3.49e-04 -4.2 4.10e-06 1.9 1.00e+00 1.00e+00h 1 87 6.1628482e-04 7.57e-12 2.03e-04 -4.0 8.36e-06 1.4 1.00e+00 1.00e+00h 1 88 6.1617988e-04 3.14e-11 1.80e-04 -4.0 2.22e-05 0.9 1.00e+00 1.00e+00h 1 89 6.1610596e-04 5.39e-11 1.59e-04 -4.0 5.90e-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 90 6.1623124e-04 8.41e-09 1.59e+02 -4.0 4.45e+01 - 1.61e-02 9.66e-06f 12 91 6.4325715e-04 3.14e-03 1.57e+02 -4.1 2.56e+01 - 2.92e-02 1.02e-02f 3 92 1.2291292e-01 4.65e-04 1.57e+02 -4.1 7.89e-02 -0.0 1.00e+00 1.00e+00h 1 93 2.1071020e-03 1.51e-06 3.94e-01 -4.1 4.32e-03 1.3 1.00e+00 1.00e+00f 1 94 9.2977326e-04 6.60e-07 1.99e-02 -4.1 3.12e-03 0.8 1.00e+00 1.00e+00h 1 95 9.5717175e-04 5.40e-07 1.58e+02 -4.1 2.78e-03 0.3 1.00e+00 1.00e+00h 1 96 8.2349973e-04 7.66e-08 1.58e+02 -4.1 1.03e-03 -0.1 1.00e+00 1.00e+00h 1 97 6.3321604e-04 8.28e-08 7.20e-02 -4.1 1.08e-03 1.2 1.00e+00 1.00e+00h 1 98 6.2795167e-04 2.41e-08 2.87e-03 -4.1 5.69e-04 0.7 1.00e+00 1.00e+00h 1 99 6.3399667e-04 1.62e-08 1.58e+02 -4.1 4.55e-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 100 6.2887788e-04 2.18e-09 1.58e+02 -4.1 1.64e-04 -0.3 1.00e+00 1.00e+00h 1 Number of Iterations....: 100 (scaled) (unscaled) Objective...............: 6.2887788495458229e-04 6.2887788495458229e-04 Dual infeasibility......: 1.5761881282646993e+02 1.5761881282646993e+02 Constraint violation....: 2.1767629920788067e-09 2.1767629920788067e-09 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 7.9423455660686823e-05 7.9423455660686823e-05 Overall NLP error.......: 1.5761881282646993e+02 1.5761881282646993e+02 Number of objective function evaluations = 168 Number of objective gradient evaluations = 101 Number of equality constraint evaluations = 168 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 = 54.642 EXIT: Maximum Number of Iterations Exceeded. constructing SmoothPulseProblem for UnitaryTrajectory... applying timesteps_all_equal constraint: Δt This is Ipopt version 3.14.19, running with linear solver MUMPS 5.8.2. 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.6733944e+01 7.55e+00 1.63e+00 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 5.0285613e+01 2.01e+00 3.33e+01 -0.6 4.26e+00 - 2.78e-01 7.34e-01f 1 2 1.8884195e+01 1.45e+00 1.20e+02 -0.4 4.05e+00 0.0 1.18e-01 2.80e-01f 1 3 3.1129280e+01 6.68e-04 1.54e+02 0.3 1.14e+00 1.3 7.39e-01 1.00e+00f 1 4 3.0431647e+01 1.36e-05 1.98e+01 -0.4 4.46e-02 1.8 9.98e-01 1.00e+00f 1 5 1.5483309e+01 3.25e-03 6.45e+01 -0.0 1.86e+00 - 1.00e+00 4.89e-01f 1 6 1.2815623e+01 4.46e-04 7.88e+00 -0.3 4.51e-01 - 9.84e-01 1.00e+00f 1 7 1.1274150e+01 3.91e-04 3.00e+02 -0.5 1.30e+00 - 9.69e-01 1.16e-01f 2 8 5.1581766e+00 1.27e-03 1.50e+02 -0.7 7.83e-01 - 1.00e+00 7.17e-01f 1 9 6.1855150e+00 1.11e-03 4.28e+02 -0.5 1.51e+00 - 1.00e+00 2.74e-01h 2 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 1.2468589e+01 3.94e-04 1.53e+03 0.2 2.96e-01 1.3 3.88e-01 1.00e+00f 1 11 1.1140764e+01 1.24e-04 2.45e+00 -0.6 1.13e-01 - 1.00e+00 1.00e+00f 1 12 8.2286606e+00 1.26e-04 1.28e+00 -1.5 8.95e-02 0.8 1.00e+00 1.00e+00f 1 13 2.6577167e+00 7.02e-04 1.83e+00 -1.9 2.07e-01 0.3 9.95e-01 1.00e+00f 1 14 1.5218417e+00 9.57e-04 5.60e+00 -1.4 2.33e+00 - 1.00e+00 3.25e-01f 1 15 3.4746245e+00 2.35e-03 3.93e+00 -0.9 1.07e+00 - 6.91e-01 9.23e-01f 1 16 2.4085653e+00 8.12e-04 3.73e-01 -1.1 5.46e-01 - 1.00e+00 1.00e+00f 1 17 1.1413248e+00 1.83e-04 4.79e-01 -1.9 2.14e-01 - 9.83e-01 1.00e+00f 1 18 5.9131732e-01 1.26e-04 1.71e-01 -3.0 2.72e-01 - 9.95e-01 1.00e+00f 1 19 4.6657444e-01 1.93e-05 3.47e-02 -4.0 9.95e-02 - 1.00e+00 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 20 4.4724860e-01 8.31e-07 1.42e-03 -4.0 2.12e-02 - 1.00e+00 1.00e+00h 1 21 4.4647263e-01 1.42e-09 4.31e-06 -4.0 8.58e-04 - 1.00e+00 1.00e+00h 1 22 4.4646514e-01 2.03e-13 9.57e-09 -4.0 2.13e-05 - 1.00e+00 1.00e+00h 1 Number of Iterations....: 22 (scaled) (unscaled) Objective...............: 4.4646514488186911e-01 4.4646514488186911e-01 Dual infeasibility......: 9.5719518366198242e-09 9.5719518366198242e-09 Constraint violation....: 2.0328920838363906e-13 2.0328920838363906e-13 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 1.0000000286201575e-04 1.0000000286201575e-04 Overall NLP error.......: 9.5719518366198242e-09 9.5719518366198242e-09 Number of objective function evaluations = 28 Number of objective gradient evaluations = 23 Number of equality constraint evaluations = 28 Number of inequality constraint evaluations = 0 Number of equality constraint Jacobian evaluations = 23 Number of inequality constraint Jacobian evaluations = 0 Number of Lagrangian Hessian evaluations = 22 Total seconds in IPOPT = 46.763 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.8.2. 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.2058485e+02 4.90e-02 2.01e+01 0.0 0.00e+00 - 0.00e+00 0.00e+00 0 1 8.8549535e+01 7.33e-02 3.56e+02 -0.4 4.04e-01 2.0 9.85e-01 9.81e-01f 1 2 8.2424777e+01 2.26e-03 6.40e+01 -0.7 1.19e-01 2.4 1.00e+00 1.00e+00f 1 3 7.6923712e+01 2.03e-03 4.89e+01 -1.4 1.70e-01 1.9 1.00e+00 2.25e-01f 1 4 1.5197680e+01 7.95e-02 5.87e+01 -1.8 5.91e-01 1.5 1.00e+00 1.00e+00f 1 5 1.5285690e+01 3.18e-02 2.46e+01 -2.0 2.78e-01 1.0 1.00e+00 5.27e-01h 1 6 9.9105895e+00 8.91e-03 4.45e+00 -4.0 1.71e-01 0.5 1.03e-01 1.00e+00f 1 7 1.0678191e+01 1.51e-03 4.39e-01 -2.3 5.37e-02 0.9 1.00e+00 1.00e+00h 1 8 5.6758825e+00 6.99e-03 2.68e+00 -1.6 2.41e-01 0.5 9.61e-01 1.00e+00f 1 9 4.7311120e+00 5.97e-04 4.93e-01 -2.2 3.90e-02 0.9 9.83e-01 1.00e+00f 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 10 2.0776065e+00 3.29e-03 5.10e-01 -2.3 8.70e-02 0.4 9.97e-01 1.00e+00f 1 11 7.1370979e-02 1.21e-02 2.30e+00 -1.9 1.76e-01 -0.1 1.00e+00 9.84e-01f 1 12 2.6754670e-01 3.44e-04 2.35e-01 -2.6 4.34e-02 0.4 1.00e+00 1.00e+00h 1 13 1.7736086e-02 6.86e-04 2.10e-01 -3.4 5.15e-02 -0.1 9.98e-01 1.00e+00f 1 14 6.7491031e-03 4.43e-04 4.38e+00 -4.0 4.27e-02 -0.6 1.00e+00 5.46e-01h 1 15 9.3624873e-03 1.05e-04 1.91e-01 -4.0 1.94e-02 -0.2 1.00e+00 1.00e+00h 1 16 6.3986930e-03 6.43e-05 3.28e+00 -4.0 1.57e-02 -0.6 1.00e+00 8.02e-01h 1 17 8.0839586e-03 6.86e-06 2.93e-01 -4.0 3.87e-03 -0.2 1.00e+00 1.00e+00h 1 18 6.5081169e-03 1.07e-05 1.89e-01 -4.0 7.06e-03 -0.7 1.00e+00 1.00e+00h 1 19 6.6464856e-03 1.16e-06 5.76e-02 -4.0 2.67e-03 -0.3 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 20 6.4740646e-03 4.55e-06 1.45e-01 -4.0 5.86e-03 -0.7 1.00e+00 1.00e+00h 1 21 6.3526579e-03 4.41e-05 2.62e-01 -4.0 1.58e-02 -1.2 1.00e+00 1.00e+00h 1 22 6.3376562e-03 3.93e-06 3.38e-02 -4.0 5.40e-03 -0.8 1.00e+00 1.00e+00h 1 23 6.2042292e-03 5.14e-05 1.66e-01 -4.0 1.44e-02 -1.3 1.00e+00 1.00e+00h 1 24 6.1662638e-03 4.32e-06 5.98e-03 -4.0 5.20e-03 -0.8 1.00e+00 1.00e+00h 1 25 6.0682750e-03 4.69e-05 3.27e-02 -4.0 1.58e-02 -1.3 1.00e+00 1.00e+00h 1 26 5.7563354e-03 4.64e-04 4.41e-01 -4.0 4.60e-02 -1.8 1.00e+00 1.00e+00h 1 27 5.7425319e-03 7.73e-07 2.77e-03 -4.0 1.85e-03 -0.5 1.00e+00 1.00e+00h 1 28 5.7333237e-03 6.71e-06 4.72e-02 -4.0 6.26e-03 -0.9 1.00e+00 1.00e+00h 1 29 5.6271051e-03 7.55e-05 5.58e-02 -4.0 1.99e-02 -1.4 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 30 5.5845927e-03 7.05e-06 2.24e-02 -4.0 6.94e-03 -1.0 1.00e+00 1.00e+00h 1 31 5.4980822e-03 6.29e-05 6.45e-02 -4.0 2.11e-02 -1.5 1.00e+00 1.00e+00h 1 32 5.4596603e-03 9.84e-06 6.56e-03 -4.0 8.05e-03 -1.1 1.00e+00 1.00e+00h 1 33 5.3584152e-03 8.22e-05 7.64e-03 -4.0 2.36e-02 -1.5 1.00e+00 1.00e+00h 1 34 5.0948278e-03 6.28e-04 2.25e-01 -4.0 6.01e-02 -2.0 1.00e+00 1.00e+00h 1 35 5.0810489e-03 1.93e-06 6.82e-04 -4.0 2.21e-03 -0.7 1.00e+00 1.00e+00h 1 36 5.0631711e-03 1.78e-05 4.48e-02 -4.0 9.37e-03 -1.2 1.00e+00 1.00e+00h 1 37 4.9518874e-03 2.11e-04 7.76e-02 -4.0 2.91e-02 -1.6 1.00e+00 1.00e+00h 1 38 4.9118134e-03 2.03e-05 2.74e-02 -4.0 8.73e-03 -1.2 1.00e+00 1.00e+00h 1 39 4.8245205e-03 1.75e-04 1.17e-01 -4.0 2.79e-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 40 4.7854053e-03 3.08e-05 5.25e-03 -4.0 1.11e-02 -1.3 1.00e+00 1.00e+00h 1 41 4.6773028e-03 2.18e-04 8.77e-02 -4.0 3.01e-02 -1.7 1.00e+00 1.00e+00h 1 42 4.6436489e-03 2.95e-05 1.32e-02 -4.0 1.04e-02 -1.3 1.00e+00 1.00e+00h 1 43 4.5479809e-03 2.88e-04 6.25e-02 -4.0 3.35e-02 -1.8 1.00e+00 1.00e+00h 1 44 4.5077841e-03 3.42e-05 2.81e-02 -4.0 9.45e-03 -1.4 1.00e+00 1.00e+00h 1 45 4.4053444e-03 2.72e-04 4.87e-03 -4.0 3.03e-02 -1.8 1.00e+00 1.00e+00h 1 46 4.2004622e-03 1.19e-03 9.94e-02 -4.0 6.74e-02 -2.3 1.00e+00 1.00e+00h 1 47 4.1860823e-03 6.65e-06 2.54e-03 -4.0 3.20e-03 -1.0 1.00e+00 1.00e+00h 1 48 4.1571146e-03 3.91e-05 4.60e-03 -4.0 1.03e-02 -1.5 1.00e+00 1.00e+00h 1 49 4.0787505e-03 2.93e-04 5.88e-03 -4.0 2.79e-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 50 3.9359790e-03 6.62e-04 6.33e-02 -4.0 5.69e-02 -2.4 1.00e+00 1.00e+00h 1 51 3.9244201e-03 4.61e-06 5.12e-04 -4.0 3.02e-03 -1.1 1.00e+00 1.00e+00h 1 52 3.8975925e-03 3.35e-05 2.51e-02 -4.0 9.22e-03 -1.6 1.00e+00 1.00e+00h 1 53 3.8364886e-03 2.68e-04 4.30e-02 -4.0 3.34e-02 -2.0 1.00e+00 1.00e+00h 1 54 3.8225014e-03 4.30e-05 1.49e-02 -4.0 9.55e-03 -1.6 1.00e+00 1.00e+00h 1 55 3.7586648e-03 2.61e-04 1.78e-01 -4.0 2.60e-02 -2.1 1.00e+00 1.00e+00h 1 56 3.7616636e-03 1.02e-06 5.89e-03 -4.0 1.46e-03 -0.8 1.00e+00 1.00e+00h 1 57 3.7625695e-03 4.22e-06 9.36e-03 -4.0 3.96e-03 -1.2 1.00e+00 1.00e+00h 1 58 3.6475788e-03 2.35e-05 1.08e-02 -4.1 1.10e-02 -1.7 1.00e+00 1.00e+00h 1 59 3.6836781e-03 3.07e-04 9.36e-03 -4.0 3.67e-02 -2.2 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 60 3.5761954e-03 2.25e-05 7.14e-03 -4.1 1.11e-02 -1.8 1.00e+00 1.00e+00h 1 61 3.6299471e-03 2.48e-04 2.38e-02 -4.0 3.22e-02 -2.2 1.00e+00 1.00e+00h 1 62 3.5170157e-03 3.39e-05 5.50e-03 -4.1 1.20e-02 -1.8 1.00e+00 1.00e+00h 1 63 3.5894642e-03 1.53e-04 1.21e-01 -4.0 3.71e-02 -2.3 1.00e+00 1.00e+00h 1 64 3.4945268e-03 9.14e-07 5.44e-03 -4.1 1.90e-03 -1.0 1.00e+00 1.00e+00h 1 65 3.5596569e-03 3.82e-06 3.48e-02 -4.0 5.54e-03 -1.4 1.00e+00 1.00e+00h 1 66 3.4585142e-03 1.03e-04 8.16e-02 -4.1 1.81e-02 -1.9 1.00e+00 1.00e+00h 1 67 3.5426122e-03 9.42e-06 1.00e-02 -4.0 6.42e-03 -1.5 1.00e+00 1.00e+00h 1 68 3.4283644e-03 2.75e-05 2.70e-02 -4.1 1.41e-02 -2.0 1.00e+00 1.00e+00h 1 69 3.5108543e-03 4.27e-06 5.80e-04 -4.0 6.80e-03 -1.5 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 70 3.4081856e-03 2.89e-05 1.05e-02 -4.1 1.56e-02 -2.0 1.00e+00 1.00e+00h 1 71 3.5423805e-03 6.22e-04 5.40e-02 -4.0 5.48e-02 -2.5 1.00e+00 1.00e+00h 1 72 3.4491670e-03 3.29e-06 2.74e-03 -4.1 2.05e-03 -1.2 1.00e+00 1.00e+00h 1 73 3.5219886e-03 5.92e-06 1.47e-02 -4.0 7.63e-03 -1.6 1.00e+00 1.00e+00h 1 74 3.4077427e-03 1.41e-04 3.14e-02 -4.1 2.13e-02 -2.1 1.00e+00 1.00e+00h 1 75 3.4913566e-03 7.01e-06 1.31e-02 -4.0 8.25e-03 -1.7 1.00e+00 1.00e+00h 1 76 3.3632816e-03 8.19e-05 6.88e-02 -4.1 1.90e-02 -2.2 1.00e+00 1.00e+00h 1 77 3.4529461e-03 1.86e-04 7.48e-03 -4.0 1.78e-02 -1.7 1.00e+00 1.00e+00h 1 78 3.3681576e-03 7.90e-05 2.26e-02 -4.1 2.04e-02 -2.2 1.00e+00 1.00e+00h 1 79 3.4576720e-03 1.00e-05 1.20e-02 -4.0 9.56e-03 -1.8 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 80 3.3316303e-03 1.35e-04 7.20e-02 -4.1 2.19e-02 -2.3 1.00e+00 1.00e+00h 1 81 3.4097206e-03 4.31e-07 5.16e-03 -4.0 1.30e-03 -0.9 1.00e+00 1.00e+00h 1 82 3.4183120e-03 2.49e-06 1.02e-02 -4.0 3.86e-03 -1.4 1.00e+00 1.00e+00h 1 83 3.3294349e-03 1.27e-05 2.13e-03 -4.1 9.12e-03 -1.9 1.00e+00 1.00e+00h 1 84 3.4448180e-03 1.54e-04 8.62e-03 -4.0 3.30e-02 -2.4 1.00e+00 1.00e+00h 1 85 3.3477360e-03 1.76e-05 1.13e-03 -4.1 9.67e-03 -2.0 1.00e+00 1.00e+00h 1 86 3.3165779e-03 1.22e-04 4.05e-02 -4.1 2.80e-02 -2.4 1.00e+00 1.00e+00h 1 87 3.4010458e-03 2.85e-07 5.99e-04 -4.0 1.64e-03 -1.1 1.00e+00 1.00e+00h 1 88 3.4062442e-03 3.20e-06 5.63e-04 -4.0 4.89e-03 -1.6 1.00e+00 1.00e+00h 1 89 3.3152936e-03 1.82e-05 3.47e-03 -4.1 1.17e-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 3.4511668e-03 1.91e-04 7.72e-02 -4.0 4.15e-02 -2.5 1.00e+00 1.00e+00h 1 91 3.3618315e-03 3.09e-06 2.09e-03 -4.1 2.03e-03 -1.2 1.00e+00 1.00e+00h 1 92 3.4359274e-03 3.18e-06 2.70e-02 -4.0 5.75e-03 -1.7 1.00e+00 1.00e+00h 1 93 3.3603168e-03 9.12e-07 9.50e-03 -4.1 1.90e-03 -1.3 1.00e+00 1.00e+00h 1 94 3.4414265e-03 3.44e-06 1.04e-02 -4.0 6.41e-03 -1.7 1.00e+00 1.00e+00h 1 95 3.3389430e-03 8.68e-05 1.14e-02 -4.1 1.78e-02 -2.2 1.00e+00 1.00e+00h 1 96 3.3323097e-03 8.22e-06 5.91e-03 -4.1 6.04e-03 -1.8 1.00e+00 1.00e+00h 1 97 3.4279504e-03 3.99e-05 4.45e-02 -4.0 2.07e-02 -2.3 1.00e+00 1.00e+00h 1 98 3.3528241e-03 3.10e-07 3.50e-03 -4.1 1.22e-03 -0.9 1.00e+00 1.00e+00h 1 99 3.4299056e-03 7.10e-07 1.02e-02 -4.0 2.89e-03 -1.4 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 100 3.3400923e-03 1.24e-05 8.97e-03 -4.1 7.94e-03 -1.9 1.00e+00 1.00e+00h 1 101 3.3086351e-03 1.25e-04 1.10e-02 -4.1 2.32e-02 -2.4 1.00e+00 1.00e+00h 1 102 3.3542731e-03 4.26e-04 5.33e+00 -4.0 7.51e-01 -2.8 3.41e-01 3.03e-02h 4 103 3.2481048e-03 5.19e-04 1.88e-02 -4.2 2.47e-02 -2.4 1.00e+00 1.00e+00h 1 104 3.2525912e-03 2.63e-03 2.54e-01 -4.0 2.32e+00 -2.9 1.16e-01 5.74e-02h 3 105 3.1848932e-03 4.45e-05 1.24e-03 -4.2 9.43e-03 -1.6 1.00e+00 1.00e+00h 1 106 3.3702419e-03 1.47e-05 2.42e-02 -4.0 1.15e-02 -2.0 1.00e+00 1.00e+00h 1 107 3.2880975e-03 7.66e-06 8.67e-03 -4.1 4.23e-03 -1.6 1.00e+00 1.00e+00h 1 108 3.3870088e-03 1.59e-05 9.09e-04 -4.0 1.33e-02 -2.1 1.00e+00 1.00e+00h 1 109 3.2781097e-03 1.85e-04 8.94e-03 -4.1 3.61e-02 -2.6 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 110 3.3185795e-03 7.18e-03 1.12e-01 -4.1 1.42e-01 -3.0 1.00e+00 1.00e+00h 1 111 3.3010322e-03 3.03e-05 8.64e-03 -4.1 1.38e-02 -0.8 1.00e+00 1.00e+00h 1 112 3.2681913e-03 2.90e-06 4.72e-02 -4.1 3.80e-03 -1.3 1.00e+00 1.00e+00h 1 113 3.3760245e-03 1.02e-08 2.46e-03 -4.0 2.09e-04 0.9 1.00e+00 1.00e+00h 1 114 3.3774110e-03 1.68e-09 6.11e-04 -4.0 4.82e-05 0.5 1.00e+00 1.00e+00h 1 115 3.3753949e-03 1.58e-08 4.08e-04 -4.0 1.55e-04 -0.0 1.00e+00 1.00e+00h 1 116 3.3749720e-03 1.43e-07 3.07e-04 -4.0 4.71e-04 -0.5 1.00e+00 1.00e+00h 1 117 3.3751712e-03 1.25e-06 1.39e-04 -4.0 1.39e-03 -1.0 1.00e+00 1.00e+00h 1 118 3.2907344e-03 7.83e-06 5.50e-03 -4.1 3.11e-03 -1.4 1.00e+00 1.00e+00h 1 119 3.2779759e-03 6.29e-05 1.37e-02 -4.1 1.03e-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 3.2576164e-03 2.70e-04 4.12e-03 -4.1 2.27e-02 -2.4 1.00e+00 1.00e+00h 1 121 3.2189076e-03 5.12e-04 5.92e-02 -4.1 5.14e-02 -2.9 1.00e+00 1.00e+00h 1 122 3.2169416e-03 3.09e-06 3.97e-03 -4.1 2.50e-03 -1.5 1.00e+00 1.00e+00h 1 123 3.3127502e-03 2.98e-05 2.07e-02 -4.0 9.53e-03 -2.0 1.00e+00 1.00e+00h 1 124 3.2331547e-03 4.84e-06 6.67e-03 -4.1 3.97e-03 -1.6 1.00e+00 1.00e+00h 1 125 3.2338704e-03 1.25e-05 2.62e-04 -4.1 8.00e-03 -2.1 1.00e+00 1.00e+00h 1 126 2.7372941e-03 1.08e-03 8.08e-02 -6.1 5.33e-02 -2.6 8.67e-01 1.00e+00h 1 127 3.1459308e-03 8.75e-06 1.82e-03 -4.0 4.12e-03 -1.2 1.00e+00 1.00e+00h 1 128 3.1304067e-03 2.09e-05 1.23e-01 -4.0 8.37e-03 -1.7 1.00e+00 1.00e+00h 1 129 3.1776943e-03 4.94e-07 4.97e-02 -4.0 8.94e-04 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 130 3.1764200e-03 1.44e-07 4.00e-03 -4.0 4.53e-04 0.1 1.00e+00 1.00e+00h 1 131 3.1679953e-03 1.14e-07 1.82e-03 -4.0 5.86e-04 -0.4 1.00e+00 1.00e+00h 1 132 3.1597584e-03 6.52e-07 6.80e-03 -4.0 1.00e-03 -0.9 1.00e+00 1.00e+00h 1 133 3.1659911e-03 5.03e-06 2.75e-03 -4.0 2.58e-03 -1.4 1.00e+00 1.00e+00h 1 134 3.0943557e-03 1.35e-05 1.17e-03 -4.1 5.13e-03 -1.9 1.00e+00 1.00e+00h 1 135 3.2322154e-03 2.12e-04 2.15e-02 -4.0 1.76e-02 -2.3 1.00e+00 1.00e+00h 1 136 3.1522338e-03 1.78e-05 2.54e-03 -4.1 6.72e-03 -1.9 1.00e+00 1.00e+00h 1 137 3.1699197e-03 1.68e-05 1.92e-02 -4.1 1.55e-02 -2.4 1.00e+00 1.00e+00h 1 138 3.1738405e-03 2.85e-06 2.40e-03 -4.1 5.80e-03 -2.0 1.00e+00 1.00e+00h 1 139 3.1841321e-03 2.38e-05 4.34e-03 -4.1 1.72e-02 -2.4 1.00e+00 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 140 3.2109770e-03 3.01e-04 5.88e-02 -4.1 4.95e-02 -2.9 1.00e+00 1.00e+00h 1 141 3.2098542e-03 2.18e-06 1.93e-04 -4.1 2.44e-03 -1.6 1.00e+00 1.00e+00h 1 142 2.8238923e-03 1.04e-04 4.64e-02 -6.1 1.84e-02 -2.1 9.58e-01 1.00e+00h 1 143 3.2494631e-03 4.00e-05 6.45e-02 -4.0 7.36e-03 -1.6 1.00e+00 1.00e+00h 1 144 3.1697313e-03 1.73e-05 2.07e-03 -4.1 7.75e-03 -2.1 1.00e+00 1.00e+00h 1 145 3.1906787e-03 5.30e-05 1.85e-02 -4.1 2.29e-02 -2.6 1.00e+00 1.00e+00h 1 146 3.1879642e-03 1.11e-07 1.58e-03 -4.1 1.07e-03 -1.3 1.00e+00 1.00e+00h 1 147 3.1921336e-03 9.57e-07 2.88e-03 -4.1 3.21e-03 -1.7 1.00e+00 1.00e+00h 1 148 3.1991654e-03 9.49e-06 7.51e-05 -4.1 9.56e-03 -2.2 1.00e+00 1.00e+00h 1 149 2.6889508e-03 1.40e-03 1.09e-01 -6.1 7.07e-02 -2.7 8.46e-01 1.00e+00h 1 iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls 150 3.0972630e-03 2.87e-05 2.93e-03 -4.0 7.32e-03 -1.4 1.00e+00 1.00e+00h 1 Number of Iterations....: 150 (scaled) (unscaled) Objective...............: 3.0972630329299377e-03 3.0972630329299377e-03 Dual infeasibility......: 2.9337899633207543e-03 2.9337899633207543e-03 Constraint violation....: 2.8718371594183001e-05 2.8718371594183001e-05 Variable bound violation: 0.0000000000000000e+00 0.0000000000000000e+00 Complementarity.........: 1.0389925318543357e-04 1.0389925318543357e-04 Overall NLP error.......: 2.9337899633207543e-03 2.9337899633207543e-03 Number of objective function evaluations = 159 Number of objective gradient evaluations = 151 Number of equality constraint evaluations = 159 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 = 90.268 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 at ensemble_trajectory.jl:109 [inlined] └ @ Core ~/.julia/packages/PiccoloQuantumObjects/vwXsF/src/quantum_trajectories/ensemble_trajectory.jl:109 applying timesteps_all_equal constraint: Δt 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 initializing optimizer... 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 This is Ipopt version 3.14.19, running with linear solver MUMPS 5.8.2. 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 = 11.603 EXIT: Maximum Number of Iterations Exceeded. initializing optimizer... 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 This is Ipopt version 3.14.19, running with linear solver MUMPS 5.8.2. 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 = 0.002 EXIT: Invalid number in NLP function or derivative detected. Test Summary: | Pass Broken Total Time Package | 221 4 225 31m30.4s Testing QuantumCollocation tests passed Testing completed after 1943.41s PkgEval succeeded after 2078.35s