updated with most recent developments

oop_imsrg/doc/generator.html

@@ -41,10 +41,15 @@
 </font></dt><dd>
 <dl>
 <dt><font face="helvetica, arial"><a href="generator.html#WegnerGenerator">WegnerGenerator</a>
+</font></dt><dd>
+<dl>
+<dt><font face="helvetica, arial"><a href="generator.html#WegnerGenerator3B">WegnerGenerator3B</a>
 </font></dt></dl>
 </dd>
 </dl>
 </dd>
+</dl>
+</dd>
 </dl>
  <p>
 <table width="100%" cellspacing=0 cellpadding=2 border=0 summary="section">
@@ -77,7 +82,8 @@ Data descriptors defined here:<br>
 <tr bgcolor="#ffc8d8"><td rowspan=2><tt>&nbsp;&nbsp;&nbsp;</tt></td>
 <td colspan=2><tt><a href="#WegnerGenerator">WegnerGenerator</a>(h,&nbsp;occ_t)<br>
 &nbsp;<br>
-Calculate&nbsp;Wegner's&nbsp;generator&nbsp;for&nbsp;a&nbsp;normal&nbsp;ordered&nbsp;Hamiltonian.<br>&nbsp;</tt></td></tr>
+Calculate&nbsp;Wegner's&nbsp;generator&nbsp;for&nbsp;a&nbsp;normal&nbsp;ordered&nbsp;Hamiltonian.<br>
+Truncate&nbsp;at&nbsp;two-body&nbsp;interactions.<br>&nbsp;</tt></td></tr>
 <tr><td>&nbsp;</td>
 <td width="100%"><dl><dt>Method resolution order:</dt>
 <dd><a href="generator.html#WegnerGenerator">WegnerGenerator</a></dd>
@@ -93,7 +99,7 @@ Arguments:<br>
 h&nbsp;--&nbsp;Hamiltonian&nbsp;<a href="builtins.html#object">object</a>&nbsp;(must&nbsp;be&nbsp;normal-ordered)<br>
 occ_t&nbsp;--&nbsp;OccupationTensor&nbsp;<a href="builtins.html#object">object</a></tt></dd></dl>
 
-<dl><dt><a name="WegnerGenerator-calc_eta"><strong>calc_eta</strong></a>(self)</dt><dd><tt>Calculate&nbsp;the&nbsp;generator.&nbsp;The&nbsp;terms&nbsp;are&nbsp;defined&nbsp;in&nbsp;An&nbsp;<br>
+<dl><dt><a name="WegnerGenerator-calc_eta"><strong>calc_eta</strong></a>(self)</dt><dd><tt>Calculate&nbsp;the&nbsp;generator.&nbsp;The&nbsp;terms&nbsp;are&nbsp;defined&nbsp;in&nbsp;An<br>
 Advanced&nbsp;Course&nbsp;in&nbsp;Computation&nbsp;Nuclear&nbsp;Physics,&nbsp;Ch.10.<br>
 See&nbsp;also&nbsp;dx.doi.org/10.1016/j.physrep.2015.12.007<br>
 &nbsp;<br>
@@ -102,8 +108,100 @@ Returns:<br>
 (eta1B,&nbsp;--&nbsp;one-body&nbsp;generator<br>
 &nbsp;eta2B)&nbsp;--&nbsp;two-body&nbsp;generator</tt></dd></dl>
 
