Add test for noise kwarg propagation

Added test case to verify that user-provided noise keyword argument
is properly propagated when creating SDEProblem from SDESystem.

The test verifies:
- User-provided noise is actually used (noise.curW is modified)
- Using the same noise via NoiseWrapper gives deterministic results
- Different seeds produce different results (default behavior still works)

🤖 Generated with [Claude Code](https://claude.com/claude-code)

Co-Authored-By: Claude <noreply@anthropic.com>

Author: ChrisRackauckas

test/sdesystem.jl

@@ -1,5 +1,6 @@
 using ModelingToolkit, StaticArrays, LinearAlgebra
 using StochasticDiffEq, OrdinaryDiffEq, SparseArrays
+using DiffEqNoiseProcess: NoiseWrapper
 using Random, Test
 using Setfield
 using Statistics
@@ -953,3 +954,47 @@ end
     @test ModelingToolkit.isbrownian(p)
     @test ModelingToolkit.isbrownian(q)
 end
+
+@testset "noise kwarg propagation (issue #3664)" begin
+    @parameters σ ρ β
+    @variables x(tt) y(tt) z(tt)
+
+    u0 = [1.0, 0.0, 0.0]
+    T = (0.0, 5.0)
+
+    eqs = [D(x) ~ σ * (y - x),
+        D(y) ~ x * (ρ - z) - y,
+        D(z) ~ x * y - β * z]
+    noiseeqs = [3.0,
+        3.0,
+        3.0]
+    @mtkbuild sde_lorentz = SDESystem(eqs, noiseeqs, tt, [x, y, z], [σ, ρ, β])
+    parammap = [σ, ρ, β] .=> [10, 28.0, 8 / 3]
+
+    # Test that user-provided noise is respected
+    Random.seed!(1)
+    noise1 = StochasticDiffEq.RealWienerProcess(0.0, 0.0, 0.0; save_everystep = true)
+    u0_dict = Dict(unknowns(sde_lorentz) .=> u0)
+    prob1 = SDEProblem(sde_lorentz, merge(u0_dict, Dict(parammap)), T; noise = noise1)
+    sol1 = solve(prob1, SRIW1())
+
+    # Verify noise was actually used (curW should be modified)
+    @test noise1.curW != 0.0
+
+    # Test that using the same noise via NoiseWrapper gives deterministic results
+    noise2 = NoiseWrapper(noise1)
+    prob2 = SDEProblem(sde_lorentz, merge(u0_dict, Dict(parammap)), T; noise = noise2)
+    sol2 = solve(prob2, SRIW1())
+
+    # Same noise should give same results
+    @test sol1.u[end] ≈ sol2.u[end]
+
+    # Test that without providing noise, different results are obtained
+    Random.seed!(1)
+    prob3 = SDEProblem(sde_lorentz, merge(u0_dict, Dict(parammap)), T)
+    Random.seed!(2)
+    prob4 = SDEProblem(sde_lorentz, merge(u0_dict, Dict(parammap)), T)
+    sol3 = solve(prob3, SRIW1(), seed = 1)
+    sol4 = solve(prob4, SRIW1(), seed = 2)
+    @test !(sol3.u[end] ≈ sol4.u[end])
+end