BugFix: remove unnecessary broadcast in jacobiany

Author: hwpang

src/Domain.jl

@@ -1627,7 +1627,7 @@ end
         if isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
             end
             @views @inbounds @fastmath jac[domain.indexes[1]:domain.indexes[2],domain.indexes[3]] .-= flow.*ns.*(-C)
         end
@@ -1665,7 +1665,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
         elseif isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
                 @inbounds @fastmath dCvavedni = cpdivR[i]*R/N
                 @fastmath dTdt = (P*V/N*flow)/(N*Cvave)
                 @fastmath ddnidTdt = -dTdt*(dCvavedni/Cvave)
@@ -1720,7 +1720,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
         elseif isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
             end
         end
     end
@@ -1747,7 +1747,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
         if isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
             end
         end
     end
@@ -1788,7 +1788,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
             flow = inter.F(t)
             @fastmath dTdt = (P*V/N*flow)/(N*Cvave)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
                 @inbounds @fastmath dCvavedni = cpdivR[i]*R/N
                 @fastmath ddnidTdt = -dTdt*(dCvavedni/Cvave)
                 @inbounds jac[domain.indexes[3],i] -= ddnidTdt
@@ -1843,7 +1843,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
         elseif isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
             end
             @views @inbounds @fastmath jac[domain.indexes[1]:domain.indexes[2],domain.indexes[4]] .-= -flow.*ns./N/V
         end
@@ -1866,7 +1866,7 @@ function jacobiany!(jac::Q,y::U,p::W,t::Z,domain::D,interfaces::Q3,colorvec::Q2=
         if isa(inter,Outlet) && domain == inter.domain
             flow = inter.F(t)
             @simd for i in domain.indexes[1]:domain.indexes[2]
-                @inbounds @fastmath jac[i,i] .-= flow/N
+                @inbounds @fastmath jac[i,i] -= flow/N
             end
         end
     end