David Cournapeau wrote: > On Wed, Nov 7, 2012 at 1:56 PM, Neal Becker <ndbeck...@gmail.com> wrote: >> David Cournapeau wrote: >> >>> On Wed, Nov 7, 2012 at 12:35 PM, Neal Becker <ndbeck...@gmail.com> wrote: >>>> I'm trying to do a bit of benchmarking to see if amd libm/acml will help >>>> me. >>>> >>>> I got an idea that instead of building all of numpy/scipy and all of my >>>> custom modules against these libraries, I could simply use: >>>> >>>> >> LD_PRELOAD=/opt/amdlibm-3.0.2/lib/dynamic/libamdlibm.so:/opt/acml5.2.0/gfortran64/lib/libacml.so >>>> <my program here> >>>> >>>> I'm hoping that both numpy and my own dll's then will take advantage of >>>> these libraries. >>>> >>>> Do you think this will work? >>> >>> Quite unlikely depending on your configuration, because those >>> libraries are rarely if ever ABI compatible (that's why it is such a >>> pain to support). >>> >>> David >> >> When you say quite unlikely (to work), you mean >> >> a) unlikely that libm/acml will be used to resolve symbols in numpy/dlls at >> runtime (e.g., exp)? >> >> or >> >> b) program may produce wrong results and/or crash ? > > Both, actually. That's not something I would use myself. Did you try > openblas ? It is open source, simple to build, and is pretty fast, > > David
In my current work, probably the largest bottlenecks are 'max*', which are log (\sum e^(x_i)) _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion