classy 3.2.3.2__tar.gz → 3.2.5.0__tar.gz

{classy-3.2.3.2/classy.egg-info → classy-3.2.5.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: classy
-Version: 3.2.3.2
+Version: 3.2.5.0
 Summary: The class code from https://lesgourg.github.io/class_public/class.html
 Home-page: http://www.class-code.net
 Author: Julien Lesgourgues & Thomas Tram

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/.github/workflows/test_nightly.yml

@@ -41,7 +41,7 @@ jobs:
     if: needs.iftest.outputs.should_run > 0
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: Create or update virtual Python environment
@@ -60,7 +60,7 @@ jobs:
         OMP_NUM_THREADS=32 COMPARE_OUTPUT_GAUGE=1 TEST_LEVEL=3 nosetests -v -a test_scenario test_class.py --nologcapture --nocapture
     - name: Upload plots 📤
       if: success() || failure()
-      uses: actions/upload-artifact@v2
+      uses: actions/upload-artifact@v4
       with:
         name: ComparePlots
         path: main_class/python/faulty_figs
@@ -72,23 +72,29 @@ jobs:
     timeout-minutes: 1440
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: Create or update virtual Python environment
       run: |
-        rm -f venv/bin/python
+        rm -rf venv
         virtualenv venv
         source venv/bin/activate
         pip install --upgrade pip
         pip install numpy scipy pandas matplotlib cython nose parameterized
+        deactivate
     - name: make
-      run: cd main_class && make OMPFLAG="" -j
+      run: |
+        source venv/bin/activate
+        cd main_class
+        make -j
+        deactivate
     - name: Generate input files
       run: |
         source venv/bin/activate
         cd main_class/python
-        TEST_LEVEL=3 CLASS_VERBOSE=1 nosetests -a dump_ini_files test_class.py
+        TEST_LEVEL=3 VALGRIND_MODE=1 CLASS_VERBOSE=1 python "$(which nosetests)" -a dump_ini_files test_class.py
+        deactivate
     - name: Valgrind 🤖
       run: |
         cd main_class/python/faulty_figs
@@ -98,11 +104,11 @@ jobs:
         find . -name "*.ini" -type f -print0 | xargs -0 -I {} -n 1 -P 32 bash -c ': \
           && printf "\nk_step_sub = 10.0\nk_step_super = 0.4\nk_per_decade_for_pk = 2\nk_per_decade_for_bao = 5\n" >> {}'
         find . -name "*.ini" -type f -print0 | xargs -0 -I {} -n 1 -P 32 bash -c ': \
-          && (valgrind --track-origins=yes --show-leak-kinds=all --leak-check=full --error-exitcode=3 ./class {} &> {}.out && echo {} "...ok") \
+          && (OMP_NUM_THREADS=1 valgrind --track-origins=yes --show-leak-kinds=all --leak-check=full --error-exitcode=3 ./class {} &> {}.out && echo {} "...ok") \
           || (echo {} "...fail" && cp {}* valgrind_output && exit 3)'
     - name: Upload errors 📤
       if: success() || failure()
-      uses: actions/upload-artifact@v2
+      uses: actions/upload-artifact@v4
       with:
         name: ValgrindOutputWideShallow
         path: main_class/python/faulty_figs/valgrind_output
@@ -114,23 +120,29 @@ jobs:
     timeout-minutes: 1440
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: Create or update virtual Python environment
       run: |
-        rm -f venv/bin/python
+        rm -rf venv
         virtualenv venv
         source venv/bin/activate
         pip install --upgrade pip
         pip install numpy scipy pandas matplotlib cython nose parameterized
+        deactivate
     - name: make
-      run: cd main_class && make OMPFLAG="" -j
+      run: |
+        source venv/bin/activate
+        cd main_class
+        make -j
+        deactivate
     - name: Generate input files
       run: |
         source venv/bin/activate
         cd main_class/python
         TEST_LEVEL=1 CLASS_VERBOSE=1 nosetests -a dump_ini_files test_class.py
+        deactivate
     - name: Valgrind 🤖
       run: |
         cd main_class/python/faulty_figs
@@ -138,11 +150,11 @@ jobs:
         rm -rf valgrind_output && mkdir valgrind_output
         cp ../../class .
         find . -name "*.ini" -type f -print0 | xargs -0 -I {} -n 1 -P 32 bash -c ': \
-          && (valgrind --track-origins=yes --show-leak-kinds=all --leak-check=full --error-exitcode=3 ./class {} &> {}.out && echo {} "...ok") \
+          && (OMP_NUM_THREADS=1 valgrind --track-origins=yes --show-leak-kinds=all --leak-check=full --error-exitcode=3 ./class {} &> {}.out && echo {} "...ok") \
           || (echo {} "...fail" && cp {}* valgrind_output && exit 3)'
     - name: Upload errors 📤
       if: success() || failure()
-      uses: actions/upload-artifact@v2
+      uses: actions/upload-artifact@v4
       with:
         name: ValgrindOutputNarrowDeep
         path: main_class/python/faulty_figs/valgrind_output

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/.github/workflows/test_on_pull_request.yml

@@ -8,11 +8,11 @@ jobs:
 
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: Checkout reference 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         ref: devel
         path: ref_class
@@ -23,8 +23,9 @@ jobs:
         source venv/bin/activate
         pip install --upgrade pip
         pip install numpy scipy pandas matplotlib cython nose parameterized
+        deactivate
     - name: make classy
-      run: source venv/bin/activate && cd main_class && make -j
+      run: source venv/bin/activate && cd main_class && make -j && deactivate
     - name: make reference classy
       run: |
         source venv/bin/activate
@@ -32,14 +33,16 @@ jobs:
         sed -i.bak "s/'classy'/'classyref'/g" python/setup.py
         sed -i.bak 's/"classy"/"classyref"/g' python/setup.py
         make -j
+        deactivate
     - name: Testing 🤖
       run: |
         source venv/bin/activate
         cd main_class/python
-        OMP_NUM_THREADS=16 COMPARE_OUTPUT_REF=1 TEST_LEVEL=2 nosetests -v -a test_scenario test_class.py --nologcapture --nocapture
+        OMP_NUM_THREADS=16 COMPARE_OUTPUT_REF=1 TEST_LEVEL=3 nosetests -v -a test_scenario test_class.py --nologcapture --nocapture
+        deactivate
     - name: Upload plots 📤
       if: success() || failure()
-      uses: actions/upload-artifact@v2
+      uses: actions/upload-artifact@v4
       with:
           name: BranchVsMasterPlots
           path: main_class/python/faulty_figs

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/.github/workflows/test_on_push.yml

@@ -8,7 +8,7 @@ jobs:
 
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: make
@@ -21,7 +21,7 @@ jobs:
 
     steps:
     - name: Checkout 🛎
-      uses: actions/checkout@v2
+      uses: actions/checkout@v4
       with:
         path: main_class
     - name: Create or update virtual Python environment
@@ -31,10 +31,12 @@ jobs:
         source venv/bin/activate
         pip install --upgrade pip
         pip install numpy scipy pandas matplotlib cython nose parameterized
+        deactivate
     - name: make classy
-      run: source venv/bin/activate && cd main_class && make -j
+      run: source venv/bin/activate && cd main_class && make -j && deactivate
     - name: Testing 🤖
       run: |
         source venv/bin/activate
         cd main_class/python
-        OMP_NUM_THREADS=16 TEST_LEVEL=1 nosetests -v -a test_scenario test_class.py --nologcapture --nocapture
+        OMP_NUM_THREADS=16 TEST_LEVEL=1 python "$(which nosetests)" -v -a test_scenario test_class.py --nologcapture --nocapture
+        deactivate

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/Makefile

@@ -195,7 +195,14 @@ tar: $(C_ALL) $(C_TEST) $(H_ALL) $(PRE_ALL) $(INI_ALL) $(MISC_FILES) $(HYREC) $(
 	tar czvf class.tar.gz $(C_ALL) $(H_ALL) $(PRE_ALL) $(INI_ALL) $(MISC_FILES) $(HYREC) $(PYTHON_FILES)
 
 classy: libclass.a python/classy.pyx python/cclassy.pxd
-	cd python; export CC=$(CC); $(PYTHON) setup.py install || $(PYTHON) setup.py install --user
+	cd python; export CC=$(CC); output=$$($(PYTHON) -m pip install . 2>&1); \
+    echo "$$output"; \
+    if echo "$$output" | grep -q "ERROR: Cannot uninstall"; then \
+        site_packages=$$($(PYTHON) -c "import distutils.sysconfig; print(distutils.sysconfig.get_python_lib())" || $(PYTHON) -c "import site; print(site.getsitepackages()[0])") && \
+		echo "Cleaning up previous installation in: $$site_packages" && \
+        rm -rf $$site_packages/classy* && \
+        $(PYTHON) -m pip install .; \
+    fi
 
 clean: .base
 	rm -rf $(WRKDIR);

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/explanatory.ini

@@ -151,6 +151,9 @@ gauge = synchronous
 #      (default: set to 'no')
 #nbody_gauge_transfer_functions = yes
 
+# 4.c) Do you want the source functions for total non-relativistic matter, delta_m and theta_m, and baryon+cdm, delta_cb and theta_cb, to be outputed in the current gauge (the one selected in 4.a), instead of being automatically expressed as a gauge-invariant (GI) quantity, likeL: delta_m^GI=delta_m + 3 a H theta_m/k2, theta_m^GI=theta_m+alpha*k2 (default: no, that is, convert to GI)
+#matter_source_in_current_gauge= no
+
 # 5) Hubble parameter : either 'H0' in km/s/Mpc or 'h' (or 'theta_s_100'), where
 #    the latter is the peak scale parameter defined exactly as 100(ds_dec/da_dec)
 #    with a decoupling time given by maximum of visibility function (quite different

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/external/RealSpaceInterface/Calc2D/CalculationClass.py

@@ -71,7 +71,7 @@ class Calculation(object):
         Valuenew = dict()
         FValue_abs = np.abs(self.FValue)
         _min, _max = FValue_abs.min(), FValue_abs.max()
-        dimensions = (self.endshape[0] / 2, self.endshape[1])
+        dimensions = (self.endshape[0] // 2, self.endshape[1])
         for quantity, FT in FValuenew.items():
             FT_abs = np.abs(FT)
             FT_normalized = cv2.resize(FT_abs, dimensions).ravel()
@@ -98,7 +98,7 @@ class Calculation(object):
         return Value.ravel(), cv2.resize(
             (np.abs(self.FValue) - np.abs(self.FValue).min()) /
             (np.abs(self.FValue).max() - np.abs(self.FValue).min()),
-            (self.endshape[0] / 2, self.endshape[1])).ravel(), (minimum,
+            (self.endshape[0] // 2, self.endshape[1])).ravel(), (minimum,
                                                                 maximum)
 
     def getTransferData(self, redshiftindex):

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/external/RealSpaceInterface/Calc2D/Database.py

@@ -14,22 +14,22 @@ class Database:
         self.db_path = os.path.join(directory, db_file)
         if not os.path.exists(self.db_path):
             logging.info("No database found; Creating one at {}.".format(self.db_path))
-            with open(self.db_path, "w") as f:
+            with open(self.db_path, "wb") as f:
                 pickle.dump(dict(), f)
 
         self.db = self.__read_database()
 
     def __read_database(self):
-        with open(self.db_path) as f:
-            return pickle.load(f)
+        with open(self.db_path, "rb") as f:
+            return pickle.load(f, encoding='latin1')
 
     def __write_database(self):
-        with open(self.db_path, "w") as f:
-            pickle.dump(self.db, f)
+        with open(self.db_path, "wb") as f:
+            pickle.dump(self.db, f, protocol=pickle.HIGHEST_PROTOCOL)
 
     def __create_file(self, data):
         filename = str(uuid.uuid4())
-        with open(os.path.join(self.directory, filename), "w") as f:
+        with open(os.path.join(self.directory, filename), "wb") as f:
             pickle.dump(data, f)
         return filename
 
@@ -40,8 +40,8 @@ class Database:
         frozen_key = self.__get_frozen_key(key)
         if frozen_key in self.db:
             filename = self.db[frozen_key]
-            with open(os.path.join(self.directory, filename)) as f:
-                return pickle.load(f)
+            with open(os.path.join(self.directory, filename), "rb") as f:
+                return pickle.load(f, encoding='latin1')
         else:
             raise KeyError("No data for key: {}".format(key))
 
@@ -55,4 +55,4 @@ class Database:
         Return whether `self` contains a record
         for the given `key`.
         """
-        return self.__get_frozen_key(key) in self.db
+        return self.__get_frozen_key(key) in self.db

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/external/RealSpaceInterface/Calc2D/TransferFunction.py

@@ -42,6 +42,7 @@ def ComputeTransferFunctionList(cosmologicalParameters, redshift, kperdecade=200
     class_settings.update({
         "output": "mTk",
         "gauge": "newtonian",
+        "matter_source_in_current_gauge": "yes",
         "evolver": "1",
         "P_k_max_h/Mpc": P_k_max,
         "k_per_decade_for_pk": kperdecade,

classy-3.2.5.0/class_public/external/RealSpaceInterface/README

@@ -0,0 +1,23 @@
+CLASS real space interface
+07.2015: started by Max Beutelspacher
+09.2018: improved by Georgios Samaras and released in class v2.7.0
+10.2025: for v3.2.5: made compatible with python 3, credits Timothy Morton and Mowen Zhao
+
+
+For installation of python packages, run
+
+    pip install -r requirements.txt
+
+To see the real space interface, go to this directory (external/RealSpaceInterface/) and launch the application with
+
+    python tornadoserver.py
+
+Then in any browser open the URL
+
+    http://localhost:7777
+
+------------------------------------------------------------
+
+Cache files are located in cache/, so to clear the cache, run
+
+    rm cache/*

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/external/RealSpaceInterface/tornadoserver.py

@@ -143,8 +143,8 @@ class DataConnection(tornado.websocket.WebSocketHandler):
         self.write_message(json.dumps({'type': 'extrema', 'extrema': extrema}))
         progress = float(redindex) / len(self.calc.redshift)
 
-        real = {quantity: base64.b64encode(data.astype(np.float32)) for quantity, data in Valuenew.iteritems()}
-        transfer = {quantity: base64.b64encode(data.astype(np.float32)) for quantity, data in TransferData.iteritems()}
+        real = {quantity: base64.b64encode(data.astype(np.float32)).decode() for quantity, data in Valuenew.items()}
+        transfer = {quantity: base64.b64encode(data.astype(np.float32)).decode() for quantity, data in TransferData.items()}
         self.write_message(
             json.dumps({
                 'type': 'data',
@@ -166,9 +166,9 @@ class DataConnection(tornado.websocket.WebSocketHandler):
         extremastring = json.dumps({'type': 'extrema', 'extrema': extrema})
         datastring = json.dumps({
             'type': 'data',
-            'real': base64.b64encode(Value.astype(np.float32)),
+            'real': base64.b64encode(Value.astype(np.float32)).decode(),
             'fourier': [],
-            'transfer': base64.b64encode(TransferData.astype(np.float32)),
+            'transfer': base64.b64encode(TransferData.astype(np.float32)).decode(),
             'k': krange.tolist()
             })
         self.write_message(extremastring)
@@ -212,6 +212,7 @@ class DataConnection(tornado.websocket.WebSocketHandler):
 
             z_of_decoupling = self.calc.z_dec
             frame_of_decoupling = np.argmin(np.abs(z_of_decoupling - self.calc.redshift))
+            frame_of_decoupling = int(frame_of_decoupling) # this is np.int64 object, which json doesn't like
             if self.calc.redshift[frame_of_decoupling] > z_of_decoupling:
                 frame_of_decoupling -= 1
             messages.append({

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/background.h

@@ -10,10 +10,6 @@
 #include "dei_rkck.h"
 #include "parser.h"
 
-/** list of possible types of spatial curvature */
-
-enum spatial_curvature {flat,open,closed};
-
 /** list of possible parametrisations of the DE equation of state */
 
 enum equation_of_state {CLP,EDE};
@@ -328,6 +324,7 @@ struct background
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
   //@}
 };
 

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/common.h

@@ -15,7 +15,7 @@
 #ifndef __COMMON__
 #define __COMMON__
 
-#define _VERSION_ "v3.2.3"
+#define _VERSION_ "v3.2.5"
 
 /* @cond INCLUDE_WITH_DOXYGEN */
 

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/distortions.h

@@ -181,6 +181,8 @@ struct distortions
 
   ErrorMsg error_message;    /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 
 };

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/fourier.h

@@ -192,6 +192,8 @@ struct fourier {
 
   ErrorMsg error_message; 	/**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 };
 
@@ -354,6 +356,7 @@ extern "C" {
 
   int fourier_get_k_list(
                          struct precision *ppr,
+                         struct primordial *ppm,
                          struct perturbations * ppt,
                          struct fourier * pfo
                          );

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/harmonic.h

@@ -134,6 +134,8 @@ struct harmonic {
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 };
 

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/input.h

@@ -203,10 +203,11 @@
  * For shooting method: definition of the possible targets
  */
 
-enum target_names {theta_s, theta_s_100, Omega_dcdmdr, omega_dcdmdr, Omega_scf, Omega_ini_dcdm, omega_ini_dcdm, sigma8, S8};
+enum target_names {theta_s, theta_s_100, Neff, Omega_dcdmdr, omega_dcdmdr, Omega_scf, Omega_ini_dcdm, omega_ini_dcdm, sigma8, S8};
 /* Important: Keep this number equal to the number of target_names (except sigma8, S8), and keep sigma8, S8 at the very end */
-#define _NUM_TARGETS_ 7
-/* Important: add one for each new target_names */
+#define _NUM_TARGETS_ 8
+
+/* Until which class stage is being computed? Pretty much fixed list, don't change. */
 enum computation_stage {cs_background, cs_thermodynamics, cs_perturbations, cs_primordial, cs_nonlinear, cs_transfer, cs_spectra};
 
 /**

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/lensing.h

@@ -89,6 +89,8 @@ struct lensing {
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 };
 

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/parser.h

@@ -100,6 +100,19 @@ extern "C" {
                            int N_options,
                            int* valid);
 
+  int parser_extend(struct file_content * pfc,
+                    int N_extend,
+                    ErrorMsg errmsg);
+
+  int parser_copy(struct file_content * pfc_source,
+                  struct file_content * pfc_destination,
+                  int index_start,
+                  int index_end);
+
+  int parser_init_from_pfc(struct file_content * pfc_source,
+                           struct file_content * pfc_destination,
+                           ErrorMsg errmsg);
+
 #ifdef __cplusplus
 }
 #endif

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/perturbations.h

@@ -198,6 +198,8 @@ struct perturbations
 
   enum possible_gauges gauge; /**< gauge in which to perform this calculation */
 
+  short has_matter_source_in_current_gauge; /**< whether to keep matter and baryon+CDM sources in current gauge, instead of automatic conversion to gauge-invariant variables */
+
   //@}
 
   /** @name - indices running on modes (scalar, vector, tensor) */
@@ -440,6 +442,8 @@ struct perturbations
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 
 };

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/primordial.h

@@ -300,6 +300,8 @@ struct primordial {
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
 };
 
 struct primordial_inflation_parameters_and_workspace {

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/thermodynamics.h

@@ -328,6 +328,7 @@ struct thermodynamics
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
   //@}
 
 };

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/include/transfer.h

@@ -5,8 +5,6 @@
 
 #include "fourier.h"
 #include "hyperspherical.h"
-#include <sys/shm.h>
-#include <sys/stat.h>
 #include "errno.h"
 
 /* macro: test if index_tt is in the range between index and index+num, while the flag is true */
@@ -207,6 +205,8 @@ struct transfer {
 
   ErrorMsg error_message; /**< zone for writing error messages */
 
+  short is_allocated; /**< flag is set to true if allocated */
+
   //@}
 };
 
@@ -675,7 +675,7 @@ extern "C" {
   int transfer_workspace_init(
                               struct transfer * ptr,
                               struct precision * ppr,
-                              struct transfer_workspace **ptw,
+                              struct transfer_workspace *ptw,
                               int perturbations_tau_size,
                               int tau_size_max,
                               double K,

{classy-3.2.3.2 → classy-3.2.5.0}/class_public/python/cclassy.pxd

@@ -53,6 +53,7 @@ cdef extern from "class.h":
         ErrorMsg error_message
 
     cdef struct background:
+        short is_allocated
         ErrorMsg error_message
         int bg_size
         int index_bg_ang_distance
@@ -106,6 +107,7 @@ cdef extern from "class.h":
         double tau_eq
 
     cdef struct thermodynamics:
+        short is_allocated
         ErrorMsg error_message
         int th_size
         int index_th_xe
@@ -141,6 +143,7 @@ cdef extern from "class.h":
         int tt_size
 
     cdef struct perturbations:
+        short is_allocated
         ErrorMsg error_message
         short has_scalars
         short has_vectors
@@ -277,9 +280,11 @@ cdef extern from "class.h":
         int ln_tau_size
 
     cdef struct transfer:
+        short is_allocated
         ErrorMsg error_message
 
     cdef struct primordial:
+        short is_allocated
         ErrorMsg error_message
         double k_pivot
         double A_s
@@ -311,6 +316,7 @@ cdef extern from "class.h":
         int lnk_size
 
     cdef struct harmonic:
+        short is_allocated
         ErrorMsg error_message
         int has_tt
         int has_te
@@ -350,6 +356,7 @@ cdef extern from "class.h":
         ErrorMsg error_message
 
     cdef struct distortions:
+        short is_allocated
         double * sd_parameter_table
         int index_type_g
         int index_type_mu
@@ -365,6 +372,7 @@ cdef extern from "class.h":
         ErrorMsg error_message
 
     cdef struct lensing:
+        short is_allocated
         int has_tt
         int has_ee
         int has_te
@@ -395,6 +403,7 @@ cdef extern from "class.h":
         ErrorMsg error_message
 
     cdef struct fourier:
+        short is_allocated
         short has_pk_matter
         int method
         int ic_size