1
1
module LAPACK2
2
2
3
+ using libblastrampoline_jll
3
4
using LinearAlgebra
4
5
using LinearAlgebra: BlasInt, chkstride1, LAPACKException
5
6
using LinearAlgebra. BLAS: @blasfunc
6
7
using LinearAlgebra. LAPACK: chkdiag, chkside, chkuplo
7
8
9
+ liblapack_name = libblastrampoline_jll. libblastrampoline
10
+
8
11
# # Standard QR/QL
9
12
function steqr! (compz:: Char ,
10
13
d:: StridedVector{Float64} ,
@@ -34,7 +37,7 @@ module LAPACK2
34
37
# Allocations
35
38
info = Vector {BlasInt} (undef, 1 )
36
39
37
- ccall ((@blasfunc (" dsteqr_" ), Base . liblapack_name),Cvoid,
40
+ ccall ((@blasfunc (" dsteqr_" ), liblapack_name),Cvoid,
38
41
(Ref{UInt8}, Ref{BlasInt}, Ptr{Float64}, Ptr{Float64},
39
42
Ptr{Float64}, Ref{BlasInt}, Ptr{Float64}, Ptr{BlasInt}),
40
43
compz, n, d, e,
@@ -56,7 +59,7 @@ module LAPACK2
56
59
# Allocations
57
60
info = BlasInt[0 ]
58
61
59
- ccall ((@blasfunc (" dsterf_" ), Base . liblapack_name), Cvoid,
62
+ ccall ((@blasfunc (" dsterf_" ), liblapack_name), Cvoid,
60
63
(Ref{BlasInt}, Ptr{Float64}, Ptr{Float64}, Ptr{BlasInt}),
61
64
n, d, e, info)
62
65
@@ -87,7 +90,7 @@ module LAPACK2
87
90
# Allocations
88
91
info = BlasInt[0 ]
89
92
90
- ccall ((@blasfunc (" dstedc_" ), Base . liblapack_name), Cvoid,
93
+ ccall ((@blasfunc (" dstedc_" ), liblapack_name), Cvoid,
91
94
(Ref{UInt8}, Ref{BlasInt}, Ptr{Float64}, Ptr{Float64},
92
95
Ptr{Float64}, Ref{BlasInt}, Ptr{Float64}, Ref{BlasInt},
93
96
Ptr{BlasInt}, Ref{BlasInt}, Ptr{BlasInt}),
@@ -153,7 +156,7 @@ module LAPACK2
153
156
tryrac = BlasInt[1 ]
154
157
info = Vector {BlasInt} (undef, 1 )
155
158
156
- ccall ((@blasfunc ($ lsymb), Base . liblapack_name), Cvoid,
159
+ ccall ((@blasfunc ($ lsymb), liblapack_name), Cvoid,
157
160
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
158
161
Ptr{$ elty}, Ref{$ elty}, Ref{$ elty}, Ref{BlasInt},
159
162
Ref{BlasInt}, Ptr{BlasInt}, Ptr{$ elty}, Ptr{$ elty},
@@ -242,7 +245,7 @@ module LAPACK2
242
245
243
246
info = Ref {BlasInt} (0 )
244
247
245
- ccall ((@blasfunc (" dlahqr_" ), Base . liblapack_name), Cvoid,
248
+ ccall ((@blasfunc (" dlahqr_" ), liblapack_name), Cvoid,
246
249
(Ref{BlasInt}, Ref{BlasInt}, Ref{BlasInt}, Ref{BlasInt},
247
250
Ref{BlasInt}, Ptr{Float64}, Ref{BlasInt}, Ptr{Float64},
248
251
Ptr{Float64}, Ref{BlasInt}, Ref{BlasInt}, Ptr{Float64},
@@ -295,7 +298,7 @@ module LAPACK2
295
298
296
299
info = Vector {BlasInt} (undef, 1 )
297
300
298
- ccall ((@blasfunc (:dpteqr_ ), Base . liblapack_name), Cvoid,
301
+ ccall ((@blasfunc (:dpteqr_ ), liblapack_name), Cvoid,
299
302
(Ref{UInt8}, Ref{BlasInt}, Ptr{Float64}, Ptr{Float64},
300
303
Ptr{Float64}, Ref{BlasInt}, Ptr{Float64}, Ptr{BlasInt}),
301
304
compz, n, d, e,
@@ -321,7 +324,7 @@ module LAPACK2
321
324
liwork = BlasInt (- 1 )
322
325
info = BlasInt[0 ]
323
326
for i = 1 : 2
324
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
327
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
325
328
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
326
329
Ref{BlasInt}, Ptr{$ elty}, Ptr{$ elty}, Ref{BlasInt},
327
330
Ptr{BlasInt}, Ref{BlasInt}, Ptr{BlasInt}),
@@ -359,7 +362,7 @@ module LAPACK2
359
362
liwork = BlasInt (- 1 )
360
363
info = BlasInt[0 ]
361
364
for i = 1 : 2
362
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
365
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
363
366
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
364
367
Ref{BlasInt}, Ptr{$ relty}, Ptr{$ elty}, Ref{BlasInt},
365
368
Ptr{$ relty}, Ref{BlasInt}, Ptr{BlasInt}, Ref{BlasInt},
@@ -443,7 +446,7 @@ for (f, elty) in ((:dtgevc_, :Float64),
443
446
m = BlasInt[0 ]
444
447
info = BlasInt[0 ]
445
448
446
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
449
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
447
450
(Ref{UInt8}, Ref{UInt8}, Ptr{BlasInt}, Ref{BlasInt},
448
451
Ptr{$ elty}, Ref{BlasInt}, Ptr{$ elty}, Ref{BlasInt},
449
452
Ptr{$ elty}, Ref{BlasInt}, Ptr{$ elty}, Ref{BlasInt},
@@ -520,7 +523,7 @@ for (f, elty, relty) in ((:dsfrk_, :Float64, :Float64),
520
523
end
521
524
lda = max (1 , stride (A, 2 ))
522
525
523
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
526
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
524
527
(Ref{UInt8}, Ref{UInt8}, Ref{UInt8}, Ref{BlasInt},
525
528
Ref{BlasInt}, Ref{$ relty}, Ptr{$ elty}, Ref{BlasInt},
526
529
Ref{$ relty}, Ptr{$ elty}),
@@ -544,7 +547,7 @@ for (f, elty) in ((:dpftrf_, :Float64),
544
547
n = round (Int,div (sqrt (8 length (A)), 2 ))
545
548
info = Vector {BlasInt} (undef, 1 )
546
549
547
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
550
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
548
551
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
549
552
Ptr{BlasInt}),
550
553
transr, uplo, n, A,
@@ -566,7 +569,7 @@ for (f, elty) in ((:dpftri_, :Float64),
566
569
n = round (Int,div (sqrt (8 length (A)), 2 ))
567
570
info = Vector {BlasInt} (undef, 1 )
568
571
569
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
572
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
570
573
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
571
574
Ptr{BlasInt}),
572
575
transr, uplo, n, A,
@@ -595,7 +598,7 @@ for (f, elty) in ((:dpftrs_, :Float64),
595
598
ldb = max (1 , stride (B, 2 ))
596
599
info = Vector {BlasInt} (undef, 1 )
597
600
598
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
601
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
599
602
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ref{BlasInt},
600
603
Ptr{$ elty}, Ptr{$ elty}, Ref{BlasInt}, Ptr{BlasInt}),
601
604
transr, uplo, n, nhrs,
@@ -631,7 +634,7 @@ for (f, elty) in ((:dtfsm_, :Float64),
631
634
end
632
635
ldb = max (1 , stride (B, 2 ))
633
636
634
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
637
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
635
638
(Ref{UInt8}, Ref{UInt8}, Ref{UInt8}, Ref{UInt8},
636
639
Ref{UInt8}, Ref{BlasInt}, Ref{BlasInt}, Ref{$ elty},
637
640
Ptr{$ elty}, Ptr{$ elty}, Ref{BlasInt}),
@@ -657,7 +660,7 @@ for (f, elty) in ((:dtftri_, :Float64),
657
660
n = round (Int,div (sqrt (8 length (A)), 2 ))
658
661
info = Vector {BlasInt} (undef, 1 )
659
662
660
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
663
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
661
664
(Ref{UInt8}, Ref{UInt8}, Ref{UInt8}, Ref{BlasInt},
662
665
Ptr{$ elty}, Ptr{BlasInt}),
663
666
transr, uplo, diag, n,
@@ -681,7 +684,7 @@ for (f, elty) in ((:dtfttr_, :Float64),
681
684
info = Vector {BlasInt} (undef, 1 )
682
685
A = similar (Arf, $ elty, n, n)
683
686
684
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
687
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
685
688
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
686
689
Ptr{$ elty}, Ref{BlasInt}, Ptr{BlasInt}),
687
690
transr, uplo, n, Arf,
@@ -707,7 +710,7 @@ for (f, elty) in ((:dtrttf_, :Float64),
707
710
info = Vector {BlasInt} (undef, 1 )
708
711
Arf = similar (A, $ elty, div (n* (n+ 1 ), 2 ))
709
712
710
- ccall ((@blasfunc ($ f), Base . liblapack_name), Cvoid,
713
+ ccall ((@blasfunc ($ f), liblapack_name), Cvoid,
711
714
(Ref{UInt8}, Ref{UInt8}, Ref{BlasInt}, Ptr{$ elty},
712
715
Ref{BlasInt}, Ptr{$ elty}, Ptr{BlasInt}),
713
716
transr, uplo, n, A,
0 commit comments