+<dl><dt><a name="WegnerGenerator-decouple_OD"><strong>decouple_OD</strong></a>(self)</dt><dd><tt>Decouple&nbsp;the&nbsp;off-/diagonal&nbsp;elements&nbsp;from&nbsp;each&nbsp;other&nbsp;in<br>
+the&nbsp;one-&nbsp;and&nbsp;two-body&nbsp;tensors.&nbsp;This&nbsp;procedure&nbsp;is&nbsp;outlined&nbsp;in<br>
+An&nbsp;Advanced&nbsp;Course&nbsp;in&nbsp;Computation&nbsp;Nuclear&nbsp;Physics,&nbsp;Ch.10.<br>
+&nbsp;<br>
+Returns:<br>
+&nbsp;<br>
+(fd,&nbsp;--&nbsp;diagonal&nbsp;part&nbsp;of&nbsp;f<br>
+&nbsp;fod,&nbsp;--&nbsp;off-diagonal&nbsp;part&nbsp;of&nbsp;f<br>
+&nbsp;Gd,&nbsp;--&nbsp;diagonal&nbsp;part&nbsp;of&nbsp;G<br>
+&nbsp;God)&nbsp;--&nbsp;off-diagonal&nbsp;part&nbsp;of&nbsp;G</tt></dd></dl>
+
+<hr>
+Data descriptors defined here:<br>
+<dl><dt><strong>G</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+f&nbsp;--&nbsp;two-body&nbsp;tensor&nbsp;elements&nbsp;(initialized&nbsp;by&nbsp;Hamiltonian&nbsp;object)</tt></dd>
+</dl>
+<dl><dt><strong>f</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+f&nbsp;--&nbsp;one-body&nbsp;tensor&nbsp;elements&nbsp;(initialized&nbsp;by&nbsp;Hamiltonian&nbsp;object)</tt></dd>
+</dl>
+<hr>
+Data descriptors inherited from <a href="generator.html#Generator">Generator</a>:<br>
+<dl><dt><strong>__dict__</strong></dt>
+<dd><tt>dictionary&nbsp;for&nbsp;instance&nbsp;variables&nbsp;(if&nbsp;defined)</tt></dd>
+</dl>
+<dl><dt><strong>__weakref__</strong></dt>
+<dd><tt>list&nbsp;of&nbsp;weak&nbsp;references&nbsp;to&nbsp;the&nbsp;object&nbsp;(if&nbsp;defined)</tt></dd>
+</dl>
+</td></tr></table> <p>
+<table width="100%" cellspacing=0 cellpadding=2 border=0 summary="section">
+<tr bgcolor="#ffc8d8">
+<td colspan=3 valign=bottom>&nbsp;<br>
+<font color="#000000" face="helvetica, arial"><a name="WegnerGenerator3B">class <strong>WegnerGenerator3B</strong></a>(<a href="generator.html#WegnerGenerator">WegnerGenerator</a>)</font></td></tr>
+    
+<tr bgcolor="#ffc8d8"><td rowspan=2><tt>&nbsp;&nbsp;&nbsp;</tt></td>
+<td colspan=2><tt><a href="#WegnerGenerator3B">WegnerGenerator3B</a>(h,&nbsp;occ_t)<br>
+&nbsp;<br>
+Calculate&nbsp;Wegner's&nbsp;generator&nbsp;for&nbsp;a&nbsp;normal&nbsp;ordered&nbsp;Hamiltonian.<br>
+Truncate&nbsp;at&nbsp;three-body&nbsp;interactions.&nbsp;Inherits&nbsp;from&nbsp;<a href="#WegnerGenerator">WegnerGenerator</a>.<br>&nbsp;</tt></td></tr>
+<tr><td>&nbsp;</td>
+<td width="100%"><dl><dt>Method resolution order:</dt>
+<dd><a href="generator.html#WegnerGenerator3B">WegnerGenerator3B</a></dd>
+<dd><a href="generator.html#WegnerGenerator">WegnerGenerator</a></dd>
+<dd><a href="generator.html#Generator">Generator</a></dd>
+<dd><a href="builtins.html#object">builtins.object</a></dd>
+</dl>
+<hr>
+Methods defined here:<br>
+<dl><dt><a name="WegnerGenerator3B-__init__"><strong>__init__</strong></a>(self, h, occ_t)</dt><dd><tt>Class&nbsp;constructor.&nbsp;Instantiate&nbsp;<a href="#WegnerGenerator">WegnerGenerator</a>&nbsp;<a href="builtins.html#object">object</a>.<br>
+&nbsp;<br>
+Arguments:<br>
+&nbsp;<br>
+h&nbsp;--&nbsp;Hamiltonian&nbsp;<a href="builtins.html#object">object</a>&nbsp;(must&nbsp;be&nbsp;normal-ordered)<br>
+occ_t&nbsp;--&nbsp;OccupationTensor&nbsp;<a href="builtins.html#object">object</a></tt></dd></dl>
+
+<dl><dt><a name="WegnerGenerator3B-calc_eta"><strong>calc_eta</strong></a>(self)</dt><dd><tt>Inherits&nbsp;from&nbsp;<a href="#WegnerGenerator">WegnerGenerator</a>.<br>
+&nbsp;<br>
+Calculate&nbsp;the&nbsp;generator.&nbsp;The&nbsp;terms&nbsp;are&nbsp;defined&nbsp;in&nbsp;An<br>
+Advanced&nbsp;Course&nbsp;in&nbsp;Computation&nbsp;Nuclear&nbsp;Physics,&nbsp;Ch.10.<br>
+See&nbsp;also&nbsp;dx.doi.org/10.1016/j.physrep.2015.12.007<br>
+&nbsp;<br>
+Returns:<br>
+&nbsp;<br>
+(eta1B,&nbsp;--&nbsp;one-body&nbsp;generator<br>
+&nbsp;eta2B,&nbsp;--&nbsp;two-body&nbsp;generator<br>
+&nbsp;eta3B)&nbsp;--&nbsp;three-body&nbsp;generator</tt></dd></dl>
+
+<dl><dt><a name="WegnerGenerator3B-decouple_OD"><strong>decouple_OD</strong></a>(self)</dt><dd><tt>Inherits&nbsp;from&nbsp;<a href="#WegnerGenerator">WegnerGenerator</a>.<br>
+&nbsp;<br>
+Decouple&nbsp;the&nbsp;off-/diagonal&nbsp;elements&nbsp;from&nbsp;each&nbsp;other&nbsp;in<br>
+the&nbsp;one-&nbsp;and&nbsp;two-body&nbsp;tensors.&nbsp;This&nbsp;procedure&nbsp;is&nbsp;outlined&nbsp;in<br>
+An&nbsp;Advanced&nbsp;Course&nbsp;in&nbsp;Computation&nbsp;Nuclear&nbsp;Physics,&nbsp;Ch.10.<br>
+&nbsp;<br>
+Returns:<br>
+&nbsp;<br>
+(fd,&nbsp;--&nbsp;diagonal&nbsp;part&nbsp;of&nbsp;f<br>
+&nbsp;fod,&nbsp;--&nbsp;off-diagonal&nbsp;part&nbsp;of&nbsp;f<br>
+&nbsp;Gd,&nbsp;--&nbsp;diagonal&nbsp;part&nbsp;of&nbsp;G<br>
+&nbsp;God,&nbsp;--&nbsp;off-diagonal&nbsp;part&nbsp;of&nbsp;G<br>
+&nbsp;Wd,&nbsp;--&nbsp;diagonal&nbsp;part&nbsp;of&nbsp;W<br>
+&nbsp;Wod)&nbsp;--&nbsp;off-diagonal&nbsp;part&nbsp;of&nbsp;W</tt></dd></dl>
+
 <hr>
 Data descriptors defined here:<br>
