-
Notifications
You must be signed in to change notification settings - Fork 2
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Trying to get more coding examples into repos.
- Loading branch information
Srinath Vadlamani
committed
Sep 1, 2011
1 parent
0423ff4
commit 2e7209e
Showing
8 changed files
with
292 additions
and
0 deletions.
There are no files selected for viewing
Binary file not shown.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,9 @@ | ||
|
||
|
||
PROGRAM hello | ||
IMPLICIT NONE | ||
|
||
WRITE(*,*) "HELLO WORLD" | ||
|
||
|
||
END PROGRAM hello |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -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() | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -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) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -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() | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -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") |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,7 @@ | ||
#! /opt/local/bin/python | ||
|
||
import newClass | ||
|
||
sp=newClass.simplePrint() | ||
sp.what("hello") | ||
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -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 <http://www.gnu.org/licenses/>. | ||
|
||
# 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: | ||
|