multithreading - Julia and Lapack: pstrf multithreaded but trtrs not -


my program written in julia not yield expected computational performance. program first computes cholesky decomposition of large matrix a using cholfact!, a = l'l. solves lx = b different b using backslash operator.

this results in straight calls lapack. function cholfact! implemented pstrf! , backslash operator uses trtrs!. these correct lapack functions use. while function pstrf! executed in parallel, function trtrs! not. profiler tells me of runtime spent on trtrs!. lines of code in program are

f = cholfact!(a, :l, pivot = true)     # precomputation, executed once 

and

x = f[:l]\b[f.piv]    # inside loop, b computed x every step 

why there difference between 2 lapack functions? how can parallel execution of pstrf!?


Comments

Popular posts from this blog

c++ - QTextObjectInterface with Qml TextEdit (QQuickTextEdit) -

javascript - angular ng-required radio button not toggling required off in firefox 33, OK in chrome -

xcode - Swift Playground - Files are not readable -