+<dl><dt><strong>W</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+W&nbsp;--&nbsp;three-body&nbsp;tensor&nbsp;elements&nbsp;(initialized&nbsp;by&nbsp;Hamiltonian&nbsp;object)</tt></dd>
+</dl>
+<hr>
+Data descriptors inherited from <a href="generator.html#WegnerGenerator">WegnerGenerator</a>:<br>
 <dl><dt><strong>G</strong></dt>
 <dd><tt>Returns:<br>
 &nbsp;<br>

oop_imsrg/doc/hamiltonian.html

@@ -97,6 +97,7 @@ n_particle_states -- number of particles states in
 &nbsp;<br>
 Keyword&nbsp;arguments:<br>
 &nbsp;<br>
+ref&nbsp;--&nbsp;the&nbsp;reference&nbsp;state.&nbsp;must&nbsp;match&nbsp;dimensions&nbsp;imposed&nbsp;by&nbsp;arugments&nbsp;(default:&nbsp;[1,1,1,1,0,0,0,0])<br>
 d&nbsp;--&nbsp;the&nbsp;energy&nbsp;level&nbsp;spacing&nbsp;(default:&nbsp;1.0)<br>
 g&nbsp;--&nbsp;the&nbsp;pairing&nbsp;strength&nbsp;(default:&nbsp;0.5)<br>
 pb&nbsp;--&nbsp;strength&nbsp;of&nbsp;the&nbsp;pair-breaking&nbsp;term&nbsp;(operates&nbsp;in&nbsp;double&nbsp;particle&nbsp;basis)&nbsp;(default:&nbsp;0.0)</tt></dd></dl>

oop_imsrg/doc/main.html

@@ -8,7 +8,7 @@
 <td valign=bottom>&nbsp;<br>
 <font color="#ffffff" face="helvetica, arial">&nbsp;<br><big><big><strong>main</strong></big></big></font></td
 ><td align=right valign=bottom
