Skip to content

Commit ab2a726

Browse files
update unstable check
1 parent 0ecdbb5 commit ab2a726

File tree

1 file changed

+4
-1
lines changed

1 file changed

+4
-1
lines changed

src/constants.jl

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -2,4 +2,7 @@
22
@inline UNITLESS_ABS2(x) = abs2(x)/(typeof(x)(one(x))*typeof(x)(one(x)))
33
@inline ODE_DEFAULT_NORM(u) = sqrt(sum(UNITLESS_ABS2,u) / length(u))
44
@inline ODE_DEFAULT_PROG_MESSAGE(dt,t,u) = "dt="*string(dt)*"\nt="*string(t)*"\nmax u="*string(maximum(abs.(u)))
5-
@inline ODE_DEFAULT_UNSTABLE_CHECK(dt,t,u) = any(isnan,u)
5+
@inline ODE_DEFAULT_UNSTABLE_CHECK(dt,t,u) = false
6+
(p::typeof(ODE_DEFAULT_UNSTABLE_CHECK))(dt,t,u::AbstractFloat) = isnan(u)
7+
(p::typeof(ODE_DEFAULT_UNSTABLE_CHECK)){T<:AbstractFloat}(dt,t,u::AbstractArray{T}) = any(isnan,u)
8+
(p::typeof(ODE_DEFAULT_UNSTABLE_CHECK))(dt,t,u::ArrayPartition) = any(isnan,chain(u.x...))

0 commit comments

Comments
 (0)