Test failure in Jacobian test

[ViralBShah]

With Julia 1.3.0 on 64-bit windows and mac. Sundials.jl v3.8.1.

(v1.3) pkg> test Sundials
   Testing Sundials
 Resolving package versions...
 Installed LaTeXStrings ──────────── v1.0.3
 Installed MatrixFactorizations ──── v0.2.1
 Installed NameResolution ────────── v0.1.3
 Installed Quadmath ──────────────── v0.5.2
 Installed DiffEqProblemLibrary ──── v4.6.0
 Installed GenericSchur ──────────── v0.3.0
 Installed ForwardDiff ───────────── v0.10.7
 Installed DiffEqOperators ───────── v4.6.0
 Installed GenericSVD ────────────── v0.2.2
 Installed FillArrays ────────────── v0.8.2
 Installed Tables ────────────────── v0.2.11
 Installed JuliaVariables ────────── v0.2.0
 Installed DiffRules ─────────────── v0.1.0
 Installed FixedPolynomials ──────── v0.4.0
 Installed BlockArrays ───────────── v0.10.2
 Installed ModelingToolkit ───────── v1.0.3
 Installed ArrayLayouts ──────────── v0.1.5
 Installed Formatting ────────────── v0.4.1
 Installed TimerOutputs ──────────── v0.5.3
 Installed ProgressMeter ─────────── v1.2.0
 Installed PoissonRandom ─────────── v0.4.0
 Installed PrettyPrint ───────────── v0.1.0
 Installed MixedSubdivisions ─────── v1.0.0
 Installed NNlib ─────────────────── v0.6.0
 Installed Crayons ───────────────── v4.0.1
 Installed DataValueInterfaces ───── v1.0.0
 Installed DiffResults ───────────── v0.0.4
 Installed DiffEqBiological ──────── v4.1.0
 Installed LazyArrays ────────────── v0.14.10
 Installed Polynomials ───────────── v0.6.0
 Installed BlockBandedMatrices ───── v0.6.0
 Installed HomotopyContinuation ──── v1.3.1
 Installed DynamicPolynomials ────── v0.3.3
 Installed NaNMath ───────────────── v0.3.3
 Installed StaticPolynomials ─────── v1.3.3
 Installed DiffEqJump ────────────── v6.3.0
 Installed CanonicalTraits ───────── v0.1.0
 Installed MultivariatePolynomials ─ v0.3.2
 Installed Latexify ──────────────── v0.12.3
 Installed BandedMatrices ────────── v0.14.2
 Installed ElasticArrays ─────────── v1.0.0
 Installed SpecialFunctions ──────── v0.9.0
 Installed ProjectiveVectors ─────── v1.1.2
 Installed CommonSubexpressions ──── v0.2.0
 Installed GeneralizedGenerated ──── v0.2.0
 Installed PrettyTables ──────────── v0.6.0
 Installed DoubleFloats ──────────── v1.0.3
 Installed RandomNumbers ─────────── v1.3.0
 Installed SymEngine ─────────────── v0.7.0
 Installed MLStyle ───────────────── v0.3.1
  Building SymEngine → `C:\Users\viral\.julia\packages\SymEngine\zSUGO\deps\build.log`
    Status `C:\Users\viral\AppData\Local\Temp\jl_IrQX6b\Manifest.toml`
  [4fba245c] ArrayInterface v2.1.0
  [4c555306] ArrayLayouts v0.1.5
  [aae01518] BandedMatrices v0.14.2
  [b99e7846] BinaryProvider v0.5.8
  [8e7c35d0] BlockArrays v0.10.2
  [ffab5731] BlockBandedMatrices v0.6.0
  [a603d957] CanonicalTraits v0.1.0
  [bbf7d656] CommonSubexpressions v0.2.0
  [34da2185] Compat v2.2.0
  [a8cc5b0e] Crayons v4.0.1
  [9a962f9c] DataAPI v1.1.0
  [864edb3b] DataStructures v0.17.6
  [e2d170a0] DataValueInterfaces v1.0.0
  [2b5f629d] DiffEqBase v6.7.0
  [eb300fae] DiffEqBiological v4.1.0
  [01453d9d] DiffEqDiffTools v1.6.0
  [c894b116] DiffEqJump v6.3.0
  [9fdde737] DiffEqOperators v4.6.0
  [a077e3f3] DiffEqProblemLibrary v4.6.0
  [163ba53b] DiffResults v0.0.4
  [b552c78f] DiffRules v0.1.0
  [ffbed154] DocStringExtensions v0.8.1
  [497a8b3b] DoubleFloats v1.0.3
  [7c1d4256] DynamicPolynomials v0.3.3
  [fdbdab4c] ElasticArrays v1.0.0
  [1a297f60] FillArrays v0.8.2
  [3dd14ad9] FixedPolynomials v0.4.0
  [59287772] Formatting v0.4.1
  [f6369f11] ForwardDiff v0.10.7
  [069b7b12] FunctionWrappers v1.0.0
  [6b9d7cbe] GeneralizedGenerated v0.2.0
  [01680d73] GenericSVD v0.2.2
  [c145ed77] GenericSchur v0.3.0
  [f213a82b] HomotopyContinuation v1.3.1
  [42fd0dbc] IterativeSolvers v0.8.1
  [82899510] IteratorInterfaceExtensions v1.0.0
  [b14d175d] JuliaVariables v0.2.0
  [b964fa9f] LaTeXStrings v1.0.3
  [23fbe1c1] Latexify v0.12.3
  [5078a376] LazyArrays v0.14.10
  [d8e11817] MLStyle v0.3.1
  [1914dd2f] MacroTools v0.5.3
  [a3b82374] MatrixFactorizations v0.2.1
  [291d046c] MixedSubdivisions v1.0.0
  [961ee093] ModelingToolkit v1.0.3
  [46d2c3a1] MuladdMacro v0.2.1
  [102ac46a] MultivariatePolynomials v0.3.2
  [872c559c] NNlib v0.6.0
  [77ba4419] NaNMath v0.3.3
  [71a1bf82] NameResolution v0.1.3
  [efe28fd5] OpenSpecFun_jll v0.5.3+1
  [bac558e1] OrderedCollections v1.1.0
  [d96e819e] Parameters v0.12.0
  [e409e4f3] PoissonRandom v0.4.0
  [f27b6e38] Polynomials v0.6.0
  [8162dcfd] PrettyPrint v0.1.0
  [08abe8d2] PrettyTables v0.6.0
  [92933f4c] ProgressMeter v1.2.0
  [01f381cc] ProjectiveVectors v1.1.2
  [be4d8f0f] Quadmath v0.5.2
  [e6cf234a] RandomNumbers v1.3.0
  [3cdcf5f2] RecipesBase v0.7.0
  [731186ca] RecursiveArrayTools v1.2.0
  [f2c3362d] RecursiveFactorization v0.1.0
  [189a3867] Reexport v0.2.0
  [ae029012] Requires v0.5.2
  [f2b01f46] Roots v0.8.4
  [276daf66] SpecialFunctions v0.9.0
  [90137ffa] StaticArrays v0.12.1
  [62e018b1] StaticPolynomials v1.3.3
  [c3572dad] Sundials v3.8.1
  [123dc426] SymEngine v0.7.0
  [3783bdb8] TableTraits v1.0.0
  [bd369af6] Tables v0.2.11
  [a759f4b9] TimerOutputs v0.5.3
  [a2a6695c] TreeViews v0.3.0
  [2a0f44e3] Base64  [`@stdlib/Base64`]
  [ade2ca70] Dates  [`@stdlib/Dates`]
  [8bb1440f] DelimitedFiles  [`@stdlib/DelimitedFiles`]
  [8ba89e20] Distributed  [`@stdlib/Distributed`]
  [b77e0a4c] InteractiveUtils  [`@stdlib/InteractiveUtils`]
  [76f85450] LibGit2  [`@stdlib/LibGit2`]
  [8f399da3] Libdl  [`@stdlib/Libdl`]
  [37e2e46d] LinearAlgebra  [`@stdlib/LinearAlgebra`]
  [56ddb016] Logging  [`@stdlib/Logging`]
  [d6f4376e] Markdown  [`@stdlib/Markdown`]
  [a63ad114] Mmap  [`@stdlib/Mmap`]
  [44cfe95a] Pkg  [`@stdlib/Pkg`]
  [de0858da] Printf  [`@stdlib/Printf`]
  [3fa0cd96] REPL  [`@stdlib/REPL`]
  [9a3f8284] Random  [`@stdlib/Random`]
  [ea8e919c] SHA  [`@stdlib/SHA`]
  [9e88b42a] Serialization  [`@stdlib/Serialization`]
  [1a1011a3] SharedArrays  [`@stdlib/SharedArrays`]
  [6462fe0b] Sockets  [`@stdlib/Sockets`]
  [2f01184e] SparseArrays  [`@stdlib/SparseArrays`]
  [10745b16] Statistics  [`@stdlib/Statistics`]
  [4607b0f0] SuiteSparse  [`@stdlib/SuiteSparse`]
  [8dfed614] Test  [`@stdlib/Test`]
  [cf7118a7] UUIDs  [`@stdlib/UUIDs`]
  [4ec0a83e] Unicode  [`@stdlib/Unicode`]
