Skip to content

Speed up pulseq_rs import by integrating gradient moments in NumPy. - #194

Open
mzaiss wants to merge 1 commit into
mainfrom
dev_mz_pulseq_import
Open

mzaiss wants to merge 1 commit into
mainfrom
dev_mz_pulseq_import

Conversation

@mzaiss

@mzaiss mzaiss commented Sep 26, 2026

Copy link
Copy Markdown
Collaborator

Resolves issue #175 The parser was already fast. The conversion called PyGradient.integrate once per event. Moments now come from one cumulative piecewise-linear integral per gradient, which covers trapezoids, extended trapezoids, and arbitrary shapes. Results match the old importer.

exact_trajectories=True (default):

File Old New
pp_mprage.seq 26 s 3.6 s
pp_radial.seq 0.80 s 0.22 s
mrseq_spiral.seq 5.6 s 0.43 s
pp_tse.seq 0.41 s 0.10 s
mrseq_radial_1_5.seq 3.0 s 0.24 s
mrseq_spiral_1.5.seq 3.0 s 0.23 s

The parser was already fast. The conversion called PyGradient.integrate once per event. Moments now come from one cumulative piecewise-linear integral per gradient, which covers trapezoids, extended trapezoids, and arbitrary shapes. Results match the old importer.

exact_trajectories=True (default):

| File | Old | New |
|---|---:|---:|
| pp_mprage.seq | 26 s | 3.6 s |
| pp_radial.seq | 0.80 s | 0.22 s |
| mrseq_spiral.seq | 5.6 s | 0.43 s |
| pp_tse.seq | 0.41 s | 0.10 s |
| mrseq_radial_1_5.seq | 3.0 s | 0.24 s |
| mrseq_spiral_1.5.seq | 3.0 s | 0.23 s |

This branch has not been deployed

No deployments
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant