@@ -1293,13 +1293,13 @@ function operator_clenshaw2D(Jx,Jy,cfs::Vector{Vector{T}},x,y) where T
1293
1293
N= length (cfs)
1294
1294
S = domainspace (x)
1295
1295
Z= ZeroOperator (S,S)
1296
- bk1= Array ( Operator{T},N+ 1 );bk1[:] = Z
1297
- bk2= Array ( Operator{T},N+ 2 );bk2[:] = Z
1296
+ bk1= Array { Operator{T}} (undef ,N+ 1 ); fill! ( bk1 , Z)
1297
+ bk2= Array { Operator{T}} (undef ,N+ 2 ); fill! ( bk2 , Z)
1298
1298
1299
- Abk1x= Array ( Operator{T},N+ 1 );Abk1x[:] = Z
1300
- Abk1y= Array ( Operator{T},N+ 1 );Abk1y[:] = Z
1301
- Abk2x= Array ( Operator{T},N+ 1 );Abk2x[:] = Z
1302
- Abk2y= Array ( Operator{T},N+ 1 );Abk2y[:] = Z
1299
+ Abk1x= Array { Operator{T}} (undef ,N+ 1 ); fill! ( Abk1x , Z)
1300
+ Abk1y= Array { Operator{T}} (undef ,N+ 1 ); fill! ( Abk1y , Z)
1301
+ Abk2x= Array { Operator{T}} (undef ,N+ 1 ); fill! ( Abk2x , Z)
1302
+ Abk2y= Array { Operator{T}} (undef ,N+ 1 ); fill! ( Abk2y , Z)
1303
1303
1304
1304
for K= Block (N): - 1 : Block (2 )
1305
1305
Bx,By= view (Jx,K,K),view (Jy,K,K)
@@ -1317,7 +1317,7 @@ function operator_clenshaw2D(Jx,Jy,cfs::Vector{Vector{T}},x,y) where T
1317
1317
resize! (Abk2x,Int (K))
1318
1318
Abk2x[:]= bk1[1 : Int (K)]
1319
1319
resize! (Abk2y,Int (K))
1320
- Abk2y[1 : Int (K)- 1 ]= Z
1320
+ Abk2y[1 : Int (K)- 1 ] . = Ref (Z)
1321
1321
Abk2y[end ]= bk1[Int (K)+ 1 ]
1322
1322
1323
1323
Abk1x,Abk2x= Abk2x,Abk1x
0 commit comments