@@ -2,7 +2,7 @@ using LinearSolve, LinearAlgebra, SparseArrays, MultiFloats, ForwardDiff
2
2
using Test
3
3
import Random
4
4
5
- const Dual64= ForwardDiff. Dual{Nothing,Float64,1 }
5
+ const Dual64 = ForwardDiff. Dual{Nothing, Float64, 1 }
6
6
7
7
n = 8
8
8
A = Matrix (I, n, n)
@@ -141,11 +141,11 @@ end
141
141
142
142
@testset " Sparspak Factorization (Float64x1)" begin
143
143
A1 = sparse (A / 1 ) .| > Float64x1
144
- b1 = rand (n) .| > Float64x1
145
- x1 = zero (b) .| > Float64x1
146
- A2 = sparse (A / 2 ) .| > Float64x1
147
- b2 = rand (n) .| > Float64x1
148
- x2 = zero (b) .| > Float64x1
144
+ b1 = rand (n) .| > Float64x1
145
+ x1 = zero (b) .| > Float64x1
146
+ A2 = sparse (A / 2 ) .| > Float64x1
147
+ b2 = rand (n) .| > Float64x1
148
+ x2 = zero (b) .| > Float64x1
149
149
150
150
prob1 = LinearProblem (A1, b1; u0 = x1)
151
151
prob2 = LinearProblem (A2, b2; u0 = x2)
@@ -154,11 +154,11 @@ end
154
154
155
155
@testset " Sparspak Factorization (Float64x2)" begin
156
156
A1 = sparse (A / 1 ) .| > Float64x2
157
- b1 = rand (n) .| > Float64x2
158
- x1 = zero (b) .| > Float64x2
159
- A2 = sparse (A / 2 ) .| > Float64x2
160
- b2 = rand (n) .| > Float64x2
161
- x2 = zero (b) .| > Float64x2
157
+ b1 = rand (n) .| > Float64x2
158
+ x1 = zero (b) .| > Float64x2
159
+ A2 = sparse (A / 2 ) .| > Float64x2
160
+ b2 = rand (n) .| > Float64x2
161
+ x2 = zero (b) .| > Float64x2
162
162
163
163
prob1 = LinearProblem (A1, b1; u0 = x1)
164
164
prob2 = LinearProblem (A2, b2; u0 = x2)
@@ -167,18 +167,17 @@ end
167
167
168
168
@testset " Sparspak Factorization (Dual64)" begin
169
169
A1 = sparse (A / 1 ) .| > Dual64
170
- b1 = rand (n) .| > Dual64
171
- x1 = zero (b) .| > Dual64
172
- A2 = sparse (A / 2 ) .| > Dual64
173
- b2 = rand (n) .| > Dual64
174
- x2 = zero (b) .| > Dual64
170
+ b1 = rand (n) .| > Dual64
171
+ x1 = zero (b) .| > Dual64
172
+ A2 = sparse (A / 2 ) .| > Dual64
173
+ b2 = rand (n) .| > Dual64
174
+ x2 = zero (b) .| > Dual64
175
175
176
176
prob1 = LinearProblem (A1, b1; u0 = x1)
177
177
prob2 = LinearProblem (A2, b2; u0 = x2)
178
178
test_interface (SparspakFactorization (), prob1, prob2)
179
179
end
180
180
181
-
182
181
@testset " FastLAPACK Factorizations" begin
183
182
A1 = A / 1
184
183
b1 = rand (n)
0 commit comments