From 8230ec63c2b7732ddea5d6197a9cef430ffa0431 Mon Sep 17 00:00:00 2001
From: GitHub CI Documentation builder <username@users.noreply.github.com>
Date: Mon, 20 Nov 2023 14:27:52 +0000
Subject: [PATCH] Update docs

---
 classSpeciesInformation_1_1Impl-members.html |   2 +-
 classSpeciesInformation_1_1Impl.html         |  30 ++---
 classSpeciesInformation_1_1Impl.js           |  10 +-
 navtreeindex3.js                             |  10 +-
 species__info_8hpp_source.html               | 127 ++++++++++---------
 5 files changed, 91 insertions(+), 88 deletions(-)

diff --git a/classSpeciesInformation_1_1Impl-members.html b/classSpeciesInformation_1_1Impl-members.html
index a0c9937a2..1c15f72e6 100644
--- a/classSpeciesInformation_1_1Impl-members.html
+++ b/classSpeciesInformation_1_1Impl-members.html
@@ -111,7 +111,7 @@
   <tr bgcolor="#f0f0f0" class="even"><td class="entry"><b>ielec</b>() const (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
   <tr bgcolor="#f0f0f0"><td class="entry"><b>Impl</b> (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">friend</span></td></tr>
   <tr bgcolor="#f0f0f0" class="even"><td class="entry"><b>Impl</b>(Impl&lt; OMemorySpace &gt; const &amp;impl) (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span><span class="mlabel">explicit</span></td></tr>
-  <tr bgcolor="#f0f0f0"><td class="entry"><b>Impl</b>(ddc::Chunk&lt; int, discrete_domain_type &gt; charge, ddc::Chunk&lt; double, discrete_domain_type &gt; mass, ddc::Chunk&lt; double, discrete_domain_type &gt; perturb_amplitude, ddc::Chunk&lt; int, discrete_domain_type &gt; perturb_mode) (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
+  <tr bgcolor="#f0f0f0"><td class="entry"><b>Impl</b>(ddc::Chunk&lt; int, discrete_domain_type, ddc::KokkosAllocator&lt; int, MemorySpace &gt;&gt; charge, ddc::Chunk&lt; double, discrete_domain_type, ddc::KokkosAllocator&lt; double, MemorySpace &gt;&gt; mass, ddc::Chunk&lt; double, discrete_domain_type, ddc::KokkosAllocator&lt; double, MemorySpace &gt;&gt; perturb_amplitude, ddc::Chunk&lt; int, discrete_domain_type, ddc::KokkosAllocator&lt; int, MemorySpace &gt;&gt; perturb_mode) (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
   <tr bgcolor="#f0f0f0" class="even"><td class="entry"><b>masses</b>() const (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
   <tr bgcolor="#f0f0f0"><td class="entry"><b>perturb_amplitudes</b>() const (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
   <tr bgcolor="#f0f0f0" class="even"><td class="entry"><b>perturb_modes</b>() const (defined in <a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a>)</td><td class="entry"><a class="el" href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl&lt; MemorySpace &gt;</a></td><td class="entry"><span class="mlabel">inline</span></td></tr>
diff --git a/classSpeciesInformation_1_1Impl.html b/classSpeciesInformation_1_1Impl.html
index 506bd34d6..c28fe43ee 100644
--- a/classSpeciesInformation_1_1Impl.html
+++ b/classSpeciesInformation_1_1Impl.html
@@ -120,24 +120,24 @@
 template&lt;class OMemorySpace &gt; </td></tr>
 <tr class="memitem:a788f54c43e25317f947aa8be9716c3bc"><td class="memTemplItemLeft" align="right" valign="top">&#160;</td><td class="memTemplItemRight" valign="bottom"><b>Impl</b> (<a class="el" href="classSpeciesInformation_1_1Impl.html">Impl</a>&lt; OMemorySpace &gt; const &amp;impl)</td></tr>
 <tr class="separator:a788f54c43e25317f947aa8be9716c3bc"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:aee061508e334dfda1ec99068b18d8909"><td class="memItemLeft" align="right" valign="top"><a id="aee061508e334dfda1ec99068b18d8909"></a>
-&#160;</td><td class="memItemRight" valign="bottom"><b>Impl</b> (ddc::Chunk&lt; int, discrete_domain_type &gt; charge, ddc::Chunk&lt; double, discrete_domain_type &gt; mass, ddc::Chunk&lt; double, discrete_domain_type &gt; perturb_amplitude, ddc::Chunk&lt; int, discrete_domain_type &gt; perturb_mode)</td></tr>
-<tr class="separator:aee061508e334dfda1ec99068b18d8909"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:a926e1d36257535d9c850fadd85f45dd2"><td class="memItemLeft" align="right" valign="top"><a id="a926e1d36257535d9c850fadd85f45dd2"></a>
+&#160;</td><td class="memItemRight" valign="bottom"><b>Impl</b> (ddc::Chunk&lt; int, discrete_domain_type, ddc::KokkosAllocator&lt; int, MemorySpace &gt;&gt; charge, ddc::Chunk&lt; double, discrete_domain_type, ddc::KokkosAllocator&lt; double, MemorySpace &gt;&gt; mass, ddc::Chunk&lt; double, discrete_domain_type, ddc::KokkosAllocator&lt; double, MemorySpace &gt;&gt; perturb_amplitude, ddc::Chunk&lt; int, discrete_domain_type, ddc::KokkosAllocator&lt; int, MemorySpace &gt;&gt; perturb_mode)</td></tr>
+<tr class="separator:a926e1d36257535d9c850fadd85f45dd2"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:ad069f5667f3d88ab36856bd941c6e937"><td class="memItemLeft" align="right" valign="top"><a id="ad069f5667f3d88ab36856bd941c6e937"></a>
 KOKKOS_FUNCTION discrete_element_type&#160;</td><td class="memItemRight" valign="bottom"><b>ielec</b> () const</td></tr>
 <tr class="separator:ad069f5667f3d88ab36856bd941c6e937"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a22075762bac76656386615ff21f71bd8"><td class="memItemLeft" align="right" valign="top"><a id="a22075762bac76656386615ff21f71bd8"></a>
-ddc::ChunkSpan&lt; const int, discrete_domain_type &gt;&#160;</td><td class="memItemRight" valign="bottom"><b>charges</b> () const</td></tr>
-<tr class="separator:a22075762bac76656386615ff21f71bd8"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a0cee25bda8e8c98cb2e4f6b1f1bdfbaa"><td class="memItemLeft" align="right" valign="top"><a id="a0cee25bda8e8c98cb2e4f6b1f1bdfbaa"></a>
-ddc::ChunkSpan&lt; const double, discrete_domain_type &gt;&#160;</td><td class="memItemRight" valign="bottom"><b>masses</b> () const</td></tr>
-<tr class="separator:a0cee25bda8e8c98cb2e4f6b1f1bdfbaa"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a447bc6c55e4a471376bfdfcbd04db742"><td class="memItemLeft" align="right" valign="top"><a id="a447bc6c55e4a471376bfdfcbd04db742"></a>
-ddc::ChunkSpan&lt; const double, discrete_domain_type &gt;&#160;</td><td class="memItemRight" valign="bottom"><b>perturb_amplitudes</b> () const</td></tr>
-<tr class="separator:a447bc6c55e4a471376bfdfcbd04db742"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a3bf40fa531902e47b82ce46109df9259"><td class="memItemLeft" align="right" valign="top"><a id="a3bf40fa531902e47b82ce46109df9259"></a>
-ddc::ChunkSpan&lt; const int, discrete_domain_type &gt;&#160;</td><td class="memItemRight" valign="bottom"><b>perturb_modes</b> () const</td></tr>
-<tr class="separator:a3bf40fa531902e47b82ce46109df9259"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:a6fa19e2bcac84772f4de38323e4430e3"><td class="memItemLeft" align="right" valign="top"><a id="a6fa19e2bcac84772f4de38323e4430e3"></a>
+auto&#160;</td><td class="memItemRight" valign="bottom"><b>charges</b> () const</td></tr>
+<tr class="separator:a6fa19e2bcac84772f4de38323e4430e3"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:ab5fea032a19715a668e86acfa0a4e93a"><td class="memItemLeft" align="right" valign="top"><a id="ab5fea032a19715a668e86acfa0a4e93a"></a>
+auto&#160;</td><td class="memItemRight" valign="bottom"><b>masses</b> () const</td></tr>
+<tr class="separator:ab5fea032a19715a668e86acfa0a4e93a"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:a0eb655bf420dc3a05d2cbb1e176cf279"><td class="memItemLeft" align="right" valign="top"><a id="a0eb655bf420dc3a05d2cbb1e176cf279"></a>
+auto&#160;</td><td class="memItemRight" valign="bottom"><b>perturb_amplitudes</b> () const</td></tr>
+<tr class="separator:a0eb655bf420dc3a05d2cbb1e176cf279"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:aeaa6f310ad69c7e977214d455bb9a039"><td class="memItemLeft" align="right" valign="top"><a id="aeaa6f310ad69c7e977214d455bb9a039"></a>
+auto&#160;</td><td class="memItemRight" valign="bottom"><b>perturb_modes</b> () const</td></tr>
+<tr class="separator:aeaa6f310ad69c7e977214d455bb9a039"><td class="memSeparator" colspan="2">&#160;</td></tr>
 </table>
 <hr/>The documentation for this class was generated from the following file:<ul>
 <li>/__w/gyselalibxx/gyselalibxx/code_branch/src/speciesinfo/<a class="el" href="species__info_8hpp_source.html">species_info.hpp</a></li>
diff --git a/classSpeciesInformation_1_1Impl.js b/classSpeciesInformation_1_1Impl.js
index d1f1afe45..d3020eae0 100644
--- a/classSpeciesInformation_1_1Impl.js
+++ b/classSpeciesInformation_1_1Impl.js
@@ -2,11 +2,11 @@ var classSpeciesInformation_1_1Impl =
 [
     [ "discrete_dimension_type", "classSpeciesInformation_1_1Impl.html#a79dd22b07e0b671a5f85b7a1d47c10a4", null ],
     [ "Impl", "classSpeciesInformation_1_1Impl.html#a788f54c43e25317f947aa8be9716c3bc", null ],
-    [ "Impl", "classSpeciesInformation_1_1Impl.html#aee061508e334dfda1ec99068b18d8909", null ],
+    [ "Impl", "classSpeciesInformation_1_1Impl.html#a926e1d36257535d9c850fadd85f45dd2", null ],
     [ "ielec", "classSpeciesInformation_1_1Impl.html#ad069f5667f3d88ab36856bd941c6e937", null ],
-    [ "charges", "classSpeciesInformation_1_1Impl.html#a22075762bac76656386615ff21f71bd8", null ],
-    [ "masses", "classSpeciesInformation_1_1Impl.html#a0cee25bda8e8c98cb2e4f6b1f1bdfbaa", null ],
-    [ "perturb_amplitudes", "classSpeciesInformation_1_1Impl.html#a447bc6c55e4a471376bfdfcbd04db742", null ],
-    [ "perturb_modes", "classSpeciesInformation_1_1Impl.html#a3bf40fa531902e47b82ce46109df9259", null ],
+    [ "charges", "classSpeciesInformation_1_1Impl.html#a6fa19e2bcac84772f4de38323e4430e3", null ],
+    [ "masses", "classSpeciesInformation_1_1Impl.html#ab5fea032a19715a668e86acfa0a4e93a", null ],
+    [ "perturb_amplitudes", "classSpeciesInformation_1_1Impl.html#a0eb655bf420dc3a05d2cbb1e176cf279", null ],
+    [ "perturb_modes", "classSpeciesInformation_1_1Impl.html#aeaa6f310ad69c7e977214d455bb9a039", null ],
     [ "Impl", "classSpeciesInformation_1_1Impl.html#a219fbcef554ce8d99265d545825c18dc", null ]
 ];
\ No newline at end of file
diff --git a/navtreeindex3.js b/navtreeindex3.js
index 5f309e6a6..3dd2c7d03 100644
--- a/navtreeindex3.js
+++ b/navtreeindex3.js
@@ -9,15 +9,15 @@ var NAVTREEINDEX3 =
 "classSpeciesInformation.html#a8e46a5b4955a5164dbc6703fd7a86eb4":[12,122,2],
 "classSpeciesInformation.html#af77dd55f79155df4dab9bf490d900a07":[12,122,1],
 "classSpeciesInformation_1_1Impl.html":[12,122,0],
-"classSpeciesInformation_1_1Impl.html#a0cee25bda8e8c98cb2e4f6b1f1bdfbaa":[12,122,0,5],
+"classSpeciesInformation_1_1Impl.html#a0eb655bf420dc3a05d2cbb1e176cf279":[12,122,0,6],
 "classSpeciesInformation_1_1Impl.html#a219fbcef554ce8d99265d545825c18dc":[12,122,0,8],
-"classSpeciesInformation_1_1Impl.html#a22075762bac76656386615ff21f71bd8":[12,122,0,4],
-"classSpeciesInformation_1_1Impl.html#a3bf40fa531902e47b82ce46109df9259":[12,122,0,7],
-"classSpeciesInformation_1_1Impl.html#a447bc6c55e4a471376bfdfcbd04db742":[12,122,0,6],
+"classSpeciesInformation_1_1Impl.html#a6fa19e2bcac84772f4de38323e4430e3":[12,122,0,4],
 "classSpeciesInformation_1_1Impl.html#a788f54c43e25317f947aa8be9716c3bc":[12,122,0,1],
 "classSpeciesInformation_1_1Impl.html#a79dd22b07e0b671a5f85b7a1d47c10a4":[12,122,0,0],
+"classSpeciesInformation_1_1Impl.html#a926e1d36257535d9c850fadd85f45dd2":[12,122,0,2],
+"classSpeciesInformation_1_1Impl.html#ab5fea032a19715a668e86acfa0a4e93a":[12,122,0,5],
 "classSpeciesInformation_1_1Impl.html#ad069f5667f3d88ab36856bd941c6e937":[12,122,0,3],
-"classSpeciesInformation_1_1Impl.html#aee061508e334dfda1ec99068b18d8909":[12,122,0,2],
+"classSpeciesInformation_1_1Impl.html#aeaa6f310ad69c7e977214d455bb9a039":[12,122,0,7],
 "classSplineBoundaryValue.html":[12,123],
 "classSplineBoundaryValue.html#a772793f39da3ff86099ca45d56d3a4b6":[12,123,1],
 "classSplineBoundaryValue.html#ab8c255bfff4881701559b44dbe568db1":[12,123,0],
diff --git a/species__info_8hpp_source.html b/species__info_8hpp_source.html
index 605565193..db526e8d0 100644
--- a/species__info_8hpp_source.html
+++ b/species__info_8hpp_source.html
@@ -166,69 +166,72 @@
 <div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;            ddc::deepcopy(m_perturb_mode, impl.m_perturb_mode);</div>
 <div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;        }</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">   64</span>&#160;        <a class="code" href="classSpeciesInformation_1_1Impl.html">Impl</a>(ddc::Chunk&lt;int, discrete_domain_type&gt; charge,</div>
-<div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;             ddc::Chunk&lt;double, discrete_domain_type&gt; mass,</div>
-<div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;             ddc::Chunk&lt;double, discrete_domain_type&gt; perturb_amplitude,</div>
-<div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;             ddc::Chunk&lt;int, discrete_domain_type&gt; perturb_mode)</div>
-<div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;            : m_charge(std::move(charge))</div>
-<div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;            , m_mass(std::move(mass))</div>
-<div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;            , m_perturb_amplitude(std::move(perturb_amplitude))</div>
-<div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;            , m_perturb_mode(std::move(perturb_mode))</div>
-<div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;        {</div>
-<div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;            assert(charge.size() &gt;= 2);</div>
-<div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;            <span class="keywordtype">bool</span> electron_found = <span class="keyword">false</span>;</div>
-<div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;            <span class="keywordflow">for</span> (discrete_element_type <span class="keyword">const</span> isp : m_charge.domain()) {</div>
-<div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;                <span class="keywordflow">if</span> (m_charge(isp) == -1) {</div>
-<div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;                    electron_found = <span class="keyword">true</span>;</div>
-<div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;                    m_ielec = isp;</div>
-<div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;                }</div>
-<div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;            }</div>
-<div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;            <span class="keywordflow">if</span> (!electron_found) {</div>
-<div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;                <span class="keywordflow">throw</span> std::runtime_error(<span class="stringliteral">&quot;electron not found&quot;</span>);</div>
+<div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;        <a class="code" href="classSpeciesInformation_1_1Impl.html">Impl</a>(ddc::Chunk&lt;<span class="keywordtype">int</span>, discrete_domain_type, ddc::KokkosAllocator&lt;int, MemorySpace&gt;&gt; charge,</div>
+<div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;             ddc::Chunk&lt;<span class="keywordtype">double</span>, discrete_domain_type, ddc::KokkosAllocator&lt;double, MemorySpace&gt;&gt;</div>
+<div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;                     mass,</div>
+<div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;             ddc::Chunk&lt;<span class="keywordtype">double</span>, discrete_domain_type, ddc::KokkosAllocator&lt;double, MemorySpace&gt;&gt;</div>
+<div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;                     perturb_amplitude,</div>
+<div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;             ddc::Chunk&lt;<span class="keywordtype">int</span>, discrete_domain_type, ddc::KokkosAllocator&lt;int, MemorySpace&gt;&gt;</div>
+<div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;                     perturb_mode)</div>
+<div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;            : m_charge(std::move(charge))</div>
+<div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;            , m_mass(std::move(mass))</div>
+<div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;            , m_perturb_amplitude(std::move(perturb_amplitude))</div>
+<div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;            , m_perturb_mode(std::move(perturb_mode))</div>
+<div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;        {</div>
+<div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;            assert(charge.size() &gt;= 2);</div>
+<div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;            <span class="keywordtype">bool</span> electron_found = <span class="keyword">false</span>;</div>
+<div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;            <span class="keywordflow">for</span> (discrete_element_type <span class="keyword">const</span> isp : m_charge.domain()) {</div>
+<div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;                <span class="keywordflow">if</span> (m_charge(isp) == -1) {</div>
+<div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;                    electron_found = <span class="keyword">true</span>;</div>
+<div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;                    m_ielec = isp;</div>
+<div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;                }</div>
 <div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;            }</div>
-<div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;        }</div>
-<div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160; </div>
-<div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;        KOKKOS_FUNCTION discrete_element_type ielec()<span class="keyword"> const</span></div>
-<div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;<span class="keyword">        </span>{</div>
-<div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;            <span class="keywordflow">return</span> m_ielec;</div>
-<div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;        }</div>
-<div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160; </div>
-<div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;        ddc::ChunkSpan&lt;const int, discrete_domain_type&gt; charges()<span class="keyword"> const</span></div>
-<div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;<span class="keyword">        </span>{</div>
-<div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;            <span class="keywordflow">return</span> m_charge;</div>
-<div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;        }</div>
-<div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160; </div>
-<div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;        ddc::ChunkSpan&lt;const double, discrete_domain_type&gt; masses()<span class="keyword"> const</span></div>
-<div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;<span class="keyword">        </span>{</div>
-<div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160;            <span class="keywordflow">return</span> m_mass;</div>
-<div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;        }</div>
-<div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160; </div>
-<div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;        ddc::ChunkSpan&lt;const double, discrete_domain_type&gt; perturb_amplitudes()<span class="keyword"> const</span></div>
-<div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;<span class="keyword">        </span>{</div>
-<div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160;            <span class="keywordflow">return</span> m_perturb_amplitude;</div>
-<div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;        }</div>
-<div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160; </div>
-<div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;        ddc::ChunkSpan&lt;const int, discrete_domain_type&gt; perturb_modes()<span class="keyword"> const</span></div>
-<div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;<span class="keyword">        </span>{</div>
-<div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;            <span class="keywordflow">return</span> m_perturb_mode;</div>
-<div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;        }</div>
-<div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;    };</div>
-<div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;};</div>
-<div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160; </div>
-<div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;KOKKOS_INLINE_FUNCTION ddc::DiscreteElement&lt;SpeciesInformation&gt; ielec()</div>
-<div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;{</div>
-<div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;    <span class="keywordflow">return</span> ddc::discrete_space&lt;SpeciesInformation&gt;().ielec();</div>
-<div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;}</div>
-<div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160; </div>
-<div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">int</span> charge(ddc::DiscreteElement&lt;SpeciesInformation&gt; <span class="keyword">const</span> isp)</div>
-<div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;{</div>
-<div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;    <span class="keywordflow">return</span> ddc::discrete_space&lt;SpeciesInformation&gt;().charges()(isp);</div>
-<div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;}</div>
-<div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160; </div>
-<div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">double</span> mass(ddc::DiscreteElement&lt;SpeciesInformation&gt; <span class="keyword">const</span> isp)</div>
-<div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;{</div>
-<div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;    <span class="keywordflow">return</span> ddc::discrete_space&lt;SpeciesInformation&gt;().masses()(isp);</div>
-<div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;}</div>
+<div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;            <span class="keywordflow">if</span> (!electron_found) {</div>
+<div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;                <span class="keywordflow">throw</span> std::runtime_error(<span class="stringliteral">&quot;electron not found&quot;</span>);</div>
+<div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;            }</div>
+<div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;        }</div>
+<div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160; </div>
+<div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;        KOKKOS_FUNCTION discrete_element_type ielec()<span class="keyword"> const</span></div>
+<div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;<span class="keyword">        </span>{</div>
+<div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;            <span class="keywordflow">return</span> m_ielec;</div>
+<div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;        }</div>
+<div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160; </div>
+<div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;        <span class="keyword">auto</span> charges()<span class="keyword"> const</span></div>
+<div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;<span class="keyword">        </span>{</div>
+<div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;            <span class="keywordflow">return</span> m_charge.span_view();</div>
+<div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;        }</div>
+<div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160; </div>
+<div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;        <span class="keyword">auto</span> masses()<span class="keyword"> const</span></div>
+<div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;<span class="keyword">        </span>{</div>
+<div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;            <span class="keywordflow">return</span> m_mass.span_view();</div>
+<div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;        }</div>
+<div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160; </div>
+<div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;        <span class="keyword">auto</span> perturb_amplitudes()<span class="keyword"> const</span></div>
+<div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;<span class="keyword">        </span>{</div>
+<div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;            <span class="keywordflow">return</span> m_perturb_amplitude.span_view();</div>
+<div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;        }</div>
+<div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160; </div>
+<div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;        <span class="keyword">auto</span> perturb_modes()<span class="keyword"> const</span></div>
+<div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;<span class="keyword">        </span>{</div>
+<div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;            <span class="keywordflow">return</span> m_perturb_mode.span_view();</div>
+<div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;        }</div>
+<div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;    };</div>
+<div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;};</div>
+<div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160; </div>
+<div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;KOKKOS_INLINE_FUNCTION ddc::DiscreteElement&lt;SpeciesInformation&gt; ielec()</div>
+<div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160;{</div>
+<div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;    <span class="keywordflow">return</span> ddc::discrete_space&lt;SpeciesInformation&gt;().ielec();</div>
+<div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;}</div>
+<div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160; </div>
+<div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">int</span> charge(ddc::DiscreteElement&lt;SpeciesInformation&gt; <span class="keyword">const</span> isp)</div>
+<div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;{</div>
+<div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;    <span class="keywordflow">return</span> ddc::host_discrete_space&lt;SpeciesInformation&gt;().charges()(isp);</div>
+<div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;}</div>
+<div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160; </div>
+<div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">double</span> mass(ddc::DiscreteElement&lt;SpeciesInformation&gt; <span class="keyword">const</span> isp)</div>
+<div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;{</div>
+<div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;    <span class="keywordflow">return</span> ddc::host_discrete_space&lt;SpeciesInformation&gt;().masses()(isp);</div>
+<div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;}</div>
 <div class="ttc" id="aclassSpeciesInformation_1_1Impl_html"><div class="ttname"><a href="classSpeciesInformation_1_1Impl.html">SpeciesInformation::Impl</a></div><div class="ttdef"><b>Definition:</b> species_info.hpp:26</div></div>
 <div class="ttc" id="aclassSpeciesInformation_html"><div class="ttname"><a href="classSpeciesInformation.html">SpeciesInformation</a></div><div class="ttdef"><b>Definition:</b> species_info.hpp:8</div></div>
 </div><!-- fragment --></div><!-- contents -->