Skip to content

Support Dual propagation in out-of-place DAE array residuals #5053

Description

@utkuyilmaz1903

Problem

#4983 array-valued implicit DAE residuals. #5049 makes the out-of-place
ArrayMaker path Dual-safe for NonlinearProblem only.

DAEProblem{false} still uses the default zeros(size) allocator. That
buffer is Float64, so ForwardDiff.Dual on du, u, parameters, or t
is dropped. In-place DAE is unaffected.

Requested change

Give DAEProblem{false} the same residual-eltype promotion already used
on the nonlinear OOP path. Do not change nonlinear residual behavior.

Required coverage

  • ForwardDiff.jacobian / Dual eltype through du, u, tunable p, and t
  • split and non-split parameter structures
  • generated expression = Val{true} path
  • existing 4983 IIP DAE tests remain green

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions