From 2e7209e2de02a8d28706233d2629f9a6a307e3f1 Mon Sep 17 00:00:00 2001 From: Srinath Vadlamani Date: Wed, 31 Aug 2011 18:44:13 -0600 Subject: [PATCH] Trying to get more coding examples into repos. --- codingTests/fortran90/hello | Bin 0 -> 9120 bytes codingTests/fortran90/hello.f90 | 9 ++ codingTests/matplotlibing/griddata.py | 44 +++++++ codingTests/matplotlibing/interp2d.py | 38 ++++++ codingTests/matplotlibing/surfacePlot.py | 28 +++++ pythonTools/ffts/simplefft.py | 15 +++ pythonTools/testNewClass.py | 7 ++ scripting/git-version-gen.sh | 151 +++++++++++++++++++++++ 8 files changed, 292 insertions(+) create mode 100755 codingTests/fortran90/hello create mode 100644 codingTests/fortran90/hello.f90 create mode 100644 codingTests/matplotlibing/griddata.py create mode 100644 codingTests/matplotlibing/interp2d.py create mode 100644 codingTests/matplotlibing/surfacePlot.py create mode 100644 pythonTools/ffts/simplefft.py create mode 100644 pythonTools/testNewClass.py create mode 100755 scripting/git-version-gen.sh diff --git a/codingTests/fortran90/hello b/codingTests/fortran90/hello new file mode 100755 index 0000000000000000000000000000000000000000..40afee08c6cfa2a73fa0b47231d62620d1814261 GIT binary patch literal 9120 zcmeHNU1(HC6rN31P5nu>C7{+;S25J~A&K>;@u9FQrruy0qls1uW!&A|Y_@K8m%De< zD1@ejlq?~T2PslY^}&a}v`?kjzO>c{u`i`i+CCNO4@9a|L_|cq{l0r=+`Tu6eJXvq z2WHMWbLO0x`DS+a-Z}gIgIkYQ2@z=&VzNO9(Tq}BDa155X(Nh(B4vN_VE3EdM-Q{r zPZ7Vmq{cYEpre$7-6sb9U0nZ|zfC9|a}z66VGk*7bKIunQ83=cbyy5ctJS2ab6t-q zDkdMx5Tq>Im2~-WQ!w6V0ED5os~%d{4oxZ^Aalu!Cr8YD$+XZCjQ698*Q@$dfAah? z#X`BZ;9$IXi#NZ7ibtJzEk{zCBXZbEjhbGMV7!wm-i!)B9cN=%;=axRASn+d2a|Q( zuz&6r=zR@v?knR)m5Y=^r}INK|K;L!=`M)(q5M1xLiWgU1sMp7W|yT>_n3>SBX_0 zOnn=<{9B_Dxnop+_$z8mU2Za__KC_pN(*NJTJ`F}EwH=Ef!?|2*WYoRbhxYD&NPTN_`(`oy_tI^>y$9G)5oJd!io275OwMsgkuyg{ zWIPs$tzFeL55sP>J%@4zu#)($=u(tk&9V0zqI^PMw=y~k1}E3(`3*w;T>0%+?P8JN zt3E&7Z4r-sBb|5gZlXBtoPmD-n8bSr^%=$I6=%PC((kpus`uLk#qTKoQ1O-O{j%OS zltTnU1VRKt1VRKt1VRKt1VRKt1VRKt1VRKt1pZqDS`$p7xk+%D&J~7`v_4{W%?rzT zIw7UbD#@~G%aoNZ3t0LWS&$RTAxll;3BS>n?^rq8q~+w27CBTbxaJX`nW?PuuxZK6 zNXklOY)lIwO05ZLj^}L1qCtwaWjkfr;*f|+@xo{-R}lSuCz4iH95V}JIjdL@SU z&|xz!M6Tn11-utw-jVuUrME+GhWE9r*wVK!kaSM^*u`( zdSe8aE+CFKhC;i(kHOWAXn`TG?pLOjuJequ&|#*lOXpG7>t6))?*jVo0X^b5Go9+v zs({XWK<56id^eAaW$V9kBp@c1v&V)19OepH{U-|VuK{wnQpm`Bv3ROd67uzAPoI=R zs*{K0<<333_6q6le@5{!WZ7$bE#=~)AoW3!?(y-?(^AJi(54T3eX~C5OSE~1yw0%1 N@%9e39|RoMe*iBNAw~cI literal 0 HcmV?d00001 diff --git a/codingTests/fortran90/hello.f90 b/codingTests/fortran90/hello.f90 new file mode 100644 index 0000000..f9c7e9e --- /dev/null +++ b/codingTests/fortran90/hello.f90 @@ -0,0 +1,9 @@ + + + PROGRAM hello + IMPLICIT NONE + + WRITE(*,*) "HELLO WORLD" + + + END PROGRAM hello diff --git a/codingTests/matplotlibing/griddata.py b/codingTests/matplotlibing/griddata.py new file mode 100644 index 0000000..f63cd42 --- /dev/null +++ b/codingTests/matplotlibing/griddata.py @@ -0,0 +1,44 @@ +import numpy as np +# Suppose we want to interpolate the 2-D function + +def func(x, y): + return x*(1-x)*np.cos(4*np.pi*x) * np.sin(4*np.pi*y**2)**2 + +# on a grid in [0, 1]x[0, 1] + +grid_x, grid_y = np.mgrid[0:1:100j, 0:1:200j] + +# but we only know its values at 1000 data points: + +points = np.random.rand(1000, 2) +values = func(points[:,0], points[:,1]) + +# This can be done with `griddata` -- below we try out all of the +# interpolation methods: + +from scipy.interpolate import griddata +grid_z0 = griddata(points, values, (grid_x, grid_y), method='nearest') +grid_z1 = griddata(points, values, (grid_x, grid_y), method='linear') +grid_z2 = griddata(points, values, (grid_x, grid_y), method='cubic') + +# One can see that the exact result is reproduced by all of the +# methods to some degree, but for this smooth function the piecewise +# cubic interpolant gives the best results: + +import matplotlib.pyplot as plt +plt.subplot(221) +plt.imshow(func(grid_x, grid_y).T, extent=(0,1,0,1), origin='lower') +plt.plot(points[:,0], points[:,1], 'k.', ms=1) +plt.title('Original') +plt.subplot(222) +plt.imshow(grid_z0.T, extent=(0,1,0,1), origin='lower') +plt.title('Nearest') +plt.subplot(223) +plt.imshow(grid_z1.T, extent=(0,1,0,1), origin='lower') +plt.title('Linear') +plt.subplot(224) +plt.imshow(grid_z2.T, extent=(0,1,0,1), origin='lower') +plt.title('Cubic') +plt.gcf().set_size_inches(6, 6) +plt.show() + diff --git a/codingTests/matplotlibing/interp2d.py b/codingTests/matplotlibing/interp2d.py new file mode 100644 index 0000000..468abfc --- /dev/null +++ b/codingTests/matplotlibing/interp2d.py @@ -0,0 +1,38 @@ + + +# Example 2-d interpolation with interp2d + + +import numpy as np +from scipy import interpolate +import pylab as py + +def func(x,y): + return (x*y)*np.exp(-5.0*(x**2+y**2)) + +x,y = np.mgrid[-1:1:15j,-1:1:15j] +fvals=func(x,y) + +newfunc= interpolate.interp2d(x,y,fvals, kind='cubic') + + +#Notice that we evaluate on the new 1-d arrays +xnew = np.linspace(-1,1,100) +ynew=xnew + +fnew=newfunc(xnew,ynew) + + +py.figure(1) +py.clf() +py.imshow(fnew, extent=[-1,1,-1,1],cmap=py.cm.jet) + + +from enthought.mayavi import mlab + +#Need fully-flesehd out gride for surface + +x2,y2=np.mgrid[-1:1:100j,-1:1:100j] + +mlab.clf() +mlab.surf(x2,y2,fnew*2) diff --git a/codingTests/matplotlibing/surfacePlot.py b/codingTests/matplotlibing/surfacePlot.py new file mode 100644 index 0000000..2a7ee85 --- /dev/null +++ b/codingTests/matplotlibing/surfacePlot.py @@ -0,0 +1,28 @@ +# By Armin Moser + +from mpl_toolkits.mplot3d import Axes3D +import matplotlib +import numpy as np +from matplotlib import cm +from matplotlib import pyplot as plt +step = 0.04 +maxval = 1.0 +fig = plt.figure() +ax = fig.add_subplot(111, projection='3d') + +# create supporting points in polar coordinates +r = np.linspace(0,1.25,50) +p = np.linspace(0,2*np.pi,50) +R,P = np.meshgrid(r,p) +# transform them to cartesian system +X,Y = R*np.cos(P),R*np.sin(P) + +Z = ((R**2 - 1)**2) +ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=cm.jet) +ax.set_zlim3d(0, 1) +ax.set_xlabel(r'$\phi_\mathrm{real}$') +ax.set_ylabel(r'$\phi_\mathrm{im}$') +ax.set_zlabel(r'$V(\phi)$') +ax.set_xticks([]) +plt.show() + diff --git a/pythonTools/ffts/simplefft.py b/pythonTools/ffts/simplefft.py new file mode 100644 index 0000000..a6c6048 --- /dev/null +++ b/pythonTools/ffts/simplefft.py @@ -0,0 +1,15 @@ +from scipy import * +a = zeros(1000) +b = fft(a) +plot(abs(b)) +from pylab import * +plot(abs(b)) +b +a[:100]=1 +a +b = fft(a) +plot(abs(b)) +_ip.magic("history ") +_ip.magic("history -help") +#?history +_ip.magic("history -n -f simplefft.py") diff --git a/pythonTools/testNewClass.py b/pythonTools/testNewClass.py new file mode 100644 index 0000000..40b423c --- /dev/null +++ b/pythonTools/testNewClass.py @@ -0,0 +1,7 @@ +#! /opt/local/bin/python + +import newClass + +sp=newClass.simplePrint() +sp.what("hello") + diff --git a/scripting/git-version-gen.sh b/scripting/git-version-gen.sh new file mode 100755 index 0000000..67d9910 --- /dev/null +++ b/scripting/git-version-gen.sh @@ -0,0 +1,151 @@ +#!/bin/sh +# Print a version string. +scriptversion=2009-09-09.22 + +# Copyright (C) 2007-2009 Free Software Foundation +# +# This program is free software: you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 3 of the License, or +# (at your option) any later version. +# +# This program is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program. If not, see . + +# This script is derived from GIT-VERSION-GEN from GIT: http://git.or.cz/. +# It may be run two ways: +# - from a git repository in which the "git describe" command below +# produces useful output (thus requiring at least one signed tag) +# - from a non-git-repo directory containing a .tarball-version file, which +# presumes this script is invoked like "./git-version-gen .tarball-version". + +# In order to use intra-version strings in your project, you will need two +# separate generated version string files: +# +# .tarball-version - present only in a distribution tarball, and not in +# a checked-out repository. Created with contents that were learned at +# the last time autoconf was run, and used by git-version-gen. Must not +# be present in either $(srcdir) or $(builddir) for git-version-gen to +# give accurate answers during normal development with a checked out tree, +# but must be present in a tarball when there is no version control system. +# Therefore, it cannot be used in any dependencies. GNUmakefile has +# hooks to force a reconfigure at distribution time to get the value +# correct, without penalizing normal development with extra reconfigures. +# +# .version - present in a checked-out repository and in a distribution +# tarball. Usable in dependencies, particularly for files that don't +# want to depend on config.h but do want to track version changes. +# Delete this file prior to any autoconf run where you want to rebuild +# files to pick up a version string change; and leave it stale to +# minimize rebuild time after unrelated changes to configure sources. +# +# It is probably wise to add these two files to .gitignore, so that you +# don't accidentally commit either generated file. +# +# Use the following line in your configure.ac, so that $(VERSION) will +# automatically be up-to-date each time configure is run (and note that +# since configure.ac no longer includes a version string, Makefile rules +# should not depend on configure.ac for version updates). +# +# AC_INIT([GNU project], +# m4_esyscmd([build-aux/git-version-gen .tarball-version]), +# [bug-project@example]) +# +# Then use the following lines in your Makefile.am, so that .version +# will be present for dependencies, and so that .tarball-version will +# exist in distribution tarballs. +# +# BUILT_SOURCES = $(top_srcdir)/.version +# $(top_srcdir)/.version: +# echo $(VERSION) > $@-t && mv $@-t $@ +# dist-hook: +# echo $(VERSION) > $(distdir)/.tarball-version + +case $# in + 1) ;; + *) echo 1>&2 "Usage: $0 \$srcdir/.tarball-version"; exit 1;; +esac + +tarball_version_file=$1 +nl=' +' + +# First see if there is a tarball-only version file. +# then try "git describe", then default. +if test -f $tarball_version_file +then + v=`cat $tarball_version_file` || exit 1 + case $v in + *$nl*) v= ;; # reject multi-line output + [0-9]*) ;; + *) v= ;; + esac + test -z "$v" \ + && echo "$0: WARNING: $tarball_version_file seems to be damaged" 1>&2 +fi + +if test -n "$v" +then + : # use $v +elif test -d .git \ + && v=`git describe --abbrev=4 --match='v*' HEAD 2>/dev/null \ + || git describe --abbrev=4 HEAD 2>/dev/null` \ + && case $v in + v[0-9]*) ;; + *) (exit 1) ;; + esac +then + # Is this a new git that lists number of commits since the last + # tag or the previous older version that did not? + # Newer: v6.10-77-g0f8faeb + # Older: v6.10-g0f8faeb + case $v in + *-*-*) : git describe is okay three part flavor ;; + *-*) + : git describe is older two part flavor + # Recreate the number of commits and rewrite such that the + # result is the same as if we were using the newer version + # of git describe. + vtag=`echo "$v" | sed 's/-.*//'` + numcommits=`git rev-list "$vtag"..HEAD | wc -l` + v=`echo "$v" | sed "s/\(.*\)-\(.*\)/\1-$numcommits-\2/"`; + ;; + esac + + # Change the first '-' to a '.', so version-comparing tools work properly. + # Remove the "g" in git describe's output string, to save a byte. + v=`echo "$v" | sed 's/-/./;s/\(.*\)-g/\1-/'`; +else + v=UNKNOWN +fi + +v=`echo "$v" |sed 's/^v//'` + +# Don't declare a version "dirty" merely because a time stamp has changed. +git status > /dev/null 2>&1 + +dirty=`sh -c 'git diff-index --name-only HEAD' 2>/dev/null` || dirty= +case "$dirty" in + '') ;; + *) # Append the suffix only if there isn't one already. + case $v in + *-dirty) ;; + *) v="$v-dirty" ;; + esac ;; +esac + +# Omit the trailing newline, so that m4_esyscmd can use the result directly. +echo "$v" | tr -d '\012' + +# Local variables: +# eval: (add-hook 'write-file-hooks 'time-stamp) +# time-stamp-start: "scriptversion=" +# time-stamp-format: "%:y-%02m-%02d.%02H" +# time-stamp-end: "$" +# End: +