Skip to main content

Mode locking and numerical methods in fortran

I made a visualization that shows how mode locking creates short pulses; the more the number of modes in the cavity, the shorter the pulse width, with the condition that they all have the same phase. You can find it here. It's not yet complete, I also need to show how the output pulse looks when the modes have random phases. The code is in python, was written using the ipython notebook and uses ipywidgets to create the interactive plot.

On the flip side, I was hitting my head against the wall trying to adapt a code to perform numerical integration in fortran. The code would compile but not execute in the first place and google told me to move the executable to a folder in the home directory. Weirdly enough, this worked. Apparently this problem is an after effect from when i was copying these files out of my home dir earlier. chmod helped fix that problem. Having figured out how to define and call functions in fortran, I was now stuck on a type mismatch error, because the function was defined to take a certain type of input and was now being fed a different one. Adding implicit none dealt with that and just when I was happy that i could execute the file, a wild segmentation fault appeared on screen upon execution of said file. Instead of learning on the go, I should probably look up what the fortran standard is and what basic rules are to be followed. I was surprised to find that using a ! to comment files was acceptable on where as using c would throw up an error, only on gfortran though. I'm sure there are more rules that I haven't bothered with so far, that I should look into. For now, i sign off, wondering why there isn't a way to batch comment lines in fortran where as ''' ... ''' does the job in python and /*...*/ does the job in C, which themselves use # and // as comments otherwise. Apparently // is also not acceptable at times in C. Rules, Rules, Rules...

Popular posts from this blog

Animation using GNUPlot

Animation using GNUPlotI've been trying to create an animation depicting a quasar spectrum moving across the 5 SDSS pass bands with respect to redshift. It is important to visualise what emission lines are moving in and out of bands to be able to understand the color-redshift plots and the changes in it.
I've tried doing this using the animate function in matplotlib, python but i wasn't able to make it work - meaning i worked on it for a couple of days and then i gave up, not having found solutions for my problems on the internet.
And then i came across this site, where the gunn-peterson trough and the lyman alpha forest have been depicted - in a beautiful manner. And this got me interested in using js and d3 to do the animations and make it dynamic - using sliders etc.
In the meanwhile, i thought i'd look up and see if there was a way to create animations in gnuplot and whoopdedoo, what do i find but nirvana!

In the image, you see 5 static curves and one dynam…

on MOOCs.

For those of you who don't know, MOOC stands for Massively Open Online Course.

The internet is an awesome thing. It's making education free for all. Well, mostly free. But it's surprising at the width and depth of courses being offered online. And it looks like they are also having an impact on students, especially those from universities that are not top ranked. Students in all parts of the world can now get a first class education experience, thanks to courses offered by Stanford, MIT, Caltech, etc.

I'm talking about MOOCs because one of my new year resolutions is to take online courses, atleast 2 per semester (6 months). And I've chosen the following two courses on edX - Analyzing Big Data with Microsoft R Server and Data Science Essentials for now. I looked at courses on Coursera but I couldn't find any which was worthy and free. There are a lot more MOOC providers out there but let's start here. And I feel like the two courses are relevant to where I …

Pandas download statistics, PyPI and Google BigQuery - Daily downloads and downloads by latest version

Inspired by this blog post :, I wanted to play around with Google BigQuery myself. And the blog post is pretty awesome because it has sample queries. I mix and matched the examples mentioned on the blog post, intent on answering two questions - 
1. How many people download the Pandas library on a daily basis? Actually, if you think about it, it's more of a question of how many times was the pandas library downloaded in a single day, because the same person could've downloaded multiple times. Or a bot could've.
This was just a fun first query/question.
2. What is the adoption rate of different versions of the Pandas library? You might have come across similar graphs which show the adoption rate of various versions of Windows.
Answering this question is actually important because the developers should have an idea of what the most popular versions are, see whether or not users are adopting new features/changes they provide…