T=0.4, Y=[0.9899652943580929, 3.470564190666435e-5, 0.01]
roots=Int32[0, 1]
T=0.4, Y=[0.9851640906722351, 3.386241699325715e-5, 0.014802046910771114]
T=4.0, Y=[0.905509670307843, 2.2403378846050117e-5, 0.09446792631331014]
T=40.0, Y=[0.7158015797469233, 9.185043328931088e-6, 0.28418923520974726]
T=400.0, Y=[0.450520872333455, 3.2228257877077353e-6, 0.5494759048407565]
T=4000.0, Y=[0.18322169164077126, 8.943516144954548e-7, 0.8167774140076135]
T=40000.0, Y=[0.03898091132304715, 1.6216688934655267e-7, 0.9610189265100634]
T=400000.0, Y=[0.004936971391576565, 1.9844502466172285e-8, 0.9950630087639206]
T=4.0e6, Y=[0.000517010320204611, 2.0690977989759934e-9, 0.9994829876106978]
T=4.0e7, Y=[9.999999999999992e-5, 4.00039528747344e-10, 0.9998999995999602]
roots=Int32[-1, 0]
T=4.0e7, Y=[5.204926760391501e-5, 2.0820776359257418e-10, 0.9999479505241882]
T=4.0e8, Y=[5.1849463049478235e-6, 2.0739891508302815e-11, 0.9999948150329551]
T=4.0e9, Y=[5.246211769412839e-7, 2.0984858346733178e-12, 0.9999994753767246]
T=4.0e10, Y=[6.043455470531721e-8, 2.41738235380483e-13, 0.9999999395652039]
Test Summary: |
CVODE         | No tests
T=0.4, Y=[0.9851721138377025, 3.386406120817996e-5, 0.014794022101089307]
T=4.0, Y=[0.9055186785328794, 2.2404781805334787e-5, 0.09445891668531428]
T=40.0, Y=[0.715827068739702, 9.185512472082487e-6, 0.28416374574782616]
T=400.0, Y=[0.45051866843284366, 3.222907159344263e-6, 0.5494781086599967]
T=4000.0, Y=[0.18320225667994575, 8.94230214472978e-7, 0.8167968490898406]
T=40000.0, Y=[0.03898337827880326, 1.621799878073902e-7, 0.96101645954121]
T=400000.0, Y=[0.004938270331418616, 1.985351550862375e-8, 0.9950617098150669]
T=4.0e6, Y=[0.0005168101324260045, 2.068434777531092e-9, 0.9994831877991401]
T=4.0e7, Y=[5.202985953290682e-5, 2.0832525035625152e-10, 0.9999479699321431]
T=4.0e8, Y=[5.207659630111956e-6, 2.0841931899519505e-11, 0.9999947923195291]
T=4.0e9, Y=[5.208321136050987e-7, 2.0828928173656967e-12, 0.9999994791658047]
T=4.0e10, Y=[5.208358902243623e-8, 1.977013868181389e-13, 0.999999947916214]
Test Summary: |
ARK           | No tests
T=0.4, Y=[0.9899652943718806, 3.4705628119662186e-5, 0.010000000000000002]
roots=Int32[0, 1]
T=0.4, Y=[0.9851707623122912, 3.386371656828505e-5, 0.01479537397114059]
T=4.0, Y=[0.905504728443, 2.24030741323827e-5, 0.09447286848286832]
T=40.0, Y=[0.71584189841088, 9.186119502141643e-6, 0.28414891546961785]
T=400.0, Y=[0.450537891284184, 3.2231511100140762e-6, 0.5494588855647063]
T=4000.0, Y=[0.1831996994939642, 8.942143972145097e-7, 0.8167994062916398]
T=40000.0, Y=[0.03898368270471694, 1.6217829385564226e-7, 0.9610161551169892]
T=400000.0, Y=[0.004938675705721228, 1.985155431045523e-8, 0.995061304442724]
T=4.0e6, Y=[0.00051688824532254, 2.0686115064950218e-9, 0.9994831096860671]
T=4.0e7, Y=[9.999999999999957e-5, 4.0003952202291907e-10, 0.9998999995999608]
roots=Int32[-1, 0]
T=4.0e7, Y=[5.204772673760556e-5, 2.0820161918090773e-10, 0.9999479520650573]
T=4.0e8, Y=[5.205357160339493e-6, 2.0821535764088784e-11, 0.9999947946220178]
T=4.0e9, Y=[5.190800205273695e-7, 2.076321146499096e-12, 0.9999994809179031]
T=4.0e10, Y=[5.11777367290652e-8, 2.047109574722673e-13, 0.9999999488220584]
Test Summary: |
IDA           | No tests
Solution: [0.7861540663366059, 0.6180344478352109]
Residual: [4.794728554635697e-6, -3.7681823696900096e-6]
Test Summary: | Pass  Total
Kinsol        |    2      2
Test Summary: | Pass  Total
Handle Tests  |   13     13
┌ Warning: Replacing docs for `DiffEqOperators.PeriodicBC :: Union{}` in module `DiffEqOperators`
└ @ Base.Docs docs\Docs.jl:223
┌ Warning: Replacing docs for `DiffEqOperators.RobinBC :: Union{}` in module `DiffEqOperators`
└ @ Base.Docs docs\Docs.jl:223

