Add multithreading support? about fastgaussquadrature.jl HOT 3 OPEN

Have you done a profiling of the current gausslegendre? Is it Bogaert's method? (A citation should be added to the code.)

For single threads, there is pretty much no reason this kind of code would be faster in C than in Julia. But there seems to be some sillyness in the current implementation, e.g., constants like -65 / 1024 which may be evaluated every loop.

from fastgaussquadrature.jl.

Note that issue is speculative since threading is still experimental in Julia. In my personal experience, Julia's experimental threading is still about a factor of 30% off from OpenMP paradigms and less flexible too.

As far as I can tell, the current implementation is a transliteration of Bogaert's method (though some terms in the series are pared down, and you're right about the constants that are fed into the macros).

For reference, here is the original C++ that appears to have been released in the public domain and here is a C version.

from fastgaussquadrature.jl.

are these threading comments still relevant? Julia's threading is pretty solid now. Anyone done any recent tests of timing?

from fastgaussquadrature.jl.