monte carlo integration and estimating the value of pi

I have talked about estimating the value of pi earlier and have even mentioned the code that I wrote to estimate the value of pi. I did not realize that I was applying the monte carlo integration method to estimate the area of the circle. I have wanted to implement it since I started working with monte carlo methods and once I realized that the method I used to estimate pi was a monte carlo method, I immediately tried to estimate the error in the estimation, which depends on the total number of monte carlo steps or iterations. Error in monte carlo integration is a function of the inverse square root of the total number of monte carlo steps, which is exactly what i saw in the log-log plot of step size vs error in estimation. After a bit of trouble, i ported the code to fortran and got the same results, with fitting and plotting being done on gnuplot. The ipynb file is available here and the fortran file here. I now need to look for more interesting functions which I can apply the monte carlo integration method to. From what I've been reading so far, it's conducive to use monte carlo integration when one is dealing with multi dimensional integration i.e 5 or more dimensions. Atleast I have something to start from now.

On the other hand, as I mentioned yesterday, I got my first segfault on fortran. I tried debugging it with gdb but couldn't understand much. I do know that the code worked perfectly fine until I started defining subroutines to perform the function evaluation and to perform numerical integration. Either way, I need to get these two sorted tomorrow. It's been a good day, let's see if I maintain the same energy tomorrow...

Popular posts from this blog

Arxiv author affiliations using Python

So, I wanted to get author affiliation information from papers on arXiv. arXiv provides with an API to bulk query their database and get information. Following that, I look for the attribute 'arxiv:affiliation' in the html data. Here's the code - import urllib from BeautifulSoup import BeautifulStoneSoup   url = 'http://export.arxiv.org/api/query?search_query=all:astro&start=0&max_results=1000'   data = urllib.urlopen(url).read() soup = BeautifulStoneSoup(data) #print(soup.prettify()) #list = soup.findAll('arxiv:affiliation') #for i in range(len(list)): #        print list[i].contents   test = [tag.string for tag in soup.findAll('arxiv:aiffiliation')] Now, the problem I'm having is that I'm getting affiliation of all authors which I want to split into sets of affiliations of authors of a paper, which I'm stuck on at the moment. Once I get that part, I can move on to the next part of this pet project, displaying these
Read more

Farewell to Enthought

I started my professional career at Enthought , at the Pune office in India, on Feb 29, 2016. Over the course of eight long years, I spent time at the Austin office in the US, the Cambridge office in the UK, and I worked remotely from India. I was let go on March 27, 2024 . I had an unbelievably great time. I had three amazing mentors - Pankaj Pandey, Senganal Thirunavukkarasu, and Mark Dickinson. Pankaj mentored me in my first year at Enthought. I started paying attention to the craft of software after working with him. Senganal started mentoring me in my second year at Enthought, and continues to mentor me, even after the both of us have moved on from Enthought. He instilled professionalism in me, ensuring that the code I wrote produced precise scientific results, and that I tackled the known unknown aspects of a project first, to ensure timely communication about project progress. I was eager to manage people at the time and Senganal made me realize that I didn't want to manage
Read more

Elementary (particle physics), my dear Watson

Ground Control to Major Tom:  Shruti Patel checking in. Since my dear friend and fellow-physicist Rahul Poruri is homeward bound, and shall remain there for a week, I'm stepping in to continue his tradition of daily posts to Astronut . Even as I write, I'm not quite sure what this post is going to be about, so please forgive me for occasional rants and digressions.  By the way Rahul, I totally resent you for the name of this blog- I'm never going to start my own blog because there's no way I'll come up with a name as cool as Astronut! For now, I shall try and be happy with being a guest writer! Before I begin, a little bit on who I am and what I do; so that you can trust that my general views on science- particle physics in particular- are not entirely baseless. I graduated with a masters in Physics from IIT Madras in July this year, and I'm currently working towards (an eventual) PhD with the Theory Group at the Deutsches Electronen-Synchrotron (DESY)
Read more