@@ -151,23 +151,23 @@ for (T, S) ∈ ((:Interval, :Interval), (:Interval, :Any), (:Any, :Interval))
151
151
C .*= β
152
152
end
153
153
else
154
- numtype = numtype (CoefType)
154
+ T = numtype (CoefType)
155
155
mC, rC = __mul (A, B)
156
156
if isone (α)
157
157
if iszero (β)
158
- C .= interval .(numtype , mC, rC; format = :midpoint )
158
+ C .= interval .(T , mC, rC; format = :midpoint )
159
159
elseif isone (β)
160
- C .+ = interval .(numtype , mC, rC; format = :midpoint )
160
+ C .+ = interval .(T , mC, rC; format = :midpoint )
161
161
else
162
- C .= interval .(numtype , mC, rC; format = :midpoint ) .+ C .* β
162
+ C .= interval .(T , mC, rC; format = :midpoint ) .+ C .* β
163
163
end
164
164
else
165
165
if iszero (β)
166
- C .= interval .(numtype , mC, rC; format = :midpoint ) .* α
166
+ C .= interval .(T , mC, rC; format = :midpoint ) .* α
167
167
elseif isone (β)
168
- C .+ = interval .(numtype , mC, rC; format = :midpoint ) .* α
168
+ C .+ = interval .(T , mC, rC; format = :midpoint ) .* α
169
169
else
170
- C .= interval .(numtype , mC, rC; format = :midpoint ) .* α .+ C .* β
170
+ C .= interval .(T , mC, rC; format = :midpoint ) .* α .+ C .* β
171
171
end
172
172
end
173
173
end
@@ -188,7 +188,7 @@ for (T, S) ∈ ((:(Complex{<:Interval}), :(Complex{<:Interval})),
188
188
C .*= β
189
189
end
190
190
else
191
- numtype = numtype (CoefType)
191
+ T = numtype (CoefType)
192
192
A_real, A_imag = reim (A)
193
193
B_real, B_imag = reim (B)
194
194
mC_1, rC_1 = __mul (A_real, B_real)
@@ -197,25 +197,25 @@ for (T, S) ∈ ((:(Complex{<:Interval}), :(Complex{<:Interval})),
197
197
mC_4, rC_4 = __mul (A_imag, B_real)
198
198
if isone (α)
199
199
if iszero (β)
200
- C .= complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
201
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint ))
200
+ C .= complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
201
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint ))
202
202
elseif isone (β)
203
- C .+ = complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
204
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint ))
203
+ C .+ = complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
204
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint ))
205
205
else
206
- C .= complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
207
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint )) .+ C .* β
206
+ C .= complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
207
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint )) .+ C .* β
208
208
end
209
209
else
210
210
if iszero (β)
211
- C .= complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
212
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint )) .* α
211
+ C .= complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
212
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint )) .* α
213
213
elseif isone (β)
214
- C .+ = complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
215
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint )) .* α
214
+ C .+ = complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
215
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint )) .* α
216
216
else
217
- C .= complex .(interval .(numtype , mC_1, rC_1; format = :midpoint ) .- interval .(numtype , mC_2, rC_2; format = :midpoint ),
218
- interval .(numtype , mC_3, rC_3; format = :midpoint ) .+ interval .(numtype , mC_4, rC_4; format = :midpoint )) .* α .+ C .* β
217
+ C .= complex .(interval .(T , mC_1, rC_1; format = :midpoint ) .- interval .(T , mC_2, rC_2; format = :midpoint ),
218
+ interval .(T , mC_3, rC_3; format = :midpoint ) .+ interval .(T , mC_4, rC_4; format = :midpoint )) .* α .+ C .* β
219
219
end
220
220
end
221
221
end
@@ -236,25 +236,25 @@ for (T, S) ∈ ((:(Complex{<:Interval}), :Interval), (:(Complex{<:Interval}), :A
236
236
C .*= β
237
237
end
238
238
else
239
- numtype = numtype (CoefType)
239
+ T = numtype (CoefType)
240
240
A_real, A_imag = reim (A)
241
241
mC_real, rC_real = __mul (A_real, B)
242
242
mC_imag, rC_imag = __mul (A_imag, B)
243
243
if isone (α)
244
244
if iszero (β)
245
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint ))
245
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint ))
246
246
elseif isone (β)
247
- C .+ = complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint ))
247
+ C .+ = complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint ))
248
248
else
249
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .+ C .* β
249
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .+ C .* β
250
250
end
251
251
else
252
252
if iszero (β)
253
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α
253
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α
254
254
elseif isone (β)
255
- C .+ = complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α
255
+ C .+ = complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α
256
256
else
257
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α .+ C .* β
257
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α .+ C .* β
258
258
end
259
259
end
260
260
end
@@ -272,25 +272,25 @@ for (T, S) ∈ ((:(Complex{<:Interval}), :Interval), (:(Complex{<:Interval}), :A
272
272
C .*= β
273
273
end
274
274
else
275
- numtype = numtype (CoefType)
275
+ T = numtype (CoefType)
276
276
B_real, B_imag = reim (B)
277
277
mC_real, rC_real = __mul (A, B_real)
278
278
mC_imag, rC_imag = __mul (A, B_imag)
279
279
if isone (α)
280
280
if iszero (β)
281
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint ))
281
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint ))
282
282
elseif isone (β)
283
- C .+ = complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint ))
283
+ C .+ = complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint ))
284
284
else
285
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .+ C .* β
285
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .+ C .* β
286
286
end
287
287
else
288
288
if iszero (β)
289
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α
289
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α
290
290
elseif isone (β)
291
- C .+ = complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α
291
+ C .+ = complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α
292
292
else
293
- C .= complex .(interval .(numtype , mC_real, rC_real; format = :midpoint ), interval .(numtype , mC_imag, rC_imag; format = :midpoint )) .* α .+ C .* β
293
+ C .= complex .(interval .(T , mC_real, rC_real; format = :midpoint ), interval .(T , mC_imag, rC_imag; format = :midpoint )) .* α .+ C .* β
294
294
end
295
295
end
296
296
end
0 commit comments