Allow different funcs for posterior

Author: achiefa

validphys2/src/validphys/config.py

@@ -1769,8 +1769,11 @@ def produce_total_phi_data(self, fitthcovmat):
             return validphys.results.total_phi_data_from_experiments
         return validphys.results.dataset_inputs_phi_data
 
+    def produce_pc_func_type(self, theorycovmatconfig=None):
+        return theorycovmatconfig.get('func_type', 'step')
+
     # @configparser.explicit_node
-    def produce_power_corr_dict(self, pc_parameters=None):
+    def produce_power_corr_dict(self, pc_parameters=None, pc_func_type=None):
         """The parameters for the power corrections are given as a list.
         This function converts this list into a dictionary with the keys
         being the names of the types of power corrections (e.g. `H2p`, `H2d`,...).
@@ -1782,7 +1785,12 @@ def produce_power_corr_dict(self, pc_parameters=None):
         # Loop over the parameterization for the power corrections in the runcard
         for par in pc_parameters:
             # Check that the length of shifts is one less than the length of nodes.
-            if len(par['yshift']) != len(par['nodes']) - 1:
+            if (len(par['yshift']) != len(par['nodes']) - 1) and pc_func_type != 'cubic':
+                raise ValueError(
+                    f"The length of nodes does not match that of the list in {par['ht']}."
+                    f"Check the runcard. Got {len(par['yshift'])} != {len(par['nodes'])}"
+                )
+            elif (len(par['yshift']) != len(par['nodes'])) and pc_func_type == 'cubic':
                 raise ValueError(
                     f"The length of nodes does not match that of the list in {par['ht']}."
                     f"Check the runcard. Got {len(par['yshift'])} != {len(par['nodes'])}"
@@ -1791,6 +1799,7 @@ def produce_power_corr_dict(self, pc_parameters=None):
             # Store parameters for each power correction
             pc_parameters_by_type[par['ht']] = {'yshift': par['yshift'], 'nodes': par['nodes']}
 
+        pc_parameters_by_type['func_type'] = pc_func_type
         return pc_parameters_by_type
 
     @configparser.explicit_node