diff --git a/docs/api/html/classyask_1_1yk__env-members.html b/docs/api/html/classyask_1_1yk__env-members.html
index acc16f7c..84a18bc0 100644
--- a/docs/api/html/classyask_1_1yk__env-members.html
+++ b/docs/api/html/classyask_1_1yk__env-members.html
@@ -74,14 +74,15 @@
 <p>This is the complete list of members for <a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a>, including all inherited members.</p>
 <table class="directory">
   <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#ab14fa168dc78346ac61b54c9a509099b">disable_debug_output</a>()</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
-  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a69eac2bb8a841f6259697a7dcc8cf386">get_debug_output</a>()</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
-  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#add8b12cd6d10f964665a41acedbb9b14">get_num_ranks</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
-  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a93deb7b62612742f5a85fd8e319df38c">get_rank_index</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
-  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca">global_barrier</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
-  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#abe162669b192b8c650ee7dbfbd62bb7d">is_trace_enabled</a>()</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
-  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#aa5cf0e8d885f4d9ebde6539d5246cda1">set_debug_output</a>(yask_output_ptr debug)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
-  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#ac6b39c25ae59e27dcdb405a9fe24e763">set_trace_enabled</a>(bool enable)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
-  <tr bgcolor="#f0f0f0" class="even"><td class="entry"><b>~yk_env</b>() (defined in <a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a>)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">inline</span><span class="mlabel">virtual</span></td></tr>
+  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">finalize</a>()=0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
+  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a69eac2bb8a841f6259697a7dcc8cf386">get_debug_output</a>()</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
+  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#add8b12cd6d10f964665a41acedbb9b14">get_num_ranks</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
+  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a93deb7b62612742f5a85fd8e319df38c">get_rank_index</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
+  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca">global_barrier</a>() const =0</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">pure virtual</span></td></tr>
+  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#abe162669b192b8c650ee7dbfbd62bb7d">is_trace_enabled</a>()</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
+  <tr><td class="entry"><a class="el" href="classyask_1_1yk__env.html#aa5cf0e8d885f4d9ebde6539d5246cda1">set_debug_output</a>(yask_output_ptr debug)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
+  <tr class="even"><td class="entry"><a class="el" href="classyask_1_1yk__env.html#ac6b39c25ae59e27dcdb405a9fe24e763">set_trace_enabled</a>(bool enable)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">static</span></td></tr>
+  <tr bgcolor="#f0f0f0"><td class="entry"><b>~yk_env</b>() (defined in <a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a>)</td><td class="entry"><a class="el" href="classyask_1_1yk__env.html">yask::yk_env</a></td><td class="entry"><span class="mlabel">inline</span><span class="mlabel">virtual</span></td></tr>
 </table></div><!-- contents -->
 <!-- start footer part -->
 <hr class="footer"/><address class="footer"><small>
diff --git a/docs/api/html/classyask_1_1yk__env.html b/docs/api/html/classyask_1_1yk__env.html
index 7785359a..0589773a 100644
--- a/docs/api/html/classyask_1_1yk__env.html
+++ b/docs/api/html/classyask_1_1yk__env.html
@@ -91,6 +91,9 @@
 <tr class="memitem:a2384185b4d346418a62db1310a3dddca"><td class="memItemLeft" align="right" valign="top">virtual void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca">global_barrier</a> () const =0</td></tr>
 <tr class="memdesc:a2384185b4d346418a62db1310a3dddca"><td class="mdescLeft">&#160;</td><td class="mdescRight">Wait until all ranks have reached this element.  <a href="#a2384185b4d346418a62db1310a3dddca">More...</a><br /></td></tr>
 <tr class="separator:a2384185b4d346418a62db1310a3dddca"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:ab2d17018fe9fadd9bfdc091d8bf32feb"><td class="memItemLeft" align="right" valign="top">virtual void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">finalize</a> ()=0</td></tr>
+<tr class="memdesc:ab2d17018fe9fadd9bfdc091d8bf32feb"><td class="mdescLeft">&#160;</td><td class="mdescRight">Finalize the environment.  <a href="#ab2d17018fe9fadd9bfdc091d8bf32feb">More...</a><br /></td></tr>
+<tr class="separator:ab2d17018fe9fadd9bfdc091d8bf32feb"><td class="memSeparator" colspan="2">&#160;</td></tr>
 </table><table class="memberdecls">
 <tr class="heading"><td colspan="2"><h2 class="groupheader"><a name="pub-static-methods"></a>
 Static Public Member Functions</h2></td></tr>
@@ -357,6 +360,34 @@ <h2 class="memtitle"><span class="permalink"><a href="#a2384185b4d346418a62db131
 <p>Wait until all ranks have reached this element. </p>
 <p>If MPI is enabled, calls <code>MPI_Barrier()</code>. Otherwise, has no effect. </p>
 
+</div>
+</div>
+<a id="ab2d17018fe9fadd9bfdc091d8bf32feb"></a>
+<h2 class="memtitle"><span class="permalink"><a href="#ab2d17018fe9fadd9bfdc091d8bf32feb">&#9670;&nbsp;</a></span>finalize()</h2>
+
+<div class="memitem">
+<div class="memproto">
+<table class="mlabels">
+  <tr>
+  <td class="mlabels-left">
+      <table class="memname">
+        <tr>
+          <td class="memname">virtual void yask::yk_env::finalize </td>
+          <td>(</td>
+          <td class="paramname"></td><td>)</td>
+          <td></td>
+        </tr>
+      </table>
+  </td>
+  <td class="mlabels-right">
+<span class="mlabels"><span class="mlabel">pure virtual</span></span>  </td>
+  </tr>
+</table>
+</div><div class="memdoc">
+
+<p>Finalize the environment. </p>
+<p>If MPI is enabled and YASK initialized the MPI communicator, calls <code>MPI_Finalize()</code>. This function is automatically called when a <a class="el" href="classyask_1_1yk__env.html" title="Kernel environment.">yk_env</a> object is destroyed. Cannot call <a class="el" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca" title="Wait until all ranks have reached this element.">global_barrier()</a> or any MPI-dependent API after calling this. </p>
+
 </div>
 </div>
 <hr/>The documentation for this class was generated from the following file:<ul>
diff --git a/docs/api/html/functions_f.html b/docs/api/html/functions_f.html
index ebd6ac62..9083f36f 100644
--- a/docs/api/html/functions_f.html
+++ b/docs/api/html/functions_f.html
@@ -65,6 +65,9 @@
 <div class="textblock">Here is a list of all documented class members with links to the class documentation for each member:</div>
 
 <h3><a id="index_f"></a>- f -</h3><ul>
+<li>finalize()
+: <a class="el" href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">yask::yk_env</a>
+</li>
 <li>first_domain_index()
 : <a class="el" href="classyask_1_1yc__solution__base.html#a56f644e7e6b4e96619245f217b3763b5">yask::yc_solution_base</a>
 </li>
@@ -72,7 +75,7 @@ <h3><a id="index_f"></a>- f -</h3><ul>
 : <a class="el" href="classyask_1_1yc__solution.html#a5f3624ded964c465724fd0221d5a2aaa">yask::yc_solution</a>
 </li>
 <li>format_indices()
-: <a class="el" href="classyask_1_1yk__var.html#ae929dea5359d07e7541abab78ca7b139">yask::yk_var</a>
+: <a class="el" href="classyask_1_1yk__var.html#a55d6c585b8384881065c01bb067c96dd">yask::yk_var</a>
 </li>
 <li>format_simple()
 : <a class="el" href="classyask_1_1yc__expr__node.html#a1af7948d0c2a977ed20c0b2d1d561052">yask::yc_expr_node</a>
diff --git a/docs/api/html/functions_func_f.html b/docs/api/html/functions_func_f.html
index b0d7aef0..a532d3ad 100644
--- a/docs/api/html/functions_func_f.html
+++ b/docs/api/html/functions_func_f.html
@@ -65,6 +65,9 @@
 &#160;
 
 <h3><a id="index_f"></a>- f -</h3><ul>
