tag:blogger.com,1999:blog-9318222.post113186340453101548..comments2017-11-11T13:14:22.152-08:00Comments on Connelly - Technical: FFT polynomial multiplicationConnelly Barneshttp://www.blogger.com/profile/02568908952592933174noreply@blogger.comBlogger10125tag:blogger.com,1999:blog-9318222.post-66203589237333916862016-06-14T16:55:17.615-07:002016-06-14T16:55:17.615-07:00وكالة نيوز
مباراة
1437
مصر
منوعات<a title="وكالة نيوز" href="http://www.news.wekalaa.com/" rel="nofollow">وكالة نيوز</a><br /><a title="مباراة" href="http://www.news.wekalaa.com/2016/" rel="nofollow">مباراة</a><br /><a title="1437" href="http://www.news.wekalaa.com/saudi-news/" rel="nofollow">1437</a><br /><a title="مصر" href="http://www.news.wekalaa.com/egypt-news/" rel="nofollow">مصر</a><br /><a title="منوعات" href="http://www.news.wekalaa.com/varieties/" rel="nofollow">منوعات</a>Khaled Solimanhttps://www.blogger.com/profile/06977474179483345994noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-9434823419498459032011-02-21T08:08:29.000-08:002011-02-21T08:08:29.000-08:00I remember that once I tried to made a new calcula...I remember that once I tried to made a new calculator for times of the employes to <a href="http://www.iservepharmacy.com" rel="nofollow">Viagra Online</a> labs and the source code it's very simple you don't need to made a lot of lines for that.Dakurohttps://www.blogger.com/profile/16460507435784917183noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-34925171076758368732010-09-15T13:49:09.221-07:002010-09-15T13:49:09.221-07:00what a good method and I wondered the same Yarosla...what a good method and I wondered the same Yaroslav, after reading this great post I asked myself if this method could also be applied to finding sums of random variables. If it could, it would be awesome!Sildenafilhttp://www.citratesildenafil.com/noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-18307362809461417422009-12-01T00:34:55.957-08:002009-12-01T00:34:55.957-08:00Oes Tsetnoc one of the ways in which we can learn ...<a rel="external" href="http://orientinspiration.com/2009/09/oes-tsetnoc.html" rel="nofollow">Oes Tsetnoc</a> one of the ways in which we can learn seo besides <a rel="external" href="http://duniasoer.com/archives/mengembalikan-jati-diri-bangsa.html" rel="nofollow">Mengembalikan Jati Diri Bangsa</a>. By participating in the <a rel="external" href="http://mncmakina.com/2009/09/oes-tsetnoc-seo-contest/" rel="nofollow">Oes Tsetnoc</a> or <a rel="external" href="http://ulaznice.info/mengembalikan-jati-diri-bangsa.html" rel="nofollow">Mengembalikan Jati Diri Bangsa</a> we can improve our seo skills. To find more information about <a rel="external" href="http://ulaznice.info/oes-tsetnoc.html" rel="nofollow">Oest Tsetnoc</a> please visit my <a rel="external" href="http://orientinspiration.com/2009/11/indonesian-bloggers-in-oes-tsetnoc-contest.html" rel="nofollow">Oes Tsetnoc</a> pages. And to find more information about <a rel="external" href="http://duniasoer.com" rel="nofollow">Mengembalikan Jati Diri Bangsa</a> please visit my <a rel="external" href="http://catatanblogger.info/464655-Mengembalikan-Jati-Diri-Bangsa.html" rel="nofollow">Mengembalikan Jati Diri Bangsa</a> pages. Thank you So much.<br> <a rel="external" href="http://orientinspiration.com" rel="nofollow">Oes Tsetnoc</a> | <a rel="external" href="http://duniasoer.com/archives/semangat-mengembalikan-jati-diri-bangsa.html" rel="nofollow">Semangat Mengembalikan Jati Diri Bangsa</a>Adihttps://www.blogger.com/profile/06355802787660974633noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-80987115778012600872009-11-17T08:21:42.263-08:002009-11-17T08:21:42.263-08:00Hi
I tried to add from Numeric import *
but it st...Hi<br /><br />I tried to add from Numeric import *<br />but it still complains that it can't find FFT.alekdarnoreply@blogger.comtag:blogger.com,1999:blog-9318222.post-90522871492319507052009-11-12T17:42:16.343-08:002009-11-12T17:42:16.343-08:00Find Internet Marketing resource hare Online Marke...Find <a rel="external" href="http://mncmakina.com" rel="nofollow">Internet Marketing</a> resource hare <a rel="external" href="http://mncmakina.com/2009/06/strategy-for-online-marketing-without-any-costs/" rel="nofollow">Online Marketing Strategy</a> | <a rel="external" href="http://mncmakina.com/2009/04/internet-marketing-tools/" rel="nofollow">Internet Marketing Tools</a> | <a rel="external" href="http://mncmakina.com/2009/06/dealing-with-unique-article-to-win-the-competition-of-online-marketing-campaign/" rel="nofollow">Online Marketing Campaign</a> | <a rel="external" href="http://mncmakina.com" rel="nofollow">Online Marketing Business</a> | <a rel="external" href="http://mncmakina.com/2009/04/finding-online-marketing-system-that-works/" rel="nofollow">Online Marketing System</a> | <a rel="external" href="http://mncmakina.com" rel="nofollow">Online Business</a> | <a rel="external" href="http://mncmakina.com/2009/06/tips-for-online-business-at-home/" rel="nofollow">Online Home Business</a> | <a rel="external" href="http://mncmakina.com/2009/06/tips-for-online-business-at-home/" rel="nofollow">Online Business Tips</a> | <a rel="external" href="http://mncmakina.com" rel="nofollow">Internet Marketing Online</a>Adihttps://www.blogger.com/profile/06355802787660974633noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-50597567832174595722008-03-23T05:33:00.000-07:002008-03-23T05:33:00.000-07:00FFT approach for finding convolution is indeed eff...FFT approach for finding convolution is indeed efficient but natural questions arises: what can be said about the error of this approximation. This is of course an approximation due to the representation of the numbers in computer arithmetic (in theory it should be exact). Since computation of the convolution using FFT is fast many you may use very long sequences but what about the round-off errors then? Do you know some sources in which one can find something about the error?Lukaszhttps://www.blogger.com/profile/18172638984641373036noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-1132016243315174972005-11-14T16:57:00.000-08:002005-11-14T16:57:00.000-08:00You can also take polynomial square roots and mult...You can also take polynomial square roots and multiply integers in this way. The NTT does it in exact integer arithmetic, or you can round the floating point results and rely on error bounds proved for IEEE-style floating point arithmetic.Connellyhttps://www.blogger.com/profile/02568908952592933174noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-1131937518295433272005-11-13T19:05:00.000-08:002005-11-13T19:05:00.000-08:00You're right about the IID distribution convolutio...You're right about the IID distribution convolutions. You can compute your distribution convolution with the fft_mul() code I gave.<BR/><BR/>>>> fft_mul([0.6, 0.3, 0.1], [0.8, 0.2])<BR/>[0.48, 0.36, 0.14, 0.02]<BR/><BR/>The reason it works is because the polynomial coefficients of c(x) = a(x)*b(x) are:<BR/><BR/>c[0] = a[0]*b[0]<BR/>c[1] = a[0]*b[1]+a[1]*b[0]<BR/>c[2] = a[0]*b[2]+a[1]*b[1]+a[2]*b[0]<BR/>...<BR/><BR/>That is, c[k] is the sum of all a[i] and b[j] such that i + j = k. Which is just P(X1+X2=k).<BR/><BR/>---<BR/><BR/>For integral transforms with an arbitrary kernel, you can just follow the proof of the <A HREF="http://mathworld.wolfram.com/ConvolutionTheorem.html" REL="nofollow">convolution theorem</A> to establish a sufficient condition on the kernel such that a star b = T(a)T(b).<BR/><BR/>I recall from Math 256 that analytic inverse Laplace transforms are usually impossible if you can't find them in some CRC type handbook. But apparently they can be done numerically.Connellyhttps://www.blogger.com/profile/02568908952592933174noreply@blogger.comtag:blogger.com,1999:blog-9318222.post-1131931513174986272005-11-13T17:25:00.000-08:002005-11-13T17:25:00.000-08:00Interesting post,I wonder if this method could als...Interesting post,<BR/><BR/>I wonder if this method could also be applied to finding sums of random variables.<BR/><BR/>For instance to get this picture(http://web.engr.oregonstate.edu/~bulatov/research/reports/vapnik/plot2.png), I had to get distributions of X1, X1+X2, X1+X2+X3+...+Xl Where X_i are IID variables.<BR/><BR/>You find this distribution using convolution. If we have P(X=k), then P(X1+X2)=k is <BR/><BR/>P(X=0)P(X=k)+P(X=1)P(X=k-1)+....<BR/><BR/>Which is really a convolution. So if a[k] gives P(X=k), then convolve(a,a)[k] will P(X1+X2)=k; def convolve(a,b)=[sum([a[j]*b[i-j] for j in range(0,i+1)]) for i in range(len(a))]<BR/><BR/>You can repeat the convolution l-1 times to get the sum of X1+...+Xl<BR/><BR/>The problem is that this gets really slow.<BR/><BR/>----<BR/><BR/>In statistics there's method for finding sums called "Method of Moments". They take two-sided Laplace transform of the probability density (which is called the Moment Generating Function http://en.wikipedia.org/wiki/Moment_generating_function),<BR/>and then convolutions of the densities become products.<BR/><BR/>So apparently Fourier transform isn't the only one where convolutions become products...<BR/><BR/>Is there some general theorem which says that convolutions turn into products for X,Y and Z transforms?<BR/><BR/>Curiously, in statistics books the inverse of the transform is done by "guess and check"...is the inverse transform really that hard?<BR/><BR/>Interestingly some proofs of the Central Limit Theorem operate the transformed space. IE, it's hard to see what'll happen if you take X1+X2+X3+....\infty if you viewing it as convolutions (ie, how exactly are you going to do the infinite dimensional integral?!?), but instead if you look at it as infinite product, it's easier. So it turns out that the infinite product converges to the same thing regardless of the thing you are multiplying, which could explain why people were so excited about Central Limit Theorem (http://www-stat.stanford.edu/~susan/courses/s116/node120.html)Yaroslavhttps://www.blogger.com/profile/17683284937008290386noreply@blogger.com