[IDAS ERROR]  IDASolve
  At t = 0.0395936 and h = 6.30116e-016, the corrector convergence failed repeatedly or with |h| = hmin.


[IDAS ERROR]  IDASolve
  At t = 0.0317944 and h = 3.07674e-019, the corrector convergence failed repeatedly or with |h| = hmin.

With consistent initial conditions:
With inconsistent initial conditions:
Jacobians: Error During Test at C:\Users\viral\.julia\packages\Sundials\MllUG\test\runtests.jl:29
  Got exception outside of a @test
  LoadError: ArgumentError: 1 == colptr[1] != 1
  Stacktrace:
   [1] sparse_check(::Int64, ::Array{Int64,1}, ::Array{Int64,1}, ::Array{Float64,1}) at D:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.3\SparseArrays\src\sparsematrix.jl:58
   [2] SparseMatrixCSC(::Int64, ::Int64, ::Array{Int64,1}, ::Array{Int64,1}, ::Array{Float64,1}) at D:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.3\SparseArrays\src\sparsematrix.jl:39
   [3] convert(::Type{SparseMatrixCSC}, ::Ptr{Sundials._generic_SUNMatrix}) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\types_and_consts_additions.jl:49
   [4] cvodejac(::Float64, ::Ptr{Sundials._generic_N_Vector}, ::Ptr{Sundials._generic_N_Vector}, ::Ptr{Sundials._generic_SUNMatrix}, ::Sundials.FunJac{ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Nothing,typeof(Lotka_jac),DiffEqBase.NullParameters,Nothing,SparseMatrixCSC{Float64,Int64},Array{Float64,1},Nothing,Nothing,Nothing}, ::Ptr{Sundials._generic_N_Vector}, ::Ptr{Sundials._generic_N_Vector}, ::Ptr{Sundials._generic_N_Vector}) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\common_interface\function_types.jl:66
   [5] __CVode at C:\Users\viral\.julia\packages\Sundials\MllUG\src\wrapped_api\cvodes.jl:312 [inlined]
   [6] CVode(::Sundials.Handle{Sundials.CVODEMem}, ::Float64, ::Array{Float64,1}, ::Array{Float64,1}, ::Int32) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\wrapped_api\cvodes.jl:317
   [7] solver_step(::Sundials.CVODEIntegrator{Array{Float64,1},DiffEqBase.NullParameters,Sundials.Handle{Sundials.CVODEMem},ODESolution{Float64,2,Array{Array{Float64,1},1},Nothing,Nothing,Array{Float64,1},Nothing,ODEProblem{Array{Float64,1},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}},DiffEqBase.StandardODEProblem},CVODE_BDF{:Newton,:KLU,Nothing,Nothing},DiffEqBase.HermiteInterpolation{Array{Float64,1},Array{Array{Float64,1},1},Array{Array{Float64,1},1}},DiffEqBase.DEStats},CVODE_BDF{:Newton,:KLU,Nothing,Nothing},ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Sundials.FunJac{ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Nothing,typeof(Lotka_jac),DiffEqBase.NullParameters,Nothing,SparseMatrixCSC{Float64,Int64},Array{Float64,1},Nothing,Nothing,Nothing},Ptr{Nothing},Sundials.DEOptions{DataStructures.BinaryHeap{Float64,DataStructures.LessThan},DataStructures.BinaryHeap{Float64,DataStructures.LessThan},CallbackSet{Tuple{},Tuple{}},Float64,Float64,typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE)},Array{Float64,1},Tuple{Int64},Array{Float64,1},Sundials.LinSolHandle{Sundials.KLU},Sundials.MatrixHandle{Sundials.SparseMatrix},Nothing}, ::Float64) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\common_interface\solve.jl:1091
   [8] solve!(::Sundials.CVODEIntegrator{Array{Float64,1},DiffEqBase.NullParameters,Sundials.Handle{Sundials.CVODEMem},ODESolution{Float64,2,Array{Array{Float64,1},1},Nothing,Nothing,Array{Float64,1},Nothing,ODEProblem{Array{Float64,1},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}},DiffEqBase.StandardODEProblem},CVODE_BDF{:Newton,:KLU,Nothing,Nothing},DiffEqBase.HermiteInterpolation{Array{Float64,1},Array{Array{Float64,1},1},Array{Array{Float64,1},1}},DiffEqBase.DEStats},CVODE_BDF{:Newton,:KLU,Nothing,Nothing},ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Sundials.FunJac{ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Nothing,typeof(Lotka_jac),DiffEqBase.NullParameters,Nothing,SparseMatrixCSC{Float64,Int64},Array{Float64,1},Nothing,Nothing,Nothing},Ptr{Nothing},Sundials.DEOptions{DataStructures.BinaryHeap{Float64,DataStructures.LessThan},DataStructures.BinaryHeap{Float64,DataStructures.LessThan},CallbackSet{Tuple{},Tuple{}},Float64,Float64,typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE)},Array{Float64,1},Tuple{Int64},Array{Float64,1},Sundials.LinSolHandle{Sundials.KLU},Sundials.MatrixHandle{Sundials.SparseMatrix},Nothing}) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\common_interface\solve.jl:1145
   [9] #__solve#22(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}, ::typeof(DiffEqBase.__solve), ::ODEProblem{Array{Float64,1},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}},DiffEqBase.StandardODEProblem}, ::CVODE_BDF{:Newton,:KLU,Nothing,Nothing}, ::Array{Any,1}, ::Array{Any,1}, ::Array{Any,1}, ::Type{Val{true}}) at C:\Users\viral\.julia\packages\Sundials\MllUG\src\common_interface\solve.jl:12
   [10] __solve at C:\Users\viral\.julia\packages\Sundials\MllUG\src\common_interface\solve.jl:10 [inlined] (repeats 5 times)
   [11] #solve_call#442(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}}, ::typeof(DiffEqBase.solve_call), ::ODEProblem{Array{Float64,1},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}},DiffEqBase.StandardODEProblem}, ::CVODE_BDF{:Newton,:KLU,Nothing,Nothing}) at C:\Users\viral\.julia\packages\DiffEqBase\E16PL\src\solve.jl:40
   [12] solve_call at C:\Users\viral\.julia\packages\DiffEqBase\E16PL\src\solve.jl:37 [inlined]
   [13] #solve#443 at C:\Users\viral\.julia\packages\DiffEqBase\E16PL\src\solve.jl:57 [inlined]
   [14] solve(::ODEProblem{Array{Float64,1},Tuple{Float64,Float64},true,DiffEqBase.NullParameters,ODEFunction{true,typeof(Lotka),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,typeof(Lotka_jac),SparseMatrixCSC{Float64,Int64},Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{(),Tuple{}}},DiffEqBase.StandardODEProblem}, ::CVODE_BDF{:Newton,:KLU,Nothing,Nothing}) at C:\Users\viral\.julia\packages\DiffEqBase\E16PL\src\solve.jl:45
   [15] top-level scope at C:\Users\viral\.julia\packages\Sundials\MllUG\test\common_interface\jacobians.jl:31
   [16] include at .\boot.jl:328 [inlined]
   [17] include_relative(::Module, ::String) at .\loading.jl:1105
   [18] include(::Module, ::String) at .\Base.jl:31
   [19] include(::String) at .\client.jl:424
   [20] top-level scope at C:\Users\viral\.julia\packages\Sundials\MllUG\test\runtests.jl:29
   [21] top-level scope at D:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.3\Test\src\Test.jl:1107
   [22] top-level scope at C:\Users\viral\.julia\packages\Sundials\MllUG\test\runtests.jl:29
   [23] top-level scope at D:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.3\Test\src\Test.jl:1107
   [24] top-level scope at C:\Users\viral\.julia\packages\Sundials\MllUG\test\runtests.jl:26
   [25] include at .\boot.jl:328 [inlined]
   [26] include_relative(::Module, ::String) at .\loading.jl:1105
   [27] include(::Module, ::String) at .\Base.jl:31
   [28] include(::String) at .\client.jl:424
   [29] top-level scope at none:6
   [30] eval(::Module, ::Any) at .\boot.jl:330
   [31] exec_options(::Base.JLOptions) at .\client.jl:263
   [32] _start() at .\client.jl:460
  in expression starting at C:\Users\viral\.julia\packages\Sundials\MllUG\test\common_interface\jacobians.jl:31

i.t = 0.04436589458068945
i.t = 0.11976370878673312
i.t = 0.19516152299277678
i.t = 0.27055933719882047
i.t = 0.34595715140486416
i.t = 0.42135496561090785
i.t = 0.5878600458168253
i.t = 0.7543651260227427
i.t = 0.9208702062286602
i.t = 1.0
Test error handling

[CVODES WARNING]  CVode
  Internal t = 0 and h = 0 are such that t + h = t on the next step. The solver will continue anyway.


[CVODES ERROR]  CVode
  At t = 0 and h = 0, the corrector convergence test failed repeatedly or with |h| = hmin.


[IDAS ERROR]  IDACalcIC
  Newton/Linesearch algorithm failed to converge.

Test Summary:    | Pass  Error  Total
Common Interface |   40      1     41
  CVODE          |   23            23
  ARKODE         |    3             3
  IDA            |    7             7
  Jacobians      |    1      1      2
  Callbacks      |    4             4
  Iterator       |              No tests
  Errors         |    1             1
  Mass Matrix    |    1             1
ERROR: LoadError: Some tests did not pass: 40 passed, 0 failed, 1 errored, 0 broken.
in expression starting at C:\Users\viral\.julia\packages\Sundials\MllUG\test\runtests.jl:25
ERROR: Package Sundials errored during testing

