Lightweight fields (#117)

* sph C

* sph python

* forgotten ultra

* ultra default

* ultra not needed

* mask_only

* mask only sph

* flat sky

* missing C flat tests

* just missing C car

* all C tests done

* python tests

* flaked

* docs

doc/build/html/functions.html

@@ -231,6 +231,10 @@ <h3><a class="anchor" id="index_i"></a>- i -</h3><ul>
 
 
 <h3><a class="anchor" id="index_l"></a>- l -</h3><ul>
+<li>lite
+: <a class="el" href="structnmt__field.html#a73caf9e696441d5a2c039d1d1d8bbc8d">nmt_field</a>
+, <a class="el" href="structnmt__field__flat.html#af26656930f7370e26559c6d9f0c3e853">nmt_field_flat</a>
+</li>
 <li>lmax
 : <a class="el" href="structnmt__covar__workspace.html#a29e3d7fe73f1573100318d21f1958875">nmt_covar_workspace</a>
 , <a class="el" href="structnmt__field.html#a8eb358139d9bddaa8b65dc36267ecad4">nmt_field</a>
@@ -264,6 +268,10 @@ <h3><a class="anchor" id="index_m"></a>- m -</h3><ul>
 : <a class="el" href="structnmt__field.html#aebdba52436ee80717acfb3ab9878ef69">nmt_field</a>
 , <a class="el" href="structnmt__field__flat.html#af88cd6bbdc9b3ce33bd51ace99275ed4">nmt_field_flat</a>
 </li>
+<li>mask_only
+: <a class="el" href="structnmt__field.html#ae31c0f09c4732df7c5bfb4ea86c292f7">nmt_field</a>
+, <a class="el" href="structnmt__field__flat.html#a742e5681a04b173a1f0ea63a381ad6bd">nmt_field_flat</a>
+</li>
 <li>matrix_M
 : <a class="el" href="structnmt__field.html#a111d5454c2cd405de1a7cbba122923ee">nmt_field</a>
 , <a class="el" href="structnmt__field__flat.html#a60bfe576ce32ebb6d5beedf5de42aebc">nmt_field_flat</a>

doc/build/html/functions_vars.html

@@ -231,6 +231,10 @@ <h3><a class="anchor" id="index_i"></a>- i -</h3><ul>
 
 
 <h3><a class="anchor" id="index_l"></a>- l -</h3><ul>
+<li>lite
+: <a class="el" href="structnmt__field.html#a73caf9e696441d5a2c039d1d1d8bbc8d">nmt_field</a>
+, <a class="el" href="structnmt__field__flat.html#af26656930f7370e26559c6d9f0c3e853">nmt_field_flat</a>
+</li>
 <li>lmax
 : <a class="el" href="structnmt__covar__workspace.html#a29e3d7fe73f1573100318d21f1958875">nmt_covar_workspace</a>
 , <a class="el" href="structnmt__field.html#a8eb358139d9bddaa8b65dc36267ecad4">nmt_field</a>
@@ -264,6 +268,10 @@ <h3><a class="anchor" id="index_m"></a>- m -</h3><ul>
 : <a class="el" href="structnmt__field.html#aebdba52436ee80717acfb3ab9878ef69">nmt_field</a>
 , <a class="el" href="structnmt__field__flat.html#af88cd6bbdc9b3ce33bd51ace99275ed4">nmt_field_flat</a>
 </li>
+<li>mask_only
+: <a class="el" href="structnmt__field.html#ae31c0f09c4732df7c5bfb4ea86c292f7">nmt_field</a>
+, <a class="el" href="structnmt__field__flat.html#a742e5681a04b173a1f0ea63a381ad6bd">nmt_field_flat</a>
+</li>
 <li>matrix_M
 : <a class="el" href="structnmt__field.html#a111d5454c2cd405de1a7cbba122923ee">nmt_field</a>
 , <a class="el" href="structnmt__field__flat.html#a60bfe576ce32ebb6d5beedf5de42aebc">nmt_field_flat</a>

doc/build/html/genindex.html

@@ -283,9 +283,13 @@ <h2 id="G">G</h2>
 </li>
       </ul></li>
       <li><a href="pymaster.html#pymaster.bins.NmtBin.get_ell_list">get_ell_list() (pymaster.bins.NmtBin method)</a>
+</li>
+      <li><a href="pymaster.html#pymaster.bins.NmtBin.get_ell_max">get_ell_max() (pymaster.bins.NmtBin method)</a>
 </li>
   </ul></td>
   <td style="width: 33%; vertical-align: top;"><ul>
+      <li><a href="pymaster.html#pymaster.bins.NmtBin.get_ell_min">get_ell_min() (pymaster.bins.NmtBin method)</a>
+</li>
       <li><a href="pymaster.html#pymaster.field.NmtField.get_maps">get_maps() (pymaster.field.NmtField method)</a>
 
       <ul>

doc/build/html/globals.html

@@ -161,6 +161,9 @@ <h3><a class="anchor" id="index_f"></a>- f -</h3><ul>
 <li>fs_map_product()
 : <a class="el" href="utils_8h.html#a566f3ea3981a35851b962ac3a07d621f">utils.h</a>
 </li>
+<li>fs_mapcpy()
+: <a class="el" href="utils_8h.html#a74eb40afaa25f83297338cf105c566b1">utils.h</a>
+</li>
 <li>fs_read_flat_map()
 : <a class="el" href="utils_8h.html#a2d5988385bbeac1dadc5dc9f68a7c8c0">utils.h</a>
 </li>
@@ -433,10 +436,10 @@ <h3><a class="anchor" id="index_n"></a>- n -</h3><ul>
 : <a class="el" href="namaster_8h.html#a7e478d292e9d6ebd73ae461e35c7a565">namaster.h</a>
 </li>
 <li>nmt_field_alloc_sph()
-: <a class="el" href="namaster_8h.html#a122a5d2d05d31706c6d7764b7c4899b9">namaster.h</a>
+: <a class="el" href="namaster_8h.html#abb297e45531dfe7aa8ba20994ff15beb">namaster.h</a>
 </li>
 <li>nmt_field_flat_alloc()
-: <a class="el" href="namaster_8h.html#a7e23cb755d4674b4c0a5a6af2913bd75">namaster.h</a>
+: <a class="el" href="namaster_8h.html#a6a03c8af1ba66dfab840fc083772fbb1">namaster.h</a>
 </li>
 <li>nmt_field_flat_free()
 : <a class="el" href="namaster_8h.html#ae485829665fb97c262c561d1b7b44d73">namaster.h</a>

doc/build/html/globals_func.html

@@ -155,6 +155,9 @@ <h3><a class="anchor" id="index_f"></a>- f -</h3><ul>
 <li>fs_map_product()
 : <a class="el" href="utils_8h.html#a566f3ea3981a35851b962ac3a07d621f">utils.h</a>
 </li>
+<li>fs_mapcpy()
+: <a class="el" href="utils_8h.html#a74eb40afaa25f83297338cf105c566b1">utils.h</a>
+</li>
 <li>fs_read_flat_map()
 : <a class="el" href="utils_8h.html#a2d5988385bbeac1dadc5dc9f68a7c8c0">utils.h</a>
 </li>
@@ -424,10 +427,10 @@ <h3><a class="anchor" id="index_n"></a>- n -</h3><ul>
 : <a class="el" href="namaster_8h.html#a7e478d292e9d6ebd73ae461e35c7a565">namaster.h</a>
 </li>
 <li>nmt_field_alloc_sph()
-: <a class="el" href="namaster_8h.html#a122a5d2d05d31706c6d7764b7c4899b9">namaster.h</a>
+: <a class="el" href="namaster_8h.html#abb297e45531dfe7aa8ba20994ff15beb">namaster.h</a>
 </li>
 <li>nmt_field_flat_alloc()
-: <a class="el" href="namaster_8h.html#a7e23cb755d4674b4c0a5a6af2913bd75">namaster.h</a>
+: <a class="el" href="namaster_8h.html#a6a03c8af1ba66dfab840fc083772fbb1">namaster.h</a>
 </li>
 <li>nmt_field_flat_free()
 : <a class="el" href="namaster_8h.html#ae485829665fb97c262c561d1b7b44d73">namaster.h</a>

doc/build/html/namaster_8h.html

@@ -208,9 +208,9 @@
 <tr class="memitem:ae485829665fb97c262c561d1b7b44d73"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#ae485829665fb97c262c561d1b7b44d73">nmt_field_flat_free</a> (<a class="el" href="structnmt__field__flat.html">nmt_field_flat</a> *fl)</td></tr>
 <tr class="memdesc:ae485829665fb97c262c561d1b7b44d73"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field__flat.html" title="Flat-sky field. ">nmt_field_flat</a> destructor  <a href="#ae485829665fb97c262c561d1b7b44d73">More...</a><br /></td></tr>
 <tr class="separator:ae485829665fb97c262c561d1b7b44d73"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a7e23cb755d4674b4c0a5a6af2913bd75"><td class="memItemLeft" align="right" valign="top"><a class="el" href="structnmt__field__flat.html">nmt_field_flat</a> *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#a7e23cb755d4674b4c0a5a6af2913bd75">nmt_field_flat_alloc</a> (int nx, int ny, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> lx, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ly, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *mask, int spin, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **maps, int ntemp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ***temp, int nl_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *l_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *beam, int pure_e, int pure_b, double tol_pinv, int masked_input)</td></tr>
-<tr class="memdesc:a7e23cb755d4674b4c0a5a6af2913bd75"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field__flat.html" title="Flat-sky field. ">nmt_field_flat</a> constructor  <a href="#a7e23cb755d4674b4c0a5a6af2913bd75">More...</a><br /></td></tr>
-<tr class="separator:a7e23cb755d4674b4c0a5a6af2913bd75"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:a6a03c8af1ba66dfab840fc083772fbb1"><td class="memItemLeft" align="right" valign="top"><a class="el" href="structnmt__field__flat.html">nmt_field_flat</a> *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#a6a03c8af1ba66dfab840fc083772fbb1">nmt_field_flat_alloc</a> (int nx, int ny, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> lx, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ly, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *mask, int spin, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **maps, int ntemp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ***temp, int nl_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *l_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *beam, int pure_e, int pure_b, double tol_pinv, int masked_input, int is_lite, int mask_only)</td></tr>
+<tr class="memdesc:a6a03c8af1ba66dfab840fc083772fbb1"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field__flat.html" title="Flat-sky field. ">nmt_field_flat</a> constructor  <a href="#a6a03c8af1ba66dfab840fc083772fbb1">More...</a><br /></td></tr>
+<tr class="separator:a6a03c8af1ba66dfab840fc083772fbb1"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:aa22c9c632e39990c2cc79561b24efba8"><td class="memItemLeft" align="right" valign="top"><a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#aa22c9c632e39990c2cc79561b24efba8">nmt_synfast_flat</a> (int nx, int ny, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> lx, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ly, int nfields, int *spin_arr, int nl_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *l_beam, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **beam_fields, int nl_cell, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *l_cell, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **cell_fields, int seed)</td></tr>
 <tr class="memdesc:aa22c9c632e39990c2cc79561b24efba8"><td class="mdescLeft">&#160;</td><td class="mdescRight">Gaussian realizations of flat-sky fields.  <a href="#aa22c9c632e39990c2cc79561b24efba8">More...</a><br /></td></tr>
 <tr class="separator:aa22c9c632e39990c2cc79561b24efba8"><td class="memSeparator" colspan="2">&#160;</td></tr>
@@ -232,9 +232,9 @@
 <tr class="memitem:a41926995999fb081de99f07473559af8"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#a41926995999fb081de99f07473559af8">nmt_field_free</a> (<a class="el" href="structnmt__field.html">nmt_field</a> *fl)</td></tr>
 <tr class="memdesc:a41926995999fb081de99f07473559af8"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field.html" title="Full-sky field. ">nmt_field</a> destructor.  <a href="#a41926995999fb081de99f07473559af8">More...</a><br /></td></tr>
 <tr class="separator:a41926995999fb081de99f07473559af8"><td class="memSeparator" colspan="2">&#160;</td></tr>
-<tr class="memitem:a122a5d2d05d31706c6d7764b7c4899b9"><td class="memItemLeft" align="right" valign="top"><a class="el" href="structnmt__field.html">nmt_field</a> *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#a122a5d2d05d31706c6d7764b7c4899b9">nmt_field_alloc_sph</a> (<a class="el" href="structnmt__curvedsky__info.html">nmt_curvedsky_info</a> *cs, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *mask, int spin, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **maps, int ntemp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ***temp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *beam, int pure_e, int pure_b, int n_iter_mask_purify, double tol_pinv, int niter, int masked_input)</td></tr>
-<tr class="memdesc:a122a5d2d05d31706c6d7764b7c4899b9"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field.html" title="Full-sky field. ">nmt_field</a> constructor  <a href="#a122a5d2d05d31706c6d7764b7c4899b9">More...</a><br /></td></tr>
-<tr class="separator:a122a5d2d05d31706c6d7764b7c4899b9"><td class="memSeparator" colspan="2">&#160;</td></tr>
+<tr class="memitem:abb297e45531dfe7aa8ba20994ff15beb"><td class="memItemLeft" align="right" valign="top"><a class="el" href="structnmt__field.html">nmt_field</a> *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#abb297e45531dfe7aa8ba20994ff15beb">nmt_field_alloc_sph</a> (<a class="el" href="structnmt__curvedsky__info.html">nmt_curvedsky_info</a> *cs, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *mask, int spin, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> **maps, int ntemp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> ***temp, <a class="el" href="namaster_8h.html#a420a9f7802706d142bec6741123f57a9">flouble</a> *beam, int pure_e, int pure_b, int n_iter_mask_purify, double tol_pinv, int niter, int masked_input, int is_lite, int mask_only)</td></tr>
+<tr class="memdesc:abb297e45531dfe7aa8ba20994ff15beb"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field.html" title="Full-sky field. ">nmt_field</a> constructor  <a href="#abb297e45531dfe7aa8ba20994ff15beb">More...</a><br /></td></tr>
+<tr class="separator:abb297e45531dfe7aa8ba20994ff15beb"><td class="memSeparator" colspan="2">&#160;</td></tr>
 <tr class="memitem:a30061f25d63662b222e814606c7a8780"><td class="memItemLeft" align="right" valign="top"><a class="el" href="structnmt__field.html">nmt_field</a> *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="namaster_8h.html#a30061f25d63662b222e814606c7a8780">nmt_field_read</a> (int is_healpix, char *fname_mask, char *fname_maps, char *fname_temp, char *fname_beam, int spin, int pure_e, int pure_b, int n_iter_mask_purify, double tol_pinv, int niter)</td></tr>
 <tr class="memdesc:a30061f25d63662b222e814606c7a8780"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="structnmt__field.html" title="Full-sky field. ">nmt_field</a> constructor from file.  <a href="#a30061f25d63662b222e814606c7a8780">More...</a><br /></td></tr>
 <tr class="separator:a30061f25d63662b222e814606c7a8780"><td class="memSeparator" colspan="2">&#160;</td></tr>
@@ -2813,7 +2813,7 @@ <h2 class="groupheader">Function Documentation</h2>
 
 </div>
 </div>
-<a class="anchor" id="a122a5d2d05d31706c6d7764b7c4899b9"></a>
+<a class="anchor" id="abb297e45531dfe7aa8ba20994ff15beb"></a>
 <div class="memitem">
 <div class="memproto">
       <table class="memname">
@@ -2893,7 +2893,19 @@ <h2 class="groupheader">Function Documentation</h2>
           <td class="paramkey"></td>
           <td></td>
           <td class="paramtype">int&#160;</td>
-          <td class="paramname"><em>masked_input</em>&#160;</td>
+          <td class="paramname"><em>masked_input</em>, </td>
+        </tr>
+        <tr>
+          <td class="paramkey"></td>
+          <td></td>
+          <td class="paramtype">int&#160;</td>
+          <td class="paramname"><em>is_lite</em>, </td>
+        </tr>
+        <tr>
+          <td class="paramkey"></td>
+          <td></td>
+          <td class="paramtype">int&#160;</td>
+          <td class="paramname"><em>mask_only</em>&#160;</td>
         </tr>
         <tr>
           <td></td>
@@ -2919,13 +2931,15 @@ <h2 class="groupheader">Function Documentation</h2>
     <tr><td class="paramname">tol_pinv</td><td>Contaminant deprojection requires the inversion of the template covariance matrix. This could be ill-defined if some templates are linearly related. In this case we use a pseudo-inverse that accounts for this possibility in a consistent way. Effectively this is a singular-value decomposition. All eigenvalues that are smaller than <code>tol_pinv</code> the largest eigenvalue will be discarded. </td></tr>
     <tr><td class="paramname">niter</td><td>number of iterations when computing alms (for all transforms other than the mask's). </td></tr>
     <tr><td class="paramname">masked_input</td><td>if not 0, input maps and templates have already been masked. This is not advisable if using purification. </td></tr>
+    <tr><td class="paramname">is_lite</td><td>if not 0, only the map alms and the mask will be stored. You can then use this field to compute the standard pseudo-C_ell with deprojection and purification, but you won't be able to compute the deprojection bias or examine any maps. </td></tr>
+    <tr><td class="paramname">mask_only</td><td>if not 0, this field will only store a mask and a beam. You will be able to use it to compute the PCL and covariance mode coupling matrices, but that's it (no actual power spectra, deprojection biases etc.). </td></tr>
   </table>
   </dd>
 </dl>
 
 </div>
 </div>
-<a class="anchor" id="a7e23cb755d4674b4c0a5a6af2913bd75"></a>
+<a class="anchor" id="a6a03c8af1ba66dfab840fc083772fbb1"></a>
 <div class="memitem">
 <div class="memproto">
       <table class="memname">
@@ -3023,7 +3037,19 @@ <h2 class="groupheader">Function Documentation</h2>
           <td class="paramkey"></td>
           <td></td>
           <td class="paramtype">int&#160;</td>
-          <td class="paramname"><em>masked_input</em>&#160;</td>
+          <td class="paramname"><em>masked_input</em>, </td>
+        </tr>
+        <tr>
+          <td class="paramkey"></td>
+          <td></td>
+          <td class="paramtype">int&#160;</td>
+          <td class="paramname"><em>is_lite</em>, </td>
+        </tr>
+        <tr>
+          <td class="paramkey"></td>
+          <td></td>
+          <td class="paramtype">int&#160;</td>
+          <td class="paramname"><em>mask_only</em>&#160;</td>
         </tr>
         <tr>
           <td></td>
@@ -3052,6 +3078,8 @@ <h2 class="groupheader">Function Documentation</h2>
     <tr><td class="paramname">pure_b</td><td>Set to &gt;0 if you want purified B-modes. </td></tr>
     <tr><td class="paramname">tol_pinv</td><td>Contaminant deprojection requires the inversion of the template covariance matrix. This could be ill-defined if some templates are linearly related. In this case we use a pseudo-inverse that accounts for this possibility in a consistent way. Effectively this is a singular-value decomposition. All eigenvalues that are smaller than <code>tol_pinv</code> the largest eigenvalue will be discarded. </td></tr>
     <tr><td class="paramname">masked_input</td><td>if not 0, input maps and templates have already been masked. This is not advisable if using purification. </td></tr>
+    <tr><td class="paramname">is_lite</td><td>if not 0, only the map alms and the mask will be stored. You can then use this field to compute the standard pseudo-C_ell with deprojection and purification, but you won't be able to compute the deprojection bias or examine any maps. </td></tr>
+    <tr><td class="paramname">mask_only</td><td>if not 0, this field will only store a mask and a beam. You will be able to use it to compute the PCL and covariance mode coupling matrices, but that's it (no actual power spectra, deprojection biases etc.). </td></tr>
   </table>
   </dd>
 </dl>