[eigen] sse asin implementation |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: [eigen] sse asin implementation*From*: Rohit Garg <rpg.314@xxxxxxxxx>*Date*: Mon, 30 Mar 2009 19:22:48 +0530*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:date:message-id:subject :from:to:content-type; bh=MjtRgbp8yZbaJjUlXZ/yNStZ7jSnU5xONZiUm5h+S4Y=; b=pe4oN31m6t4baf+SWqdfwBHFm2Tc3UlImUanZZGhXux7DfbZDyLCdOyZnFhp6pev17 ENsfMkdvMUXokAVyYpZnCqlW2lCWjMTcCwm/EnD2mByhajLLDFmUoJqt5NkoKuXdOP8w CoudpPf/c5257QWHH8C1wf3JVMLAaQ+mBrh+U=*Domainkey-signature*: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=YgrKPds7nIw53cVuDL7Rl0MR9dEf7pU8nxFNsENhrcNa4hxs5qUPU7AfQ3zUhTKOwP dEP1U4g5qvEYASnIoyw8wGIoW8E9YKY8DLFq/IXrdo5uOmrnDJz3vfLsB6CoqeFRqaWW UTU8fZ9rsuyXDBLlcxz8p1ZQV0XcshmCkCsKw=

This file has the patch for sse implementation of asin. I have written it along the lines of the cephes library. I have dropped one branch which cephes takes that for x<1e-4 asin(x) is the same as x. Also I don't know how to handle out of range inputs, so I dropped the error checking. Perhaps, the user should be more careful before calling this function. I dunno if it will return NaN for out of range inputs. I wrote this in the MathFunctions.h file using as much of eigen's abstraction so someone porting it to Altivec won't curse me too much. :) Also I noted that Gael while committing my input had to write a fair amount of top level code since I was writing pure intrinsic functions. So I cut some slack here and decided to write it as close to eigen's final implementation as I could. Unfortunately, it means that since I am not sure of what else do I need to patch, I couldn't test it. :( But it is a very simple one, so perhaps no errors made it through. :) On a related note, perhaps we should have an ei_psel function since it is going to be used a lot (in transcendental functions atleast) so would help with porting to altivec. -- Rohit Garg http://rpg-314.blogspot.com/ Senior Undergraduate Department of Physics Indian Institute of Technology Bombay

**Attachment:
asin_patch**

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Another LDLt issue** - Next by Date:
**[eigen] discussion on special matrices** - Previous by thread:
**Re: [eigen] Fwd: problem with compile on AIX system** - Next by thread:
**Re: [eigen] sse asin implementation**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |