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

Reply via email to