[Swift-commit] r8221 - www/ATPESC_2014
wozniak at ci.uchicago.edu
wozniak at ci.uchicago.edu
Thu Aug 14 16:01:14 CDT 2014
Author: wozniak
Date: 2014-08-14 16:01:14 -0500 (Thu, 14 Aug 2014)
New Revision: 8221
Modified:
www/ATPESC_2014/tutorial.html
Log:
Add use case
Modified: www/ATPESC_2014/tutorial.html
===================================================================
--- www/ATPESC_2014/tutorial.html 2014-08-14 20:47:54 UTC (rev 8220)
+++ www/ATPESC_2014/tutorial.html 2014-08-14 21:01:14 UTC (rev 8221)
@@ -3,26 +3,29 @@
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en">
<head>
<meta http-equiv="Content-Type" content="application/xhtml+xml; charset=UTF-8" />
-<meta name="generator" content="AsciiDoc 8.6.9" />
-<title>ATPESC 2014 Swift Tutorial</title>
+<meta name="generator" content="AsciiDoc 8.6.3" />
+<title>Tutorial: Swift/T</title>
<style type="text/css">
-/* Shared CSS for AsciiDoc xhtml11 and html5 backends */
-
-/* Default font. */
-body {
- font-family: Georgia,serif;
-}
-
-/* Title font. */
+/* Sans-serif font. */
h1, h2, h3, h4, h5, h6,
div.title, caption.title,
thead, p.table.header,
-#toctitle,
-#author, #revnumber, #revdate, #revremark,
-#footer {
+div#toctitle,
+span#author, span#revnumber, span#revdate, span#revremark,
+div#footer {
font-family: Arial,Helvetica,sans-serif;
}
+/* Serif font. */
+div.sectionbody {
+ font-family: Georgia,"Times New Roman",Times,serif;
+}
+
+/* Monospace font. */
+tt {
+ font-size: inherit;
+}
+
body {
margin: 1em 5% 1em 5%;
}
@@ -45,6 +48,11 @@
color: #083194;
}
+tt {
+ font-size: inherit;
+ color: navy;
+}
+
h1, h2, h3, h4, h5, h6 {
color: #527bbd;
margin-top: 1.2em;
@@ -64,9 +72,6 @@
h3 + * {
clear: left;
}
-h5 {
- font-size: 1.0em;
-}
div.sectionbody {
margin-left: 0;
@@ -87,47 +92,41 @@
ul > li { color: #aaa; }
ul > li > * { color: black; }
-.monospaced, code, pre {
- font-family: "Courier New", Courier, monospace;
- font-size: inherit;
- color: navy;
+pre {
padding: 0;
margin: 0;
}
-pre {
- white-space: pre-wrap;
-}
-#author {
+span#author {
color: #527bbd;
font-weight: bold;
font-size: 1.1em;
}
-#email {
+span#email {
}
-#revnumber, #revdate, #revremark {
+span#revnumber, span#revdate, span#revremark {
}
-#footer {
+div#footer {
font-size: small;
border-top: 2px solid silver;
padding-top: 0.5em;
margin-top: 4.0em;
}
-#footer-text {
+div#footer-text {
float: left;
padding-bottom: 0.5em;
}
-#footer-badges {
+div#footer-badges {
float: right;
padding-bottom: 0.5em;
}
-#preamble {
+div#preamble {
margin-top: 1.5em;
margin-bottom: 1.5em;
}
-div.imageblock, div.exampleblock, div.verseblock,
+div.tableblock, div.imageblock, div.exampleblock, div.verseblock,
div.quoteblock, div.literalblock, div.listingblock, div.sidebarblock,
div.admonitionblock {
margin-top: 1.0em;
@@ -185,7 +184,7 @@
margin-left: 1.0em;
margin-right: 10%;
border-left: 5px solid #f0f0f0;
- color: #888;
+ color: #777777;
}
div.quoteblock > div.attribution {
@@ -225,7 +224,7 @@
}
div.imageblock div.content { padding-left: 0; }
-span.image img { border-style: none; vertical-align: text-bottom; }
+span.image img { border-style: none; }
a.image:visited { color: white; }
dl {
@@ -268,13 +267,36 @@
margin-bottom: 0.1em;
}
+div.tableblock > table {
+ border: 3px solid #527bbd;
+}
+thead, p.table.header {
+ font-weight: bold;
+ color: #527bbd;
+}
tfoot {
font-weight: bold;
}
td > div.verse {
white-space: pre;
}
+p.table {
+ margin-top: 0;
+}
+/* Because the table frame attribute is overriden by CSS in most browsers. */
+div.tableblock > table[frame="void"] {
+ border-style: none;
+}
+div.tableblock > table[frame="hsides"] {
+ border-left-style: none;
+ border-right-style: none;
+}
+div.tableblock > table[frame="vsides"] {
+ border-top-style: none;
+ border-bottom-style: none;
+}
+
div.hdlist {
margin-top: 0.8em;
margin-bottom: 0.8em;
@@ -340,14 +362,14 @@
}
@media print {
- #footer-badges { display: none; }
+ div#footer-badges { display: none; }
}
-#toc {
+div#toc {
margin-bottom: 2.5em;
}
-#toctitle {
+div#toctitle {
color: #527bbd;
font-size: 1.1em;
font-weight: bold;
@@ -355,7 +377,7 @@
margin-bottom: 0.1em;
}
-div.toclevel0, div.toclevel1, div.toclevel2, div.toclevel3, div.toclevel4 {
+div.toclevel1, div.toclevel2, div.toclevel3, div.toclevel4 {
margin-top: 0;
margin-bottom: 0;
}
@@ -372,170 +394,7 @@
font-size: 0.9em;
}
-span.aqua { color: aqua; }
-span.black { color: black; }
-span.blue { color: blue; }
-span.fuchsia { color: fuchsia; }
-span.gray { color: gray; }
-span.green { color: green; }
-span.lime { color: lime; }
-span.maroon { color: maroon; }
-span.navy { color: navy; }
-span.olive { color: olive; }
-span.purple { color: purple; }
-span.red { color: red; }
-span.silver { color: silver; }
-span.teal { color: teal; }
-span.white { color: white; }
-span.yellow { color: yellow; }
-span.aqua-background { background: aqua; }
-span.black-background { background: black; }
-span.blue-background { background: blue; }
-span.fuchsia-background { background: fuchsia; }
-span.gray-background { background: gray; }
-span.green-background { background: green; }
-span.lime-background { background: lime; }
-span.maroon-background { background: maroon; }
-span.navy-background { background: navy; }
-span.olive-background { background: olive; }
-span.purple-background { background: purple; }
-span.red-background { background: red; }
-span.silver-background { background: silver; }
-span.teal-background { background: teal; }
-span.white-background { background: white; }
-span.yellow-background { background: yellow; }
-
-span.big { font-size: 2em; }
-span.small { font-size: 0.6em; }
-
-span.underline { text-decoration: underline; }
-span.overline { text-decoration: overline; }
-span.line-through { text-decoration: line-through; }
-
-div.unbreakable { page-break-inside: avoid; }
-
-
-/*
- * xhtml11 specific
- *
- * */
-
-div.tableblock {
- margin-top: 1.0em;
- margin-bottom: 1.5em;
-}
-div.tableblock > table {
- border: 3px solid #527bbd;
-}
-thead, p.table.header {
- font-weight: bold;
- color: #527bbd;
-}
-p.table {
- margin-top: 0;
-}
-/* Because the table frame attribute is overriden by CSS in most browsers. */
-div.tableblock > table[frame="void"] {
- border-style: none;
-}
-div.tableblock > table[frame="hsides"] {
- border-left-style: none;
- border-right-style: none;
-}
-div.tableblock > table[frame="vsides"] {
- border-top-style: none;
- border-bottom-style: none;
-}
-
-
-/*
- * html5 specific
- *
- * */
-
-table.tableblock {
- margin-top: 1.0em;
- margin-bottom: 1.5em;
-}
-thead, p.tableblock.header {
- font-weight: bold;
- color: #527bbd;
-}
-p.tableblock {
- margin-top: 0;
-}
-table.tableblock {
- border-width: 3px;
- border-spacing: 0px;
- border-style: solid;
- border-color: #527bbd;
- border-collapse: collapse;
-}
-th.tableblock, td.tableblock {
- border-width: 1px;
- padding: 4px;
- border-style: solid;
- border-color: #527bbd;
-}
-
-table.tableblock.frame-topbot {
- border-left-style: hidden;
- border-right-style: hidden;
-}
-table.tableblock.frame-sides {
- border-top-style: hidden;
- border-bottom-style: hidden;
-}
-table.tableblock.frame-none {
- border-style: hidden;
-}
-
-th.tableblock.halign-left, td.tableblock.halign-left {
- text-align: left;
-}
-th.tableblock.halign-center, td.tableblock.halign-center {
- text-align: center;
-}
-th.tableblock.halign-right, td.tableblock.halign-right {
- text-align: right;
-}
-
-th.tableblock.valign-top, td.tableblock.valign-top {
- vertical-align: top;
-}
-th.tableblock.valign-middle, td.tableblock.valign-middle {
- vertical-align: middle;
-}
-th.tableblock.valign-bottom, td.tableblock.valign-bottom {
- vertical-align: bottom;
-}
-
-
-/*
- * manpage specific
- *
- * */
-
-body.manpage h1 {
- padding-top: 0.5em;
- padding-bottom: 0.5em;
- border-top: 2px solid silver;
- border-bottom: 2px solid silver;
-}
-body.manpage h2 {
- border-style: none;
-}
-body.manpage div.sectionbody {
- margin-left: 3em;
-}
-
- at media print {
- body.manpage div#toc { display: none; }
-}
-
-
-
/* SWIFT/T GUIDE CUSTOMIZATIONS */
a:visited {
@@ -548,6 +407,7 @@
</style>
<script type="text/javascript">
/*<+'])');
+ var re = new RegExp('[hH]([2-'+(toclevels+1)+'])');
// Function that scans the DOM tree for header elements (the DOM2
// nodeIterator API would be a better technique but not supported by all
// browsers).
@@ -609,25 +469,6 @@
}
var toc = document.getElementById("toc");
- if (!toc) {
- return;
- }
-
- // Delete existing TOC entries in case we're reloading the TOC.
- var tocEntriesToRemove = [];
- var i;
- for (i = 0; i < toc.childNodes.length; i++) {
- var entry = toc.childNodes[i];
- if (entry.nodeName.toLowerCase() == 'div'
- && entry.getAttribute("class")
- && entry.getAttribute("class").match(/^toclevel/))
- tocEntriesToRemove.push(entry);
- }
- for (i = 0; i < tocEntriesToRemove.length; i++) {
- toc.removeChild(tocEntriesToRemove[i]);
- }
-
- // Rebuild TOC entries.
var entries = tocEntries(document.getElementById("content"), toclevels);
for (var i = 0; i < entries.length; ++i) {
var entry = entries[i];
@@ -655,44 +496,24 @@
*/
footnotes: function () {
- // Delete existing footnote entries in case we're reloading the footnodes.
- var i;
+ var cont = document.getElementById("content");
var noteholder = document.getElementById("footnotes");
- if (!noteholder) {
- return;
- }
- var entriesToRemove = [];
- for (i = 0; i < noteholder.childNodes.length; i++) {
- var entry = noteholder.childNodes[i];
- if (entry.nodeName.toLowerCase() == 'div' && entry.getAttribute("class") == "footnote")
- entriesToRemove.push(entry);
- }
- for (i = 0; i < entriesToRemove.length; i++) {
- noteholder.removeChild(entriesToRemove[i]);
- }
-
- // Rebuild footnote entries.
- var cont = document.getElementById("content");
var spans = cont.getElementsByTagName("span");
var refs = {};
var n = 0;
for (i=0; i<spans.length; i++) {
if (spans[i].className == "footnote") {
n++;
- var note = spans[i].getAttribute("data-note");
- if (!note) {
- // Use [\s\S] in place of . so multi-line matches work.
- // Because JavaScript has no s (dotall) regex flag.
- note = spans[i].innerHTML.match(/\s*\[([\s\S]*)]\s*/)[1];
- spans[i].innerHTML =
- "[<a id='_footnoteref_" + n + "' href='#_footnote_" + n +
- "' title='View footnote' class='footnote'>" + n + "</a>]";
- spans[i].setAttribute("data-note", note);
- }
+ // Use [\s\S] in place of . so multi-line matches work.
+ // Because JavaScript has no s (dotall) regex flag.
+ note = spans[i].innerHTML.match(/\s*\[([\s\S]*)]\s*/)[1];
noteholder.innerHTML +=
"<div class='footnote' id='_footnote_" + n + "'>" +
"<a href='#_footnoteref_" + n + "' title='Return to text'>" +
n + "</a>. " + note + "</div>";
+ spans[i].innerHTML =
+ "[<a id='_footnoteref_" + n + "' href='#_footnote_" + n +
+ "' title='View footnote' class='footnote'>" + n + "</a>]";
var id =spans[i].getAttribute("id");
if (id != null) refs["#"+id] = n;
}
@@ -712,62 +533,433 @@
}
}
}
-},
-
-install: function(toclevels) {
- var timerId;
-
- function reinstall() {
- asciidoc.footnotes();
- if (toclevels) {
- asciidoc.toc(toclevels);
- }
- }
-
- function reinstallAndRemoveTimer() {
- clearInterval(timerId);
- reinstall();
- }
-
- timerId = setInterval(reinstall, 500);
- if (document.addEventListener)
- document.addEventListener("DOMContentLoaded", reinstallAndRemoveTimer, false);
- else
- window.onload = reinstallAndRemoveTimer;
}
}
-asciidoc.install();
/*]]>*/
</script>
</head>
<body class="article" style="max-width:750px">
<div id="header">
-<h1>ATPESC 2014 Swift Tutorial</h1>
+<h1>Tutorial: Swift/T</h1>
+<div id="toc">
+ <div id="toctitle">Table of Contents</div>
+ <noscript><p><b>JavaScript must be enabled in your browser to display the table of contents.</b></p></noscript>
</div>
+</div>
<div id="content">
<div id="preamble">
<div class="sectionbody">
-<div class="paragraph"><p>This is the Swift tutorial for <a href="http://extremecomputingtraining.anl.gov">ATPESC 2014</a>.</p></div>
+<div class="paragraph"><p>This tutorial will walk through 4 examples of the use of Swift/T, in
+addition to a molecular dynamics example.</p></div>
+</div>
+</div>
+<div class="sect1">
+<h2 id="_swift_t_materials">Swift/T materials</h2>
+<div class="sectionbody">
+<div class="paragraph"><p>For more information on Swift/T, see:</p></div>
<div class="ulist"><ul>
<li>
<p>
-<a href="./swift-k/README">Swift/K tutorial</a>
+<a href="http://swift-lang.org/Swift-T">Swift/T home</a>
</p>
</li>
<li>
<p>
-<a href="./swift-t/tutorial.html">Swift/T tutorial</a>
+<a href="http://swift-lang.org/Swift-T/guide.html">Swift/T user guide</a>
</p>
</li>
</ul></div>
</div>
</div>
+<div class="sect1">
+<h2 id="_files">Files</h2>
+<div class="sectionbody">
+<div class="paragraph"><p>The examples are in the <tt>examples/</tt> directory.</p></div>
+<div class="paragraph"><p>Each example directory is numbered from <tt>01</tt> to <tt>05</tt>.</p></div>
+<div class="paragraph"><p>Configuration files are in <tt>cfg/</tt>. These refer to public
+installations of Swift/T on Cetus/Mira.</p></div>
+<div class="paragraph"><p>Each example contains a <tt>*.swift</tt> Swift script that is the Swift
+program. The Swift/T compiler <tt>stc</tt> translates the Swift script into
+an MPI program encoded in a <tt>*.tcl</tt> file. This is launched by using
+the <tt>turbine</tt> program or a scheduler-specific script, such as the
+provided <tt>turbine-cobalt-run.zsh</tt> script.</p></div>
+<div class="paragraph"><p>There is also a <tt>run-local</tt> script to run the Swift script
+on the login node and a <tt>run-compute</tt> script to submit the Swift
+program to Cobalt. These scripts add <tt>stc</tt> and <tt>turbine</tt> entries to
+your <tt>PATH</tt> automatically.</p></div>
</div>
+</div>
+<div class="sect1">
+<h2 id="_swift_t_installations">Swift/T installations</h2>
+<div class="sectionbody">
+<div class="paragraph"><p>These instructions assume you will be running a preinstalled version
+of Swift/T on Cetus or Mira.</p></div>
+<div class="paragraph"><p>Two Swift/T installations are provided- one that runs on the login
+node and one that runs on the compute nodes (via Cobalt). The login
+node installation should be used with caution- expensive computation
+will impact other users on the login nodes.</p></div>
+</div>
+</div>
+<div class="sect1">
+<h2 id="_examples">Examples</h2>
+<div class="sectionbody">
+<div class="sect2">
+<h3 id="_01_hello_world">01 Hello world</h3>
+<div class="paragraph"><p>This is a very simple Swift/T program that illustrates Swift/T usage.</p></div>
+<div class="olist arabic"><ol class="arabic">
+<li>
+<p>
+<tt>cd</tt> into the directory <tt>01-hello</tt>
+</p>
+</li>
+<li>
+<p>
+Run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>./run-local.sh</tt></pre>
+</div></div>
+</li>
+<li>
+<p>
+This program is a simple Hello World example.
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>/**
+ Example 1 - HELLO.SWIFT
+*/
+
+import io;
+
+main
+{
+ printf("Hello world!");
+}</tt></pre>
+</div></div>
+</li>
+<li>
+<p>
+Run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>./run-compute.sh</tt></pre>
+</div></div>
+<div class="paragraph"><p>This will report something like:</p></div>
+<div class="listingblock">
+<div class="content">
+<pre><tt>SETTINGS=../../cfg/settings-compute.sh
+TURBINE_OUTPUT=/home/user/turbine-output/2014/08/14/10/06/20
+SCRIPT_NAME=hello.tcl
+NODES=2
+project: ATPESC2014
+JOB_ID=312170</tt></pre>
+</div></div>
+<div class="paragraph"><p>Turbine will copy <tt>hello.tcl</tt> to <tt>TURBINE_OUTPUT</tt> and run from that
+directory. The Cobalt log may be found there. Standard output is
+directed to <tt>TURBINE_OUTPUT/output.txt</tt>.</p></div>
+<div class="paragraph"><p>The contents of <tt>output.txt</tt> are the Blue Gene runtime output plus the
+Swift-generated line:</p></div>
+<div class="listingblock">
+<div class="content">
+<pre><tt>stdout[0]: Hello world!</tt></pre>
+</div></div>
+</li>
+</ol></div>
+</div>
+<div class="sect2">
+<h3 id="_02_loop">02 Loop</h3>
+<div class="paragraph"><p>This example shows one way to obtain concurrency from Swift/T.</p></div>
+<div class="listingblock">
+<div class="content">
+<pre><tt>/**
+ Example 2 - LOOP.SWIFT
+*/
+
+import io;
+
+ at dispatch=WORKER
+(int j) compute(int i) "turbine" "0.0"
+[
+----
+ after [ expr <<i>> * 1000 ]
+ puts "compute: <<i>> [ expr [clock seconds] % 1000]"
+ set <<j>> [ expr <<i>> * 2 ]
+----
+];
+
+main
+{
+ printf("Running loop...");
+ foreach i in [0:9]
+ {
+ j = compute(i);
+ printf("compute(%i) = %i", i, j);
+ }
+}</tt></pre>
+</div></div>
+<div class="paragraph"><p>The program executes <tt>compute()</tt> 10 times concurrently, because the
+Swift <tt>foreach</tt> loop is concurrent. The <tt>compute()</tt> function is
+implemented as a Tcl fragment that does 3 things:</p></div>
+<div class="olist arabic"><ol class="arabic">
+<li>
+<p>
+Delay progress with <tt>after</tt> (sleeps for <tt>i*1000</tt> milliseconds)
+</p>
+</li>
+<li>
+<p>
+Print a timestamp for reference
+</p>
+</li>
+<li>
+<p>
+Set a result value
+</p>
+</li>
+</ol></div>
+<div class="paragraph"><p>The result value is printed by Swift’s <tt>printf()</tt>.</p></div>
+<div class="paragraph"><p>Note the use of the <tt><<<tt>x</tt>>></tt> syntax, which maps Swift variable <tt>x</tt> to a
+Tcl variable.</p></div>
+<div class="paragraph"><p>This script may be run as the Hello World case, but the run scripts
+accept Turbine runtime arguments:</p></div>
+<div class="dlist"><dl>
+<dt class="hdlist1">
+<tt>-l</tt>
+</dt>
+<dd>
+<p>
+Prepend outputs with the MPI rank number (like <tt>mpiexec -l</tt>)
+</p>
+</dd>
+<dt class="hdlist1">
+<tt>-n N</tt>
+</dt>
+<dd>
+<p>
+Run <tt>N</tt> MPI processes (like <tt>mpiexec -n N</tt>). Note that one
+process is dedicated for Swift logic, so choose <tt>N>1</tt>.
+</p>
+</dd>
+</dl></div>
+<div class="paragraph"><p>Add more processes and see how that improves time to job
+completion. Check <tt>output.txt</tt> for the <tt>ADLB Total Elapsed Time:</tt>,
+which is a good measure of total execution time (discounts queue and
+allocation startup/shutdown time).</p></div>
+<div class="admonitionblock">
+<table><tr>
+<td class="icon">
+<div class="title">Note</div>
+</td>
+<td class="content">Ease of calling to Tcl makes for ease of calling to native code,
+as shown in the next example.</td>
+</tr></table>
+</div>
+</div>
+<div class="sect2">
+<h3 id="_03_blas">03 BLAS</h3>
+<div class="paragraph"><p>In this example, we use Swift to call a BLAS routine for dot product.</p></div>
+<div class="listingblock">
+<div class="content">
+<pre><tt>/**
+ Example 3 - DOT.SWIFT
+*/
+
+// Swift/T libraries:
+import blob;
+import io;
+import matrix;
+
+// BLAS library for this example:
+import blas;
+
+main
+{
+ // Swift floats are 64-bit
+ float A[];
+ float B[];
+ int N = 3;
+ foreach i in [0:N-1]
+ {
+ A[i] = itof(i); // integer-to-float
+ B[i] = itof(i+1);
+ }
+
+ // Use => for output ordering
+ printf("A") =>
+ vector_print(A) =>
+ printf("\nB") =>
+ vector_print(B)
+ =>
+ float z = blas_ddot(N, A, B);
+
+ printf("\nA*B=%0.1f", z);
+}</tt></pre>
+</div></div>
+<div class="paragraph"><p>In this example, we created a stub of a Swift BLAS library in
+<tt>blas.swift</tt>. This calls the CBLAS <tt>cblas_ddot()</tt> function through a
+Tcl interface.</p></div>
+<div class="olist arabic"><ol class="arabic">
+<li>
+<p>
+Run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>./build.sh</tt></pre>
+</div></div>
+<div class="paragraph"><p>This creates a Tcl package for BLAS (using a hand-coded method- we
+could also use <a href="http://swig.org">SWIG</a> to automate this).</p></div>
+</li>
+<li>
+<p>
+Run <tt>./run-*.sh</tt>.
+</p>
+</li>
+</ol></div>
+<div class="paragraph"><p>In <tt>dot.swift</tt>, the Swift statement <tt>import blas;</tt> brings <tt>blas.swift</tt>
+into the program. The Swift call to <tt>blas_ddot()</tt> operates on the
+Swift arrays <tt>A</tt> and <tt>B</tt>. In <tt>blas.swift</tt>, the Swift arrays are
+converted to C-compatible arrays of <tt>double</tt>, which are represented in
+Swift/T with the <tt>blob</tt> type. These are passed to the BLAS package.</p></div>
+<div class="paragraph"><p>Thus, calling native code from Swift/T is comparable in complexity to
+calling native code from other scripting languages such as Tcl or
+Python.</p></div>
+</div>
+<div class="sect2">
+<h3 id="_04_python">04 Python</h3>
+<div class="paragraph"><p>Another good way to call native code and bind it together in a script
+is through Python. Swift/T has convenient interfaces for Python, Tcl,
+R, Julia, and Qt Script (and the shell). In this example, we call
+into Python and its numerical package Numpy, which uses BLAS
+internally for numerical operations.</p></div>
+<div class="listingblock">
+<div class="content">
+<pre><tt>/**
+ Example 4 - ADD.SWIFT
+*/
+
+// Swift/T libraries:
+import io;
+import python;
+import string;
+
+// This example: See ./numpy.swift
+import numpy;
+
+f()
+{
+ python(
+----
+print("Python works!")
+repr(0)
+----);
+}
+
+main
+{
+ f();
+ matrix A1 = eye(3);
+ matrix A2 = eye(3);
+ matrix sum = add(A1, A2);
+ printf("2*eye(3)=%s", sum);
+}</tt></pre>
+</div></div>
+<div class="paragraph"><p>In this example, we also created a stub Swift numerical library in
+<tt>numpy.swift</tt>. This is imported into the main Swift script
+<tt>add.swift</tt>.</p></div>
+<div class="paragraph"><p>First, this program calls <tt>f()</tt>, which is a simple use of the
+Swift/T-Python interface. The Python code fragment must return a
+value in its last line as a string- in this case, we return <tt>"0"</tt>.</p></div>
+<div class="paragraph"><p>Second, we use the <tt>numpy.swift</tt> features. We allocate 2 Numpy arrays
+<em>I<sub>3</sub></em> and add them together, producing <em>2I<sub>3</sub></em>.</p></div>
+<div class="paragraph"><p>Swift/T can call into arbitrary Python modules, just set <tt>PYTHONPATH</tt>.
+C/C++ code can be easily wrapped for Python with
+<a href="http://swig.org">SWIG</a>, Fortran code with SciPy’s
+<a href="http://wiki.scipy.org/Cookbook/F2Py">F2PY</a>.</p></div>
+<div class="paragraph"><p><strong><span class="red">★</span></strong> Thus, a variety of scripts and
+libraries can be easily integrated with Swift/T and run concurrently
+on massive computers like the Blue Gene/Q!</p></div>
+</div>
+<div class="sect2">
+<h3 id="_05_molecular_dynamics">05 Molecular dynamics</h3>
+<div class="paragraph"><p>In our most realistic example, we bundle a small molecular dynamics
+simulator <a href="http://people.sc.fsu.edu/~jburkardt/c_src/md/md.html">MD</a> for
+Swift/T and call it concurrently.</p></div>
+<div class="olist arabic"><ol class="arabic">
+<li>
+<p>
+<tt>cd</tt> to <tt>examples/md-src/</tt> and run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>make</tt></pre>
+</div></div>
+<div class="paragraph"><p>This builds the stand-alone <tt>md</tt> program. You can test it with <tt>test-md-*.sh</tt>.</p></div>
+</li>
+<li>
+<p>
+Run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>make swift-pkg</tt></pre>
+</div></div>
+<div class="paragraph"><p>This builds the Tcl package by:</p></div>
+<div class="olist loweralpha"><ol class="loweralpha">
+<li>
+<p>
+Running SWIG to expose the core <tt>simulate()</tt> function to Tcl;
+</p>
+</li>
+<li>
+<p>
+Compiling the SWIG-generated wrapper code and linking it as a Tcl
+package.
+</p>
+<div class="paragraph"><p>You can test the Tcl function with <tt>./test-md.tcl</tt>, which simply loads
+the package and calls <tt>simulate</tt> once.</p></div>
+</li>
+</ol></div>
+</li>
+<li>
+<p>
+<tt>cd</tt> back to <tt>examples/05-md/</tt>.
+</p>
+</li>
+<li>
+<p>
+The <tt>md.swift</tt> library shows the calling sequence for MD.
+</p>
+</li>
+<li>
+<p>
+The <tt>run-md.swift</tt> script shows how to call MD once. The output files are indicated in the script, in addition to the usual <tt>output.txt</tt>. Run this with <tt>./run-*-md.sh</tt>.
+</p>
+</li>
+<li>
+<p>
+A concurrent ensemble of simulations is launched by
+<tt>run-mds.swift</tt>. This script accepts many arguments on the command
+line: see <tt>run-*-mds.sh</tt>. For example, you can run:
+</p>
+<div class="listingblock">
+<div class="content">
+<pre><tt>./run-local-mds.sh 1 3 2 10 10 0.1 1 "0.1 1.0 0.2 0.05 50.0 0.1" 2.5 2.0 42</tt></pre>
+</div></div>
+</li>
+</ol></div>
+</div>
+</div>
+</div>
+</div>
<div id="footnotes"><hr /></div>
<div id="footer">
<div id="footer-text">
-Last updated 2014-08-14 12:44:15 CDT
+Last updated 2014-08-14 22:59:40 RDT
</div>
</div>
</body>
More information about the Swift-commit
mailing list