Computing Low-Rank Approximation of a Dense Matrix on Multicore CPUs with a GPU and Its Application to Solving a Hierarchically Semiseparable Linear System of Equations

<div>Performance of StruMF solving 3D Poisson equation using the GPU (<svg height="12.1266pt" id="M301" style="vertical-align:-3.42784pt" version="1.1" viewbox="-0.0498162 -8.69876 47.9081 12.1266" width="47.9081pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.0135,0,0,-0.0135,0,0)"><path d="M495 86L479 114C446 82 419 66 409 66C401 66 401 72 406 97C420 166 436 231 453 297C489 435 454 448 428 448C406 448 384 439 354 422C305 394 222 327 161 247H159L183 345C200 415 194 448 173 448C143 448 82 410 23 351L38 325C64 349 95 371 105 371C111 371 116 365 109 336L25 -4L31 -12C50 -4 77 3 107 9C119 69 132 122 145 168C197 254 321 381 370 381C387 381 393 374 378 305L329 95C309 17 320 -12 345 -12C372 -12 430 19 495 86Z" id="g113-111"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="502" vert-adv-y="502"></glyph.data></g><g transform="matrix(.0095,0,0,-0.0095,6.733,3.264)"><path d="M545 403C545 425 527 451 498 451C448 451 400 404 312 286L290 341C262 412 246 451 219 451C182 451 138 404 92 345L112 323C152 368 169 378 181 378C192 378 201 359 216 321L257 214C184 115 142 69 113 69C102 69 93 73 88 79S78 89 68 89C47 89 24 61 24 40C24 8 49 -12 75 -12C125 -12 175 31 271 175L312 64C330 15 357 -12 382 -12C421 -12 475 35 515 98L496 121C466 88 439 62 420 62C402 62 384 92 363 147L325 247C349 280 374 309 397 333C422 361 446 381 462 381C471 381 481 374 488 366C493 360 499 357 505 357C521 357 545 380 545 403Z" id="g50-121"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="569" vert-adv-y="569"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,16.433,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,28.171,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,34.674,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,41.177,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>).</div>

Scientific Programming

fig9

Figure 9

Figure 9: Computing Low-Rank Approximation of a Dense Matrix on Multicore CPUs with a GPU and Its Application to Solving a Hierarchically Semiseparable Linear System of Equations 

Figure 9 | Computing Low-Rank Approximation of a Dense Matrix on Multicore CPUs with a GPU and Its Application to Solving a Hierarchically Semiseparable Linear System of Equations 