+<li>finalize()
+: <a class="el" href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">yask::yk_env</a>
+</li>
 <li>first_domain_index()
 : <a class="el" href="classyask_1_1yc__solution__base.html#a56f644e7e6b4e96619245f217b3763b5">yask::yc_solution_base</a>
 </li>
@@ -72,7 +75,7 @@ <h3><a id="index_f"></a>- f -</h3><ul>
 : <a class="el" href="classyask_1_1yc__solution.html#a5f3624ded964c465724fd0221d5a2aaa">yask::yc_solution</a>
 </li>
 <li>format_indices()
-: <a class="el" href="classyask_1_1yk__var.html#ae929dea5359d07e7541abab78ca7b139">yask::yk_var</a>
+: <a class="el" href="classyask_1_1yk__var.html#a55d6c585b8384881065c01bb067c96dd">yask::yk_var</a>
 </li>
 <li>format_simple()
 : <a class="el" href="classyask_1_1yc__expr__node.html#a1af7948d0c2a977ed20c0b2d1d561052">yask::yc_expr_node</a>
diff --git a/docs/api/html/search/all_6.js b/docs/api/html/search/all_6.js
index 8a6cb6d5..f9deef96 100644
--- a/docs/api/html/search/all_6.js
+++ b/docs/api/html/search/all_6.js
@@ -1,5 +1,6 @@
 var searchData=
 [
+  ['finalize',['finalize',['../classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb',1,'yask::yk_env']]],
   ['first_5fdomain_5findex',['first_domain_index',['../classyask_1_1yc__solution__base.html#a56f644e7e6b4e96619245f217b3763b5',1,'yask::yc_solution_base']]],
   ['format',['format',['../classyask_1_1yc__solution.html#a5f3624ded964c465724fd0221d5a2aaa',1,'yask::yc_solution']]],
   ['format_5findices',['format_indices',['../classyask_1_1yk__var.html#a55d6c585b8384881065c01bb067c96dd',1,'yask::yk_var::format_indices(const idx_t_vec &amp;indices) const =0'],['../classyask_1_1yk__var.html#ae929dea5359d07e7541abab78ca7b139',1,'yask::yk_var::format_indices(const idx_t_init_list &amp;indices) const =0']]],
diff --git a/docs/api/html/search/functions_6.js b/docs/api/html/search/functions_6.js
index 8a6cb6d5..f9deef96 100644
--- a/docs/api/html/search/functions_6.js
+++ b/docs/api/html/search/functions_6.js
@@ -1,5 +1,6 @@
 var searchData=
 [
+  ['finalize',['finalize',['../classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb',1,'yask::yk_env']]],
   ['first_5fdomain_5findex',['first_domain_index',['../classyask_1_1yc__solution__base.html#a56f644e7e6b4e96619245f217b3763b5',1,'yask::yc_solution_base']]],
   ['format',['format',['../classyask_1_1yc__solution.html#a5f3624ded964c465724fd0221d5a2aaa',1,'yask::yc_solution']]],
   ['format_5findices',['format_indices',['../classyask_1_1yk__var.html#a55d6c585b8384881065c01bb067c96dd',1,'yask::yk_var::format_indices(const idx_t_vec &amp;indices) const =0'],['../classyask_1_1yk__var.html#ae929dea5359d07e7541abab78ca7b139',1,'yask::yk_var::format_indices(const idx_t_init_list &amp;indices) const =0']]],
diff --git a/docs/api/html/yask__kernel__api_8hpp_source.html b/docs/api/html/yask__kernel__api_8hpp_source.html
index 5a6e4d06..1b700cc4 100644
--- a/docs/api/html/yask__kernel__api_8hpp_source.html
+++ b/docs/api/html/yask__kernel__api_8hpp_source.html
@@ -70,17 +70,18 @@
 <div class="title">yask_kernel_api.hpp</div>  </div>
 </div><!--header-->
 <div class="contents">
-<a href="yask__kernel__api_8hpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="comment">/*****************************************************************************</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment">YASK: Yet Another Stencil Kit</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="comment">Copyright (c) 2014-2022, Intel Corporation</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="comment">Permission is hereby granted, free of charge, to any person obtaining a copy</span></div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="comment">of this software and associated documentation files (the &quot;Software&quot;), to</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="comment">deal in the Software without restriction, including without limitation the</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="comment">rights to use, copy, modify, merge, publish, distribute, sublicense, and/or</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;<span class="comment">sell copies of the Software, and to permit persons to whom the Software is</span></div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="comment">furnished to do so, subject to the following conditions:</span></div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="comment">* The above copyright notice and this permission notice shall be included in</span></div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="comment">  all copies or substantial portions of the Software.</span></div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="comment">THE SOFTWARE IS PROVIDED &quot;AS IS&quot;, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR</span></div><div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="comment">IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,</span></div><div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;<span class="comment">FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE</span></div><div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="comment">AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER</span></div><div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;<span class="comment">LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING</span></div><div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;<span class="comment">FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS</span></div><div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment">IN THE SOFTWARE.</span></div><div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;<span class="comment">*****************************************************************************/</span></div><div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;</div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;</div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;<span class="comment">// This file uses Doxygen 1.8 markup for API documentation-generation.</span></div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;<span class="comment">// See http://www.stack.nl/~dimitri/doxygen.</span></div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;<span class="comment"></span><span class="preprocessor">#pragma once</span></div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;</div><div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yask__common__api_8hpp.html">yask_common_api.hpp</a>&quot;</span></div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;<span class="preprocessor">#include &lt;vector&gt;</span></div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;<span class="preprocessor">#include &lt;cinttypes&gt;</span></div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;</div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;<span class="preprocessor">#ifndef MPI_VERSION</span></div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;<span class="keyword">typedef</span> <span class="keywordtype">int</span> MPI_Comm;</div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;<span class="preprocessor">#endif</span></div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;</div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;<span class="keyword">namespace </span>yask {</div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;</div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;    <span class="comment">// Forward declarations of classes and pointers.</span></div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;</div><div class="line"><a name="l00052"></a><span class="lineno"><a class="line" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">   52</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__env.html">yk_env</a>;</div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_env&gt; <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a>;</div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;</div><div class="line"><a name="l00056"></a><span class="lineno"><a class="line" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">   56</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__solution.html">yk_solution</a>;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_solution&gt; <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a>;</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;</div><div class="line"><a name="l00060"></a><span class="lineno"><a class="line" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">   60</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__var.html">yk_var</a>;</div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_var&gt; <a class="code" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">yk_var_ptr</a>;</div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;</div><div class="line"><a name="l00064"></a><span class="lineno"><a class="line" href="group__yk.html#ga12d1d46aeb01bd7509a8dc3251657f75">   64</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__stats.html">yk_stats</a>;</div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_stats&gt; <a class="code" href="group__yk.html#ga12d1d46aeb01bd7509a8dc3251657f75">yk_stats_ptr</a>;</div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;} <span class="comment">// namespace yask.</span></div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yk__solution__api_8hpp.html">aux/yk_solution_api.hpp</a>&quot;</span></div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yk__var__api_8hpp.html">aux/yk_var_api.hpp</a>&quot;</span></div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;</div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;<span class="keyword">namespace </span>yask {</div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;</div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;    <span class="keyword">class </span>yk_factory {</div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;    <span class="keyword">public</span>:</div><div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;        yk_factory();</div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;        <span class="keyword">virtual</span> ~yk_factory() {}</div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;</div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;</div><div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;        <span class="keyword">virtual</span> std::string</div><div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a049d5803312219e78797cced8aca0b9a">get_version_string</a>();</div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;</div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;</div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a></div><div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#afccc05b0095d99f89b7164b603b53d00">new_env</a>() <span class="keyword">const</span>;</div><div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;</div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;</div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a></div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#afccc05b0095d99f89b7164b603b53d00">new_env</a>(MPI_Comm comm) <span class="keyword">const</span>;</div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;</div><div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;</div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a></div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a755b1bfc0dd9bfddfe80d924a188b350">new_solution</a>(<a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a> env ) <span class="keyword">const</span>;</div><div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;</div><div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;</div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a></div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a755b1bfc0dd9bfddfe80d924a188b350">new_solution</a>(<a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a> env ,</div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;                     <span class="keyword">const</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a> source ) <span class="keyword">const</span>;</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;    }; <span class="comment">// yk_factory.</span></div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;</div><div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;</div><div class="line"><a name="l00167"></a><span class="lineno"><a class="line" href="classyask_1_1yk__env.html">  167</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__env.html">yk_env</a> {</div><div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;    <span class="keyword">public</span>:</div><div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;        <span class="keyword">virtual</span> ~<a class="code" href="classyask_1_1yk__env.html">yk_env</a>() {}</div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;</div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;</div><div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#aa5cf0e8d885f4d9ebde6539d5246cda1">set_debug_output</a>(<a class="code" href="group__yask.html#ga605185252a3f3f917593c83b7dde4b66">yask_output_ptr</a> debug );</div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;</div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;</div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#ab14fa168dc78346ac61b54c9a509099b">disable_debug_output</a>();</div><div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;</div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;</div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;        <span class="keyword">static</span> <a class="code" href="group__yask.html#ga605185252a3f3f917593c83b7dde4b66">yask_output_ptr</a></div><div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#a69eac2bb8a841f6259697a7dcc8cf386">get_debug_output</a>();</div><div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;</div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;</div><div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#ac6b39c25ae59e27dcdb405a9fe24e763">set_trace_enabled</a>(<span class="keywordtype">bool</span> enable);</div><div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;</div><div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;</div><div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">bool</span></div><div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#abe162669b192b8c650ee7dbfbd62bb7d">is_trace_enabled</a>();</div><div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;</div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">int</span> <a class="code" href="classyask_1_1yk__env.html#add8b12cd6d10f964665a41acedbb9b14">get_num_ranks</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;</div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;</div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">int</span> <a class="code" href="classyask_1_1yk__env.html#a93deb7b62612742f5a85fd8e319df38c">get_rank_index</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;</div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;</div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca">global_barrier</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;</div><div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;    }; <span class="comment">// yk_env.</span></div><div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;</div><div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;    <a class="code" href="yask__common__api_8hpp.html#af7d3d837169568cf38a2efc3e7b04123">YASK_DEPRECATED</a></div><div class="line"><a name="l00245"></a><span class="lineno"><a class="line" href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">  245</a></span>&#160;    <span class="keyword">typedef</span> <a class="code" href="classyask_1_1yk__var.html">yk_var</a> <a class="code" href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yk_grid</a>;</div><div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;    <a class="code" href="yask__common__api_8hpp.html#af7d3d837169568cf38a2efc3e7b04123">YASK_DEPRECATED</a></div><div class="line"><a name="l00248"></a><span class="lineno"><a class="line" href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">  248</a></span>&#160;    <span class="keyword">typedef</span> <a class="code" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">yk_var_ptr</a> <a class="code" href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">yk_grid_ptr</a>;</div><div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;</div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;} <span class="comment">// namespace yask.</span></div><div class="ttc" id="aclassyask_1_1yk__env_html_ab14fa168dc78346ac61b54c9a509099b"><div class="ttname"><a href="classyask_1_1yk__env.html#ab14fa168dc78346ac61b54c9a509099b">yask::yk_env::disable_debug_output</a></div><div class="ttdeci">static void disable_debug_output()</div><div class="ttdoc">Disable the debug output.</div></div>
+<a href="yask__kernel__api_8hpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="comment">/*****************************************************************************</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment">YASK: Yet Another Stencil Kit</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="comment">Copyright (c) 2014-2022, Intel Corporation</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="comment">Permission is hereby granted, free of charge, to any person obtaining a copy</span></div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="comment">of this software and associated documentation files (the &quot;Software&quot;), to</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="comment">deal in the Software without restriction, including without limitation the</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="comment">rights to use, copy, modify, merge, publish, distribute, sublicense, and/or</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;<span class="comment">sell copies of the Software, and to permit persons to whom the Software is</span></div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="comment">furnished to do so, subject to the following conditions:</span></div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="comment">* The above copyright notice and this permission notice shall be included in</span></div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="comment">  all copies or substantial portions of the Software.</span></div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="comment">THE SOFTWARE IS PROVIDED &quot;AS IS&quot;, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR</span></div><div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="comment">IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,</span></div><div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;<span class="comment">FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE</span></div><div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="comment">AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER</span></div><div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;<span class="comment">LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING</span></div><div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;<span class="comment">FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS</span></div><div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment">IN THE SOFTWARE.</span></div><div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;<span class="comment"></span></div><div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;<span class="comment">*****************************************************************************/</span></div><div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;</div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;</div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;<span class="comment">// This file uses Doxygen 1.8 markup for API documentation-generation.</span></div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;<span class="comment">// See http://www.stack.nl/~dimitri/doxygen.</span></div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;<span class="comment"></span><span class="preprocessor">#pragma once</span></div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;</div><div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yask__common__api_8hpp.html">yask_common_api.hpp</a>&quot;</span></div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;<span class="preprocessor">#include &lt;vector&gt;</span></div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;<span class="preprocessor">#include &lt;cinttypes&gt;</span></div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;</div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;<span class="preprocessor">#ifndef MPI_VERSION</span></div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;<span class="keyword">typedef</span> <span class="keywordtype">int</span> MPI_Comm;</div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;<span class="preprocessor">#endif</span></div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;</div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;<span class="keyword">namespace </span>yask {</div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;</div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;    <span class="comment">// Forward declarations of classes and pointers.</span></div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;</div><div class="line"><a name="l00052"></a><span class="lineno"><a class="line" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">   52</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__env.html">yk_env</a>;</div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_env&gt; <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a>;</div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;</div><div class="line"><a name="l00056"></a><span class="lineno"><a class="line" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">   56</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__solution.html">yk_solution</a>;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_solution&gt; <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a>;</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;</div><div class="line"><a name="l00060"></a><span class="lineno"><a class="line" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">   60</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__var.html">yk_var</a>;</div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_var&gt; <a class="code" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">yk_var_ptr</a>;</div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;</div><div class="line"><a name="l00064"></a><span class="lineno"><a class="line" href="group__yk.html#ga12d1d46aeb01bd7509a8dc3251657f75">   64</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__stats.html">yk_stats</a>;</div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;    <span class="keyword">typedef</span> std::shared_ptr&lt;yk_stats&gt; <a class="code" href="group__yk.html#ga12d1d46aeb01bd7509a8dc3251657f75">yk_stats_ptr</a>;</div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;} <span class="comment">// namespace yask.</span></div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yk__solution__api_8hpp.html">aux/yk_solution_api.hpp</a>&quot;</span></div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="yk__var__api_8hpp.html">aux/yk_var_api.hpp</a>&quot;</span></div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;</div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;<span class="keyword">namespace </span>yask {</div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;</div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;    <span class="keyword">class </span>yk_factory {</div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;    <span class="keyword">public</span>:</div><div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;        yk_factory();</div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;        <span class="keyword">virtual</span> ~yk_factory() {}</div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;</div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;</div><div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;        <span class="keyword">virtual</span> std::string</div><div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a049d5803312219e78797cced8aca0b9a">get_version_string</a>();</div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;</div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;</div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a></div><div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#afccc05b0095d99f89b7164b603b53d00">new_env</a>() <span class="keyword">const</span>;</div><div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;</div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;</div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a></div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#afccc05b0095d99f89b7164b603b53d00">new_env</a>(MPI_Comm comm) <span class="keyword">const</span>;</div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;</div><div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;</div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a></div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a755b1bfc0dd9bfddfe80d924a188b350">new_solution</a>(<a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a> env ) <span class="keyword">const</span>;</div><div class="line"><a name="l00147"></a><span class="lineno">  147</span>&#160;</div><div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;</div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;        <span class="keyword">virtual</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a></div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;        <a class="code" href="classyask_1_1yk__factory.html#a755b1bfc0dd9bfddfe80d924a188b350">new_solution</a>(<a class="code" href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yk_env_ptr</a> env ,</div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;                     <span class="keyword">const</span> <a class="code" href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yk_solution_ptr</a> source ) <span class="keyword">const</span>;</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;    }; <span class="comment">// yk_factory.</span></div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;</div><div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;</div><div class="line"><a name="l00167"></a><span class="lineno"><a class="line" href="classyask_1_1yk__env.html">  167</a></span>&#160;    <span class="keyword">class </span><a class="code" href="classyask_1_1yk__env.html">yk_env</a> {</div><div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;    <span class="keyword">public</span>:</div><div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;        <span class="keyword">virtual</span> ~<a class="code" href="classyask_1_1yk__env.html">yk_env</a>() {}</div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;</div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;</div><div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#aa5cf0e8d885f4d9ebde6539d5246cda1">set_debug_output</a>(<a class="code" href="group__yask.html#ga605185252a3f3f917593c83b7dde4b66">yask_output_ptr</a> debug );</div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;</div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;</div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#ab14fa168dc78346ac61b54c9a509099b">disable_debug_output</a>();</div><div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;</div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;</div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;        <span class="keyword">static</span> <a class="code" href="group__yask.html#ga605185252a3f3f917593c83b7dde4b66">yask_output_ptr</a></div><div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#a69eac2bb8a841f6259697a7dcc8cf386">get_debug_output</a>();</div><div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;</div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;</div><div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#ac6b39c25ae59e27dcdb405a9fe24e763">set_trace_enabled</a>(<span class="keywordtype">bool</span> enable);</div><div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;</div><div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;</div><div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;        <span class="keyword">static</span> <span class="keywordtype">bool</span></div><div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#abe162669b192b8c650ee7dbfbd62bb7d">is_trace_enabled</a>();</div><div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;</div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">int</span> <a class="code" href="classyask_1_1yk__env.html#add8b12cd6d10f964665a41acedbb9b14">get_num_ranks</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;</div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;</div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">int</span> <a class="code" href="classyask_1_1yk__env.html#a93deb7b62612742f5a85fd8e319df38c">get_rank_index</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;</div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;</div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#a2384185b4d346418a62db1310a3dddca">global_barrier</a>() <span class="keyword">const</span> =0;</div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;</div><div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;</div><div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;        <span class="keyword">virtual</span> <span class="keywordtype">void</span></div><div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;        <a class="code" href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">finalize</a>() =0;</div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;</div><div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;    }; <span class="comment">// yk_env.</span></div><div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;</div><div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;    <a class="code" href="yask__common__api_8hpp.html#af7d3d837169568cf38a2efc3e7b04123">YASK_DEPRECATED</a></div><div class="line"><a name="l00257"></a><span class="lineno"><a class="line" href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">  257</a></span>&#160;    <span class="keyword">typedef</span> <a class="code" href="classyask_1_1yk__var.html">yk_var</a> <a class="code" href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yk_grid</a>;</div><div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;    <a class="code" href="yask__common__api_8hpp.html#af7d3d837169568cf38a2efc3e7b04123">YASK_DEPRECATED</a></div><div class="line"><a name="l00260"></a><span class="lineno"><a class="line" href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">  260</a></span>&#160;    <span class="keyword">typedef</span> <a class="code" href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">yk_var_ptr</a> <a class="code" href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">yk_grid_ptr</a>;</div><div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;</div><div class="line"><a name="l00264"></a><span class="lineno">  264</span>&#160;} <span class="comment">// namespace yask.</span></div><div class="ttc" id="aclassyask_1_1yk__env_html_ab14fa168dc78346ac61b54c9a509099b"><div class="ttname"><a href="classyask_1_1yk__env.html#ab14fa168dc78346ac61b54c9a509099b">yask::yk_env::disable_debug_output</a></div><div class="ttdeci">static void disable_debug_output()</div><div class="ttdoc">Disable the debug output.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__stats_html"><div class="ttname"><a href="classyask_1_1yk__stats.html">yask::yk_stats</a></div><div class="ttdoc">Statistics from calls to run_solution().</div><div class="ttdef"><b>Definition:</b> yk_solution_api.hpp:1219</div></div>
 <div class="ttc" id="agroup__yk_html_ga12d1d46aeb01bd7509a8dc3251657f75"><div class="ttname"><a href="group__yk.html#ga12d1d46aeb01bd7509a8dc3251657f75">yask::yk_stats_ptr</a></div><div class="ttdeci">std::shared_ptr&lt; yk_stats &gt; yk_stats_ptr</div><div class="ttdoc">Shared pointer to yk_stats.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:64</div></div>
 <div class="ttc" id="aclassyask_1_1yk__env_html"><div class="ttname"><a href="classyask_1_1yk__env.html">yask::yk_env</a></div><div class="ttdoc">Kernel environment.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:167</div></div>
 <div class="ttc" id="aclassyask_1_1yk__factory_html_afccc05b0095d99f89b7164b603b53d00"><div class="ttname"><a href="classyask_1_1yk__factory.html#afccc05b0095d99f89b7164b603b53d00">yask::yk_factory::new_env</a></div><div class="ttdeci">virtual yk_env_ptr new_env() const</div><div class="ttdoc">Create an object to hold environment information.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__env_html_abe162669b192b8c650ee7dbfbd62bb7d"><div class="ttname"><a href="classyask_1_1yk__env.html#abe162669b192b8c650ee7dbfbd62bb7d">yask::yk_env::is_trace_enabled</a></div><div class="ttdeci">static bool is_trace_enabled()</div><div class="ttdoc">Get whether tracing is enabled.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__env_html_add8b12cd6d10f964665a41acedbb9b14"><div class="ttname"><a href="classyask_1_1yk__env.html#add8b12cd6d10f964665a41acedbb9b14">yask::yk_env::get_num_ranks</a></div><div class="ttdeci">virtual int get_num_ranks() const =0</div><div class="ttdoc">Get number of MPI ranks.</div></div>
+<div class="ttc" id="aclassyask_1_1yk__env_html_ab2d17018fe9fadd9bfdc091d8bf32feb"><div class="ttname"><a href="classyask_1_1yk__env.html#ab2d17018fe9fadd9bfdc091d8bf32feb">yask::yk_env::finalize</a></div><div class="ttdeci">virtual void finalize()=0</div><div class="ttdoc">Finalize the environment.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__env_html_a69eac2bb8a841f6259697a7dcc8cf386"><div class="ttname"><a href="classyask_1_1yk__env.html#a69eac2bb8a841f6259697a7dcc8cf386">yask::yk_env::get_debug_output</a></div><div class="ttdeci">static yask_output_ptr get_debug_output()</div><div class="ttdoc">Get object to receive debug output.</div></div>
 <div class="ttc" id="agroup__yk_html_ga8dc62f5599d5c5eb9f7583d7d6a63df1"><div class="ttname"><a href="group__yk.html#ga8dc62f5599d5c5eb9f7583d7d6a63df1">yask::yk_env_ptr</a></div><div class="ttdeci">std::shared_ptr&lt; yk_env &gt; yk_env_ptr</div><div class="ttdoc">Shared pointer to yk_env.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:52</div></div>
 <div class="ttc" id="ayk__solution__api_8hpp_html"><div class="ttname"><a href="yk__solution__api_8hpp.html">yk_solution_api.hpp</a></div></div>
-<div class="ttc" id="agroup__yk_html_gaf6e19ac605b32b47d4edc5a8985b3c5d"><div class="ttname"><a href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">yask::yk_grid_ptr</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var_ptr yk_grid_ptr</div><div class="ttdoc">[Deprecated] Use yk_var_ptr.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:248</div></div>
+<div class="ttc" id="agroup__yk_html_gaf6e19ac605b32b47d4edc5a8985b3c5d"><div class="ttname"><a href="group__yk.html#gaf6e19ac605b32b47d4edc5a8985b3c5d">yask::yk_grid_ptr</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var_ptr yk_grid_ptr</div><div class="ttdoc">[Deprecated] Use yk_var_ptr.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:260</div></div>
 <div class="ttc" id="aclassyask_1_1yk__factory_html_a049d5803312219e78797cced8aca0b9a"><div class="ttname"><a href="classyask_1_1yk__factory.html#a049d5803312219e78797cced8aca0b9a">yask::yk_factory::get_version_string</a></div><div class="ttdeci">virtual std::string get_version_string()</div><div class="ttdoc">Version information.</div></div>
 <div class="ttc" id="agroup__yask_html_ga605185252a3f3f917593c83b7dde4b66"><div class="ttname"><a href="group__yask.html#ga605185252a3f3f917593c83b7dde4b66">yask::yask_output_ptr</a></div><div class="ttdeci">std::shared_ptr&lt; yask_output &gt; yask_output_ptr</div><div class="ttdoc">Shared pointer to yask_output.</div><div class="ttdef"><b>Definition:</b> yask_common_api.hpp:94</div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html"><div class="ttname"><a href="classyask_1_1yk__var.html">yask::yk_var</a></div><div class="ttdoc">A run-time YASK data container.</div><div class="ttdef"><b>Definition:</b> yk_var_api.hpp:185</div></div>
@@ -93,7 +94,7 @@
 <div class="ttc" id="agroup__yk_html_ga2debaa7135bb46dfc295ca623bee2876"><div class="ttname"><a href="group__yk.html#ga2debaa7135bb46dfc295ca623bee2876">yask::yk_solution_ptr</a></div><div class="ttdeci">std::shared_ptr&lt; yk_solution &gt; yk_solution_ptr</div><div class="ttdoc">Shared pointer to yk_solution.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:56</div></div>
 <div class="ttc" id="aclassyask_1_1yk__env_html_aa5cf0e8d885f4d9ebde6539d5246cda1"><div class="ttname"><a href="classyask_1_1yk__env.html#aa5cf0e8d885f4d9ebde6539d5246cda1">yask::yk_env::set_debug_output</a></div><div class="ttdeci">static void set_debug_output(yask_output_ptr debug)</div><div class="ttdoc">Set object to receive debug output.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__solution_html"><div class="ttname"><a href="classyask_1_1yk__solution.html">yask::yk_solution</a></div><div class="ttdoc">Stencil solution as defined by the generated code from the YASK stencil compiler.</div><div class="ttdef"><b>Definition:</b> yk_solution_api.hpp:82</div></div>
-<div class="ttc" id="agroup__yk_html_gab1c5abbc86c9fdde32def4217482cc63"><div class="ttname"><a href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yask::yk_grid</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var yk_grid</div><div class="ttdoc">[Deprecated] Use yk_var.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:245</div></div>
+<div class="ttc" id="agroup__yk_html_gab1c5abbc86c9fdde32def4217482cc63"><div class="ttname"><a href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yask::yk_grid</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var yk_grid</div><div class="ttdoc">[Deprecated] Use yk_var.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:257</div></div>
 <div class="ttc" id="ayk__var__api_8hpp_html"><div class="ttname"><a href="yk__var__api_8hpp.html">yk_var_api.hpp</a></div></div>
 <div class="ttc" id="agroup__yk_html_ga95f75e2b515e5455b570ae705115696a"><div class="ttname"><a href="group__yk.html#ga95f75e2b515e5455b570ae705115696a">yask::yk_var_ptr</a></div><div class="ttdeci">std::shared_ptr&lt; yk_var &gt; yk_var_ptr</div><div class="ttdoc">Shared pointer to yk_var.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:60</div></div>
 </div><!-- fragment --></div><!-- contents -->
diff --git a/docs/api/html/yk__var__api_8hpp_source.html b/docs/api/html/yk__var__api_8hpp_source.html
index 0d091f95..1080a972 100644
--- a/docs/api/html/yk__var__api_8hpp_source.html
+++ b/docs/api/html/yk__var__api_8hpp_source.html
@@ -132,7 +132,7 @@
 <div class="ttc" id="ayask__kernel__api_8hpp_html"><div class="ttname"><a href="yask__kernel__api_8hpp.html">yask_kernel_api.hpp</a></div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html_a420e0dde2114ba663069b16d579072eb"><div class="ttname"><a href="classyask_1_1yk__var.html#a420e0dde2114ba663069b16d579072eb">yask::yk_var::is_fixed_size</a></div><div class="ttdeci">virtual bool is_fixed_size() const =0</div><div class="ttdoc">Determine whether this var is not automatically resized based on the solution.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html_a4a4e42c4cfc5b3b6a5d998611d0d8602"><div class="ttname"><a href="classyask_1_1yk__var.html#a4a4e42c4cfc5b3b6a5d998611d0d8602">yask::yk_var::get_raw_storage_buffer</a></div><div class="ttdeci">virtual void * get_raw_storage_buffer()=0</div><div class="ttdoc">[Advanced] Get pointer to raw data storage buffer.</div></div>
-<div class="ttc" id="agroup__yk_html_gab1c5abbc86c9fdde32def4217482cc63"><div class="ttname"><a href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yask::yk_grid</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var yk_grid</div><div class="ttdoc">[Deprecated] Use yk_var.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:245</div></div>
+<div class="ttc" id="agroup__yk_html_gab1c5abbc86c9fdde32def4217482cc63"><div class="ttname"><a href="group__yk.html#gab1c5abbc86c9fdde32def4217482cc63">yask::yk_grid</a></div><div class="ttdeci">YASK_DEPRECATED typedef yk_var yk_grid</div><div class="ttdoc">[Deprecated] Use yk_var.</div><div class="ttdef"><b>Definition:</b> yask_kernel_api.hpp:257</div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html_ae62308be69ca643714879b0a1362de0e"><div class="ttname"><a href="classyask_1_1yk__var.html#ae62308be69ca643714879b0a1362de0e">yask::yk_var::get_left_extra_pad_size</a></div><div class="ttdeci">virtual idx_t get_left_extra_pad_size(const std::string &amp;dim) const =0</div><div class="ttdoc">Get the actual extra left padding in the specified dimension.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html_a3611ad6d130cb5ee7e030e6ad0c24f5b"><div class="ttname"><a href="classyask_1_1yk__var.html#a3611ad6d130cb5ee7e030e6ad0c24f5b">yask::yk_var::is_storage_layout_identical</a></div><div class="ttdeci">virtual bool is_storage_layout_identical(const yk_var_ptr other) const =0</div><div class="ttdoc">[Advanced] Determines whether storage layout is the same as another var.</div></div>
 <div class="ttc" id="aclassyask_1_1yk__var_html_aaf5b3f06b832bbf77b8f722fda8a2998"><div class="ttname"><a href="classyask_1_1yk__var.html#aaf5b3f06b832bbf77b8f722fda8a2998">yask::yk_var::get_first_local_index</a></div><div class="ttdeci">virtual idx_t get_first_local_index(const std::string &amp;dim) const =0</div><div class="ttdoc">Get the first valid index in this rank in the specified dimension.</div></div>
diff --git a/include/yask_kernel_api.hpp b/include/yask_kernel_api.hpp
index 42e78ba4..f0d493ba 100644
--- a/include/yask_kernel_api.hpp
+++ b/include/yask_kernel_api.hpp
@@ -238,6 +238,18 @@ namespace yask {
         virtual void
         global_barrier() const =0;
 
+        /// Finalize the environment.
+        /**
+           If MPI is enabled and YASK initialized the MPI communicator,
+           calls `MPI_Finalize()`.
+           This function is automatically called when a yk_env object
+           is destroyed.
+           Cannot call global_barrier() or any MPI-dependent API after
+           calling this.
+        */
+        virtual void
+        finalize() =0;
+
     }; // yk_env.
 
     /// **[Deprecated]** Use yk_var.
diff --git a/src/common/common_utils.cpp b/src/common/common_utils.cpp
index 65d7085f..7ce586ad 100644
--- a/src/common/common_utils.cpp
+++ b/src/common/common_utils.cpp
@@ -43,7 +43,7 @@ namespace yask {
     // for numbers above 9 (at least up to 99).
 
     // Format: "major.minor.patch[-alpha|-beta]".
-    const string version = "4.00.02";
+    const string version = "4.00.03";
 
     string yask_get_version_string() {
         return version;
diff --git a/src/kernel/Makefile b/src/kernel/Makefile
index 05f801b6..b6f66673 100644
--- a/src/kernel/Makefile
+++ b/src/kernel/Makefile
@@ -965,15 +965,19 @@ $(YK_API_TEST_EXEC): $(YK_TEST_SRC_DIR)/yask_kernel_api_test.cpp $(YK_LIB)
 	@ls -l $@
 
 # Build and run C++ tests.
+# First run on 1 rank, then multiple ranks if ranks>1.
 cxx-yk-api-test:
 	$(TEST_MAKE) $(YK_API_TEST_EXEC) allow_new_var_types=1
 	@echo '*** Running the C++ YASK kernel API test...'
 	$(RUN_PREFIX) $(YK_API_TEST_EXEC)
+	if (( $(ranks) > 1 )); then $(RUN_PREFIX) mpirun -np $(ranks) $(YK_API_TEST_EXEC); fi
 
 # Run Python kernel API test.
+# First run on 1 rank, then multiple ranks if ranks>1.
 py-yk-api-test: $(YK_TEST_SRC_DIR)/yask_kernel_api_test.py $(YK_PY_LIB)
 	@echo '*** Running the Python YASK kernel API test...'
 	$(RUN_PYTHON) $<
+	if (( $(ranks) > 1 )); then mpirun -np $(ranks) $(RUN_PYTHON) $<; fi
 
 # Build C++ kernel tests with exceptions.
 $(YK_API_TEST_EXEC_WITH_EXCEPTION): $(YK_TEST_SRC_DIR)/yask_kernel_api_exception_test.cpp $(YK_LIB)
@@ -1063,10 +1067,10 @@ yk-test-no-yc: kernel-no-yc $(YK_SCRIPT)
 
 # Run the kernel API tests for C++ and Python with and w/o expected exceptions.
 api-tests:
-	$(MAKE) clean; $(MAKE) cxx-yk-api-test real_bytes=8 stencil=iso3dfd
-	$(MAKE) clean; $(MAKE) py-yk-api-test stencil=iso3dfd
-	$(MAKE) clean; $(MAKE) cxx-yk-api-test-with-exception real_bytes=8 stencil=iso3dfd
-	$(MAKE) clean; $(MAKE) py-yk-api-test-with-exception stencil=iso3dfd
+	$(MAKE) clean; $(MAKE) cxx-yk-api-test real_bytes=8
+	$(MAKE) clean; $(MAKE) py-yk-api-test
+	$(MAKE) clean; $(MAKE) cxx-yk-api-test-with-exception real_bytes=8
+	$(MAKE) clean; $(MAKE) py-yk-api-test-with-exception
 
 # Run several stencils using built-in validation.
 # NB: set arch var as applicable.
diff --git a/src/kernel/lib/settings.hpp b/src/kernel/lib/settings.hpp
index 498e3b1f..fdf7c495 100644
--- a/src/kernel/lib/settings.hpp
+++ b/src/kernel/lib/settings.hpp
@@ -70,6 +70,7 @@ namespace yask {
         MPI_Group group = MPI_GROUP_NULL;
         int num_ranks = 1;        // total number of ranks.
         int my_rank = 0;          // MPI-assigned index.
+        bool finalize_needed = false;
 
         // Vars for shared-mem ranks.
         MPI_Comm shm_comm = MPI_COMM_NULL; // shm communicator.
@@ -81,12 +82,24 @@ namespace yask {
         int max_threads=0;      // initial value from OMP.
 
         KernelEnv() { }
-        virtual ~KernelEnv() { }
+        virtual ~KernelEnv() {
+            finalize();
+        }
 
         // Init MPI, OMP, etc.
         // This is normally called very early in the program.
         virtual void init_env(int* argc, char*** argv, MPI_Comm comm);
 
+        virtual void finalize() {
+            TRACE_MSG("finalize_needed = " << finalize_needed);
+            if (comm != MPI_COMM_NULL && finalize_needed) {
+                MPI_Finalize();
+                comm = MPI_COMM_NULL;
+                shm_comm = MPI_COMM_NULL;
+            }
+            finalize_needed = false;
+        }
+
         // Lock.
         static void set_debug_lock() {
             if (!_debug_lock_init_done) {
diff --git a/src/kernel/lib/setup.cpp b/src/kernel/lib/setup.cpp
index 5ed9e1c2..68ce154f 100644
--- a/src/kernel/lib/setup.cpp
+++ b/src/kernel/lib/setup.cpp
@@ -61,6 +61,7 @@ namespace yask {
                     THROW_YASK_EXCEPTION("error: MPI_THREAD_SERIALIZED or MPI_THREAD_MULTIPLE not provided");
                 }
                 is_init = true;
+                finalize_needed = true;
             }
             comm = MPI_COMM_WORLD;
         }
diff --git a/src/kernel/lib/yk_var.hpp b/src/kernel/lib/yk_var.hpp
index 0582e5da..45733b1d 100644
--- a/src/kernel/lib/yk_var.hpp
+++ b/src/kernel/lib/yk_var.hpp
@@ -654,7 +654,11 @@ namespace yask {
         // Find sizes needed for slicing.
         inline Indices get_slice_range(const Indices& first_indices,
                                        const Indices& last_indices) const {
-            return last_indices.add_const(1).sub_elements(first_indices);
+            host_assert(first_indices.get_num_dims() == last_indices.get_num_dims());
+            Indices range = first_indices;
+            for (int i = 0; i < first_indices.get_num_dims(); i++)
+                range[i] = std::max(idx_t(0), last_indices[i] + 1 - first_indices[i]);
+            return range;
         }
 
         // Accessors to GenericVar.
@@ -955,8 +959,8 @@ namespace yask {
         template<typename Visitor>
         idx_t _visit_elements_in_slice(bool strict_indices,
                                        void* buffer_ptr,
-                                       const Indices& first_indices,
-                                       const Indices& last_indices,
+                                       const Indices& in_first_indices,
+                                       const Indices& in_last_indices,
                                        bool on_device) {
             STATE_VARS(this);
             if (get_storage() == 0) {
@@ -969,11 +973,17 @@ namespace yask {
             }
 
             TRACE_MSG(Visitor::fname() << ": " << make_info_string() << " [" <<
-                      make_index_string(first_indices) << " ... " <<
-                      make_index_string(last_indices) << "] with buffer at " <<
+                      make_index_string(in_first_indices) << " ... " <<
+                      make_index_string(in_last_indices) << "] with buffer at " <<
                       buffer_ptr << " on " << (on_device ? "OMP device" : "host"));
-            check_indices(first_indices, Visitor::fname(), strict_indices, true, false);
-            check_indices(last_indices, Visitor::fname(), strict_indices, true, false);
+            Indices first_indices, last_indices;
+            check_indices(in_first_indices, Visitor::fname(), strict_indices,
+                          true, false, &first_indices);
+            check_indices(in_last_indices, Visitor::fname(), strict_indices,
+                          true, false, &last_indices);
+            TRACE_MSG(Visitor::fname() << ": clipped to  [" <<
+                      make_index_string(first_indices) << " ... " <<
+                      make_index_string(last_indices) << "]");
 
             // Find range.
             auto range = get_slice_range(first_indices, last_indices);
@@ -1927,31 +1937,31 @@ namespace yask {
         }
         virtual idx_t set_element(double val,
                                   const Indices& indices,
-                                  bool strict_indices = false);
+                                  bool strict_indices);
         virtual idx_t set_element(double val,
                                   const VarIndices& indices,
-                                  bool strict_indices = false) {
+                                  bool strict_indices) {
             const Indices indices2(indices);
             return set_element(val, indices2, strict_indices);
         }
         virtual idx_t set_element(double val,
                                   const idx_t_init_list& indices,
-                                  bool strict_indices = false) {
+                                  bool strict_indices) {
             const Indices indices2(indices);
             return set_element(val, indices2, strict_indices);
         }
         virtual idx_t add_to_element(double val,
                                      const Indices& indices,
-                                     bool strict_indices = false);
+                                     bool strict_indices);
         virtual idx_t add_to_element(double val,
                                      const VarIndices& indices,
-                                     bool strict_indices = false) {
+                                     bool strict_indices) {
             const Indices indices2(indices);
             return add_to_element(val, indices2, strict_indices);
         }
         virtual idx_t add_to_element(double val,
                                      const idx_t_init_list& indices,
-                                     bool strict_indices = false) {
+                                     bool strict_indices) {
             const Indices indices2(indices);
             return add_to_element(val, indices2, strict_indices);
         }
diff --git a/src/kernel/tests/var_test.cpp b/src/kernel/tests/var_test.cpp
index 300d8c5c..ee599164 100644
--- a/src/kernel/tests/var_test.cpp
+++ b/src/kernel/tests/var_test.cpp
@@ -87,8 +87,8 @@ void run_tests(int argc, char* argv[]) {
 
         double val = 3.14;
         os << "Testing with " << val << endl;
-        g0->set_element(val, {});
-        g1->set_element(val, {});
+        g0->set_element(val, {}, true);
+        g1->set_element(val, {}, true);
         auto v0 = g0->get_element({});
         auto v1 = g1->get_element({});
         assert(v0 == v1);
@@ -133,8 +133,8 @@ void run_tests(int argc, char* argv[]) {
 
         IdxTuple first, last;
         for (auto dname : gdims) {
-            first.add_dim_back(dname, g3->get_first_rank_alloc_index(dname));
-            last.add_dim_back(dname, g3->get_last_rank_alloc_index(dname));
+            first.add_dim_back(dname, g3->get_first_local_index(dname));
+            last.add_dim_back(dname, g3->get_last_local_index(dname));
         }
         Indices firsti(first), lasti(last);
 
@@ -263,6 +263,7 @@ void run_tests(int argc, char* argv[]) {
         offload_map_free(buf, sz);
         os << "Exiting 3-D test\n";
     }
+    kenv->finalize();
 }
 
 int main(int argc, char* argv[]) {
diff --git a/src/kernel/tests/yask_kernel_api_exception_test.cpp b/src/kernel/tests/yask_kernel_api_exception_test.cpp
index dbea63a3..e8f1cddb 100644
--- a/src/kernel/tests/yask_kernel_api_exception_test.cpp
+++ b/src/kernel/tests/yask_kernel_api_exception_test.cpp
@@ -99,6 +99,7 @@ int main() {
     soln->run_solution(1, 10);
     soln->end_solution();
     soln->get_stats();
+    env->finalize();
 
     // TODO: add exception tests for the methods below:
     // StencilContext::calc_region
@@ -128,7 +129,7 @@ int main() {
         exit(1);
     }
     else
-        cout << "End of YASK kernel API test with exception.\n";
+        cout << "End of YASK kernel C++ API test with exception.\n";
 
     return 0;
 }
diff --git a/src/kernel/tests/yask_kernel_api_exception_test.py b/src/kernel/tests/yask_kernel_api_exception_test.py
index 639c3fee..d8a135ce 100755
--- a/src/kernel/tests/yask_kernel_api_exception_test.py
+++ b/src/kernel/tests/yask_kernel_api_exception_test.py
@@ -105,9 +105,10 @@
     soln.end_solution()
     soln.get_stats()
     #print("Debug output captured:\n", debug_output.get_string())
-
+    env.finalize()
+    
     if num_exception != 2:
         print("There is a problem in exception test.")
         exit(1)
     else:
-        print("End of YASK kernel API test with exception.")
+        print("End of YASK Python kernel API test with exception.")
diff --git a/src/kernel/tests/yask_kernel_api_test.cpp b/src/kernel/tests/yask_kernel_api_test.cpp
index 101cf43c..eac9ffb4 100644
--- a/src/kernel/tests/yask_kernel_api_test.cpp
+++ b/src/kernel/tests/yask_kernel_api_test.cpp
@@ -36,23 +36,30 @@ IN THE SOFTWARE.
 using namespace std;
 using namespace yask;
 
-int main() {
-
-    // The factory from which all other kernel objects are made.
-    yk_factory kfac;
-
-    // Initalize MPI, etc.
-    auto env = kfac.new_env();
+int main(int argc, char** argv) {
 
+    int rank_num = -1;
     try {
 
+        // The factory from which all other kernel objects are made.
+        yk_factory kfac;
+        
+        // Initalize MPI, etc.
+        auto env = kfac.new_env();
+        rank_num = env->get_rank_index();
+
+        // Local options.
+        for (int i = 1; i < argc; i++)
+            if (string(argv[i]) == "-trace")
+                env->set_trace_enabled(true);
+        
         // Create solution.
         auto soln = kfac.new_solution(env);
 
+        // Apply any YASK command-line options.
+        soln->apply_command_line_options(argc, argv);
+
         // Show output only from last rank.
-        // This is an example of using the rank APIs,
-        // the yask_output_factory, and set_debug_output().
-        int rank_num = env->get_rank_index();
         if (rank_num < env->get_num_ranks() - 1) {
             yk_env::disable_debug_output();
             cout << "Suppressing output on rank " << rank_num << ".\n";
@@ -187,7 +194,7 @@ int main() {
                 continue;
 
             // Create indices describing a subset of the overall domain.
-            vector<idx_t> first_indices, last_indices;
+            vector<idx_t> first_cube_indices, last_cube_indices;
             for (auto dname : var->get_dim_names()) {
                 idx_t first_idx = 0, last_idx = 0;
 
@@ -222,61 +229,69 @@ int main() {
                 }
 
                 // Add indices to index vectors.
-                first_indices.push_back(first_idx);
-                last_indices.push_back(last_idx);
+                first_cube_indices.push_back(first_idx);
+                last_cube_indices.push_back(last_idx);
             }
-            assert(first_indices.size() == ndims);
-            assert(last_indices.size() == ndims);
+            assert(first_cube_indices.size() == ndims);
+            assert(last_cube_indices.size() == ndims);
 
             // Print range of inner cube.
             if (ndims == 0)
                 os << "      inner-cube is a single point for this scalar var.\n";
             else
                 os << "      range of inner-cube: " <<
-                var->format_indices(first_indices) << " to " <<
-                var->format_indices(last_indices) << ".\n";
+                var->format_indices(first_cube_indices) << " to " <<
+                var->format_indices(last_cube_indices) << ".\n";
 
             // Check 2 corners of inner-cube to make sure they were set properly.
-            if (var->are_indices_local(first_indices)) {
-                auto val2 = var->get_element(first_indices);
+            if (var->are_indices_local(first_cube_indices)) {
+                auto val2 = var->get_element(first_cube_indices);
                 os << "      first element in inner cube == " << val2 << ".\n";
                 assert(val2 == outer_val);
             }
             else
-                os << "      first element NOT in rank.\n";
-            if (var->are_indices_local(last_indices)) {
-                auto val2 = var->get_element(last_indices);
+                os << "      first element NOT in this rank.\n";
+            if (var->are_indices_local(last_cube_indices)) {
+                auto val2 = var->get_element(last_cube_indices);
                 os << "      last element in inner cube == " << val2 << ".\n";
                 assert(val2 == outer_val);
             }
             else
-                os << "      last element NOT in rank.\n";
+                os << "      last element NOT in this rank.\n";
 
             // Re-init the values in the inner cube using the indices created above.
+            // Use 'strict_indices = false' because first/last_cube_indices are global,
+            // so part or all of cube may be outside this rank.
             double inner_val = 2.0;
-            bool strict_indices = false; // because first/last_indices are global.
-            idx_t nset = var->set_elements_in_slice_same(inner_val, first_indices, last_indices, strict_indices);
+            bool strict_indices = false;
+            idx_t nset = var->set_elements_in_slice_same(inner_val, first_cube_indices,
+                                                         last_cube_indices, strict_indices);
             os << "      " << nset << " element(s) in inner-cube set to " << inner_val << ".\n";
             
-            // Check values on 2 corners of the inner cube to make sure they were set properly.
-            if (var->are_indices_local(first_indices)) {
-                auto val2 = var->get_element(first_indices);
+            // Check values on 2 corners of the inner cube to make sure they
+            // were set properly.  Only check if element is in this rank.
+            if (var->are_indices_local(first_cube_indices)) {
+                auto val2 = var->get_element(first_cube_indices);
                 os << "      first element in inner cube == " << val2 << ".\n";
                 assert(val2 == inner_val);
             }
-            if (var->are_indices_local(last_indices)) {
-                auto val2 = var->get_element(last_indices);
+            if (var->are_indices_local(last_cube_indices)) {
+                auto val2 = var->get_element(last_cube_indices);
                 os << "      last element in inner cube == " << val2 << ".\n";
                 assert(val2 == inner_val);
             }
 
-            // Add to value in first corner of the inner cube if it's in this rank.
-            nset = var->add_to_element(1.0, first_indices);
-            os << "      " << nset << " element(s) updated.\n";
-
-            // Check to make sure they were updated.
-            if (var->are_indices_local(first_indices)) {
-                auto val2 = var->get_element(first_indices);
+            // Add to value in first corner of the inner cube if it's in
+            // this rank.  Here, we can use the default 'strict_indices =
+            // true' parameter to add_to_element() because the call is
+            // inside the check for local indices.
+            if (var->are_indices_local(first_cube_indices)) {
+                nset = var->add_to_element(1.0, first_cube_indices);
+                os << "      " << nset << " element(s) updated.\n";
+                assert(nset == 1);
+
+                // Check to make sure they were updated.
+                auto val2 = var->get_element(first_cube_indices);
                 os << "      first element in inner cube == " << val2 << ".\n";
                 assert(val2 == inner_val + 1.0);
             }
@@ -295,19 +310,21 @@ int main() {
         // Apply the stencil solution to the data.
         soln->reset_auto_tuner(false);
         env->global_barrier();
-        os << "Running the solution for 1 step...\n";
+        os << "Running the solution on " << env->get_num_ranks() << " rank(s)...\n";
+        os << "Running for 1 step...\n";
         soln->run_solution(0);
-        os << "Running the solution for 4 more steps...\n";
+        os << "Running for 4 more steps...\n";
         soln->run_solution(1, 4);
 
         soln->end_solution();
         soln->get_stats();
-        os << "End of YASK kernel API test.\n";
+        env->finalize();
+        os << "End of YASK C++ kernel API test.\n";
         return 0;
     }
     catch (yask_exception e) {
         cerr << "YASK kernel API test: " << e.get_message() <<
-            " on rank " << env->get_rank_index() << ".\n";
+            " on rank " << rank_num << ".\n";
         return 1;
     }
 }
diff --git a/src/kernel/tests/yask_kernel_api_test.py b/src/kernel/tests/yask_kernel_api_test.py
index 32f796ee..01f92074 100755
--- a/src/kernel/tests/yask_kernel_api_test.py
+++ b/src/kernel/tests/yask_kernel_api_test.py
@@ -86,11 +86,13 @@ def read_var(var, timestep) :
     print("Testing reading var '", var.get_name(), "' at time ", timestep, "...")
     ndarray, first_indices, last_indices, point = make_ndarray(var, timestep)
 
-    print("Reading 1 element...")
+    print("Reading corner elements...")
     val1 = var.get_element(first_indices)
     print("Read value ", val1)
+    val2 = var.get_element(last_indices)
+    print("Read value ", val2)
 
-    print("Reading all element(s) in ndarray...")
+    print("Reading element(s) in ndarray from ", first_indices, " to ", last_indices, "...")
     nread = var.get_elements_in_slice(ndarray.data, first_indices, last_indices)
     print(ndarray)
 
@@ -159,6 +161,10 @@ def init_var(var, timestep) :
 
     # Initalize MPI, etc.
     env = kfac.new_env()
+    if env.get_rank_index() > 0:
+        env.disable_debug_output();
+    else:
+        env.set_trace_enabled(True);
 
     # Create solution.
     soln = kfac.new_solution(env)
@@ -208,24 +214,27 @@ def init_var(var, timestep) :
 
     # Print some info about the solution.
     print("Stencil-solution '", name, "':")
-    print("  Step dimension: ", soln.get_step_dim_name())
-    print("  Domain dimensions: ", soln.get_domain_dim_names())
+    print("  Step dimension:", soln.get_step_dim_name())
+    print("  Domain dimensions:", soln.get_domain_dim_names())
     print("  Vars:")
     for var in soln.get_vars() :
         print("    ", var.get_name(), var.get_dim_names())
         for dname in var.get_dim_names() :
             if dname in soln.get_domain_dim_names() :
                 print("      '", dname, "' allowed domain index range in this rank: ",
-                      var.get_first_rank_alloc_index(dname), " ... ",
+                      var.get_first_rank_alloc_index(dname), "...",
                       var.get_last_rank_alloc_index(dname))
             elif dname == soln.get_step_dim_name() :
                 print("      '", dname, "' allowed step index range: ",
-                      var.get_first_valid_step_index(), " ... ",
+                      var.get_first_valid_step_index(), "...",
                       var.get_last_valid_step_index())
             else :
                 print("      '", dname, "' allowed misc index range: ",
-                      var.get_first_misc_index(dname), " ... ",
+                      var.get_first_misc_index(dname), "...",
                       var.get_last_misc_index(dname))
+        print("      allowed range in all dims: ",
+              var.get_first_local_index_vec(), "...",
+              var.get_last_local_index_vec())
 
     # Init the vars.
     for var in soln.get_vars() :
@@ -292,7 +301,7 @@ def init_var(var, timestep) :
                 last_indices += [var.get_last_misc_index(dname)]
 
         # Init value at one point.
-        nset = var.set_element(15.0, one_indices)
+        nset = var.set_element(15.0, one_indices, False)
         print("Set ", nset, " element(s) in rank ", env.get_rank_index())
 
         # Init the values within the small cube.
@@ -304,22 +313,24 @@ def init_var(var, timestep) :
 
     # Apply the stencil solution to the data.
     env.global_barrier()
-    print("Running the solution for 1 step...")
+    print("Running the solution on", env.get_num_ranks(), "rank(s)");
+    print("Running for 1 step...")
     soln.run_solution(0)
 
     # Print result at timestep 1.
     for var in soln.get_vars() :
         read_var(var, 1)
 
-    print("Running the solution for 10 more steps...")
-    soln.run_solution(1, 10)
+    print("Running for 4 more steps...")
+    soln.run_solution(1, 4)
 
-    # Print final result at timestep 11, assuming update was to t+1.
+    # Print final result at timestep 5, assuming last update was to t+1.
     for var in soln.get_vars() :
-        read_var(var, 11)
+        read_var(var, 5)
 
     soln.end_solution()
     soln.get_stats()
+    env.finalize()
 
     #print("Debug output captured:\n", debug_output.get_string())
     print("End of YASK Python kernel API test.")
diff --git a/src/kernel/yask_main.cpp b/src/kernel/yask_main.cpp
index b48b0585..4515d771 100644
--- a/src/kernel/yask_main.cpp
+++ b/src/kernel/yask_main.cpp
@@ -652,7 +652,7 @@ int main(int argc, char** argv)
         if (!ok)
             exit_yask(1);
 
-        MPI_Finalize();
+        kenv->finalize();
         os << "YASK DONE." << endl << div_line << flush;
     }
     catch (yask_exception& e) {