Jeff Fessler
Jeff Fessler
The README interface table https://github.com/JuliaLinearAlgebra/ToeplitzMatrices.jl#implemented-interface is missing `*` and `mul!` and probably more like `\` and `ldiv!`. Maybe they "go without saying" but I'm not 100% sure they all apply...
Currently `npzread` is "hardwired" to read from a file: https://github.com/fhs/NPZ.jl/blob/838df340de225c9bc523372c07cd8efb291469ea/src/NPZ.jl#L271 It would be more flexible to expose a version of it that allows reading from an IOStream, e.g., for data...
If I understand correctly, the current approach is that user decides between cpu and gpu versions by either * CPU: `p = plan_nfft(x, N)` or (currently equivalently I think) `p...
The result of a 1D NFFT of length N and a (N,1)-dimensional NFFT are quite inconsistent. Here is a MWE: ``` using NFFT M = 4 x1 = collect((-Int(M/2)):(Int(M/2)-1))/M N1...
spares us from more dependabot requests later... This should be preferable to #90
- 3d radial with simple DCF for `mri_example_3d` - CI e.g. per https://github.com/mathworks/ci-configuration-examples - mexPrintf so that mex errors go to console not terminal
MH reports that putting images from J-SENSE recon code output https://sigpy.readthedocs.io/en/latest/generated/sigpy.mri.app.JsenseRecon.html#sigpy.mri.app.JsenseRecon into `mri_field_map_reg.m` causes amplification of pixels outside the head in this line: https://github.com/JeffFessler/mirt/blob/43b20b5ea04b5bb61979aea7574ca73e3ea1b291/mri/fieldmap/mri_field_map_reg.m#L96 https://github.com/JeffFessler/mirt/blob/43b20b5ea04b5bb61979aea7574ca73e3ea1b291/mri/fieldmap/mri_field_map_reg.m#L148 Masking the fieldmaps (using ESPIRT...
The current code for `lldl` fails for Unitful sparse (or dense) arrays. MWE: ```julia using SparseArrays: spdiagm using LimitedLDLFactorizations: lldl using Unitful: s A = spdiagm(Float32.(1:3)*s^2) # unitful sparse array...
The README implies that this package is in the general registry, but I do not find it there: https://github.com/JuliaRegistries/General/tree/master/S Any idea why not? BTW, it would still be nice to...