{"action":"created","issue":{"url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/291","repository_url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl","labels_url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/291/labels{/name}","comments_url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/291/comments","events_url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/291/events","html_url":"https://github.com/JuliaDiffEq/DifferentialEquations.jl/issues/291","id":324874410,"node_id":"MDU6SXNzdWUzMjQ4NzQ0MTA=","number":291,"title":"Jacobians for arrays of static arrays","user":{"login":"ChrisRackauckas","id":1814174,"node_id":"MDQ6VXNlcjE4MTQxNzQ=","avatar_url":"https://avatars2.githubusercontent.com/u/1814174?v=4","gravatar_id":"","url":"https://api.github.com/users/ChrisRackauckas","html_url":"https://github.com/ChrisRackauckas","followers_url":"https://api.github.com/users/ChrisRackauckas/followers","following_url":"https://api.github.com/users/ChrisRackauckas/following{/other_user}","gists_url":"https://api.github.com/users/ChrisRackauckas/gists{/gist_id}","starred_url":"https://api.github.com/users/ChrisRackauckas/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/ChrisRackauckas/subscriptions","organizations_url":"https://api.github.com/users/ChrisRackauckas/orgs","repos_url":"https://api.github.com/users/ChrisRackauckas/repos","events_url":"https://api.github.com/users/ChrisRackauckas/events{/privacy}","received_events_url":"https://api.github.com/users/ChrisRackauckas/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2018-05-21T10:58:16Z","updated_at":"2018-07-04T10:02:43Z","closed_at":null,"author_association":"MEMBER","body":"```julia\r\nusing StaticArrays, OrdinaryDiffEq\r\n\r\n@inline @inbounds function loop(u, p, t)\r\n    σ = p[1]; ρ = p[2]; β = p[3]\r\n    du1 = σ*(u[2]-u[1])\r\n    du2 = u[1]*(ρ-u[3]) - u[2]\r\n    du3 = u[1]*u[2] - β*u[3]\r\n    return SVector{3}(du1, du2, du3)\r\nend\r\n\r\np = [10, 28, 8/3]\r\n\r\nstate = [0, 10.0, 0]\r\n\r\nstates = [state, state + 1e-9rand(3)]\r\n\r\nst = [SVector{3}(s) for s in states]\r\nL = length(st)\r\n# The following may be inneficient\r\nparalleleom = (du, u, p, t) -> begin\r\n    for i in 1:L\r\n        @inbounds du[i] = loop(u[i], p, t)\r\n    end\r\nend\r\n\r\n\r\npprob = ODEProblem(paralleleom, st, (0.0, Inf), p)\r\n\r\ninit(pprob, Vern9(); save_everystep = false)\r\ninit(pprob, Tsit5(); save_everystep = false)\r\ninit(pprob, KenCarp5(); save_everystep = false)\r\n```\r\n\r\nThis might just need generic linalg tools to really be handled appropriately."},"comment":{"url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/comments/402431123","html_url":"https://github.com/JuliaDiffEq/DifferentialEquations.jl/issues/291#issuecomment-402431123","issue_url":"https://api.github.com/repos/JuliaDiffEq/DifferentialEquations.jl/issues/291","id":402431123,"node_id":"MDEyOklzc3VlQ29tbWVudDQwMjQzMTEyMw==","user":{"login":"ChrisRackauckas","id":1814174,"node_id":"MDQ6VXNlcjE4MTQxNzQ=","avatar_url":"https://avatars2.githubusercontent.com/u/1814174?v=4","gravatar_id":"","url":"https://api.github.com/users/ChrisRackauckas","html_url":"https://github.com/ChrisRackauckas","followers_url":"https://api.github.com/users/ChrisRackauckas/followers","following_url":"https://api.github.com/users/ChrisRackauckas/following{/other_user}","gists_url":"https://api.github.com/users/ChrisRackauckas/gists{/gist_id}","starred_url":"https://api.github.com/users/ChrisRackauckas/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/ChrisRackauckas/subscriptions","organizations_url":"https://api.github.com/users/ChrisRackauckas/orgs","repos_url":"https://api.github.com/users/ChrisRackauckas/repos","events_url":"https://api.github.com/users/ChrisRackauckas/events{/privacy}","received_events_url":"https://api.github.com/users/ChrisRackauckas/received_events","type":"User","site_admin":false},"created_at":"2018-07-04T10:02:43Z","updated_at":"2018-07-04T10:02:43Z","author_association":"MEMBER","body":"Using reinterpret can be an option to solve this."}}
{"actor":{"display_login":"ChrisRackauckas"}}