Candidates interested in pursuing a PhD and passionate about working in the area of computational metal plasticity may send me an email

Useful tools

Gnuplot

A very useful opensource plotting tool with lot of control.

Gnuplot script to obtain the COVID-19 update in India

The script downloads the data from covid19india.org's repository and plots the daily cases till now. Cases for individual cities can also be plotted. You need to have cygwin installed if you are using windows.

Kolkata.png

Delhi.png

Howrah.png

Here is the code to generate this plot

if(dwn==1){system "wget https://api.covid19india.org/csv/latest/case_time_series.csv -O ./case_time_series.csv"} resetf(x)=a*exp(-0.5*((x-b)**2)/c**2)a=100;b=40;c=20;array A[10];A[2]="Daily Confirmed";A[4]="Daily Recovered";A[6]="Daily Deceased";set datafile separator ","set grid ytics front lw 3set boxwidth set isosamples 50set style fill solid 0.5 border -1set term pngcairo size 1920,1080 enhanced font "Arial Narrow,18"set xlabel "Days"set ylabel "New Cases/Day"set output 'India_'.A[coln].'.png';print 'Processing 'stats'case_time_series.csv' u 1:coln nooutput if(fityes==1){fit f(x) 'case_time_series.csv' u (strptime("%d %B",strcol(1))/3600/24-strptime("%d %B","31 January")/3600/24):coln via a,b,c;g(x)=f((x-strptime("%d %m %Y","31 01 2020"))/3600/24);totproj = 1./2*(sqrt(2.)*sqrt(pi)*c*erf(1./2*sqrt(2.)*(b + 10000.)/c) - sqrt(2.)*sqrt(pi)*c*erf(1./2*sqrt(2.)*(b - 10000.)/c))*a}else {totproj = 0};set xdata time;set timefmt "%d-%m-%Y"set format x "%b %d"set xtics format "%b %d"if(fityes == 0){maxy=ceil(STATS_max_y/100.)*100;set yrange[0:maxy];}set xrange[strptime("%d-%m-%Y",startday):strptime("%d-%m-%Y",endday)]set xtics startday,dayint*24*3600,enddayset mxtics dayintset key left topset label "{/:Bold Total till today : }".int(STATS_sum_y) tc rgb "blue" at graph 0.45,0.95 front;set label "{/:Bold Projected Total cases : }".int(totproj) tc rgb "red" at graph 0.02,0.65 front left;p 'case_time_series.csv' u (strptime("%d %B",strcol(1)) - strptime("%d %m %Y","31 01 1970") + strptime("%d %m %Y","31 01 2020")):coln w boxes lc coln lw 1.5 t "India ".columnhead(coln), 'case_time_series.csv' u (strptime("%d %B",strcol(1)) - strptime("%d %m %Y","31 01 1970") + strptime("%d %m %Y","31 01 2020")):coln every ::STATS_records-1::STATS_records w boxes lc coln+2 lw 1.5 notitle, '' u (strptime("%d %B",strcol(1)) - strptime("%d %m %Y","31 01 1970") + strptime("%d %m %Y","31 01 2020")):coln:(strcol(coln)) every ::STATS_records::STATS_records w labels boxed offset 0.0,0.5 font "Arial Narrow-Bold,15" tc rgb "red" notitle,(fityes==1) ? g(x) :1/0 lc rgb "black" lw 4 t "Fit";set output;reset

Use the following command to run the script from terminal -

> gnuplot -e "fityes=1;dwn=1;dayint=30;coln=2;endday=31-12-2020"

Asymptote

It is a very useful vector graphics language which comes bundled with the texlive package. One can create very useful 2D or 3D interactive vector graphics using asymptote. It can be used as a standalone to generate a PDF of the vector graphics, or can be used within a latex code to embed the graphics in a latex document. The standalone asymptote files can be saved with an extension ".asy" and can be compiled using the following command -

> asy filename.asy

An example of such vector graphics created with asymptote is shown below.

Please download the file and open it with Adobe Acrobat Reader to view the content.

indent.pdf

The code snippet for the above vector graphics: (get the actual file here)

Octave

An open source equivalent of Matlab. It's available here - https://www.gnu.org/software/octave/

Page updated

Google Sites

Report abuse