@@ -151,9 +151,9 @@ function spmatmul{Tv,Ti}(A::SparseMatrixCSC{Tv,Ti}, B::SparseMatrixCSC{Tv,Ti};
151
151
colptrB = B. colptr; rowvalB = B. rowval; nzvalB = B. nzval
152
152
# TODO : Need better estimation of result space
153
153
nnzC = min (mA* nB, length (nzvalA) + length (nzvalB))
154
- colptrC = Array (Ti, nB+ 1 )
155
- rowvalC = Array (Ti, nnzC)
156
- nzvalC = Array (Tv, nnzC)
154
+ colptrC = Array {Ti} ( nB+ 1 )
155
+ rowvalC = Array {Ti} ( nnzC)
156
+ nzvalC = Array {Tv} ( nnzC)
157
157
158
158
@inbounds begin
159
159
ip = 1
@@ -302,7 +302,7 @@ function triu{Tv,Ti}(S::SparseMatrixCSC{Tv,Ti}, k::Integer=0)
302
302
if (k > 0 && k > n) || (k < 0 && - k > m)
303
303
throw (BoundsError ())
304
304
end
305
- colptr = Array (Ti, n+ 1 )
305
+ colptr = Array {Ti} ( n+ 1 )
306
306
nnz = 0
307
307
for col = 1 : min (max (k+ 1 ,1 ), n+ 1 )
308
308
colptr[col] = 1
@@ -314,8 +314,8 @@ function triu{Tv,Ti}(S::SparseMatrixCSC{Tv,Ti}, k::Integer=0)
314
314
end
315
315
colptr[col+ 1 ] = nnz+ 1
316
316
end
317
- rowval = Array (Ti, nnz)
318
- nzval = Array (Tv, nnz)
317
+ rowval = Array {Ti} ( nnz)
318
+ nzval = Array {Tv} ( nnz)
319
319
A = SparseMatrixCSC (m, n, colptr, rowval, nzval)
320
320
for col = max (k+ 1 ,1 ) : n
321
321
c1 = S. colptr[col]
@@ -333,7 +333,7 @@ function tril{Tv,Ti}(S::SparseMatrixCSC{Tv,Ti}, k::Integer=0)
333
333
if (k > 0 && k > n) || (k < 0 && - k > m)
334
334
throw (BoundsError ())
335
335
end
336
- colptr = Array (Ti, n+ 1 )
336
+ colptr = Array {Ti} ( n+ 1 )
337
337
nnz = 0
338
338
colptr[1 ] = 1
339
339
for col = 1 : min (n, m+ k)
@@ -347,8 +347,8 @@ function tril{Tv,Ti}(S::SparseMatrixCSC{Tv,Ti}, k::Integer=0)
347
347
for col = max (min (n, m+ k)+ 2 ,1 ) : n+ 1
348
348
colptr[col] = nnz+ 1
349
349
end
350
- rowval = Array (Ti, nnz)
351
- nzval = Array (Tv, nnz)
350
+ rowval = Array {Ti} ( nnz)
351
+ nzval = Array {Tv} ( nnz)
352
352
A = SparseMatrixCSC (m, n, colptr, rowval, nzval)
353
353
for col = 1 : min (n, m+ k)
354
354
c1 = S. colptr[col+ 1 ]- 1
@@ -367,10 +367,10 @@ end
367
367
function sparse_diff1 {Tv,Ti} (S:: SparseMatrixCSC{Tv,Ti} )
368
368
m,n = size (S)
369
369
m > 1 || return SparseMatrixCSC (0 , n, ones (Ti,n+ 1 ), Ti[], Tv[])
370
- colptr = Array (Ti, n+ 1 )
370
+ colptr = Array {Ti} ( n+ 1 )
371
371
numnz = 2 * nnz (S) # upper bound; will shrink later
372
- rowval = Array (Ti, numnz)
373
- nzval = Array (Tv, numnz)
372
+ rowval = Array {Ti} ( numnz)
373
+ nzval = Array {Tv} ( numnz)
374
374
numnz = 0
375
375
colptr[1 ] = 1
376
376
for col = 1 : n
@@ -406,10 +406,10 @@ end
406
406
407
407
function sparse_diff2 {Tv,Ti} (a:: SparseMatrixCSC{Tv,Ti} )
408
408
m,n = size (a)
409
- colptr = Array (Ti, max (n,1 ))
409
+ colptr = Array {Ti} ( max (n,1 ))
410
410
numnz = 2 * nnz (a) # upper bound; will shrink later
411
- rowval = Array (Ti, numnz)
412
- nzval = Array (Tv, numnz)
411
+ rowval = Array {Ti} ( numnz)
412
+ nzval = Array {Tv} ( numnz)
413
413
414
414
z = zero (Tv)
415
415
@@ -561,8 +561,8 @@ function normestinv{T}(A::SparseMatrixCSC{T}, t::Integer = min(2,maximum(size(A)
561
561
if t > n
562
562
throw (ArgumentError (" number of blocks must not be greater than $n " ))
563
563
end
564
- ind = Array ( Int64, n)
565
- ind_hist = Array ( Int64, maxiter * t)
564
+ ind = Array { Int64} ( n)
565
+ ind_hist = Array { Int64} ( maxiter * t)
566
566
567
567
Ti = typeof (float (zero (T)))
568
568
@@ -584,7 +584,7 @@ function normestinv{T}(A::SparseMatrixCSC{T}, t::Integer = min(2,maximum(size(A)
584
584
end
585
585
586
586
# Generate the block matrix
587
- X = Array (Ti, n, t)
587
+ X = Array {Ti} ( n, t)
588
588
X[1 : n,1 ] = 1
589
589
for j = 2 : t
590
590
while true
@@ -734,9 +734,9 @@ function kron{Tv,Ti}(a::SparseMatrixCSC{Tv,Ti}, b::SparseMatrixCSC{Tv,Ti})
734
734
735
735
m,n = mA* mB, nA* nB
736
736
737
- colptr = Array (Ti, n+ 1 )
738
- rowval = Array (Ti, numnz)
739
- nzval = Array (Tv, numnz)
737
+ colptr = Array {Ti} ( n+ 1 )
738
+ rowval = Array {Ti} ( numnz)
739
+ nzval = Array {Tv} ( numnz)
740
740
741
741
colptr[1 ] = 1
742
742
0 commit comments