Accelerated Computing in Magnetic Resonance Imaging: Real-Time Imaging Using Nonlinear Inverse Reconstruction

<div>Run time (i.e., minimal wall-clock time in ms) measurement of a single-precision 2D squared complex-to-complex out-of-place Fourier transform of FFTW, clFFT, and cuFFT libraries. clFFT and cuFFT benchmarks were obtained for NVIDIA GeForce Titan Black; the FFTW benchmark was obtained for Intel Xeon E5-2650 utilizing 8 threads. FFTW is initialized with <span class="nowrap"><svg height="8.98583pt" id="M1" style="vertical-align:-0.2324905pt" version="1.1" viewbox="-0.0498162 -8.75334 98.7613 8.98583" width="98.7613pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M493 503C489 551 484 614 483 650H43V622C120 616 128 611 128 525V126C128 40 120 34 40 28V0H312V28C221 34 213 40 213 126V307H316C407 307 412 296 424 227H453V420H424C412 355 407 346 316 346H213V584C213 613 216 616 246 616H322C398 616 419 607 436 579C449 559 455 539 464 499L493 503Z" id="g190-71"></path></g><g transform="matrix(.013,0,0,-0.013,6.877,0)"><path d="M493 503C489 551 484 614 483 650H43V622C120 616 128 611 128 525V126C128 40 120 34 40 28V0H312V28C221 34 213 40 213 126V307H316C407 307 412 296 424 227H453V420H424C412 355 407 346 316 346H213V584C213 613 216 616 246 616H322C398 616 419 607 436 579C449 559 455 539 464 499L493 503Z" id="g190-71"></path></g><g transform="matrix(.013,0,0,-0.013,13.754,0)"><path d="M592 498C586 556 582 629 581 675H560C545 656 537 650 510 650H115C87 650 74 652 60 675H40C38 620 33 555 28 495H57C69 539 77 568 89 585C102 607 122 616 205 616H267V123C267 44 259 34 168 28V0H454V28C360 34 352 44 352 123V616H422C498 616 514 608 531 584C542 568 552 543 562 495L592 498Z" id="g190-85"></path></g><g transform="matrix(.013,0,0,-0.013,22.127,0)"><path d="M956 650H736V622C814 613 818 602 804 542C789 471 733 275 691 141H689C628 314 570 476 514 641H485C434 478 377 312 323 142H321C277 284 238 419 199 554C184 605 186 613 259 622V650H16V622C77 615 92 606 117 524C172 346 228 170 281 -11H312C363 151 419 312 475 488H477C528 343 597 150 652 -11H683C765 269 826 458 848 522C878 605 887 613 956 622V650Z" id="g190-88"></path></g><rect height="0.3985" width="3.91466" x="35.5333" y="-0.3985"></rect><g transform="matrix(.013,0,0,-0.013,39.448,0)"><path d="M861 0V28C774 35 771 41 768 147L759 509C756 612 762 614 851 622V650H681L449 149L221 650H57V622C148 613 153 609 144 479L130 271C123 166 117 123 111 88C104 46 85 34 26 28V0H259V28C192 35 169 42 167 90C166 130 166 173 170 256L185 541H187L411 7H431L675 555H679L683 147C683 41 680 35 598 28V0H861Z" id="g190-78"></path></g><g transform="matrix(.013,0,0,-0.013,51.031,0)"><path d="M517 162C503 123 484 88 467 68C445 42 417 34 341 34C291 34 256 34 237 47C219 59 213 81 213 128V317H308C395 317 402 311 415 240H444V431H415C403 364 398 356 307 356H213V584C213 613 215 616 246 616H322C394 616 421 609 435 587C448 566 458 544 467 502L496 506C493 557 488 625 488 650H42V622C120 616 128 612 128 523V125C128 43 120 35 29 28V0H511C520 31 540 125 546 158L517 162Z" id="g190-70"></path></g><g transform="matrix(.013,0,0,-0.013,58.415,0)"><path d="M673 0V28C608 34 594 43 563 129C499 303 432 494 370 665L339 656L132 132C98 43 84 37 20 28V0H238V28C162 35 154 45 175 108C188 150 203 192 219 237H432C454 176 474 122 487 83S491 36 428 28V0H673ZM418 280H234C265 362 296 450 328 535H330L418 280Z" id="g190-66"></path></g><g transform="matrix(.013,0,0,-0.013,67.359,0)"><path d="M409 504C401 567 396 607 392 642C354 654 312 665 266 665C137 665 60 583 60 487C60 374 161 325 225 290C300 250 355 215 355 141C355 68 311 21 235 21C131 21 86 122 71 183L41 176C48 128 61 42 68 21C78 16 93 8 118 0C142 -7 175 -15 216 -15C349 -15 438 69 438 174C438 287 344 333 265 374C186 414 138 449 138 522C138 576 172 631 249 631C336 631 363 562 380 499L409 504Z" id="g190-84"></path></g><g transform="matrix(.013,0,0,-0.013,73.286,0)"><path d="M716 650H478V622C549 617 570 605 575 566C579 535 582 493 582 401V326C582 161 537 40 393 40C269 40 205 133 205 294V521C205 611 213 617 291 622V650H33V622C111 615 120 611 120 521V272C120 154 148 80 208 34C252 0 309 -14 368 -14C434 -14 501 7 548 55C605 116 623 201 623 319V401C623 493 626 536 629 563C633 605 656 618 716 622V650Z" id="g190-86"></path></g><g transform="matrix(.013,0,0,-0.013,82.854,0)"><path d="M631 18C609 24 585 35 559 65C534 91 514 117 478 169C448 214 406 281 389 313C462 346 516 399 516 485C516 545 490 590 449 616C412 641 363 650 290 650H42V622C120 615 128 612 128 527V125C128 40 120 34 38 28V0H300V28C221 34 212 40 212 125V284H244C295 284 312 272 329 244C359 195 395 133 430 84C475 19 516 -3 592 -7C603 -8 615 -8 627 -8L631 18ZM212 316V563C212 591 215 602 223 607C231 613 248 617 277 617C352 617 423 577 423 469C423 415 407 375 368 345C343 324 310 316 260 316H212Z" id="g190-83"></path></g><g transform="matrix(.013,0,0,-0.013,90.927,0)"><path d="M517 162C503 123 484 88 467 68C445 42 417 34 341 34C291 34 256 34 237 47C219 59 213 81 213 128V317H308C395 317 402 311 415 240H444V431H415C403 364 398 356 307 356H213V584C213 613 215 616 246 616H322C394 616 421 609 435 587C448 566 458 544 467 502L496 506C493 557 488 625 488 650H42V622C120 616 128 612 128 523V125C128 43 120 35 29 28V0H511C520 31 540 125 546 158L517 162Z" id="g190-70"></path></g></svg>.</span> The clFFT library has only limited support for mixed radix FFTs and achieves good performance for power of 2 FFTs only.</div>

Computational and Mathematical Methods in Medicine

fig1

Figure 1

Figure 1: Accelerated Computing in Magnetic Resonance Imaging: Real-Time Imaging Using Nonlinear Inverse Reconstruction 