Skip to main content

Pairwise Summation algorithm implementation with python

Pairwise summation algorithm is a summation algorithm like Kahan Algorithm but does perform better when compared to the Kahan Algorithm. This is because of the Divide and Conquer approach followed by this algorithm. According to Wikipedia, this is the official summation algorithm used in Numpy.


If you don't know what divide and conquer is all about, then you should see the following video: lecture 3: Divide and Conquer by MITOCW.

To be simple, this algorithm takes a list of numbers whose length is at least greater than 0 and also a sufficiently small number N < n(length of the list).

Now the list is broken down into sublists of length N, we find the sum of these sublists and finally add all the summations together to get the result. As simple as that.

Up next is the program for this algorithm which I have written using python. But if you want further reference then you can see the Wikipedia pairwise summation algorithm. Remember that you can only understand the algorithm better when you will take an example. For example, you can plug in [1, 2, 3] as the list. Take a paper and write down the output of each and every statement.

Program

Program is self explanative and I don't think this one needs a detailed explanation.

Let us test if the algorithm is performing as expected. I am just testing if there is precision or not. You can also use timeit to check the speed of the algorithm
>>> aList = [0.1]*10
>>> aList
[0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1]
>>> pairwise(aList)
1.0
>>> pairwise(aList) == 1.0
True
>>> 
As you can see the algorithm performs as expected.  But remember that the condition N < n is taken for granted. If you change the value of N to be less than the length of the list then you will get an error.

You can replace the lines 10-12 with sum function.

But guys I don't know why, this algorithm is giving a wrong result when the input is [0.1, 0.1, 0.1]. Remaining results are satisfactory, at-least the test cases which I have done! If you are aware of why this is happening, then please, please comment in the comment box.

If you have any doubt or haven't understood anything, then comment in the comment box below and I will be glad to help you.

If you have any suggestions or feedback or if you have found a typo then please comment in the comment box below and I will be very happy to see your post.

You can also contact me.

Thank you. Have a nice day😃.

Popular posts from this blog

Making a quiz web app with python and flask

Edit : When you are creating a web app with h tml templates, then y ou will have to sa ve the html file in templates folder in the Current Wor ki ng Directory( CWD). If you save the file in the C W D directl y you will get a TemplateNotFound error. Thank you Udhay for pointing it out.   In this post we will create a quiz website using python . I will be using the flask framework . After reading this tutorial you will learn form submission , flask templates , python code in flask templates , shuffling the questions and options with the random module and few others.  Please note that this tutorial is not big as it seems to be. Some of the code has been rewritten to maintain consistency and also font size is somewhat big so that your eyes won't get stressed reading this tutorial. Also the content has not occupied the full width of the page. In this tutorial I am assuming that you are having a very basic understanding of the flask framework . Please refer the documentation

Problem 11 Project Euler Solution with python

Largest product in a grid In the 20×20 grid below, four numbers along a diagonal line have been marked in red. 08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08 49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00 81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65 52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91 22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80 24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50 32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70 67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21 24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72 21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95 78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92 16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57 86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58 19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40 04 52 08 83 97 35 99 16 07

Problem 60 Project Euler Solution with python

Prime pair sets The primes 3, 7, 109, and 673, are quite remarkable. By taking any two primes and concatenating them in any order the result will always be prime. For example, taking 7 and 109, both 7109 and 1097 are prime. The sum of these four primes, 792, represents the lowest sum for a set of four primes with this property. Find the lowest sum for a set of five primes for which any two primes concatenate to produce another prime. This problem is j u st a brute force problem. If you have come here because you don't know the limit upto which you will h ave to gener ate the prime numbers t hen go ahe ad and t r y with 10,000 . When I first start ed solving the problem I chose 1 million(beca use most of the problem s on project E uler have this limit ), but it took very long for the computer to fin d the solution. After searching on the internet then I found many people choosing 10, 000 so I have changed my in put f rom 1 million to 10000 and the output was f ast. He