(v1.3) pkg> st
    Status `C:\Users\viral\.julia\environments\v1.3\Project.toml`
  [c3572dad] Sundials v3.8.1

julia> versioninfo()
Julia Version 1.3.0
Commit 46ce4d7933 (2019-11-26 06:09 UTC)
Platform Info:
  OS: Windows (x86_64-w64-mingw32)
  CPU: Intel(R) Core(TM) i5-3340M CPU @ 2.70GHz
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-6.0.1 (ORCJIT, ivybridge)

julia>                                                                                                                  

[ChrisRackauckas]

It looks like the v1.3 change found a bug in the sparse Jacobian handling. I don't know why it never errored before, but

(J.rowval, J.colptr, J.nzval, jac_prototype.nzval) = ([0, 0, 0, 0], [0, 0, 0], [0.0, 0.0, 0.0, 0.0], [0.025161239785942757, 0.9748387602140572, -1.0, -2.0])
(J.rowval, J.colptr, J.nzval, jac_prototype.nzval) = ([0, 0, 0, 0], [0, 0, 0], [0.0, 0.0, 0.0, 0.0], [-1.6967896443241721, 2.696789644324172, -1.4739916695950448, -1.5260083304049552])
(J.rowval, J.colptr, J.nzval, jac_prototype.nzval) = ([0, 0, 0, 0], [0, 0, 0], [0.0, 0.0, 0.0, 0.0], [0.3008406498733516, 0.6991593501266484, -1.0319599755403681, -1.9680400244596319])

on previous versions, so it wasn't getting wrapped properly. On v1.3 Julia catches this and errors. If that is avoided, the nzval changes sizes from 4 to 2 after 1 step, so something really weird is going on. Sounds like it is grabbing the wrong pieces of memory or something... but anyways, it's a miracle that test was passing before. So I'm just setting it to broken for now.

[ViralBShah]

IIRC, we started zeroing out the sparse matrix data structures, which forces these sort of bugs.

[ChrisRackauckas]

oh wait...

https://github.com/JuliaDiffEq/Sundials.jl/blob/master/src/types_and_consts_additions.jl#L38-L50

so I grab the stuff from Sundials and call SparseMatrixCSC, but then the constructor zeros everything now?

[ViralBShah]

No, I don't think that constructor will zero anything. The zeroing is only in the case of creating empty structure (which would in the past have remained uninitialized). We are probably checking the consistency of the sparse data structure more than we used to though.