66
77
88![ JSON 0.21.4] ( https://img.shields.io/badge/JSON-0.21.4-000?style=flat-square&labelColor=999 )
9- [ ![ LinearOperators 2.5.2 ] ( https://img.shields.io/badge/LinearOperators-2.5.2 -4b0082?style=flat-square&labelColor=9558b2 )] ( https://jso.dev/LinearOperators.jl/stable/ )
10- [ ![ Krylov 0.9.4 ] ( https://img.shields.io/badge/Krylov-0.9.4 -4b0082?style=flat-square&labelColor=9558b2 )] ( https://jso.dev/Krylov.jl/stable/ )
11- ![ FFTW 1.7.1 ] ( https://img.shields.io/badge/FFTW-1.7.1 -000?style=flat-square&labelColor=999 )
9+ [ ![ LinearOperators 2.11.0 ] ( https://img.shields.io/badge/LinearOperators-2.11.0 -4b0082?style=flat-square&labelColor=9558b2 )] ( https://jso.dev/LinearOperators.jl/stable/ )
10+ [ ![ Krylov 0.9.10 ] ( https://img.shields.io/badge/Krylov-0.9.10 -4b0082?style=flat-square&labelColor=9558b2 )] ( https://jso.dev/Krylov.jl/stable/ )
11+ ![ FFTW 1.10.0 ] ( https://img.shields.io/badge/FFTW-1.10.0 -000?style=flat-square&labelColor=999 )
1212
1313
1414
@@ -38,8 +38,8 @@ y = op * x
3838``` plaintext
39395-element Vector{Float64}:
4040 0.24766385713448147
41- 0.34471427066907173
42- 0.6399025892338569
41+ 0.3447142706690718
42+ 0.6399025892338568
4343 0.246664873173519
4444 0.5258420478788777
4545```
@@ -62,7 +62,7 @@ norm(A \ v - op * v) / norm(v)
6262```
6363
6464``` plaintext
65- 1.148733996146629e-13
65+ 3.3842227055482783e-15
6666```
6767
6868
@@ -138,7 +138,7 @@ norm(dft' * y - x) # DFT is a unitary operator
138138```
139139
140140``` plaintext
141- 3.3612925285989963e -16
141+ 1.1320399892723098e -16
142142```
143143
144144
@@ -149,16 +149,16 @@ transpose(dft) * y
149149
150150``` plaintext
15115110-element Vector{ComplexF64}:
152- 0.7113168325963566 - 0.0im
153- 0.2148297961020408 - 0.0im
154- 0.9453427110792467 - 0.0im
155- 0.7309321298104021 - 0.0im
156- 0.23194328498409336 - 0.0im
157- 0.9501874396162999 - 0.0im
158- 0.5123847829172379 - 0.0im
159- 0.9037088931078092 - 0.0im
160- 0.005729552514365821 + 0.0im
161- 0.2588088046728426 - 0.0im
152+ 0.003820639240554513 - 0.0im
153+ 0.8695984692566511 - 0.0im
154+ 0.6447374729383316 - 0.0im
155+ 0.8365781167218365 - 0.0im
156+ 0.7119679627204016 - 0.0im
157+ 0.4093590753517967 - 0.0im
158+ 0.8669244192719149 - 0.0im
159+ 0.4754090095137926 - 0.0im
160+ 0.5455204733218788 + 0.0im
161+ 0.042978926103728554 - 0.0im
162162```
163163
164164
280280```
281281
282282``` plaintext
283- ex = InexactError(:Float64, Float64, 0.47507530806279163 - 0.3832301880068187im )
283+ ex = InexactError(:Float64, ( Float64, 0.5164858914848688 + 0.24488071994631527im) )
284284```
285285
286286
296296```
297297
298298``` plaintext
299- ex = InexactError(:Float64, Float64, 0.8090169943749475 - 0.5877852522924731im)
299+ ex = InexactError(:Float64, ( Float64, 0.8090169943749475 - 0.5877852522924731im) )
300300```
301301
302302
@@ -317,7 +317,7 @@ norm(b - opAAT * x)
317317```
318318
319319``` plaintext
320- 2.3624810325472306e -13
320+ 1.9772018739020697e -13
321321```
322322
323323
344344```
345345
346346``` plaintext
347- 3.919077465054229e -13
347+ 3.815809259303818e -13
348348```
349349
350350
@@ -355,7 +355,6 @@ There is also a Limited-Memory SR1 operator for which only the forward form is i
355355Note that the SR1 operator can be indefinite; therefore, its inverse form is less relevant than for the BFGS approximation.
356356For this reason, the inverse form is not implemented for the SR1 operator.
357357
358-
359358## Restriction, extension and slices
360359
361360The restriction operator restricts a vector to a set of indices.
@@ -425,9 +424,9 @@ A[I,J] * ones(3)
425424
426425``` plaintext
4274263-element Vector{Float64}:
428- 1.043315409692812
429- 1.126646200633167
430- 1.2315932943898407
427+ 1.855452693501958
428+ 1.2659861763815514
429+ 1.5303428526550302
431430```
432431
433432
@@ -438,9 +437,9 @@ opRestriction(I, 5) * opA * opExtension(J, 5) * ones(3)
438437
439438``` plaintext
4404393-element Vector{Float64}:
441- 1.043315409692812
442- 1.126646200633167
443- 1.2315932943898407
440+ 1.8554526935019582
441+ 1.2659861763815514
442+ 1.5303428526550302
444443```
445444
446445
@@ -455,7 +454,7 @@ opA[1,:] * ones(5)
455454
456455``` plaintext
4574561-element Vector{Float64}:
458- 1.9634567072673046
457+ 2.9159013209470097
459458```
460459
461460
@@ -466,11 +465,11 @@ opA[:,1] * ones(1)
466465
467466``` plaintext
4684675-element Vector{Float64}:
469- 0.02203447865143171
470- 0.18781833339603016
471- 0.1692361439290926
472- 0.20276666600590854
473- 0.24969824326802326
468+ 0.7011787151612684
469+ 0.4203121328637791
470+ 0.19620605317757256
471+ 0.6722900514676972
472+ 0.7640075887373438
474473```
475474
476475
@@ -481,6 +480,6 @@ opA[1,1] * ones(1)
481480
482481``` plaintext
4834821-element Vector{Float64}:
484- 0.02203447865143171
483+ 0.7011787151612684
485484```
486485
0 commit comments