From fbba7395130ba0cf609a6b53766e597bb7c90fe2 Mon Sep 17 00:00:00 2001 From: Pepijn de Vos Date: Thu, 16 Nov 2023 17:28:17 +0100 Subject: [PATCH] apparently type inference doesn't like the broadcast macro thing --- src/dense/generic_dense.jl | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/dense/generic_dense.jl b/src/dense/generic_dense.jl index cdef8d227f..0299f689b7 100644 --- a/src/dense/generic_dense.jl +++ b/src/dense/generic_dense.jl @@ -641,8 +641,8 @@ Herimte Interpolation, chosen if no other dispatch for ode_interpolant @muladd function hermite_interpolant(Θ, dt, y₀, y₁, k, ::Type{Val{false}}, idxs::Nothing, T::Type{Val{0}}, dv=1) # Default interpolant is Hermite #@.. broadcast=false (1-Θ)*y₀+Θ*y₁+Θ*(Θ-1)*((1-2Θ)*(y₁-y₀)+(Θ-1)*dt*k[1] + Θ*dt*k[2]) - @inbounds @.. (1 - Θ) * y₀ + Θ * y₁ + - dv * Θ * (Θ - 1) * ((1 - 2Θ) * (y₁ - y₀) + (Θ - 1) * dt * k[1] + Θ * dt * k[2]) + @inbounds (1 - Θ) * y₀ + Θ * y₁ + + dv .* (Θ * (Θ - 1) * ((1 - 2Θ) * (y₁ - y₀) + (Θ - 1) * dt * k[1] + Θ * dt * k[2])) end @muladd function hermite_interpolant(Θ, dt, y₀, y₁, k, ::Type{Val{true}}, idxs::Nothing, @@ -664,10 +664,10 @@ end @muladd function hermite_interpolant(Θ, dt, y₀, y₁, k, cache, idxs, T::Type{Val{0}}, dv=1) # Default interpolant is Hermite # return @.. broadcast=false (1-Θ)*y₀[idxs]+Θ*y₁[idxs]+Θ*(Θ-1)*((1-2Θ)*(y₁[idxs]-y₀[idxs])+(Θ-1)*dt*k[1][idxs] + Θ*dt*k[2][idxs]) - return @.. (1 - Θ) * y₀[idxs] + Θ * y₁[idxs] + - dv * Θ * (Θ - 1) * + return (1 - Θ) * y₀[idxs] + Θ * y₁[idxs] + + dv .* (Θ * (Θ - 1) * ((1 - 2Θ) * (y₁[idxs] - y₀[idxs]) + (Θ - 1) * dt * k[1][idxs] + - Θ * dt * k[2][idxs]) + Θ * dt * k[2][idxs])) end @muladd function hermite_interpolant!(out, Θ, dt, y₀, y₁, k, idxs::Nothing, T::Type{Val{0}}, dv=1) # Default interpolant is Hermite