diff options
| -rw-r--r-- | day1/data/plot1.png | bin | 0 -> 12830 bytes | |||
| -rw-r--r-- | day1/data/plot10.png | bin | 0 -> 25192 bytes | |||
| -rw-r--r-- | day1/data/plot11.png | bin | 0 -> 24874 bytes | |||
| -rw-r--r-- | day1/data/plot2.png | bin | 0 -> 16342 bytes | |||
| -rw-r--r-- | day1/data/plot3.png | bin | 0 -> 18972 bytes | |||
| -rw-r--r-- | day1/data/plot4.png | bin | 0 -> 21119 bytes | |||
| -rw-r--r-- | day1/data/plot5.png | bin | 0 -> 22972 bytes | |||
| -rw-r--r-- | day1/data/plot6.png | bin | 0 -> 24653 bytes | |||
| -rw-r--r-- | day1/data/plot7.png | bin | 0 -> 25381 bytes | |||
| -rw-r--r-- | day1/data/plot8.png | bin | 0 -> 25165 bytes | |||
| -rw-r--r-- | day1/data/plot9.png | bin | 0 -> 25054 bytes | |||
| -rw-r--r-- | day1/exercises.tex | 187 | ||||
| -rw-r--r-- | day2/session3.tex | 697 |
13 files changed, 594 insertions, 290 deletions
diff --git a/day1/data/plot1.png b/day1/data/plot1.png Binary files differnew file mode 100644 index 0000000..6ebbd3d --- /dev/null +++ b/day1/data/plot1.png diff --git a/day1/data/plot10.png b/day1/data/plot10.png Binary files differnew file mode 100644 index 0000000..b1571f7 --- /dev/null +++ b/day1/data/plot10.png diff --git a/day1/data/plot11.png b/day1/data/plot11.png Binary files differnew file mode 100644 index 0000000..caa1348 --- /dev/null +++ b/day1/data/plot11.png diff --git a/day1/data/plot2.png b/day1/data/plot2.png Binary files differnew file mode 100644 index 0000000..34e0ebb --- /dev/null +++ b/day1/data/plot2.png diff --git a/day1/data/plot3.png b/day1/data/plot3.png Binary files differnew file mode 100644 index 0000000..35ed5fc --- /dev/null +++ b/day1/data/plot3.png diff --git a/day1/data/plot4.png b/day1/data/plot4.png Binary files differnew file mode 100644 index 0000000..11007a1 --- /dev/null +++ b/day1/data/plot4.png diff --git a/day1/data/plot5.png b/day1/data/plot5.png Binary files differnew file mode 100644 index 0000000..19dde13 --- /dev/null +++ b/day1/data/plot5.png diff --git a/day1/data/plot6.png b/day1/data/plot6.png Binary files differnew file mode 100644 index 0000000..297822c --- /dev/null +++ b/day1/data/plot6.png diff --git a/day1/data/plot7.png b/day1/data/plot7.png Binary files differnew file mode 100644 index 0000000..eb0d68a --- /dev/null +++ b/day1/data/plot7.png diff --git a/day1/data/plot8.png b/day1/data/plot8.png Binary files differnew file mode 100644 index 0000000..e144853 --- /dev/null +++ b/day1/data/plot8.png diff --git a/day1/data/plot9.png b/day1/data/plot9.png Binary files differnew file mode 100644 index 0000000..dd64ef1 --- /dev/null +++ b/day1/data/plot9.png diff --git a/day1/exercises.tex b/day1/exercises.tex new file mode 100644 index 0000000..da88e05 --- /dev/null +++ b/day1/exercises.tex @@ -0,0 +1,187 @@ +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%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[] {7 November, 2009\\Day 1, Session 5} +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +%\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} + \frametitle{Problem 1} + \begin{itemize} + \item Open file 'pos.txt', it has X and Y Coordinate of a particle under motion + \item Plot X vs Y Graph. + \item Label both the axes. + \item What kind of motion is this? + \item Title the graph accordingly. + \item Annotate the position where vertical velocity is zero. + \end{itemize} +\end{frame} + +\begin{frame} + \frametitle{Problem 2} +Write a Program that plots a regular n-gon(Let n = 5). +\end{frame} + +\begin{frame}[fragile] + \frametitle{Problem 3} +Create a sequence of images in which the damped oscillator($e^{x/10}sin(x)$) slowly evolves over time. +\begin{columns} +\column{0.35\textwidth} +\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot2} +\column{0.35\textwidth} +\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot4} +\column{0.35\textwidth} +\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot6} +\end{columns} +\begin{block}{Hint} +\small + \begin{lstlisting} +savefig('plot'+str(i)+'.png') #i is int variable + \end{lstlisting} +\end{block} +\end{frame} + +\begin{frame} + \frametitle{Problem 4} + Legendre polynomials $P_n(x)$ are defined by the following recurrence relation + +\center{$(n+1)P_{n+1}(x) - (2n+1)xP_n(x) + nP_{n-1}(x) = 0$}\\ + +with $P_0(x) = 1$, $P_1(x) = x$ and $P_2(x) = (3x^2 - 1)/2$. Compute the next three + Legendre polynomials and plot all 6 over the interval [-1,1]. +\end{frame} +\end{document} + +%% \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} diff --git a/day2/session3.tex b/day2/session3.tex index d224ef3..af621e4 100644 --- a/day2/session3.tex +++ b/day2/session3.tex @@ -1,33 +1,48 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%Tutorial slides on Python. +% Tutorial slides on Python. % -% Author: FOSSEE -% Copyright (c) 2009, FOSSEE, IIT Bombay +% 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[Basic Python]{Python language: Data structures and functions} +\title[3D Plotting]{3D data Visualization} -\author[FOSSEE Team] {The FOSSEE Group} +\author[FOSSEE] {FOSSEE} \institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay} -\date[] {8 November, 2009\\Day 2, Session 3} +\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: @@ -102,19 +132,12 @@ \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 + \maketitle \end{frame} \begin{frame} @@ -123,325 +146,419 @@ % You might wish to add the option [pausesections] \end{frame} -\section{Functions} -\subsection{Default arguments} -\begin{frame}[fragile] - \frametitle{Functions: default arguments} - \begin{lstlisting} -In []: greet = 'hello world' +\section{3D Data Visualization} -In []: greet.split() -Out[]: ['hello', 'world'] +\begin{frame} + \frametitle{What is visualization?} + \Large + \begin{center} + Visual representation of data + \end{center} +\end{frame} -In []: line = 'Rossum, Guido, 54, 46, 55' -In []: line.split(',') -Out[]: ['Rossum', ' Guido', ' 54', - ' 46', ' 55'] - \end{lstlisting} +%% \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} + \frametitle{3D visualization} + \Large + \begin{center} + Harder but important + \end{center} \end{frame} -\begin{frame}[fragile] - \frametitle{Functions: default arguments \ldots} - \small - \begin{lstlisting} -def ask_ok(prompt, complaint='Yes or no!'): - while True: - ok = raw_input(prompt) - if ok in ('y', 'ye', 'yes'): - return True - if ok in ('n', 'no', 'nop', - 'nope'): - return False - print complaint - -ask_ok('?') -ask_ok('?', '[Y/N]') - \end{lstlisting} -\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} + +\section{Tools available} + +\subsection{mlab} + +\begin{frame} + {Overview} + \Large + \begin{itemize} + \item Simple + \item Convenient + \item Full-featured + \end{itemize} +\end{frame} -\subsection{Keyword arguments} \begin{frame}[fragile] - \frametitle{Functions: Keyword arguments} -We have seen the following - \begin{lstlisting} -In []: legend(['sin(2y)'], - loc='center') -In []: plot(y, sin(y), 'g', - linewidth=2) -In []: annotate('local max', - xy=(1.5, 1)) -In []: pie(science.values(), - labels=science.keys()) - \end{lstlisting} + + \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{Functions: keyword arguments \ldots} - \small - \begin{lstlisting} -def ask_ok(prompt, complaint='Yes or no!'): - while True: - ok = raw_input(prompt) - if ok in ('y', 'ye', 'yes'): - return True - if ok in ('n', 'no', 'nop', - 'nope'): - return False - print complaint - -ask_ok(prompt='?') -ask_ok(prompt='?', complaint='[y/n]') -ask_ok(complaint='[y/n]', prompt='?') -\end{lstlisting} + \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} -\subsection{Built-in functions} \begin{frame} - {Before writing a function} - \begin{itemize} - \item Variety of builtin 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} + {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} -\subsection{Exercises} -\begin{frame}{Problem set 3: Problem 3.1} - Write a function to return the gcd of two numbers. +\begin{frame}[fragile] + \frametitle{Clearing the plot area} + \PythonCode{In []: mlab.clf()} \end{frame} -\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. +\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}{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} +\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 + \end{frame} -\section{Modules} \begin{frame}[fragile] - \frametitle{\typ{from} \ldots \typ{import} magic} - \begin{lstlisting} -from scipy.interpolate import splrep -from scipy.interpolate import splev + \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} -from scipy.integrate import quad -from scipy.integrate import odeint + \emphbar{\PythonCode{In []: mlab.surf(x, y, z)}} + + \alert{Assumes the points are rectilinear} -from scipy.optimize import fsolve - \end{lstlisting} -\emphbar{All the above statements import one function into your namespace} \end{frame} \begin{frame}[fragile] - \frametitle{Running scripts from command line} - \small - \begin{itemize} - \item Start cmd - \item cd to Desktop - \item python sine\_plot.py - \end{itemize} - \pause + \frametitle{mgrid} \begin{lstlisting} -Traceback (most recent call last): - File "sine_plot.py", line 1, in <module> - x = linspace(-5*pi, 5*pi, 500) -NameError: name 'linspace' is not defined - \end{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} \end{frame} \begin{frame}[fragile] - \frametitle{Remedy} - \emphbar{Adding what lines to sine\_plot.py makes this program work?} - \begin{lstlisting} -from scipy import * - \end{lstlisting} -\alert{Now run python sine\_plot.py again!} - \pause + \frametitle{Example} \begin{lstlisting} -Traceback (most recent call last): - File "sine_plot.py", line 4, in <module> - plot(x, x, 'b') -NameError: name 'plot' is not defined - \end{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] - \frametitle{Remedy \ldots} - \emphbar{What should we add now?} - \begin{lstlisting} -from pylab import * - \end{lstlisting} -\alert{Now run python sine\_plot.py again!!} + \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] - \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} + + \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] - \frametitle{Package hierarchies} - \begin{lstlisting} -from scipy.interpolate import splev -from scipy.integrate import quad + \myemph{\Large 3D vector data: \PythonCode{mlab.quiver3d}} + \vspace*{0.25in} -from scipy.optimize import fsolve - \end{lstlisting} -\end{frame} + \pgfimage[width=2in]{MEDIA/m2/mlab/quiver3d_ex}\\ + +\begin{lstlisting} +In []: mlab.test_quiver3d() +\end{lstlisting} -\begin{frame}[fragile] - \frametitle{from \ldots import in a 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} +\emphbar{\PythonCode{obj = mlab.quiver3d(x, y, z, u, v, w)}} +\inctime{20} \end{frame} -\begin{frame}[fragile] - \frametitle{from \ldots import in a 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} + +\subsection{Mayavi2} \begin{frame} - \frametitle{Modules: Standard library} + \frametitle{Introduction to Mayavi} \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/} + \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} -\inctime{5} + \pause + \begin{block}{The Goal} + Provide a \alert{flexible} library/app for all of these needs! + \end{block} \end{frame} -\begin{frame}[fragile] - \frametitle{Modules of special interest} - \begin{description}[matplotlibfor2d] - \item[\typ{pylab}] Easy, interactive, 2D plotting +\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} - \item[\typ{scipy}] arrays, statistics, optimization, integration, linear - algebra, Fourier transforms, signal and image processing, - genetic algorithms, ODE solvers, special functions, and more - \item[\typ{Mayavi}] Easy, interactive, 3D plotting - \end{description} +\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} -\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} + \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{Using Objects} - \begin{itemize} - \item Creating Objects - \begin{itemize} - \item Initialization - \end{itemize} - \begin{lstlisting} -In []: a = str() + \frametitle{Simple solution} -In []: b = "Hello World" + \begin{lstlisting} +lat = linspace(-89, 89, 37) +lon = linspace(0, 360, 37) +z = linspace(0, 100, 11) \end{lstlisting} - \item Object Manipulation - \begin{itemize} - \item Object methods - \item ``.'' operator - \end{itemize} - \begin{lstlisting} -In []: "Hello World".split() -Out[]: ['Hello', 'World'] +\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{itemize} \end{frame} \begin{frame}[fragile] - \frametitle{Objects provide consistency} - \small + \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} -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 +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} - \inctime{10} \end{frame} \begin{frame} - \frametitle{What did we learn?} + \frametitle{We have covered:} \begin{itemize} - \item Functions: Default and Keyword arguments - \item Modules - \item Objects + \item Need of visualization. + \item Using mlab to create 3 D plots. + \item Mayavi Toolkit. \end{itemize} \end{frame} -\end{document}
\ No newline at end of file +\end{document} + |