-><font color="#ffffff" face="helvetica, arial"><a href=".">index</a><br><a href="file:c%3A%5Cusers%5Cdavison%5Cresearch%5Cim-srg_tensorflow%5Coop_imsrg%5Cmain.py">c:\users\davison\research\im-srg_tensorflow\oop_imsrg\main.py</a></font></td></tr></table>
+><font color="#ffffff" face="helvetica, arial"><a href=".">index</a><br><a href="file:c%3A%5Cusers%5Cdavison%5Cresearch%5Cim-srg_tensorflow%5Cmain.py">c:\users\davison\research\im-srg_tensorflow\main.py</a></font></td></tr></table>
     <p><tt>#&nbsp;Main&nbsp;program&nbsp;for&nbsp;IM-SRG.</tt></p>
 <p>
 <table width="100%" cellspacing=0 cellpadding=2 border=0 summary="section">
@@ -18,7 +18,7 @@
     
 <tr><td bgcolor="#aa55cc"><tt>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</tt></td><td>&nbsp;</td>
 <td width="100%"><table width="100%" summary="list"><tr><td width="25%" valign=top><a href="tensornetwork.backends.html">tensornetwork.backends</a><br>
-<a href="ci_pairing.cipy_pairing_plus_ph.html">ci_pairing.cipy_pairing_plus_ph</a><br>
+<a href="oop_imsrg.ci_pairing.cipy_pairing_plus_ph.html">oop_imsrg.ci_pairing.cipy_pairing_plus_ph</a><br>
 <a href="tensornetwork.config.html">tensornetwork.config</a><br>
 <a href="tensornetwork.contractors.html">tensornetwork.contractors</a><br>
 <a href="glob.html">glob</a><br>

oop_imsrg/doc/occupation_tensors.html

@@ -66,12 +66,12 @@ Data descriptors defined here:<br>
 <dd><tt>list&nbsp;of&nbsp;weak&nbsp;references&nbsp;to&nbsp;the&nbsp;object&nbsp;(if&nbsp;defined)</tt></dd>
 </dl>
 <dl><dt><strong>occA</strong></dt>
-<dd><tt>Returns:&nbsp;<br>
+<dd><tt>Returns:<br>
 &nbsp;<br>
 occA&nbsp;--&nbsp;represents&nbsp;n_a&nbsp;-&nbsp;n_b.</tt></dd>
 </dl>
 <dl><dt><strong>occB</strong></dt>
-<dd><tt>Returns:&nbsp;<br>
+<dd><tt>Returns:<br>
 &nbsp;<br>
 occB&nbsp;--&nbsp;represents&nbsp;1&nbsp;-&nbsp;n_a&nbsp;-&nbsp;n_b.</tt></dd>
 </dl>
@@ -85,5 +85,36 @@ occC&nbsp;--&nbsp;represents&nbsp;n_a*n_b&nbsp;+&nbsp;(1-n_a-n_b)*n_c</tt></dd>
 &nbsp;<br>
 occD&nbsp;--&nbsp;represents&nbsp;na*nb*(1-nc-nd)&nbsp;+&nbsp;na*nb*nc*nd</tt></dd>
 </dl>
+<dl><dt><strong>occE</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occE&nbsp;--&nbsp;represents&nbsp;n_a*n_b*n_c*(1-n_d)*(1-n_e)*(1-n_f)</tt></dd>
+</dl>
+<dl><dt><strong>occF</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occF&nbsp;--&nbsp;represents&nbsp;n_a*n_b*(1-n_c)*(1-n_d)*(1-n_e)&nbsp;+<br>
+&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(1-n_a)*(1-n_b)*n_c*n_d*n_e</tt></dd>
+</dl>
+<dl><dt><strong>occG</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occG&nbsp;--&nbsp;represents&nbsp;n_a*(1-n_b)*(1-n_c)&nbsp;+&nbsp;(1-n_a)*n_b*n_c</tt></dd>
+</dl>
+<dl><dt><strong>occH</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occH&nbsp;--&nbsp;represents&nbsp;n_a*(1-n_b)*(1-n_c)*(1-n_d)&nbsp;-&nbsp;(1-n_a)*n_b*n_c*n_d</tt></dd>
+</dl>
+<dl><dt><strong>occI</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occI&nbsp;--&nbsp;represents&nbsp;(1-n_a)*(1-n_b)*n_c*n_d&nbsp;-&nbsp;n_a*n_b*(1-n_c)*(1-n_d)</tt></dd>
+</dl>
+<dl><dt><strong>occJ</strong></dt>
+<dd><tt>Returns:<br>
+&nbsp;<br>
+occJ&nbsp;--&nbsp;represents&nbsp;n_a*n_b*n_c&nbsp;+&nbsp;(1-n_a)*(1-n_b)*(1-n_c)</tt></dd>
+</dl>
 </td></tr></table></td></tr></table>
 </body></html>
\ No newline at end of file