René's URL Explorer Experiment

{"@context":"https://schema.org","@type":"DiscussionForumPosting","headline":" solve() inside gfor() for a large number of independent matrices ","articleBody":"Right now, `solve()` is not optimized to solve for a large number of independent linear systems `A_i x_i = b_i`, where `i = 1,..., N`\n\n`solve()` uses `magma_getrf_gpu()` (https://github.com/maxhutch/magma/blob/master/src/sgetrs_gpu.cpp) from Magma as its backend. Magma also has `magma_sgetrs_batched()` (https://github.com/maxhutch/magma/blob/master/src/sgetrs_batched.cpp) which could be used for batch computations (like those inside gfor?)\n\nPossibly connected issue: https://github.com/arrayfire/arrayfire/issues/483\n\nRelevant post on forums: https://groups.google.com/forum/#!topic/arrayfire-users/RZLKGjsJR3g\n\nDetails of the batched implementation in Magma:\n1) http://icl.cs.utk.edu/news_pub/submissions/batchedlu.pdf\n2) http://on-demand.gputechconf.com/gtc/2015/presentation/S5476-Stanimire-Tomov.pdf (starts from slide 36)\n","author":{"url":"https://github.com/mchandra","@type":"Person","name":"mchandra"},"datePublished":"2015-10-26T19:25:24.000Z","interactionStatistic":{"@type":"InteractionCounter","interactionType":"https://schema.org/CommentAction","userInteractionCount":2},"url":"https://github.com/1075/arrayfire/issues/1075"}