diff options
Diffstat (limited to 'day2')
| -rw-r--r-- | day2/exercises.tex | 201 | ||||
| -rw-r--r-- | day2/session1.tex | 45 | ||||
| -rw-r--r-- | day2/session2.tex | 79 | ||||
| -rw-r--r-- | day2/session3.tex | 674 | ||||
| -rw-r--r-- | day2/session4.tex | 332 | ||||
| -rw-r--r-- | day2/session5.tex | 563 |
6 files changed, 1168 insertions, 726 deletions
diff --git a/day2/exercises.tex b/day2/exercises.tex new file mode 100644 index 0000000..6c3fd19 --- /dev/null +++ b/day2/exercises.tex @@ -0,0 +1,201 @@ +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%Tutorial slides on Python. +% +% Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in> +% Copyright (c) 2005-2009, Prabhu Ramachandran +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +\documentclass[14pt,compress]{beamer} +%\documentclass[draft]{beamer} +%\documentclass[compress,handout]{beamer} +%\usepackage{pgfpages} +%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm] + +% Modified from: generic-ornate-15min-45min.de.tex +\mode<presentation> +{ + \usetheme{Warsaw} + \useoutertheme{infolines} + \setbeamercovered{transparent} +} + +\usepackage[english]{babel} +\usepackage[latin1]{inputenc} +%\usepackage{times} +\usepackage[T1]{fontenc} + +% Taken from Fernando's slides. +\usepackage{ae,aecompl} +\usepackage{mathpazo,courier,euler} +\usepackage[scaled=.95]{helvet} + +\definecolor{darkgreen}{rgb}{0,0.5,0} + +\usepackage{listings} +\lstset{language=Python, + basicstyle=\ttfamily\bfseries, + commentstyle=\color{red}\itshape, + stringstyle=\color{darkgreen}, + showstringspaces=false, + keywordstyle=\color{blue}\bfseries} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Macros +\setbeamercolor{emphbar}{bg=blue!20, fg=black} +\newcommand{\emphbar}[1] +{\begin{beamercolorbox}[rounded=true]{emphbar} + {#1} + \end{beamercolorbox} +} +\newcounter{time} +\setcounter{time}{0} +\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}} + +\newcommand{\typ}[1]{\texttt{#1}} + +\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}} } + +%%% This is from Fernando's setup. +% \usepackage{color} +% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2} +% % Use and configure listings package for nicely formatted code +% \usepackage{listings} +% \lstset{ +% language=Python, +% basicstyle=\small\ttfamily, +% commentstyle=\ttfamily\color{blue}, +% stringstyle=\ttfamily\color{orange}, +% showstringspaces=false, +% breaklines=true, +% postbreak = \space\dots +% } + + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Title page +\title[Exercises]{Exercises} + +\author[FOSSEE] {FOSSEE} + +\institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay} +\date[] {8 November, 2009\\Day 2} +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo} +%\logo{\pgfuseimage{iitmlogo}} + + +%% Delete this, if you do not want the table of contents to pop up at +%% the beginning of each subsection: +\AtBeginSubsection[] +{ + \begin{frame}<beamer> + \frametitle{Outline} + \tableofcontents[currentsection,currentsubsection] + \end{frame} +} + + +% If you wish to uncover everything in a step-wise fashion, uncomment +% the following command: +%\beamerdefaultoverlayspecification{<+->} + +%\includeonlyframes{current,current1,current2,current3,current4,current5,current6} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% DOCUMENT STARTS +\begin{document} + +\begin{frame} + \titlepage +\end{frame} + +\begin{frame}{Problem 1.1} + The aliquot of a number is defined as: the sum of the \emph{proper} divisors of the number. For example:\\ +aliquot(12) = 1 + 2 + 3 + 4 + 6 = 16.\\ + Write a function that returns the aliquot number of a given number. +\end{frame} + +\begin{frame}{Problem 1.2} + Pair of numbers (a, b) is said to be \alert{amicable} if aliquot number of a is b and aliquot number of b is a.\\ + Example: \texttt{220, 284}\\ + Write a program that prints all four digit amicable pairs. + +\inctime{20} +\end{frame} + +%% \begin{frame}{Problem 2} +%% Given an empty chessboard and one Bishop placed in any s%quare, say (r, c), generate the list of all squares the Bi%shop could move to. + +%% \end{frame} + +\begin{frame}[fragile] + \frametitle{Problem 2} + Given a string of numbers like, ``1, 3-7, 12, 15, 18-21'', produce the following list \\ + \begin{lstlisting} + [1,3,4,5,6,7,12,15,18,19,20,21] + \end{lstlisting} +\inctime{10} +\end{frame} + +\begin{frame} + \frametitle{Problem 3} + Count frequencies of words in a file. +\inctime{25} +\end{frame} + +\begin{frame}[fragile] + \frametitle{Problem set 4} + Central difference + \begin{equation*} + \frac{sin(x+h)-sin(x-h)}{2h} + \end{equation*} + \begin{lstlisting} + In []: x = linspace(0, 2*pi, 100) + In []: y = sin(x) + In []: deltax = x[1] - x[0] + \end{lstlisting} + \pause + \begin{enumerate} + \item Given this, get the finite difference of sin in the range 0 to 2*pi + \end{enumerate} +\end{frame} + +\begin{frame} + \frametitle{Problem Set 5} + \begin{itemize} + \item[5.1] Write a function that plots any regular n-gon given \typ{n}. + \item[5.2] Consider the logistic map, $f(x) = kx(1-x)$, plot it for + $k=2.5, 3.5$ and $4$ in the same plot. +\end{itemize} +\end{frame} + +\begin{frame}[fragile] +\frametitle{Problem Set 5} + \begin{columns} + \column{0.6\textwidth} + \small{ + \begin{itemize} + \item[3] Consider the iteration $x_{n+1} = f(x_n)$ where $f(x) = kx(1-x)$. Plot the successive iterates of this process as explained below. + \end{itemize}} + \column{0.35\textwidth} + \hspace*{-0.5in} + \includegraphics[height=1.6in, interpolate=true]{data/cobweb} +\end{columns} +\end{frame} + +\begin{frame} + \frametitle{Problem Set 5.3} + Plot the cobweb plot as follows: + \begin{enumerate} + \item Start at $(x_0, 0)$ ($\implies$ i=0) + \item Draw a line to $(x_i, f(x_i))$ + \item Set $x_{i+1} = f(x_i)$ + \item Draw a line to $(x_{i+1}, x_{i+1})$ + \item $(i\implies i+1)$ + \item Repeat from 2 for as long as you want + \end{enumerate} +\inctime{20} +\end{frame} + +\end{document} diff --git a/day2/session1.tex b/day2/session1.tex index 2d5288b..df67b5c 100644 --- a/day2/session1.tex +++ b/day2/session1.tex @@ -51,7 +51,8 @@ \setcounter{time}{0} \newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}} -\newcommand{\typ}[1]{\texttt{#1}} +\newcommand{\typ}[1]{\lstinline{#1}} + \newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}} } @@ -130,7 +131,7 @@ \begin{itemize} \item Numbers: float, int, complex \item Strings - \item Boolean + \item Booleans \end{itemize} \end{frame} @@ -138,7 +139,7 @@ \begin{frame}[fragile] \frametitle{Numbers} \begin{itemize} - \item \kwrd{int}\\ \kwrd{int} = whole number, no matter what the size! + \item \kwrd{int}\\ whole number, no matter what the size! \begin{lstlisting} In []: a = 13 @@ -167,9 +168,9 @@ Out[]: 3.0 \end{lstlisting} \end{frame} -\subsection{Boolean} +\subsection{Booleans} \begin{frame}[fragile] - \frametitle{Boolean} + \frametitle{Booleans} \begin{lstlisting} In []: t = True @@ -211,15 +212,19 @@ In []: w = "hello" In []: print w[0] + w[2] + w[-1] Out[]: hlo -In []: len(w) # guess what +In []: len(w) Out[]: 5 \end{lstlisting} \end{frame} \begin{frame}[fragile] \frametitle{Strings \ldots} + \emphbar{Strings are immutable} + \begin{lstlisting} +In []: w[0] = 'H' + \end{lstlisting} + \pause \begin{lstlisting} -In []: w[0] = 'H' # Can't do that! -------------------------------------------- TypeError Traceback (most recent call last) @@ -248,7 +253,7 @@ Out[]: 'hello world' \end{lstlisting} \end{frame} -\begin{frame} +\begin{frame}[fragile] \frametitle{A bit about IPython} \begin{itemize} \item IPython provides better help @@ -265,7 +270,7 @@ In []: a.s<Tab> \end{frame} \begin{frame}[fragile] -\frametitle{Still with strings} + \frametitle{Still with strings} \begin{itemize} \item We saw split() yesterday \item join() is the opposite of split() @@ -294,6 +299,7 @@ Out[]: 'x is 1, y is 1.234' \section{Operators} \begin{frame}[fragile] \frametitle{Arithmetic operators} + \small \begin{lstlisting} In []: 1786 % 12 Out[]: 10 @@ -351,15 +357,15 @@ In []: a /= 5 \begin{frame}[fragile] \frametitle{String operations} \begin{lstlisting} -In []: s = 'Hello ' +In []: s = 'Hello' In []: p = 'World' In []: s + p -Out[]: 'Hello World' +Out[]: 'HelloWorld' In []: s * 4 -Out[]: 'Hello Hello Hello Hello' +Out[]: 'HelloHelloHelloHello' \end{lstlisting} \end{frame} @@ -367,13 +373,16 @@ Out[]: 'Hello Hello Hello Hello' \frametitle{String operations \ldots} \begin{lstlisting} In []: s * s + \end{lstlisting} + \pause + \begin{lstlisting} -------------------------------------------- TypeError Traceback (most recent call last) /<ipython console> in <module>() -TypeError: can't multiply sequence by - non-int of type 'str' +TypeError: can`t multiply sequence by + non-int of type `str` \end{lstlisting} \end{frame} @@ -401,7 +410,7 @@ Out[]: False In []: int(17 / 2.0) Out[]: 8 -In []: float(17 / 2) # Recall +In []: float(17 / 2) Out[]: 8.0 In []: str(17 / 2.0) @@ -455,7 +464,7 @@ Enter a value: 5 \frametitle{Simple IO: Console output} \begin{itemize} \item \texttt{print} is straight forward - \item Put the following code snippet in a file + \item Put the following code snippet in a file \emph{hello1.py} \end{itemize} \begin{lstlisting} print "Hello" @@ -470,7 +479,7 @@ World \begin{frame}[fragile] \frametitle{Simple IO: Console output \ldots} -Put the following code snippet in a file +Put the following code snippet in a file \emph{hello2.py} \begin{lstlisting} print "Hello", print "World" @@ -522,7 +531,7 @@ In []: if x < 0: \begin{lstlisting} ... a = raw_input('Enter number(Q to quit):') -num = int(a) if a != 'Q' else break +num = int(a) if a != 'Q' else 0 ... \end{lstlisting} \end{frame} diff --git a/day2/session2.tex b/day2/session2.tex index c20a2bb..14fe857 100644 --- a/day2/session2.tex +++ b/day2/session2.tex @@ -51,7 +51,7 @@ \setcounter{time}{0} \newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}} -\newcommand{\typ}[1]{\texttt{#1}} +\newcommand{\typ}[1]{\lstinline{#1}} \newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}} } @@ -127,10 +127,10 @@ \subsection{Basic Looping} \begin{frame}[fragile] \frametitle{\typ{while}} -Example: Fibonacci series +\begin{block}{Example: Fibonacci series} + Sum of previous two elements defines the next +\end{block} \begin{lstlisting} -# the sum of two elements -# defines the next In []: a, b = 0, 1 In []: while b < 10: ...: print b, @@ -145,23 +145,24 @@ In []: while b < 10: \frametitle{\typ{range()}} \kwrd{range([start,] stop[, step])}\\ \begin{itemize} - \item range() returns a list of integers - \item The \emph{start} and the \emph{step} arguments are optional - \item \emph{stop} argument is not included in the list + \item \typ{range()} returns a list of integers + \item The \typ{start} and the \typ{step} arguments are optional + \item \typ{stop} is not included in the list \end{itemize} \vspace*{.5in} -\begin{itemize} - \item \alert{Anything within \typ{[]} is optional} +\begin{block}{Documentation convention} \begin{itemize} - \item Nothing to do with Python. + \item \alert{Anything within \typ{[]} is optional} + \begin{itemize} + \item Nothing to do with Python. + \end{itemize} \end{itemize} -\end{itemize} - +\end{block} \end{frame} \begin{frame}[fragile] - \frametitle{\typ{for} \ldots \typ{range()}} -Example: print squares of first \typ{n} numbers + \frametitle{\texttt{for} \ldots \typ{range()}} +Example: print squares of first \typ{5} numbers \begin{lstlisting} In []: for i in range(5): ....: print i, i * i @@ -176,7 +177,7 @@ In []: for i in range(5): \end{frame} \begin{frame}[fragile] - \frametitle{\typ{for} \ldots \typ{range()}} + \frametitle{\texttt{for} \ldots \typ{range()}} Example: print squares of odd numbers from 3 to 9 \begin{lstlisting} In []: for i in range(3, 10, 2): @@ -240,6 +241,9 @@ num = [1, 2, 3, 4] \begin{lstlisting} In []: num = [1, 2, 3, 4] +In []: num + [9, 10, 11] +Out[]: [1, 2, 3, 4, 9, 10, 11] + In []: num.append([9, 10, 11]) In []: num @@ -325,7 +329,6 @@ In []: t = (1, 2, 3, 4, 5, 6, 7, 8) In []: t[0] + t[3] + t[-1] Out[]: 13 -# Try the following! In []: t[4] = 7 \end{lstlisting} \pause @@ -361,8 +364,8 @@ In []: player['Avg'] Out[]: 52.530000000000001 \end{lstlisting} \begin{block}{Note!} - Duplicate keys are not allowed!\\ - Dictionaries are iterable through keys. + Duplicate keys $\Rightarrow$ overwritten!\\ + You can iterate through a dictionary using keys. \end{block} \end{frame} @@ -375,6 +378,12 @@ Out[]: True In []: 'Econ' in player Out[]: False \end{lstlisting} + \begin{block}{Note} + \begin{itemize} + \item We can check for the containership of keys only + \item Not values + \end{itemize} + \end{block} \end{frame} \begin{frame}[fragile] @@ -384,17 +393,18 @@ In []: player.keys() Out[]: ['Runs', 'Inn', 'Avg', 'Mat'] In []: player.values() -Out[]: [10823, 233, 52.530000000000001, 134] +Out[]: [10823, 233, + 52.530000000000001, 134] \end{lstlisting} \end{frame} \begin{frame} {Problem Set 2.1: Problem 2.1.1} -You are given date strings of the form ``29, Jul 2009'', or ``4 January 2008''. In other words a number, a string and another number, with a comma sometimes separating the items.\\Write a function that takes such a string and returns a tuple (yyyy, mm, dd) where all three elements are ints. +You are given date strings of the form ``29 Jul, 2009'', or ``4 January 2008''. In other words a number, a string and another number, with a comma sometimes separating the items.\\Write a function that takes such a string and returns a tuple (yyyy, mm, dd) where all three elements are ints. \end{frame} -\subsection{Set} +\subsection{Sets} \begin{frame}[fragile] - \frametitle{Set} + \frametitle{Sets} \begin{itemize} \item Simplest container, mutable \item No ordering, no duplicates @@ -447,18 +457,14 @@ Out[]: 5 \end{frame} \begin{frame} - \frametitle{Problem set 2.2} - \begin{description} - \item[2.2.1] Given a dictionary of the names of students and their marks, identify how many duplicate marks are there? and what are these? -\end{description} -\inctime{15} + \frametitle{Problem set 2.2: Problem 2.2.1} +Given a dictionary of the names of students and their marks, identify how many duplicate marks are there? and what are these? \end{frame} \begin{frame} - \frametitle{Problem set 2.2} - \begin{description} - \item[2.2.2] Given a list of words, find all the anagrams in the list -\end{description} + \frametitle{Problem 2.2.2} +Given a list of words, find all the anagrams in the list. + \inctime{15} \end{frame} @@ -468,7 +474,6 @@ Out[]: 5 \begin{itemize} \item \kwrd{def} - keyword to define a function \item Arguments are local to a function - \item Docstrings are important! \item Functions can return multiple values \end{itemize} \end{frame} @@ -487,6 +492,7 @@ def signum( r ): else: return 0 \end{lstlisting} + \emphbar{Note docstrings} \end{frame} \begin{frame}[fragile] @@ -517,8 +523,15 @@ def what( n ): \frametitle{What did we learn?} \begin{itemize} \item Loops: \kwrd{while}, \kwrd{for} - \item Advanced Data structures + \item Advanced Data structures: + \begin{itemize} + \item Lists + \item Tuples + \item Dictionaries + \item Sets + \end{itemize} \item Functions + \item Docstrings \end{itemize} \end{frame} diff --git a/day2/session3.tex b/day2/session3.tex index af621e4..068f17e 100644 --- a/day2/session3.tex +++ b/day2/session3.tex @@ -1,48 +1,33 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% Tutorial slides on Python. +%Tutorial slides on Python. % -% Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in> -% Copyright (c) 2005-2009, Prabhu Ramachandran +% Author: FOSSEE +% Copyright (c) 2009, FOSSEE, IIT Bombay %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\documentclass[compress,14pt]{beamer} -% \documentclass[handout]{beamer} -% \usepackage{pgfpages} -% \pgfpagesuselayout{4 on 1}[a4paper,border, shrink=5mm,landscape] -\usepackage{tikz} -\newcommand{\hyperlinkmovie}{} -%\usepackage{movie15} - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% Note that in presentation mode -% \paperwidth 364.19536pt -% \paperheight 273.14662pt -% h/w = 0.888 - +\documentclass[14pt,compress]{beamer} +%\documentclass[draft]{beamer} +%\documentclass[compress,handout]{beamer} +%\usepackage{pgfpages} +%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm] +% Modified from: generic-ornate-15min-45min.de.tex \mode<presentation> { \usetheme{Warsaw} - %\usetheme{Boadilla} - %\usetheme{default} \useoutertheme{infolines} \setbeamercovered{transparent} } -% To remove navigation symbols -\setbeamertemplate{navigation symbols}{} - -\usepackage{amsmath} \usepackage[english]{babel} \usepackage[latin1]{inputenc} -\usepackage{times} +%\usepackage{times} \usepackage[T1]{fontenc} % Taken from Fernando's slides. \usepackage{ae,aecompl} \usepackage{mathpazo,courier,euler} \usepackage[scaled=.95]{helvet} -\usepackage{pgf} \definecolor{darkgreen}{rgb}{0,0.5,0} @@ -55,65 +40,50 @@ keywordstyle=\color{blue}\bfseries} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% My Macros -\setbeamercolor{postit}{bg=yellow,fg=black} +% Macros \setbeamercolor{emphbar}{bg=blue!20, fg=black} \newcommand{\emphbar}[1] {\begin{beamercolorbox}[rounded=true]{emphbar} {#1} \end{beamercolorbox} } -%{\centerline{\fcolorbox{gray!50} {blue!10}{ -%\begin{minipage}{0.9\linewidth} -% {#1} -%\end{minipage} -% }}} - -\newcommand{\myemph}[1]{\structure{\emph{#1}}} -\newcommand{\PythonCode}[1]{\lstinline{#1}} - -\newcommand{\tvtk}{\texttt{tvtk}} -\newcommand{\mlab}{\texttt{mlab}} - \newcounter{time} \setcounter{time}{0} -\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\vspace*{0.1in}\tiny \thetime\ m}} +\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}} -\newcommand\BackgroundPicture[1]{% - \setbeamertemplate{background}{% - \parbox[c][\paperheight]{\paperwidth}{% - \vfill \hfill - \hfill \vfill -}}} +\newcommand{\typ}[1]{\lstinline{#1}} -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% Configuring the theme -%\setbeamercolor{normal text}{fg=white} -%\setbeamercolor{background canvas}{bg=black} +\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}} } +%%% This is from Fernando's setup. +% \usepackage{color} +% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2} +% % Use and configure listings package for nicely formatted code +% \usepackage{listings} +% \lstset{ +% language=Python, +% basicstyle=\small\ttfamily, +% commentstyle=\ttfamily\color{blue}, +% stringstyle=\ttfamily\color{orange}, +% showstringspaces=false, +% breaklines=true, +% postbreak = \space\dots +% } %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Title page -\title[3D Plotting]{3D data Visualization} +\title[Basic Python]{Python language: Functions, modules and objects} -\author[FOSSEE] {FOSSEE} +\author[FOSSEE Team] {The FOSSEE Group} \institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay} -\date[] {8 November, 2009\\Day 2, Session 5} - +\date[] {8 November, 2009\\Day 2, Session 3} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%\pgfdeclareimage[height=0.75cm]{iitblogo}{iitblogo} -%\logo{\pgfuseimage{iitblogo}} +%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo} +%\logo{\pgfuseimage{iitmlogo}} -\AtBeginSection[] -{ - \begin{frame}<beamer> - \frametitle{Outline} - \tableofcontents[currentsection,currentsubsection] - \end{frame} -} %% Delete this, if you do not want the table of contents to pop up at %% the beginning of each subsection: @@ -132,12 +102,19 @@ \tableofcontents[currentsection,currentsubsection] \end{frame} } + +% If you wish to uncover everything in a step-wise fashion, uncomment +% the following command: +%\beamerdefaultoverlayspecification{<+->} + +%\includeonlyframes{current,current1,current2,current3,current4,current5,current6} + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % DOCUMENT STARTS \begin{document} \begin{frame} - \maketitle + \titlepage \end{frame} \begin{frame} @@ -146,419 +123,312 @@ % You might wish to add the option [pausesections] \end{frame} -\section{3D Data Visualization} - -\begin{frame} - \frametitle{What is visualization?} - \Large - \begin{center} - Visual representation of data - \end{center} -\end{frame} +\section{Functions} +\subsection{Default arguments} +\begin{frame}[fragile] + \frametitle{Functions: default arguments} + \begin{lstlisting} +In []: greet = 'hello world' +In []: greet.split() +Out[]: ['hello', 'world'] -%% \begin{frame} -%% \frametitle{Is this new?} -%% \begin{center} -%% We have moved from: -%% \end{center} -%% \begin{columns} -%% \column{} -%% \hspace*{-1in} -%% \includegraphics[width=1.75in,height=1.75in, interpolate=true]{data/3832} -%% \column{}\hspace*{-0.25in} -%% To -%% \column{} -%% \hspace*{-1in} -%% \includegraphics[width=1.75in, height=1.75in, interpolate=true]{data/torus} -%% \end{columns} -%% \end{frame} +In []: line = 'Rossum, Guido, 54, 46, 55' -\begin{frame} - \frametitle{3D visualization} - \Large - \begin{center} - Harder but important - \end{center} -\end{frame} - -\begin{frame} - \frametitle{Is this Graphics?} - \Large - \begin{center} - Visualization is about data! - \end{center} +In []: line.split(',') +Out[]: ['Rossum', ' Guido', ' 54', + ' 46', ' 55'] + \end{lstlisting} \end{frame} -\begin{frame} - \frametitle{Examples: trajectory in space} - \Large - \begin{center} - \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex} - \end{center} -\end{frame} +\begin{frame}[fragile] + \frametitle{Functions: default arguments \ldots} + \begin{lstlisting} +In []: def welcome(greet, name="World"): + .... print greet, name -\begin{frame} - \frametitle{Examples: Fire in a room} - \Large - \begin{center} - Demo of data - \end{center} -\inctime{10} -\end{frame} +In []: welcome("Hello") +Hello World -\section{Tools available} +In []: welcome("Hi", "Guido") +Hi Guido + \end{lstlisting} +\end{frame} -\subsection{mlab} +\subsection{Keyword arguments} +\begin{frame}[fragile] + \frametitle{Functions: Keyword arguments} +We have seen the following +\begin{lstlisting} +In []: legend(['sin(2y)'], + loc = 'center') -\begin{frame} - {Overview} - \Large - \begin{itemize} - \item Simple - \item Convenient - \item Full-featured - \end{itemize} -\end{frame} +In []: plot(y, sin(y), 'g', + linewidth = 2) -\begin{frame}[fragile] +In []: annotate('local max', + xy = (1.5, 1)) - \frametitle{Getting started} - \myemph{\Large Vanilla:} - \begin{lstlisting}[language=bash] - $ ipython -wthread - \end{lstlisting} - \myemph{\Large with Pylab:} - \begin{lstlisting}[language=bash] - $ ipython -pylab -wthread - \end{lstlisting} +In []: pie(science.values(), + labels = science.keys()) + \end{lstlisting} \end{frame} \begin{frame}[fragile] - \frametitle{Using mlab} - - \begin{lstlisting} -In []:from enthought.mayavi import mlab - \end{lstlisting} - - \vspace*{0.5in} + \frametitle{Functions: keyword arguments \ldots} + \begin{lstlisting} +In []: def welcome(greet, name="World"): + .... print greet, name - \myemph{\Large Try these} +In []: welcome("Hello", "James") +Hello James - \vspace*{0.25in} +In []: welcome("Hi", name="Guido") +Hi Guido - \begin{lstlisting} -In []: mlab.test_<TAB> -In []: mlab.test_contour3d() -In []: mlab.test_contour3d?? - \end{lstlisting} +In []: welcome(name="Guido", greet="Hey") +Hey Guido + \end{lstlisting} \end{frame} +\subsection{Built-in functions} \begin{frame} - {Exploring the view} - \begin{columns} - \column{0.6\textwidth} - \pgfimage[width=3in]{MEDIA/m2/contour3d} - \column{0.4\textwidth} - \begin{itemize} - \item Mouse - \item Keyboard - \item Toolbar - \item Mayavi icon\pgfimage[width=0.2in]{MEDIA/m2/m2_icon} - \end{itemize} - \end{columns} + {Before writing a function} + \begin{itemize} + \item Variety of built-in functions are available + \item \typ{abs, any, all, len, max, min} + \item \typ{pow, range, sum, type} + \item Refer here: + \url{http://docs.python.org/library/functions.html} + \end{itemize} + \inctime{10} \end{frame} -\begin{frame}[fragile] - \frametitle{Clearing the plot area} - \PythonCode{In []: mlab.clf()} +\subsection{Exercises} +\begin{frame}{Problem set 3: Problem 3.1} + Write a function to return the gcd of two numbers. \end{frame} -\begin{frame}[fragile] - \frametitle{\mlab\ plotting functions} - \begin{columns} - \column{0.25\textwidth} - \myemph{\Large 0D data} - \column{0.5\textwidth} - \pgfimage[width=2in]{MEDIA/m2/mlab/points3d_ex} - \end{columns} - - \begin{lstlisting} -In []: t = linspace(0, 2*pi, 50) -In []: u = cos(t) * pi -In []: x, y, z = sin(u), cos(u), sin(t) - \end{lstlisting} - \emphbar{\PythonCode{In []: mlab.points3d(x, y, z)}} +\begin{frame}{Problem 3.2} +Write a program to print all primitive pythagorean triads (a, b, c) where a, b are in the range 1---100 \\ +A pythagorean triad $(a,b,c)$ has the property $a^2 + b^2 = c^2$.\\By primitive we mean triads that do not `depend' on others. For example, (4,3,5) is a variant of (3,4,5) and hence is not primitive. And (10,24,26) is easily derived from (5,12,13) and is also not primitive. \end{frame} -\begin{frame} - \begin{columns} - \column{0.25\textwidth} - \myemph{\Large 1D data} - \column{0.5\textwidth} - \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex} - \end{columns} - \PythonCode{In []: mlab.clf()} - \emphbar{\PythonCode{In []: mlab.plot3d(x, y, z, t)}} - - Plots lines between the points - +\begin{frame}{Problem 3.3} + Write a program that generates a list of all four digit numbers that have all their digits even and are perfect squares.\newline\\\emph{For example, the output should include 6400 but not 8100 (one digit is odd) or 4248 (not a perfect square).} + +\inctime{15} \end{frame} +\section{Modules} \begin{frame}[fragile] - \begin{columns} - \column{0.25\textwidth} - \myemph{\Large 2D data} - \column{0.5\textwidth} - \pgfimage[width=2in]{MEDIA/m2/mlab/surf_ex} - \end{columns} - \begin{lstlisting} -In []: x, y = mgrid[-3:3:100j,-3:3:100j] -In []: z = sin(x*x + y*y) - \end{lstlisting} - - \emphbar{\PythonCode{In []: mlab.surf(x, y, z)}} - - \alert{Assumes the points are rectilinear} + \frametitle{\texttt{from} \ldots \texttt{import} magic} + \begin{lstlisting} +from scipy.integrate import odeint +from scipy.optimize import fsolve + \end{lstlisting} +\emphbar{Above statements import a function to our namespace} \end{frame} \begin{frame}[fragile] - \frametitle{mgrid} + \frametitle{Running scripts from command line} + \small + \begin{itemize} + \item Fire up a terminal + \item python four\_plot.py + \end{itemize} + \pause \begin{lstlisting} -In []: mgrid[0:3,0:3] -Out[]: -array([[[0, 0, 0], - [1, 1, 1], - [2, 2, 2]], - - [[0, 1, 2], - [0, 1, 2], - [0, 1, 2]]]) - -In []: mgrid[-1:1:5j] -Out[]: array([-1., -0.5, 0., 0.5, 1.]) -\end{lstlisting} +Traceback (most recent call last): + File "four_plot.py", line 1, in <module> + x = linspace(-5*pi, 5*pi, 500) +NameError: name 'linspace' is not defined + \end{lstlisting} \end{frame} \begin{frame}[fragile] - \frametitle{Example} + \frametitle{Remedy} + \begin{lstlisting} +from scipy import * + \end{lstlisting} +\alert{Now run python four\_plot.py again!} + \pause \begin{lstlisting} -In []: x, y = mgrid[-1:1:5j, -1:1:5j] -In []: z = x*x + y*y - -In []: z -Out[]: -array([[ 2. , 1.25, 1. , 1.25, 2. ], - [ 1.25, 0.5 , 0.25, 0.5 , 1.25], - [ 1. , 0.25, 0. , 0.25, 1. ], - [ 1.25, 0.5 , 0.25, 0.5 , 1.25], - [ 2. , 1.25, 1. , 1.25, 2. ]]) -\end{lstlisting} +Traceback (most recent call last): + File "four_plot.py", line 4, in <module> + plot(x, x, 'b') +NameError: name 'plot' is not defined + \end{lstlisting} \end{frame} \begin{frame}[fragile] - \myemph{\Large 2D data: \texttt{mlab.mesh}} - \vspace*{0.25in} - - \emphbar{\PythonCode{In []: mlab.mesh(x, y, z)}} - - \alert{Points needn't be regular} - - \vspace*{0.25in} -\begin{lstlisting} -In []: phi, theta = mgrid[0:pi:20j, -... 0:2*pi:20j] -In []: x = sin(phi)*cos(theta) -In []: y = sin(phi)*sin(theta) -In []: z = cos(phi) -In []: mlab.mesh(x, y, z, -... representation= -... 'wireframe') -\end{lstlisting} - + \frametitle{Remedy \ldots} + \begin{lstlisting} +from pylab import * + \end{lstlisting} +\alert{Now run python four\_plot.py again!!} \end{frame} \begin{frame}[fragile] - - \begin{columns} - \column{0.25\textwidth} - \myemph{\Large 3D data} - \column{0.5\textwidth} - \pgfimage[width=1.5in]{MEDIA/m2/mlab/contour3d}\\ - \end{columns} -\begin{lstlisting} -In []: x, y, z = mgrid[-5:5:64j, -... -5:5:64j, -... -5:5:64j] -In []: mlab.contour3d(x*x*0.5 + y*y + - z*z*2) -\end{lstlisting} + \frametitle{Modules} + \begin{itemize} + \item The \kwrd{import} keyword ``loads'' a module + \item One can also use: + \begin{lstlisting} +In []: from scipy import * +In []: from scipy import linspace + \end{lstlisting} + \item What is the difference? + \item \alert{Use the former only in interactive mode} + \end{itemize} \end{frame} \begin{frame}[fragile] + \frametitle{Package hierarchies} + \begin{lstlisting} +from scipy.integrate import odeint - \myemph{\Large 3D vector data: \PythonCode{mlab.quiver3d}} - \vspace*{0.25in} - - \pgfimage[width=2in]{MEDIA/m2/mlab/quiver3d_ex}\\ - -\begin{lstlisting} -In []: mlab.test_quiver3d() -\end{lstlisting} - -\emphbar{\PythonCode{obj = mlab.quiver3d(x, y, z, u, v, w)}} -\inctime{20} -\end{frame} - - -\subsection{Mayavi2} - -\begin{frame} - \frametitle{Introduction to Mayavi} - \begin{itemize} - \item Most scientists not interested in details of visualization - \item Visualization of data files with a nice UI - \item Interactive visualization of data (think Matlab) - \item Embedding visualizations in applications - \item Customization - \end{itemize} - \pause - \begin{block}{The Goal} - Provide a \alert{flexible} library/app for all of these needs! - \end{block} -\end{frame} - -\begin{frame} - {Overview of features} - \vspace*{-0.3in} - \begin{center} - \hspace*{-0.2in}\pgfimage[width=5in]{MEDIA/m2/m2_app3_3} - \end{center} +from scipy.optimize import fsolve + \end{lstlisting} \end{frame} - -\begin{frame} - \frametitle{Mayavi in applications} - \vspace*{-0.3in} - \begin{center} - \hspace*{-0.2in}\pgfimage[width=4.5in]{MEDIA/m2/m2_envisage} - \end{center} +\begin{frame}[fragile] + \frametitle{\texttt{from} \ldots \texttt{import} - conventional way!} + \small + \begin{lstlisting} +from scipy import linspace, pi, sin +from pylab import plot, legend, annotate +from pylab import xlim, ylim + +x = linspace(-5*pi, 5*pi, 500) +plot(x, x, 'b') +plot(x, -x, 'b') +plot(x, sin(x), 'g', linewidth=2) +plot(x, x*sin(x), 'r', linewidth=3) +legend(['x', '-x', 'sin(x)', 'xsin(x)']) +annotate('origin', xy = (0, 0)) +xlim(-5*pi, 5*pi) +ylim(-5*pi, 5*pi) + \end{lstlisting} \end{frame} -\begin{frame} - \frametitle{Live in your dialogs} - \vspace*{0.1in} - \begin{center} - \hspace*{-0.2in}\pgfimage[width=2.5in]{MEDIA/m2/mlab_tui} - \end{center} +\begin{frame}[fragile] + \frametitle{\texttt{from} \ldots \texttt{import} - conventional way!} + \small + \begin{lstlisting} +import scipy +import pylab + +x = scipy.linspace(-5*scipy.pi, 5*scipy.pi, 500) +pylab.plot(x, x, 'b') +pylab.plot(x, -x, 'b') +pylab.plot(x, scipy.sin(x), 'g', linewidth=2) +pylab.plot(x, x*scipy.sin(x), 'r', linewidth=3) +pylab.legend(['x', '-x', 'sin(x)', 'xsin(x)']) +pylab.annotate('origin', xy = (0, 0)) +pylab.xlim(-5*scipy.pi, 5*scipy.pi) +pylab.ylim(-5*scipy.pi, 5*scipy.pi) + \end{lstlisting} \end{frame} \begin{frame} - {Exploring the documentation} - \begin{center} - \pgfimage[width=4in]{MEDIA/m2/m2_ug_doc} - \end{center} + \frametitle{Modules: Standard library} + \begin{itemize} + \item Very powerful, ``Batteries included'' + \item Some standard modules: + \begin{itemize} + \item Math: \typ{math}, \typ{random} + \item Internet access: \typ{urllib2}, \typ{smtplib} + \item System, Command line arguments: \typ{sys} + \item Operating system interface: \typ{os} + \item Regular expressions: \typ{re} + \item Compression: \typ{gzip}, \typ{zipfile}, and \typ{tarfile} + \item And a whole lot more! + \end{itemize} + \item Check out the Python Library reference: + \url{http://docs.python.org/library/} + \end{itemize} +\inctime{5} \end{frame} +\begin{frame}[fragile] + \frametitle{Modules of special interest} + \begin{description}[matplotlibfor2d] + \item[\texttt{pylab}] Easy, interactive, 2D plotting -\begin{frame} - \frametitle{Summary} - \begin{itemize} - \item \url{http://code.enthought.com/projects/mayavi} - \item Uses VTK (\url{www.vtk.org}) - \item BSD license - \item Linux, win32 and Mac OS X - \item Highly scriptable - \item Embed in Traits UIs (wxPython and PyQt4) - \item Envisage Plugins - \item Debian/Ubuntu/Fedora - \item \alert{Pythonic} - \end{itemize} - - \inctime{10} + \item[\texttt{scipy}] arrays, statistics, optimization, integration, linear + algebra, Fourier transforms, signal and image processing, + genetic algorithms, ODE solvers, special functions, and more + \item[\texttt{Mayavi}] Easy, interactive, 3D plotting + \end{description} \end{frame} -\begin{frame} - {Getting hands dirty!} - - \begin{block}{Motivational problem} - Atmospheric data of temperature over the surface of the earth. - Let temperature ($T$) vary linearly with height ($z$): - \begin{center} - $T = 288.15 - 6.5z$ - \end{center} - \end{block} -\end{frame} +\section{Objects} +\begin{frame}{Everything is an Object!} + \begin{itemize} + \item \typ{int} + \item \typ{float} + \item \typ{str} + \item \typ{list} + \item \typ{tuple} + \item \typ{string} + \item \typ{dictionary} + \item \typ{function} + \item User defined class is also an object! + \end{itemize} +\end {frame} \begin{frame}[fragile] - \frametitle{Simple solution} - +\frametitle{Using Objects} + \begin{itemize} + \item Creating Objects + \begin{itemize} + \item Initialization + \end{itemize} \begin{lstlisting} -lat = linspace(-89, 89, 37) -lon = linspace(0, 360, 37) -z = linspace(0, 100, 11) +In []: a = str() + +In []: b = "Hello World" \end{lstlisting} -\pause - \begin{lstlisting} -x, y, z = mgrid[0:360:37j,-89:89:37j, - 0:100:11j] -t = 288.15 - 6.5*z -mlab.contour3d(x, y, z, t) -mlab.outline() -mlab.colorbar() + \item Object Manipulation + \begin{itemize} + \item Object methods + \item ``.'' operator + \end{itemize} + \begin{lstlisting} +In []: "Hello World".split() +Out[]: ['Hello', 'World'] \end{lstlisting} + \end{itemize} \end{frame} \begin{frame}[fragile] - \frametitle{Exercise: Lorenz equation} - \begin{columns} - \column{0.25\textwidth} - \begin{eqnarray*} - \frac{d x}{dt} &=& s (y-x)\\ - \frac{d y}{d t} &=& rx -y -xz\\ - \frac{d z}{d t} &=& xy - bz\\ - \end{eqnarray*} - \column{0.25\textwidth} - Let $s=10,$ - $r=28,$ - $b=8./3.$ - \end{columns} - \structure{\Large Region of interest} - \begin{lstlisting} -x, y, z = mgrid[-50:50:20j,-50:50:20j, - -10:60:20j] - \end{lstlisting} -\inctime{20} - -\end{frame} -\begin{frame}[fragile] - \frametitle{Solution} + \frametitle{Objects provide consistency} + \small \begin{lstlisting} -def lorenz(x,y,z,s=10.,r=28.,b=8./3.): - u = s*(y-x) - v = r*x-y-x*z - w = x*y-b*z - return u,v,w -x,y,z = mgrid [-50:50:20j,-50:50:20j, - -10:60:20j ] -u,v,w = lorenz( x , y , z ) -# Your plot here -# -mlab.show() - +for element in (1, 2, 3): + print element +for key in {'one':1, 'two':2}: + print key +for char in "123": + print char +for line in open("myfile.txt"): + print line +for line in urllib2.urlopen('http://site.com'): + print line \end{lstlisting} + \inctime{10} \end{frame} \begin{frame} - \frametitle{We have covered:} + \frametitle{What did we learn?} \begin{itemize} - \item Need of visualization. - \item Using mlab to create 3 D plots. - \item Mayavi Toolkit. + \item Functions: Default and Keyword arguments + \item Modules + \item Objects \end{itemize} \end{frame} -\end{document} - +\end{document}
\ No newline at end of file diff --git a/day2/session4.tex b/day2/session4.tex index 4f286c1..f401d8d 100644 --- a/day2/session4.tex +++ b/day2/session4.tex @@ -144,7 +144,7 @@ \begin{frame}[fragile] \frametitle{gcd revisited!} \begin{itemize} - \item Open gcd.py + \item Open \texttt{gcd.py} \end{itemize} \begin{lstlisting} def gcd(a, b): @@ -156,15 +156,15 @@ print gcd(16, 76) \end{lstlisting} \begin{itemize} - \item python gcd.py + \item \texttt{python gcd.py} \end{itemize} \end{frame} \begin{frame}[fragile] \frametitle{Find lcm using our gcd module} \begin{itemize} - \item Open lcm.py - \item $lcm = \frac{a*b}{gcd(a,b)}$ + \item Open \texttt{lcm.py} + \item $lcm = \frac{a * b}{gcd(a,b)}$ \end{itemize} \begin{lstlisting} from gcd import gcd @@ -174,7 +174,7 @@ print lcm(14, 56) \end{lstlisting} \begin{itemize} - \item python lcm.py + \item \texttt{python lcm.py} \end{itemize} \begin{lstlisting} 5 @@ -185,7 +185,7 @@ \begin{frame}[fragile] \frametitle{Writing stand-alone module} -Edit gcd.py file to: +Edit \texttt{gcd.py} file to: \begin{lstlisting} def gcd(a, b): if a % b == 0: @@ -197,22 +197,23 @@ Edit gcd.py file to: print gcd(16, 76) \end{lstlisting} \begin{itemize} - \item python gcd.py - \item python lcm.py + \item \texttt{python gcd.py} + \item \texttt{python lcm.py} \end{itemize} \end{frame} \begin{frame}[fragile] - \frametitle{More use of main} - For automating tests. + \frametitle{Automating tests} \begin{lstlisting} if __name__ == '__main__': for line in open('numbers.txt'): numbers = line.split() x = int(numbers[0]) y = int(numbers[1]) - result = (int(numbers[2])) - assert gcd(x, y) == result + result = int(numbers[2]) + if gcd(x, y) != result: + print "Failed gcd test + for", x, y \end{lstlisting} \end{frame} @@ -242,7 +243,7 @@ if __name__ == '__main__': \frametitle{Code Layout} \begin{itemize} \item Indentation - \item Tabs or Spaces?? + \item Tabs or Spaces? \item Maximum Line Length \item Blank Lines \item Encodings @@ -251,8 +252,8 @@ if __name__ == '__main__': \begin{frame}{Whitespaces in Expressions} \begin{itemize} - \item When to use extraneous whitespaces?? - \item When to avoid extra whitespaces?? + \item When to use extraneous whitespaces? + \item When to avoid extra whitespaces? \item Use one statement per line \end{itemize} \end{frame} @@ -319,43 +320,36 @@ or modulo by zero \end{frame} \begin{frame}[fragile] + \frametitle{Processing user input} + \begin{lstlisting} +prompt = 'Enter a number(Q to quit): ' + +a = raw_input(prompt) + +num = int(a) if a != 'Q' else 0 + \end{lstlisting} + \emphbar{What if the user enters some other alphabet?} +\end{frame} + + +\begin{frame}[fragile] \frametitle{Handling Exceptions} - Python uses \typ{try} and \typ{except} clause. - %%Revisiting the raw\_input + Python provides \typ{try} and \typ{except} clause. \begin{lstlisting} -a = raw_input('Enter number(Q to quit):') +prompt = 'Enter a number(Q to quit): ' + +a = raw_input(prompt) try: num = int(a) print num except: if a == 'Q': - print 'Exiting...' + print "Exiting ..." else: - print 'Wrong input!' - \end{lstlisting} - - + print "Wrong input ..." + \end{lstlisting} \end{frame} -%% \begin{frame}[fragile] -%% \frametitle{Solving it with \typ{try} and \typ{except}} -%% \vspace{-0.2in} -%% \begin{lstlisting} -%% highest = 0 -%% for record in open('sslc1.txt'): -%% fields = record.split(';') -%% try: -%% total = 0 -%% for score_str in fields[3:8]: -%% score = int(score_str) -%% total += score -%% if total > highest: -%% highest = total -%% except: -%% pass -%% print highest -%% \end{lstlisting} -%% \end{frame} \subsection{Strategy} \begin{frame}[fragile] \frametitle{Debugging effectively} @@ -421,6 +415,20 @@ savefig('science.png') \inctime{10} \end{frame} +\begin{frame} + \frametitle{Summary} +We have covered: + \begin{itemize} + \item Following and Resolving Error Messages. + \item Exceptions. + \item Handling exceptions + \item Approach for Debugging. +% \item Writting and running tests. + \end{itemize} +\end{frame} + +\end{document} + %% \begin{frame} %% \frametitle{Testing} @@ -433,130 +441,116 @@ savefig('science.png') %% \end{itemize} %% \end{frame} -\section{Test Driven Approach} -\begin{frame} - \frametitle{Need for Testing!} - - \begin{itemize} - \item Quality - \item Regression - \item Documentation - \end{itemize} - %% \vspace*{0.25in} - %% \emphbar{It is to assure that section of code is working as it is supposed to work} -\end{frame} - -\begin{frame}[fragile] - \frametitle{Example} - \begin{block}{Problem Statement} - Write a function to check whether a given input - string is a palindrome. - \end{block} -\end{frame} - -\begin{frame}[fragile] - \frametitle{Function: palindrome.py} -\begin{lstlisting} -def is_palindrome(input_str): - return input_str == input_str[::-1] -\end{lstlisting} -\end{frame} - -\begin{frame}[fragile] - \frametitle{Test for the palindrome: palindrome.py} -\begin{lstlisting} -def test_function_normal_words(): - input = "noon" - assert is_palindrome(input) == True - -if __name__ == "main'': - test_function_normal_words() -\end{lstlisting} -\end{frame} - -\begin{frame}[fragile] - \frametitle{Running the tests.} -\begin{lstlisting} -$ nosetests palindrome.py -. ----------------------------------------------- -Ran 1 test in 0.001s - -OK -\end{lstlisting} -\end{frame} - -\begin{frame}[fragile] - \frametitle{Exercise: Including new tests.} -\begin{lstlisting} -def test_function_ignore_cases_words(): - input = "Noon" - assert is_palindrome(input) == True -\end{lstlisting} - \vspace*{0.25in} - Check\\ - \PythonCode{$ nosetests palindrome.py} \\ - \begin{block}{Task} - Tweak the code to pass this test. - \end{block} -\end{frame} - -%\begin{frame}[fragile] -% \frametitle{Lets write some test!} -%\begin{lstlisting} -%#for form of equation y=mx+c -%#given m and c for two equation, -%#finding the intersection point. -%def intersect(m1,c1,m2,c2): -% x = (c2-c1)/(m1-m2) -% y = m1*x+c1 -% return (x,y) -%\end{lstlisting} -% -%Create a simple test for this -% -%function which will make it fail. -% -%\inctime{15} -%\end{frame} -% - -%% \begin{frame}[fragile] -%% \frametitle{Exercise} -%% Based on Euclid's algorithm: -%% \begin{center} -%% $gcd(a,b)=gcd(b,b\%a)$ -%% \end{center} -%% gcd function can be written as: -%% \begin{lstlisting} -%% def gcd(a, b): -%% if a%b == 0: return b -%% return gcd(b, a%b) -%% \end{lstlisting} -%% \vspace*{-0.15in} -%% \begin{block}{Task} -%% \begin{itemize} -%% \item Write at least -%% two tests for above mentioned function. -%% \item Write a non recursive implementation -%% of gcd(), and test it using already -%% written tests. -%% \end{itemize} -%% \end{block} - -%% \inctime{15} -%% \end{frame} - -\begin{frame} - \frametitle{Summary} -We have coverd: - \begin{itemize} - \item Following and Resolving Error Messages. - \item Exceptions. - \item Handling exceptions - \item Approach for Debugging. - \item Writting and running tests. - \end{itemize} -\end{frame} - -\end{document} +% \section{Test Driven Approach} +% \begin{frame} +% \frametitle{Need for Testing!} +% +% \begin{itemize} +% \item Quality +% \item Regression +% \item Documentation +% \end{itemize} +% %% \vspace*{0.25in} +% %% \emphbar{It is to assure that section of code is working as it is supposed to work} +% \end{frame} +% +% \begin{frame}[fragile] +% \frametitle{Example} +% \begin{block}{Problem Statement} +% Write a function to check whether a given input +% string is a palindrome. +% \end{block} +% \end{frame} +% +% \begin{frame}[fragile] +% \frametitle{Function: palindrome.py} +% \begin{lstlisting} +% def is_palindrome(input_str): +% return input_str == input_str[::-1] +% \end{lstlisting} +% \end{frame} +% +% \begin{frame}[fragile] +% \frametitle{Test for the palindrome: palindrome.py} +% \begin{lstlisting} +% def test_function_normal_words(): +% input = "noon" +% assert is_palindrome(input) == True +% +% if __name__ == "main'': +% test_function_normal_words() +% \end{lstlisting} +% \end{frame} +% +% \begin{frame}[fragile] +% \frametitle{Running the tests.} +% \begin{lstlisting} +% $ nosetests palindrome.py +% . +% ---------------------------------------------- +% Ran 1 test in 0.001s +% +% OK +% \end{lstlisting} +% \end{frame} +% +% \begin{frame}[fragile] +% \frametitle{Exercise: Including new tests.} +% \begin{lstlisting} +% def test_function_ignore_cases_words(): +% input = "Noon" +% assert is_palindrome(input) == True +% \end{lstlisting} +% \vspace*{0.25in} +% Check\\ +% \PythonCode{$ nosetests palindrome.py} \\ +% \begin{block}{Task} +% Tweak the code to pass this test. +% \end{block} +% \end{frame} +% +% %\begin{frame}[fragile] +% % \frametitle{Lets write some test!} +% %\begin{lstlisting} +% %#for form of equation y=mx+c +% %#given m and c for two equation, +% %#finding the intersection point. +% %def intersect(m1,c1,m2,c2): +% % x = (c2-c1)/(m1-m2) +% % y = m1*x+c1 +% % return (x,y) +% %\end{lstlisting} +% % +% %Create a simple test for this +% % +% %function which will make it fail. +% % +% %\inctime{15} +% %\end{frame} +% % +% +% %% \begin{frame}[fragile] +% %% \frametitle{Exercise} +% %% Based on Euclid's algorithm: +% %% \begin{center} +% %% $gcd(a,b)=gcd(b,b\%a)$ +% %% \end{center} +% %% gcd function can be written as: +% %% \begin{lstlisting} +% %% def gcd(a, b): +% %% if a%b == 0: return b +% %% return gcd(b, a%b) +% %% \end{lstlisting} +% %% \vspace*{-0.15in} +% %% \begin{block}{Task} +% %% \begin{itemize} +% %% \item Write at least +% %% two tests for above mentioned function. +% %% \item Write a non recursive implementation +% %% of gcd(), and test it using already +% %% written tests. +% %% \end{itemize} +% %% \end{block} +% +% %% \inctime{15} +% %% \end{frame}
\ No newline at end of file diff --git a/day2/session5.tex b/day2/session5.tex index 0b13e9a..a3e63e8 100644 --- a/day2/session5.tex +++ b/day2/session5.tex @@ -1,33 +1,48 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%Tutorial slides on Python. +% Tutorial slides on Python. % % Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in> % Copyright (c) 2005-2009, Prabhu Ramachandran %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\documentclass[14pt,compress]{beamer} -%\documentclass[draft]{beamer} -%\documentclass[compress,handout]{beamer} -%\usepackage{pgfpages} -%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm] +\documentclass[compress,14pt]{beamer} +% \documentclass[handout]{beamer} +% \usepackage{pgfpages} +% \pgfpagesuselayout{4 on 1}[a4paper,border, shrink=5mm,landscape] +\usepackage{tikz} +\newcommand{\hyperlinkmovie}{} +%\usepackage{movie15} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Note that in presentation mode +% \paperwidth 364.19536pt +% \paperheight 273.14662pt +% h/w = 0.888 + -% Modified from: generic-ornate-15min-45min.de.tex \mode<presentation> { \usetheme{Warsaw} + %\usetheme{Boadilla} + %\usetheme{default} \useoutertheme{infolines} \setbeamercovered{transparent} } +% To remove navigation symbols +\setbeamertemplate{navigation symbols}{} + +\usepackage{amsmath} \usepackage[english]{babel} \usepackage[latin1]{inputenc} -%\usepackage{times} +\usepackage{times} \usepackage[T1]{fontenc} % Taken from Fernando's slides. \usepackage{ae,aecompl} \usepackage{mathpazo,courier,euler} \usepackage[scaled=.95]{helvet} +\usepackage{pgf} \definecolor{darkgreen}{rgb}{0,0.5,0} @@ -40,50 +55,65 @@ keywordstyle=\color{blue}\bfseries} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% Macros +% My Macros +\setbeamercolor{postit}{bg=yellow,fg=black} \setbeamercolor{emphbar}{bg=blue!20, fg=black} \newcommand{\emphbar}[1] {\begin{beamercolorbox}[rounded=true]{emphbar} {#1} \end{beamercolorbox} } +%{\centerline{\fcolorbox{gray!50} {blue!10}{ +%\begin{minipage}{0.9\linewidth} +% {#1} +%\end{minipage} +% }}} + +\newcommand{\myemph}[1]{\structure{\emph{#1}}} +\newcommand{\PythonCode}[1]{\lstinline{#1}} + +\newcommand{\tvtk}{\texttt{tvtk}} +\newcommand{\mlab}{\texttt{mlab}} + \newcounter{time} \setcounter{time}{0} -\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}} +\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\vspace*{0.1in}\tiny \thetime\ m}} -\newcommand{\typ}[1]{\texttt{#1}} +\newcommand\BackgroundPicture[1]{% + \setbeamertemplate{background}{% + \parbox[c][\paperheight]{\paperwidth}{% + \vfill \hfill + \hfill \vfill +}}} -\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}} } +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +% Configuring the theme +%\setbeamercolor{normal text}{fg=white} +%\setbeamercolor{background canvas}{bg=black} -%%% This is from Fernando's setup. -% \usepackage{color} -% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2} -% % Use and configure listings package for nicely formatted code -% \usepackage{listings} -% \lstset{ -% language=Python, -% basicstyle=\small\ttfamily, -% commentstyle=\ttfamily\color{blue}, -% stringstyle=\ttfamily\color{orange}, -% showstringspaces=false, -% breaklines=true, -% postbreak = \space\dots -% } %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Title page -\title[Exercises]{Exercises} +\title[3D Plotting]{3D data Visualization} \author[FOSSEE] {FOSSEE} \institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay} \date[] {8 November, 2009\\Day 2, Session 5} + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo} -%\logo{\pgfuseimage{iitmlogo}} +%\pgfdeclareimage[height=0.75cm]{iitblogo}{iitblogo} +%\logo{\pgfuseimage{iitblogo}} +\AtBeginSection[] +{ + \begin{frame}<beamer> + \frametitle{Outline} + \tableofcontents[currentsection,currentsubsection] + \end{frame} +} %% Delete this, if you do not want the table of contents to pop up at %% the beginning of each subsection: @@ -95,109 +125,434 @@ \end{frame} } - -% If you wish to uncover everything in a step-wise fashion, uncomment -% the following command: -%\beamerdefaultoverlayspecification{<+->} - -%\includeonlyframes{current,current1,current2,current3,current4,current5,current6} - +\AtBeginSection[] +{ + \begin{frame}<beamer> + \frametitle{Outline} + \tableofcontents[currentsection,currentsubsection] + \end{frame} +} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % DOCUMENT STARTS \begin{document} \begin{frame} - \titlepage + \maketitle \end{frame} -\begin{frame}{Problem 1.1} - The aliquot of a number is defined as: the sum of the \emph{proper} divisors of the number. \\For example: -\center{aliquot(12) = 1 + 2 + 3 + 4 + 6 = 16.}\\ - Write a function that returns the aliquot number of a given number. +\begin{frame} + \frametitle{Outline} + \tableofcontents + % You might wish to add the option [pausesections] \end{frame} -\begin{frame}{Problem 1.2} - Pair of numbers (a, b) is said to be \alert{amicable} if aliquot number of a is b and aliquot number of b is a.\\ - Example: \texttt{220, 284}\\ - Write a program that prints all four digit amicable pairs. - -\inctime{20} +\section{3D Data Visualization} + +\begin{frame} + \frametitle{What is visualization?} + \Large + \begin{center} + Visual representation of data + \end{center} \end{frame} -%% \begin{frame}{Problem 2} -%% Given an empty chessboard and one Bishop placed in any s%quare, say (r, c), generate the list of all squares the Bi%shop could move to. +%% \begin{frame} +%% \frametitle{Is this new?} +%% \begin{center} +%% We have moved from: +%% \end{center} +%% \begin{columns} +%% \column{} +%% \hspace*{-1in} +%% \includegraphics[width=1.75in,height=1.75in, interpolate=true]{data/3832} +%% \column{}\hspace*{-0.25in} +%% To +%% \column{} +%% \hspace*{-1in} +%% \includegraphics[width=1.75in, height=1.75in, interpolate=true]{data/torus} +%% \end{columns} %% \end{frame} -\begin{frame}[fragile] - \frametitle{Problem Set 2} - Given a string like, ``1, 3-7, 12, 15, 18-21'', produce the list \\ - \begin{lstlisting} - [1,3,4,5,6,7,12,15,18,19,20,21] - \end{lstlisting} +\begin{frame} + \frametitle{3D visualization} + \Large + \begin{center} + Harder but important + \end{center} +\end{frame} + +\begin{frame} + \frametitle{Is this Graphics?} + \Large + \begin{center} + Visualization is about data! + \end{center} +\end{frame} + +\begin{frame} + \frametitle{Examples: trajectory in space} + \Large + \begin{center} + \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex} + \end{center} +\end{frame} + +\begin{frame} + \frametitle{Examples: Fire in a room} + \Large + \begin{center} + Demo of data + \end{center} \inctime{10} \end{frame} -\begin{frame} - \frametitle{Problem Set 3} - \begin{description} - \item[3.1] Count word frequencies in a file. -\end{description} -\inctime{5} +\section{Tools available} + +\subsection{mlab} + +\begin{frame} + {Overview} + \Large + \begin{itemize} + \item Simple + \item Convenient + \item Full-featured + \end{itemize} +\end{frame} + +\begin{frame}[fragile] + + \frametitle{Getting started} + \myemph{\Large Vanilla:} + \begin{lstlisting}[language=bash] + $ ipython -wthread + \end{lstlisting} + \myemph{\Large with Pylab:} + \begin{lstlisting}[language=bash] + $ ipython -pylab -wthread + \end{lstlisting} \end{frame} \begin{frame}[fragile] - \frametitle{Problem set 4} - Central difference - \begin{equation*} - \frac{sin(x+h)-sin(x-h)}{2h} - \end{equation*} + \frametitle{Using mlab} + + \begin{lstlisting} +In []:from enthought.mayavi import mlab + \end{lstlisting} + + \vspace*{0.5in} + + \myemph{\Large Try these} + + \vspace*{0.25in} + + \begin{lstlisting} +In []: mlab.test_<TAB> +In []: mlab.test_contour3d() +In []: mlab.test_contour3d?? + \end{lstlisting} +\end{frame} + +\begin{frame} + {Exploring the view} + \begin{columns} + \column{0.6\textwidth} + \pgfimage[width=3in]{MEDIA/m2/contour3d} + \column{0.4\textwidth} + \begin{itemize} + \item Mouse + \item Keyboard + \item Toolbar + \item Mayavi icon\pgfimage[width=0.2in]{MEDIA/m2/m2_icon} + \end{itemize} + \end{columns} +\end{frame} + +\begin{frame}[fragile] + \frametitle{\mlab\ plotting functions} + \begin{columns} + \column{0.25\textwidth} + \myemph{\Large 0D data} + \column{0.5\textwidth} + \pgfimage[width=2in]{MEDIA/m2/mlab/points3d_ex} + \end{columns} + + \begin{lstlisting} +In []: t = linspace(0, 2*pi, 50) +In []: u = cos(t) * pi +In []: x, y, z = sin(u), cos(u), sin(t) + \end{lstlisting} + \emphbar{\PythonCode{In []: mlab.points3d(x, y, z)}} +\end{frame} + +\begin{frame} + \begin{columns} + \column{0.25\textwidth} + \myemph{\Large 1D data} + \column{0.5\textwidth} + \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex} + \end{columns} + \emphbar{\PythonCode{In []: mlab.plot3d(x, y, z, t)}} + + Plots lines between the points + +\end{frame} + +\begin{frame}[fragile] + \begin{columns} + \column{0.25\textwidth} + \myemph{\Large 2D data} + \column{0.5\textwidth} + \pgfimage[width=2in]{MEDIA/m2/mlab/surf_ex} + \end{columns} + \begin{lstlisting} +In []: x, y = mgrid[-3:3:100j,-3:3:100j] +In []: z = sin(x*x + y*y) + \end{lstlisting} + + \emphbar{\PythonCode{In []: mlab.surf(x, y, z)}} + + \alert{Assumes the points are rectilinear} + +\end{frame} + +\begin{frame}[fragile] + \frametitle{mgrid} \begin{lstlisting} - In []: x = linspace(0, 2*pi, 100) - In []: y = sin(x) - In []: deltax = x[1] - x[0] - \end{lstlisting} - \pause - \begin{enumerate} - \item Given this, get the finite difference of sin in the range 0 to 2*pi - \end{enumerate} +In []: mgrid[0:3,0:3] +Out[]: +array([[[0, 0, 0], + [1, 1, 1], + [2, 2, 2]], + + [[0, 1, 2], + [0, 1, 2], + [0, 1, 2]]]) + +In []: mgrid[-1:1:5j] +Out[]: array([-1., -0.5, 0., 0.5, 1.]) +\end{lstlisting} +\end{frame} + +\begin{frame}[fragile] + \frametitle{Example} + \begin{lstlisting} +In []: x, y = mgrid[-1:1:5j, -1:1:5j] +In []: z = x*x + y*y + +In []: z +Out[]: +array([[ 2. , 1.25, 1. , 1.25, 2. ], + [ 1.25, 0.5 , 0.25, 0.5 , 1.25], + [ 1. , 0.25, 0. , 0.25, 1. ], + [ 1.25, 0.5 , 0.25, 0.5 , 1.25], + [ 2. , 1.25, 1. , 1.25, 2. ]]) +\end{lstlisting} +\end{frame} + +\begin{frame}[fragile] + \myemph{\Large 2D data: \texttt{mlab.mesh}} + \vspace*{0.25in} + + \emphbar{\PythonCode{In []: mlab.mesh(x, y, z)}} + + \alert{Points needn't be regular} + + \vspace*{0.25in} +\begin{lstlisting} +In []: phi, theta = mgrid[0:pi:20j, +... 0:2*pi:20j] +In []: x = sin(phi)*cos(theta) +In []: y = sin(phi)*sin(theta) +In []: z = cos(phi) +In []: mlab.mesh(x, y, z, +... representation= +... 'wireframe') +\end{lstlisting} + \end{frame} +\begin{frame}[fragile] + + \begin{columns} + \column{0.25\textwidth} + \myemph{\Large 3D data} + \column{0.5\textwidth} + \pgfimage[width=1.5in]{MEDIA/m2/mlab/contour3d}\\ + \end{columns} +\begin{lstlisting} +In []: x, y, z = mgrid[-5:5:64j, +... -5:5:64j, +... -5:5:64j] +In []: mlab.contour3d(x*x*0.5 + y*y + + z*z*2) +\end{lstlisting} +\end{frame} + +\begin{frame}[fragile] + + \myemph{\Large 3D vector data: \PythonCode{mlab.quiver3d}} + \vspace*{0.25in} + + \pgfimage[width=2in]{MEDIA/m2/mlab/quiver3d_ex}\\ + +\begin{lstlisting} +In []: mlab.test_quiver3d() +\end{lstlisting} + +\emphbar{\PythonCode{obj = mlab.quiver3d(x, y, z, u, v, w)}} +\inctime{20} +\end{frame} + + +\subsection{Mayavi2} + \begin{frame} - \frametitle{Problem Set 5} + \frametitle{Introduction to Mayavi} \begin{itemize} - \item[5.1] Write a function that plots any regular n-gon given \typ{n}. - \item[5.2] Consider the logistic map, $f(x) = kx(1-x)$, plot it for - $k=2.5, 3.5$ and $4$ in the same plot. -\end{itemize} + \item Most scientists not interested in details of visualization + \item Visualization of data files with a nice UI + \item Interactive visualization of data (think Matlab) + \item Embedding visualizations in applications + \item Customization + \end{itemize} + \pause + \begin{block}{The Goal} + Provide a \alert{flexible} library/app for all of these needs! + \end{block} \end{frame} -\begin{frame}[fragile] -\frametitle{Problem Set 5} - \begin{columns} - \column{0.6\textwidth} - \small{ - \begin{itemize} - \item[3] Consider the iteration $x_{n+1} = f(x_n)$ where $f(x) = kx(1-x)$. Plot the successive iterates of this process as explained below. - \end{itemize}} - \column{0.35\textwidth} - \hspace*{-0.5in} - \includegraphics[height=1.6in, interpolate=true]{data/cobweb} -\end{columns} -\end{frame} - -\begin{frame} - \frametitle{Problem Set 5.3} - Plot the cobweb plot as follows: - \begin{enumerate} - \item Start at $(x_0, 0)$ ($\implies$ i=0) - \item Draw a line to $(x_i, f(x_i))$ - \item Set $x_{i+1} = f(x_i)$ - \item Draw a line to $(x_{i+1}, x_{i+1})$ - \item $(i\implies i+1)$ - \item Repeat from 2 for as long as you want - \end{enumerate} +\begin{frame} + {Overview of features} + \vspace*{-0.3in} + \begin{center} + \hspace*{-0.2in}\pgfimage[width=5in]{MEDIA/m2/m2_app3_3} + \end{center} +\end{frame} + + +\begin{frame} + \frametitle{Mayavi in applications} + \vspace*{-0.3in} + \begin{center} + \hspace*{-0.2in}\pgfimage[width=4.5in]{MEDIA/m2/m2_envisage} + \end{center} +\end{frame} + +\begin{frame} + \frametitle{Live in your dialogs} + \vspace*{0.1in} + \begin{center} + \hspace*{-0.2in}\pgfimage[width=2.5in]{MEDIA/m2/mlab_tui} + \end{center} +\end{frame} + +\begin{frame} + {Exploring the documentation} + \begin{center} + \pgfimage[width=4in]{MEDIA/m2/m2_ug_doc} + \end{center} +\end{frame} + + +\begin{frame} + \frametitle{Summary} + \begin{itemize} + \item \url{http://code.enthought.com/projects/mayavi} + \item Uses VTK (\url{www.vtk.org}) + \item BSD license + \item Linux, win32 and Mac OS X + \item Highly scriptable + \item Embed in Traits UIs (wxPython and PyQt4) + \item Envisage Plugins + \item Debian/Ubuntu/Fedora + \item \alert{Pythonic} + \end{itemize} + + \inctime{10} + +\end{frame} + +\begin{frame} + {Getting hands dirty!} + + \begin{block}{Motivational problem} + Atmospheric data of temperature over the surface of the earth. + Let temperature ($T$) vary linearly with height ($z$): + \begin{center} + $T = 288.15 - 6.5z$ + \end{center} + \end{block} +\end{frame} + +\begin{frame}[fragile] + \frametitle{Simple solution} + + \begin{lstlisting} +lat = linspace(-89, 89, 37) +lon = linspace(0, 360, 37) +z = linspace(0, 100, 11) + \end{lstlisting} +\pause + \begin{lstlisting} +x, y, z = mgrid[0:360:37j,-89:89:37j, + 0:100:11j] +t = 288.15 - 6.5*z +mlab.contour3d(x, y, z, t) +mlab.outline() +mlab.colorbar() + \end{lstlisting} +\end{frame} + +\begin{frame}[fragile] + \frametitle{Exercise: Lorenz equation} + \begin{columns} + \column{0.25\textwidth} + \begin{eqnarray*} + \frac{d x}{dt} &=& s (y-x)\\ + \frac{d y}{d t} &=& rx -y -xz\\ + \frac{d z}{d t} &=& xy - bz\\ + \end{eqnarray*} + \column{0.25\textwidth} + Let $s=10,$ + $r=28,$ + $b=8./3.$ + \end{columns} + \structure{\Large Region of interest} + \begin{lstlisting} +x, y, z = mgrid[-50:50:20j,-50:50:20j, + -10:60:20j] + \end{lstlisting} \inctime{20} + +\end{frame} +\begin{frame}[fragile] + \frametitle{Solution} + \begin{lstlisting} +def lorenz(x,y,z,s=10.,r=28.,b=8./3.): + u = s*(y-x) + v = r*x-y-x*z + w = x*y-b*z + return u,v,w +x,y,z = mgrid [-50:50:20j,-50:50:20j, + -10:60:20j ] +u,v,w = lorenz( x , y , z ) +# Your plot here +# +mlab.show() + + \end{lstlisting} +\end{frame} + +\begin{frame} + \frametitle{We have covered:} + \begin{itemize} + \item Need of visualization. + \item Using mlab to create 3 D plots. + \item Mayavi Toolkit. + \end{itemize} \end{frame} \end{document} + |