Introduction to Python
Heavily based on presentations by Matt Huenerfauth (Penn State) Guido van Rossum (Google) Richard P. Muller (Caltech) ...
Monday, October 19, 2009
Python
Open source general-purpose language. Object Oriented, Procedural, Functional Easy to interface with C/ObjC/Java/Fortran Easy-ish to interface with C++ (via SWIG) Great interactive environment
Downloads: http://www.python.org Documentation: http://www.python.org/doc/ Free book: http://www.diveintopython.org
Monday, October 19, 2009
2.5.x / 2.6.x / 3.x ???
Current version is 2.6.x Mainstream version is 2.5.x The new kid on the block is 3.x
You probably want 2.5.x unless you are starting from scratch. Then maybe 3.x
Monday, October 19, 2009
Technical Issues
Installing & Running Python
Monday, October 19, 2009
Binaries
Python comes pre-installed with Mac OS X and Linux. Windows binaries from http://python.org/ You might not have to do anything!
Monday, October 19, 2009
The Python Interpreter
Interactive interface to Python
% python Python 2.5 (r25:51908, May 25 2007, 16:14:04) [GCC 4.1.2 20061115 (prerelease) (SUSE Linux)] on linux2 Type "help", "copyright", "credits" or "license" for more information. >>>
Python interpreter evaluates inputs:
>>> 3*(7+2) 27
Python prompts with >>>. To exit Python:
CTRL-D
Monday, October 19, 2009
Running Programs on UNIX
% python filename.py
You could make the *.py file executable and add the following #!/usr/bin/env python to the top to make it runnable.
Monday, October 19, 2009
Batteries Included
Large collection of proven modules included in the standard distribution.
http://docs.python.org/modindex.html
Monday, October 19, 2009
numpy
Offers Matlab-ish capabilities within Python Fast array operations 2D arrays, multi-D arrays, linear algebra etc.
Downloads: http://numpy.scipy.org/ Tutorial: http://www.scipy.org/ Tentative_NumPy_Tutorial
Monday, October 19, 2009
matplotlib
High quality plotting library.
#!/usr/bin/env python import numpy as np import matplotlib.mlab as mlab import matplotlib.pyplot as plt mu, sigma = 100, 15 x = mu + sigma*np.random.randn(10000) # the histogram of the data n, bins, patches = plt.hist(x, 50, normed=1, facecolor='green', alpha=0.75) # add a 'best fit' line y = mlab.normpdf( bins, mu, sigma) l = plt.plot(bins, y, 'r--', linewidth=1) plt.xlabel('Smarts') plt.ylabel('Probability') plt.title(r'$\mathrm{Histogram\ of\ IQ:}\ \mu=100,\ \sigma=15$') plt.axis([40, 160, 0, 0.03]) plt.grid(True) plt.show()
Downloads: http://matplotlib.sourceforge.net/
Monday, October 19, 2009
PyFITS
FITS I/O made simple:
>>> import pyfits >>> hdulist = pyfits.open(input.fits) >>> hdulist.info() Filename: test1.fits No. Name Type Cards Dimensions Format 0 PRIMARY PrimaryHDU 220 () Int16 1 SCI ImageHDU 61 (800, 800) Float32 2 SCI ImageHDU 61 (800, 800) Float32 3 SCI ImageHDU 61 (800, 800) Float32 4 SCI ImageHDU 61 (800, 800) Float32 >>> hdulist[0].header[targname] NGC121 >>> scidata = hdulist[1].data >>> scidata.shape (800, 800) >>> scidata.dtype.name float32 >>> scidata[30:40,10:20] = scidata[1,4] = 999
Downloads: http://www.stsci.edu/resources/ software_hardware/pyfits
Monday, October 19, 2009
pyds9 / python-sao
Interaction with DS9 Display Python 1-D and 2-D arrays in DS9 Display FITS files in DS9
Downloads: Ask Eric Mandel :-) Downloads: http://code.google.com/p/python-sao/
Monday, October 19, 2009
Wrappers for Astronomical Packages
CasaPy (Casa) PYGILDAS (GILDAS) ParselTongue (AIPS) PyRAF (IRAF) PyMIDAS (MIDAS) PyIMSL (IMSL)
Monday, October 19, 2009
Custom Distributions
Python(x,y): http://www.pythonxy.com/
Python(x,y) is a free scientific and engineering development software for numerical computations, data analysis and data visualization
Sage: http://www.sagemath.org/
Sage is a free open-source mathematics software system licensed under the GPL. It combines the power of many existing open-source packages into a common Python-based interface.
Monday, October 19, 2009
Extra Astronomy Links
iPython (better shell, distributed computing): http://ipython.scipy.org/ SciPy (collection of science tools): http:// www.scipy.org/ Python Astronomy Modules: http:// astlib.sourceforge.net/ Python Astronomer Wiki: http://macsingularity.org/ astrowiki/tiki-index.php?page=python AstroPy: http://www.astro.washington.edu/users/ rowen/AstroPy.html Python for Astronomers: http://www.iac.es/ sieinvens/siepedia/pmwiki.php? n=HOWTOs.EmpezandoPython
Monday, October 19, 2009
The Basics
Monday, October 19, 2009
A Code Sample
x = 34 - 23 y = Hello z = 3.45 if z == 3.45 or y == Hello: x = x + 1 y = y + World print x print y # String concat. # A comment. # Another one.
Monday, October 19, 2009
Enough to Understand the Code
Assignment uses = and comparison uses ==. For numbers + - * / % are as expected. Special use of + for string concatenation. Special use of % for string formatting (as with printf in C) Logical operators are words (and, or, not) not symbols The basic printing command is print. The first assignment to a variable creates it.
Variable types dont need to be declared. Python figures out the variable types on its own.
Monday, October 19, 2009
Basic Datatypes
Integers (default for numbers)
z = 5 / 2 # Answer is 2, integer division.
Floats
x = 3.456
Strings
Can use or to specify. abc abc (Same thing.) Unmatched can occur within the string. matts Use triple double-quotes for multi-line strings or strings than contain both and inside of them: abc
Monday, October 19, 2009
Whitespace
Whitespace is meaningful in Python: especially indentation and placement of newlines. Use a newline to end a line of code.
Use \ when must go to next line prematurely.
No braces { } to mark blocks of code in Python Use consistent indentation instead.
The first line with less indentation is outside of the block. The first line with more indentation starts a nested block
Often a colon appears at the start of a new block. (E.g. for function and class definitions.)
Monday, October 19, 2009
Comments
Start comments with # the rest of line is ignored. Can include a documentation string as the first line of any new function or class that you define. The development environment, debugger, and other tools use it: its good style to include one.
def my_function(x, y): This is the docstring. This function does blah blah blah. # The code would go here...
Monday, October 19, 2009
Assignment
Binding a variable in Python means setting a name to hold a reference to some object.
Assignment creates references, not copies
Names in Python do not have an intrinsic type. Objects have types.
Python determines the type of the reference automatically based on the data object assigned to it.
You create a name the first time it appears on the left side of an assignment expression:
! x = 3
A reference is deleted via garbage collection after any names bound to it have passed out of scope.
Monday, October 19, 2009
Accessing Non-Existent Names
If you try to access a name before its been properly created (by placing it on the left side of an assignment), youll get an error.
>>> y Traceback (most recent call last): File "<pyshell#16>", line 1, in -toplevely NameError: name y' is not defined >>> y = 3 >>> y 3
Monday, October 19, 2009
Multiple Assignment
You can also assign to multiple names at the same time.
>>> x, y = 2, 3 >>> x 2 >>> y 3
Monday, October 19, 2009
Naming Rules
Names are case sensitive and cannot start with a number. They can contain letters, numbers, and underscores.
bob Bob _bob _2_bob_ bob_2 BoB
There are some reserved words:
and, assert, break, class, continue, def, del, elif, else, except, exec, finally, for, from, global, if, import, in, is, lambda, not, or, pass, print, raise, return, try, while
Monday, October 19, 2009
Understanding Reference Semantics in Python
Monday, October 19, 2009
Understanding Reference Semantics
Assignment manipulates references
x = y does not make a copy of the object y references x = y makes x reference the object y references
Very useful; but beware! Example:
>>> a = [1, 2, 3] >>> b = a >>> a.append(4) >>> print b [1, 2, 3, 4] # a now references the list [1, 2, 3] # b now references what a references # this changes the list a references # if we print what b references, # SURPRISE! It has changed
Why??
Monday, October 19, 2009
Understanding Reference Semantics II
There is a lot going on when we type:
x = 3 First, an integer 3 is created and stored in memory A name x is created An reference to the memory location storing the 3 is then assigned to the name x So: When we say that the value of x is 3 we mean that x now refers to the integer 3 Name: x Ref: <address1> Type: Integer Data: 3
name list
memory
Monday, October 19, 2009
Understanding Reference Semantics III
The data 3 we created is of type integer. In Python, the datatypes integer, float, and string (and tuple) are immutable. This doesnt mean we cant change the value of x, i.e. change what x refers to For example, we could increment x:
>>> x = 3 >>> x = x + 1 >>> print x 4
Monday, October 19, 2009
Understanding Reference Semantics IV
If we increment x, then whats really happening is: 1. The reference of name x is looked up. 2. The value at that reference is retrieved. >>> x = x + 1
Name: x Ref: <address1>
Type: Integer Data: 3
Monday, October 19, 2009
Understanding Reference Semantics IV
If we increment x, then whats really happening is: 1. The reference of name x is looked up. 2. The value at that reference is retrieved. 3. The 3+1 calculation occurs, producing a new data element 4 which is assigned to a fresh memory location with a new reference. >>> x = x + 1
Name: x Ref: <address1>
Type: Integer Data: 3 Type: Integer Data: 4
Monday, October 19, 2009
Understanding Reference Semantics IV
If we increment x, then whats really happening is: 1. The reference of name x is looked up. 2. The value at that reference is retrieved. 3. The 3+1 calculation occurs, producing a new data element 4 which is assigned to a fresh memory location with a new reference. 4. The name x is changed to point to this new reference. >>> x = x + 1
Name: x Ref: <address1>
Type: Integer Data: 3 Type: Integer Data: 4
Monday, October 19, 2009
Understanding Reference Semantics IV
If we increment x, then whats really happening is: 1. The reference of name x is looked up. 2. The value at that reference is retrieved. 3. The 3+1 calculation occurs, producing a new data element 4 which is assigned to a fresh memory location with a new reference. 4. The name x is changed to point to this new reference. 5. The old data 3 is garbage collected if no name still refers to it. >>> x = x + 1
Name: x Ref: <address1> Type: Integer Data: 4
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Name: x Ref: <address1>
Type: Integer Data: 3
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Name: x Ref: <address1> Name: y Ref: <address1>
Type: Integer Data: 3
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Name: x Ref: <address1> Name: y Ref: <address1>
Type: Integer Data: 3 Type: Integer Data: 4
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Name: x Ref: <address1> Name: y Ref: <address2>
Type: Integer Data: 3 Type: Integer Data: 4
Monday, October 19, 2009
Assignment 1
So, for simple built-in datatypes (integers, floats, strings), assignment behaves as you would expect:
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x # # # # Creates 3, name Creates name y, Creates ref for No effect on x, x refers to 3 refers to 3. 4. Changes y. still ref 3.
Name: x Ref: <address1> Name: y Ref: <address2>
Type: Integer Data: 3 Type: Integer Data: 4
Monday, October 19, 2009
Assignment 2
For other data types (lists, dictionaries, user-defined types), assignment works differently.
These datatypes are mutable. When we change these data, we do it in place. We dont copy them into a new memory address each time. If we type y=x and then modify y, both x and y are changed.
immutable
>>> >>> >>> >>> 3 x = 3 y = x y = 4 print x
mutable
x = some mutable object y = x make a change to y look at x x will be changed as well
Monday, October 19, 2009
Why? Changing a Shared List
a = [1, 2, 3] a a b=a b a a.append(4) b 1 2 3 4 1 2 3 1 2 3
Monday, October 19, 2009
Our surprising example surprising no more...
So now, heres our code:
>>> a = [1, 2, 3] >>> b = a >>> a.append(4) >>> print b [1, 2, 3, 4]
# a now references the list [1, 2, 3] # b now references what a references # this changes the list a references # if we print what b references, # SURPRISE! It has changed
Monday, October 19, 2009
Sequence types: Tuples, Lists, and Strings
Monday, October 19, 2009
Sequence Types
1. Tuple
A simple immutable ordered sequence of items Items can be of mixed types, including collection types
2. Strings Immutable Conceptually very much like a tuple 3. List Mutable ordered sequence of items of mixed types
Monday, October 19, 2009
Similar Syntax
All three sequence types (tuples, strings, and lists) share much of the same syntax and functionality.
Key difference:
Tuples and strings are immutable Lists are mutable
The operations shown in this section can be applied to all sequence types
most examples will just show the operation performed on one
Monday, October 19, 2009
Sequence Types 1
Tuples are defined using parentheses (and commas).
>>> tu = (23, abc, 4.56, (2,3), def)
Lists are defined using square brackets (and commas).
>>> li = [abc, 34, 4.34, 23]
Strings are defined using quotes (, , or ).
>>> st >>> st >>> st string = Hello World = Hello World = This is a multi-line that uses triple quotes.
Monday, October 19, 2009
Sequence Types 2
We can access individual members of a tuple, list, or string using square bracket array notation. Note that all are 0 based
>>> tu = (23, abc, 4.56, (2,3), def) >>> tu[1] # Second item in the tuple. abc >>> li = [abc, 34, 4.34, 23] >>> li[1] # Second item in the list. 34 >>> st = Hello World >>> st[1] # Second character in string. e
Monday, October 19, 2009
Positive and negative indices
>>> t = (23, abc, 4.56, (2,3), def)
Positive index: count from the left, starting with 0.
>>> t[1] abc
Negative lookup: count from right, starting with 1.
>>> t[-3] 4.56
Monday, October 19, 2009
Slicing: Return Copy of a Subset 1
>>> t = (23, abc, 4.56, (2,3), def)
Return a copy of the container with a subset of the original members. Start copying at the first index, and stop copying before the second index.
>>> t[1:4] (abc, 4.56, (2,3))
You can also use negative indices when slicing.
>>> t[1:-1] (abc, 4.56, (2,3))
Monday, October 19, 2009
Slicing: Return Copy of a Subset 2
>>> t = (23, abc, 4.56, (2,3), def)
Omit the first index to make a copy starting from the beginning of the container.
>>> t[:2] (23, abc)
Omit the second index to make a copy starting at the first index and going to the end of the container.
>>> t[2:] (4.56, (2,3), def)
Monday, October 19, 2009
Copying the Whole Sequence
To make a copy of an entire sequence, you can use [:].
>>> t[:] (23, abc, 4.56, (2,3), def)
Note the difference between these two lines for mutable sequences:
>>> list2 = list1 # 2 names refer to 1 ref # Changing one affects both
>>> list2 = list1[:] # Two independent copies, two refs
Monday, October 19, 2009
The in Operator
Boolean test whether a value is inside a container:
>>> t >>> 3 False >>> 4 True >>> 4 False = [1, 2, 4, 5] in t in t not in t
For strings, tests for substrings
>>> a = 'abcde'
>>> 'c' in a True >>> 'cd' in a True >>> 'ac' in a False
Be careful: the in keyword is also used in the syntax of for loops and list comprehensions.
Monday, October 19, 2009
The + Operator
The + operator produces a new tuple, list, or string whose value is the concatenation of its arguments.
>>> (1, 2, 3) + (4, 5, 6) (1, 2, 3, 4, 5, 6) >>> [1, 2, 3] + [4, 5, 6] [1, 2, 3, 4, 5, 6] >>> Hello + + World Hello World
Monday, October 19, 2009
The * Operator
The * operator produces a new tuple, list, or string that repeats the original content.
>>> (1, 2, 3) * 3 (1, 2, 3, 1, 2, 3, 1, 2, 3) >>> [1, 2, 3] * 3 [1, 2, 3, 1, 2, 3, 1, 2, 3] >>> Hello * 3 HelloHelloHello
Monday, October 19, 2009
Mutability: Tuples vs. Lists
Monday, October 19, 2009
Tuples: Immutable
>>> t = (23, abc, 4.56, (2,3), def) >>> t[2] = 3.14 Traceback (most recent call last): File "<pyshell#75>", line 1, in -topleveltu[2] = 3.14 TypeError: object doesn't support item assignment
You cant change a tuple. You can make a fresh tuple and assign its reference to a previously used name.
>>> t = (23, abc, 3.14, (2,3), def)
Monday, October 19, 2009
Lists: Mutable
>>> li = [abc, 23, 4.34, 23] >>> li[1] = 45 >>> li [abc, 45, 4.34, 23] We can change lists in place. Name li still points to the same memory reference when were done. The mutability of lists means that they arent as fast as tuples.
Monday, October 19, 2009
Operations on Lists Only 1
>>> li = [1, 11, 3, 4, 5] >>> li.append(a) # Our first exposure to method syntax >>> li [1, 11, 3, 4, 5, a] >>> li.insert(2, i) >>>li [1, 11, i, 3, 4, 5, a]
Monday, October 19, 2009
The extend method vs the + operator.
+ creates a fresh list (with a new memory reference) extend operates on list li in place.
>>> li.extend([9, 8, 7]) >>>li [1, 2, i, 3, 4, 5, a, 9, 8, 7]
Confusing: Extend takes a list as an argument. Append takes a singleton as an argument.
>>> li.append([10, 11, 12]) >>> li [1, 2, i, 3, 4, 5, a, 9, 8, 7, [10, 11, 12]]
Monday, October 19, 2009
Operations on Lists Only 3
>>> li = [a, b, c, b] >>> li.index(b) 1 >>> li.count(b) 2 >>> li.remove(b) >>> li [a, c, b] # index of first occurrence
# number of occurrences
# remove first occurrence
Monday, October 19, 2009
Operations on Lists Only 4
>>> li = [5, 2, 6, 8] >>> li.reverse() >>> li [8, 6, 2, 5] >>> li.sort() >>> li [2, 5, 6, 8] # reverse the list *in place*
# sort the list *in place*
>>> li.sort(some_function) # sort in place using user-defined comparison
Monday, October 19, 2009
Tuples vs. Lists
Lists slower but more powerful than tuples.
Lists can be modified, and they have lots of handy operations we can perform on them. Tuples are immutable and have fewer features.
To convert between tuples and lists use the list() and tuple() functions:
li = list(tu) tu = tuple(li)
Monday, October 19, 2009
Dictionaries
63
Monday, October 19, 2009
Dictionaries: A Mapping type
Dictionaries store a mapping between a set of keys and a set of values.
Keys can be any immutable type. Values can be any type A single dictionary can store values of different types
You can define, modify, view, lookup, and delete the key-value pairs in the dictionary.
Monday, October 19, 2009
Using dictionaries
>>> d = {user:bozo, pswd:1234} >>> d[user] bozo >>> d[pswd] 1234 >>> d[bozo] Traceback (innermost last): File <interactive input> line 1, in ? KeyError: bozo >>> d = {user:bozo, pswd:1234} >>> d[user] = clown >>> d {user:clown, pswd:1234} >>> d[id] = 45 >>> d {user:clown, id:45, pswd:1234}
>>> d = {user:bozo, p:1234, i:34} >>> del d[user] # Remove one. >>> d {p:1234, i:34} >>> d.clear() # Remove all. >>> d {}
>>> d = {user:bozo, p:1234, i:34} >>> d.keys() # List of keys. [user, p, i] >>> d.values() # List of values. [bozo, 1234, 34] >>> d.items() # List of item tuples. [(user,bozo), (p,1234), (i,34)]
Monday, October 19, 2009
Functions
Monday, October 19, 2009
Functions
def creates a function and assigns it a name return sends a result back to the caller Arguments are passed by assignment Arguments and return types are not declared
def <name>(arg1, arg2, ..., argN): ! <statements> ! return <value> def times(x,y): ! return x*y
Monday, October 19, 2009
Passing Arguments to Functions
Arguments are passed by assignment Passed arguments are assigned to local names Assignment to argument names don't affect the caller Changing a mutable argument may affect the caller
def changer (x,y): ! x = 2! ! ! y[0] = 'hi'! ! ! # changes local value of x only # changes shared object
Monday, October 19, 2009
Optional Arguments
Can define defaults for arguments that need not be passed
def func(a, b, c=10, d=100): ! print a, b, c, d >>> func(1,2) 1 2 10 100 >>> func(1,2,3,4) 1,2,3,4
Monday, October 19, 2009
Gotchas
All functions in Python have a return value
even if no return line inside the code.
Functions without a return return the special value None. There is no function overloading in Python.
Two different functions cant have the same name, even if they have different arguments.
Functions can be used as any other data type. They can be:
Arguments to function Return values of functions Assigned to variables Parts of tuples, lists, etc
Monday, October 19, 2009
Control of Flow
71
Monday, October 19, 2009
Examples
if x == 3: print X equals 3. elif x == 2: print X equals 2. else: print X equals something else. print This is outside the if. assert(number_of_players < 5)
x = 3 while x < 10: if x > 7: x += 2 continue x = x + 1 print Still in the loop. if x == 8: break print Outside of the loop.
for x in range(10): if x > 7: x += 2 continue x = x + 1 print Still in the loop. if x == 8: break print Outside of the loop.
Monday, October 19, 2009
Modules
Monday, October 19, 2009
Why Use Modules?
Code reuse
Routines can be called multiple times within a program Routines can be used from multiple programs
Namespace partitioning
Group data together with functions used for that data
Implementing shared services or data
Can provide global data structure that is accessed by multiple subprograms
Monday, October 19, 2009
Modules
Modules are functions and variables defined in separate files Items are imported using from or import
from module import function function() import module module.function()
Modules are namespaces
Can be used to organize variable names, i.e.
atom.position = atom.position - molecule.position
Monday, October 19, 2009
Classes and Objects
Monday, October 19, 2009
What is an Object?
A software item that contains variables and methods Object Oriented Design focuses on
Encapsulation:
dividing the code into a public interface, and a private implementation of that interface
Polymorphism:
the ability to overload standard operators so that they have appropriate behavior based on their context
Inheritance:
the ability to create subclasses that contain specializations of their parents
Monday, October 19, 2009
Example
class ! def !! !! ! def !! ! def !! !! !! atom(object): __init__(self,atno,x,y,z): self.atno = atno self.position = (x,y,z) symbol(self): # a class method return Atno_to_Symbol[atno] __repr__(self): # overloads printing return '%d %10.4f %10.4f %10.4f' % ! (self.atno, self.position[0], ! self.position[1],self.position[2])
>>> at = atom(6,0.0,1.0,2.0) >>> print at 6 0.0000 1.0000 2.0000 >>> at.symbol() 'C'
Monday, October 19, 2009
Atom Class
Overloaded the default constructor Defined class variables (atno,position) that are persistent and local to the atom object Good way to manage shared memory:
instead of passing long lists of arguments, encapsulate some of this data into an object, and pass the object. much cleaner programs result
Overloaded the print operator We now want to use the atom class to build molecules...
Monday, October 19, 2009
Molecule Class
class ! def !! !! ! def !! ! def !! !! !! !! !! molecule: __init__(self,name='Generic'): self.name = name self.atomlist = [] addatom(self,atom): self.atomlist.append(atom) __repr__(self): str = 'This is a molecule named %s\n' % self.name str = str+'It has %d atoms\n' % len(self.atomlist) for atom in self.atomlist: ! str = str + `atom` + '\n' return str
Monday, October 19, 2009
Using Molecule Class
>>> mol = molecule('Water') >>> at = atom(8,0.,0.,0.) >>> mol.addatom(at) >>> mol.addatom(atom(1,0.,0.,1.)) >>> mol.addatom(atom(1,0.,1.,0.)) >>> print mol This is a molecule named Water It has 3 atoms 8 0.000 0.000 0.000 1 0.000 0.000 1.000 1 0.000 1.000 0.000
Note that the print function calls the atoms print function
Code reuse: only have to type the code that prints an atom once; this means that if you change the atom specification, you only have one place to update.
Monday, October 19, 2009
Inheritance
class qm_molecule(molecule): def addbasis(self): self.basis = [] for atom in self.atomlist: self.basis = add_bf(atom,self.basis)
__init__, __repr__, and __addatom__ are taken from the parent class (molecule) Added a new function addbasis() to add a basis set Another example of code reuse
Basic functions don't have to be retyped, just inherited Less to rewrite when specifications change
Monday, October 19, 2009
Overloading
! ! ! ! ! ! ! ! class qm_molecule(molecule): def __repr__(self): str = 'QM Rules!\n' for atom in self.atomlist: ! str = str + `atom` + '\n' return str
Now we only inherit __init__ and addatom from the parent We define a new version of __repr__ specially for QM
Monday, October 19, 2009
Adding to Parent Functions
Sometimes you want to extend, rather than replace, the parent functions.
! !! !! class qm_molecule(molecule): def __init__(self,name="Generic",basis="6-31G**"): self.basis = basis super(qm_molecule, self).__init__(name)
Monday, October 19, 2009
Public and Private Data
In Python anything with two leading underscores is private
__a, __my_variable
Anything with one leading underscore is semiprivate, and you should feel guilty accessing this data directly.
_b
Sometimes useful as an intermediate step to making data private
Monday, October 19, 2009
The Extra Stuff...
86
Monday, October 19, 2009
File I/O, Strings, Exceptions...
>>> try: ... 1 / 0 ... except: ... print('That was silly!') ... finally: ... print('This gets executed no matter what') ... That was silly! This gets executed no matter what fileptr = open(filename) somestring = fileptr.read() for line in fileptr: print line fileptr.close()
>>> a = 1 >>> b = 2.4 >>> c = 'Tom' >>> '%s has %d coins worth a total of $%.02f' % (c, a, b) 'Tom has 1 coins worth a total of $2.40'
Monday, October 19, 2009
