OpenFisca-France 161.0.3__py3-none-any.whl → 169.14.7__py3-none-any.whl

openfisca_france/model/prelevements_obligatoires/impot_revenu/ir.py

@@ -416,10 +416,10 @@ class revenu_assimile_salaire_apres_abattements(Variable):
         revenu_assimile_salaire = individu('revenu_assimile_salaire', period)
         chomeur_longue_duree = individu('chomeur_longue_duree', period)
         frais_reels = individu('frais_reels', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
+        parameters_deductions = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
 
-        abattement_minimum = where(chomeur_longue_duree, P.abatpro.min2, P.abatpro.min)
-        abatfor = round_(min_(max_(P.taux_salaires_pensions * revenu_assimile_salaire, abattement_minimum), P.abatpro.max))
+        abattement_minimum = where(chomeur_longue_duree, parameters_deductions.abatpro.min2, parameters_deductions.abatpro.min)
+        abatfor = round_(min_(max_(parameters_deductions.abatpro.taux * revenu_assimile_salaire, abattement_minimum), parameters_deductions.abatpro.max))
         return (
             (frais_reels > abatfor)
             * (revenu_assimile_salaire - frais_reels)
@@ -430,9 +430,9 @@ class revenu_assimile_salaire_apres_abattements(Variable):
     def formula_2018_01_01(individu, period, parameters):
         revenu_assimile_salaire = individu('revenu_assimile_salaire', period)
         frais_reels = individu('frais_reels', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
+        parameters_deductions = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
 
-        abatfor = round_(min_(max_(P.taux_salaires_pensions * revenu_assimile_salaire, P.abatpro.min), P.abatpro.max))
+        abatfor = round_(min_(max_(parameters_deductions.abatpro.taux * revenu_assimile_salaire, parameters_deductions.abatpro.min), parameters_deductions.abatpro.max))
         return (
             (frais_reels > abatfor)
             * (revenu_assimile_salaire - frais_reels)
@@ -464,17 +464,17 @@ class revenu_assimile_pension_apres_abattements(Variable):
 
     def formula(individu, period, parameters):
         revenu_assimile_pension = individu('revenu_assimile_pension', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
+        parameters_deductions = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
 
         #    TODO: problème car les pensions sont majorées au niveau du foyer
     #    d11 = ( AS + BS + CS + DS + ES +
     #            AO + BO + CO + DO + EO )
     #    penv2 = (d11-f11> abatpen.max)*(penv + (d11-f11-abatpen.max)) + (d11-f11<= abatpen.max)*penv
     #    Plus d'abatement de 20% en 2006
-        return max_(0, revenu_assimile_pension - round_(max_(P.taux_salaires_pensions * revenu_assimile_pension, P.abatpen.min)))
+        return max_(0, revenu_assimile_pension - round_(max_(parameters_deductions.abatpen.taux * revenu_assimile_pension, parameters_deductions.abatpen.min)))
 
 
-#    return max_(0, revenu_assimile_pension - min_(round_(max_(P.taux_salaires_pensions*revenu_assimile_pension , P.abatpen.min)), P.abatpen.max))  le max se met au niveau du foyer
+#    return max_(0, revenu_assimile_pension - min_(round_(max_(parameters_deductions.abatpen.taux*revenu_assimile_pension , parameters_deductions.abatpen.min)), parameters_deductions.abatpen.max))  le max se met au niveau du foyer
 
 class indu_plaf_abat_pen(Variable):
     value_type = float
@@ -485,13 +485,13 @@ class indu_plaf_abat_pen(Variable):
     def formula(foyer_fiscal, period, parameters):
         rev_pen_i = foyer_fiscal.members('revenu_assimile_pension', period)
         pen_net_i = foyer_fiscal.members('revenu_assimile_pension_apres_abattements', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
+        parameters_deductions = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
 
         revenu_assimile_pension_apres_abattements = foyer_fiscal.sum(pen_net_i)
         revenu_assimile_pension = foyer_fiscal.sum(rev_pen_i)
 
         abat = revenu_assimile_pension - revenu_assimile_pension_apres_abattements
-        return abat - min_(abat, P.abatpen.max)
+        return abat - min_(abat, parameters_deductions.abatpen.max)
 
 
 class abattement_salaires_pensions(Variable):
@@ -504,9 +504,9 @@ class abattement_salaires_pensions(Variable):
     def formula(individu, period, parameters):
         revenu_assimile_salaire_apres_abattements = individu('revenu_assimile_salaire_apres_abattements', period)
         revenu_assimile_pension_apres_abattements = individu('revenu_assimile_pension_apres_abattements', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
+        parameters_deductions = parameters(period).impot_revenu.calcul_revenus_imposables.deductions
 
-        return min_(P.abat_supp.taux * max_(revenu_assimile_salaire_apres_abattements + revenu_assimile_pension_apres_abattements, 0), P.abat_supp.max)
+        return min_(parameters_deductions.abat_supp.taux * max_(revenu_assimile_salaire_apres_abattements + revenu_assimile_pension_apres_abattements, 0), parameters_deductions.abat_supp.max)
 
 
 class rente_viagere_titre_onereux(Variable):
@@ -691,6 +691,9 @@ class revenu_categoriel_capital(Variable):
         Revenus des valeurs et capitaux mobiliers
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
         f2ch = foyer_fiscal('f2ch', period)
         f2dc = foyer_fiscal('f2dc', period)
@@ -700,32 +703,34 @@ class revenu_categoriel_capital(Variable):
         f2go = foyer_fiscal('f2go', period)
         f2gr = foyer_fiscal('f2gr', period)
         f2tr = foyer_fiscal('f2tr', period)
-        rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
 
         f2dc_bis = f2dc
         f2tr_bis = f2tr
         # # Calcul du revenu catégoriel
         # 1.2 Revenus des valeurs et capitaux mobiliers
-        b12 = min_(f2ch, rvcm.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses))
+        b12 = min_(f2ch, abattement_assurance_vie)
         TOT1 = f2ch - b12  # c12
         # Part des frais s'imputant sur les revenus déclarés case DC
         den = ((f2dc_bis + f2ts) != 0) * (f2dc_bis + f2ts) + ((f2dc_bis + f2ts) == 0)
         F1 = f2ca / den * f2dc_bis  # f12
         # Revenus de capitaux mobiliers nets de frais, ouvrant droit à abattement
         # partie négative (à déduire des autres revenus nets de frais d'abattements
-        g12a = -min_(f2dc_bis * (1 - rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) - F1, 0)
+        g12a = -min_(f2dc_bis * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) - F1, 0)
         # partie positive
-        g12b = max_(f2dc_bis * (1 - rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) - F1, 0)
-        rev = g12b + f2gr + f2fu * (1 - rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement)
+        g12b = max_(f2dc_bis * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) - F1, 0)
+        rev = g12b + f2gr + f2fu * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement)
 
         # Abattements, limité au revenu
-        h12 = rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses)
+        h12 = parameters_rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses)
         TOT2 = max_(0, rev - h12)
         # i121= -min_(0,rev - h12)
 
         # Part des frais s'imputant sur les revenus déclarés ligne TS
         F2 = f2ca - F1
-        TOT3 = (f2ts - F2) + f2go * rvcm.majoration_revenus_reputes_distribues + f2tr_bis - g12a
+        TOT3 = (f2ts - F2) + f2go * parameters_rvcm.majoration_revenus_reputes_distribues + f2tr_bis - g12a
 
         DEF = deficit_rcm
         return max_(TOT1 + TOT2 + TOT3 - DEF, 0)
@@ -735,6 +740,9 @@ class revenu_categoriel_capital(Variable):
         Revenus des valeurs et capitaux mobiliers
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
         f2ch = foyer_fiscal('f2ch', period)
         f2da = foyer_fiscal('f2da', period)
@@ -745,23 +753,25 @@ class revenu_categoriel_capital(Variable):
         f2go = foyer_fiscal('f2go', period)
         f2gr = foyer_fiscal('f2gr', period)
         f2tr = foyer_fiscal('f2tr', period)
-        parameter_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
 
         part_frais_imputes_sur_f2dc = f2ca / max_(1, f2dc + f2ts) * f2dc
-        part_frais_restant_a_imputer = -min_(f2dc * (1 - parameter_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - part_frais_imputes_sur_f2dc, 0)
+        part_frais_restant_a_imputer = -min_(f2dc * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - part_frais_imputes_sur_f2dc, 0)
 
-        dividendes_apres_abattements = max_(f2dc * (1 - parameter_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - part_frais_imputes_sur_f2dc, 0)
-        revenus_assurance_vie_apres_abattements = f2ch - min_(f2ch, parameter_rvcm.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses))
+        dividendes_apres_abattements = max_(f2dc * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - part_frais_imputes_sur_f2dc, 0)
+        revenus_assurance_vie_apres_abattements = f2ch - min_(f2ch, abattement_assurance_vie)
         rvcm_apres_abattements_proportionnels = (
             revenus_assurance_vie_apres_abattements
             + dividendes_apres_abattements
             + f2gr
-            + f2fu * (1 - parameter_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0))
+            + f2fu * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0))
             )
-        rvcm_apres_abattements_proportionnels_et_fixes = max_(0, rvcm_apres_abattements_proportionnels - parameter_rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses))
+        rvcm_apres_abattements_proportionnels_et_fixes = max_(0, rvcm_apres_abattements_proportionnels - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses))
         autres_rvcm_sans_abattements = (
             f2ts - (f2ca - part_frais_imputes_sur_f2dc)
-            + f2go * parameter_rvcm.majoration_revenus_reputes_distribues
+            + f2go * parameters_rvcm.majoration_revenus_reputes_distribues
             + f2tr - part_frais_restant_a_imputer
             )
 
@@ -772,6 +782,9 @@ class revenu_categoriel_capital(Variable):
         Revenus des valeurs et capitaux mobiliers
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
         f2ca = foyer_fiscal('f2ca', period)
         f2ch = foyer_fiscal('f2ch', period)
@@ -780,15 +793,16 @@ class revenu_categoriel_capital(Variable):
         f2go = foyer_fiscal('f2go', period)
         f2tr = foyer_fiscal('f2tr', period)
         f2ts = foyer_fiscal('f2ts', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
+
         rvcm_apres_abattement = (
             f2fu + f2dc - abattement_dividende
             + f2ch - min_(f2ch, abattement_assurance_vie)
-            + f2ts + f2tr + f2go * P.majoration_revenus_reputes_distribues
+            + f2ts + f2tr + f2go * parameters_rvcm.majoration_revenus_reputes_distribues
             )
 
         return max_(0, rvcm_apres_abattement - f2ca - deficit_rcm)
@@ -798,6 +812,9 @@ class revenu_categoriel_capital(Variable):
         Revenus des valeurs et capitaux mobiliers
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
         f2ca = foyer_fiscal('f2ca', period)
         f2ch = foyer_fiscal('f2ch', period)
@@ -808,15 +825,16 @@ class revenu_categoriel_capital(Variable):
         f2ts = foyer_fiscal('f2ts', period)
         f2tt_2016 = foyer_fiscal('f2tt_2016', period)
         f2tu_2016 = foyer_fiscal('f2tu_2016', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
+
         rvcm_apres_abattement = (
             f2fu + f2dc - abattement_dividende
             + f2ch - min_(f2ch, abattement_assurance_vie)
-            + f2ts + f2tr + max_(0, f2tt_2016 - f2tu_2016) + f2go * P.majoration_revenus_reputes_distribues
+            + f2ts + f2tr + max_(0, f2tt_2016 - f2tu_2016) + f2go * parameters_rvcm.majoration_revenus_reputes_distribues
             )
 
         return max_(0, rvcm_apres_abattement - f2ca - deficit_rcm)
@@ -834,6 +852,9 @@ class revenu_categoriel_capital(Variable):
               - Brochure pratique de l'IR 2017 sur revenus 2016 : https://www.impots.gouv.fr/portail/www2/fichiers/documentation/brochure/ir_2017/files/assets/common/downloads/publication.pdf
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
         f2ca = foyer_fiscal('f2ca', period)
         f2ch = foyer_fiscal('f2ch', period)
@@ -843,15 +864,15 @@ class revenu_categoriel_capital(Variable):
         f2tr = foyer_fiscal('f2tr', period)
         f2ts = foyer_fiscal('f2ts', period)
         f2tt = foyer_fiscal('f2tt', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
         rvcm_apres_abattement = (
             f2fu + f2dc - abattement_dividende
             + f2ch - min_(f2ch, abattement_assurance_vie)
-            + f2ts + f2tr + f2tt + f2go * P.majoration_revenus_reputes_distribues
+            + f2ts + f2tr + f2tt + f2go * parameters_rvcm.majoration_revenus_reputes_distribues
             )
 
         return max_(0, rvcm_apres_abattement - f2ca - deficit_rcm)
@@ -866,8 +887,11 @@ class revenu_categoriel_capital(Variable):
         Dans ce cas, ils ne sortent pas de la variable `revenu_categoriel_capital`.
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
         imposition_au_bareme = foyer_fiscal('f2op', period)
 
         # Revenus à prendre en compte dans les deux cas: pfu ou imposition au barème
@@ -887,12 +911,13 @@ class revenu_categoriel_capital(Variable):
         f2zz = foyer_fiscal('f2zz', period)
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
         abattement_residuel = max_(abattement_assurance_vie - f2ch, 0)
         abattement_residuel2 = max_(abattement_residuel - f2vv, 0)
+
         pre_result = where(imposition_au_bareme, f2zz + max_(f2vv - abattement_residuel, 0) + max_(f2ww - abattement_residuel2, 0) + f2fu + f2dc - abattement_dividende
-            + f2ts + f2tr + f2tt + f2go * P.majoration_revenus_reputes_distribues, 0)
+            + f2ts + f2tr + f2tt + f2go * parameters_rvcm.majoration_revenus_reputes_distribues, 0)
         rvcm_apres_abattement = (
             f2yy
             + f2ch - min_(f2ch, abattement_assurance_vie)
@@ -911,8 +936,11 @@ class revenu_categoriel_capital(Variable):
             Source : Brochure pratique revenus 2019 page 123 et 340: https://www.impots.gouv.fr/www2/fichiers/documentation/brochure/ir_2020/accueil.htm
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
         imposition_au_bareme = foyer_fiscal('f2op', period)
 
         # Revenus à prendre en compte dans les deux cas: pfu ou imposition au barème
@@ -933,12 +961,13 @@ class revenu_categoriel_capital(Variable):
         f2tq = foyer_fiscal('f2tq', period)
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
         abattement_residuel = max_(abattement_assurance_vie - f2ch, 0)
         abattement_residuel2 = max_(abattement_residuel - f2vv, 0)
+
         pre_result = where(imposition_au_bareme, f2zz + max_(f2vv - abattement_residuel, 0) + max_(f2ww - abattement_residuel2, 0) + f2fu + f2dc - abattement_dividende
-            + f2ts + f2tr + f2tt + f2go * P.majoration_revenus_reputes_distribues + f2tq, 0)
+            + f2ts + f2tr + f2tt + f2go * parameters_rvcm.majoration_revenus_reputes_distribues + f2tq, 0)
         rvcm_apres_abattement = (
             f2yy
             + f2ch - min_(f2ch, abattement_assurance_vie)
@@ -957,8 +986,11 @@ class revenu_categoriel_capital(Variable):
             Source : Brochure pratique revenus 2020 page 119, 132 et 364: https://www.impots.gouv.fr/www2/fichiers/documentation/brochure/ir_2021/accueil.htm
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         deficit_rcm = foyer_fiscal('deficit_rcm', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
         imposition_au_bareme = foyer_fiscal('f2op', period)
 
         # Revenus à prendre en compte dans les deux cas: pfu ou imposition au barème
@@ -980,12 +1012,13 @@ class revenu_categoriel_capital(Variable):
         f2tz = foyer_fiscal('f2tz', period)
 
         # Revenus après abatemment
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
-        abattement_assurance_vie = P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
         abattement_residuel = max_(abattement_assurance_vie - f2ch, 0)
         abattement_residuel2 = max_(abattement_residuel - f2vv, 0)
+
         pre_result = where(imposition_au_bareme, f2zz + max_(f2vv - abattement_residuel, 0) + max_(f2ww - abattement_residuel2, 0) + f2fu + f2dc - abattement_dividende
-            + f2ts + f2tr + f2tt + f2go * P.majoration_revenus_reputes_distribues + f2tq + f2tz, 0)
+            + f2ts + f2tr + f2tt + f2go * parameters_rvcm.majoration_revenus_reputes_distribues + f2tq + f2tz, 0)
         rvcm_apres_abattement = (
             f2yy
             + f2ch - min_(f2ch, abattement_assurance_vie)
@@ -1010,7 +1043,7 @@ class rfr_rvcm_abattements_a_reintegrer(Variable):
         f2gr = foyer_fiscal('f2gr', period)
         f2fu = foyer_fiscal('f2fu', period)
         f2da = foyer_fiscal('f2da', period)  # noqa F841
-        rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
         # Calcul de i121
         # Part des frais s'imputant sur les revenus déclarés case DC
@@ -1018,20 +1051,20 @@ class rfr_rvcm_abattements_a_reintegrer(Variable):
         F1 = f2ca / den * f2dc  # f12
         # Revenus de capitaux mobiliers nets de frais, ouvrant droit à abattement
         # partie positive
-        g12b = max_(f2dc * (1 - rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - F1, 0)
-        rev = g12b + f2gr + f2fu * (1 - rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0))
+        g12b = max_(f2dc * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0)) - F1, 0)
+        rev = g12b + f2gr + f2fu * (1 - parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement * (f2da == 0))
 
         # Abattements, limité au revenu
-        h12 = rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses)
+        h12 = parameters_rvcm.revenus_capitaux_mobiliers_dividendes.abattement_forfaitaire * (1 + maries_ou_pacses)
         i121 = - min_(0, rev - h12)
-        return max_((rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) * (f2dc + f2fu) * (f2da == 0) - i121, 0)
+        return max_((parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement) * (f2dc + f2fu) * (f2da == 0) - i121, 0)
 
     def formula_2013_01_01(foyer_fiscal, period, parameters):
         f2dc = foyer_fiscal('f2dc', period)
         f2fu = foyer_fiscal('f2fu', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
-        abattement_dividende = (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement
+        abattement_dividende = (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement
 
         return abattement_dividende
 
@@ -1043,6 +1076,9 @@ class rfr_rvcm_abattements_a_reintegrer(Variable):
               Si le foyer a choisi l'imposition au barème pour les revenus éligibles au pfu, les revenus de l'assurance-vie entrent dans le calcul du RFR via `revenus_categoriel` net d'abattement.
         '''
         maries_ou_pacses = foyer_fiscal('maries_ou_pacses', period)
+        celibataire_ou_divorce = foyer_fiscal('celibataire_ou_divorce', period)
+        veuf = foyer_fiscal('veuf', period)
+        jeune_veuf = foyer_fiscal('jeune_veuf', period)
         imposition_au_bareme = foyer_fiscal('f2op', period)
         f2ch = foyer_fiscal('f2ch', period)
         f2dh = foyer_fiscal('f2dh', period)
@@ -1050,14 +1086,15 @@ class rfr_rvcm_abattements_a_reintegrer(Variable):
         f2ww = foyer_fiscal('f2ww', period)
         f2dc = foyer_fiscal('f2dc', period)
         f2fu = foyer_fiscal('f2fu', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
-        abattement_assurance_vie = where(imposition_au_bareme, 0,
-            (f2ch < P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses)) * max_(0, min_(f2vv + f2ww, P.produits_assurances_vies_assimiles.abattement * (1 + maries_ou_pacses) - f2ch - f2dh))
+        abattement_assurance_vie = parameters_rvcm.produits_assurances_vies_assimiles.abattement_couple * maries_ou_pacses + parameters_rvcm.produits_assurances_vies_assimiles.abattement_celib * (celibataire_ou_divorce | veuf | jeune_veuf)
+        abattement_assu_vie = where(imposition_au_bareme, 0,
+            (f2ch < abattement_assurance_vie) * max_(0, min_(f2vv + f2ww, abattement_assurance_vie - f2ch - f2dh))
             )
-        abattement_dividende = where(imposition_au_bareme, (f2fu + f2dc) * P.revenus_capitaux_mobiliers_dividendes.taux_abattement, 0)
+        abattement_dividende = where(imposition_au_bareme, (f2fu + f2dc) * parameters_rvcm.revenus_capitaux_mobiliers_dividendes.taux_abattement, 0)
 
-        return - abattement_assurance_vie + abattement_dividende
+        return - abattement_assu_vie + abattement_dividende
 
 
 class revenu_categoriel_foncier(Variable):
@@ -1549,8 +1586,8 @@ class ir_plaf_qf(Variable):
         residence_dom = (residence_guadeloupe_martinique_reunion | residence_guyane_mayotte)
 
         abattement_dom = (
-            residence_guadeloupe_martinique_reunion * min_(plafond_qf.abat_dom.plaf_GuadMarReu, plafond_qf.abat_dom.taux_GuadMarReu * impot_apres_reduction_complementaire)
-            + residence_guyane_mayotte * min_(plafond_qf.abat_dom.plaf_GuyMay, plafond_qf.abat_dom.taux_GuyMay * impot_apres_reduction_complementaire)
+            residence_guadeloupe_martinique_reunion * min_(plafond_qf.abat_dom.guadeloupe_martinique_reunion.plaf, plafond_qf.abat_dom.guadeloupe_martinique_reunion.taux * impot_apres_reduction_complementaire)
+            + residence_guyane_mayotte * min_(plafond_qf.abat_dom.guyane_mayotte.plaf, plafond_qf.abat_dom.guyane_mayotte.taux * impot_apres_reduction_complementaire)
             )
 
         impot_apres_abattement_dom = max_(0, impot_apres_reduction_complementaire - abattement_dom)
@@ -1637,13 +1674,13 @@ class reduction_ss_condition_revenus(Variable):
         nb_adult = foyer_fiscal('nb_adult', period)
         nb_parts = foyer_fiscal('nbptr', period)
         rfr = foyer_fiscal('rfr', period)
-        P = parameters(period).impot_revenu.calcul_impot_revenu.plaf_qf.reduction_ss_condition_revenus
+        parameters = parameters(period).impot_revenu.calcul_impot_revenu.plaf_qf.reduction_ss_condition_revenus
 
         ir_apres_plaf_qf_et_decote = ir_plaf_qf - decote
-        plafond1 = P.plafond_rfr_celib * nb_adult + P.majoration_plafond_par_demi_parts_supp * 2 * (nb_parts - nb_adult)
-        plafond2 = P.plafond_rfr_couple * nb_adult + P.majoration_plafond_par_demi_parts_supp * 2 * (nb_parts - nb_adult)
-        reduction1 = P.taux * ir_apres_plaf_qf_et_decote
-        reduction2 = P.taux * ir_apres_plaf_qf_et_decote * (plafond2 - rfr) / (plafond2 - plafond1)
+        plafond1 = parameters.plafond_rfr_celib * nb_adult + parameters.majoration_plafond_par_demi_parts_supp * 2 * (nb_parts - nb_adult)
+        plafond2 = parameters.plafond_rfr_couple * nb_adult + parameters.majoration_plafond_par_demi_parts_supp * 2 * (nb_parts - nb_adult)
+        reduction1 = parameters.taux * ir_apres_plaf_qf_et_decote
+        reduction2 = parameters.taux * ir_apres_plaf_qf_et_decote * (plafond2 - rfr) / (plafond2 - plafond1)
 
         reduction_sous_condition_de_ressources = (
             (rfr < plafond1) * reduction1
@@ -1667,8 +1704,6 @@ class nat_imp(Variable):
         credits_impot = foyer_fiscal('credits_impot', period)
         contribution_exceptionnelle_hauts_revenus = foyer_fiscal('contribution_exceptionnelle_hauts_revenus', period)
 
-        # def _nat_imp(rni, nbptr, non_imposable = law.impot_revenu.calcul_impot_revenu.non_imposable):
-        # seuil = non_imposable.seuil + (nbptr - 1)*non_imposable.supp
         return (iai - credits_impot + contribution_exceptionnelle_hauts_revenus) > 0
 
 
@@ -1689,6 +1724,15 @@ class ip_net(Variable):
 
         return around(max_(0, ir_plaf_qf - decote - reduction_ss_condition_revenus))
 
+    def formula_2020_01_01(foyer_fiscal, period, parameters):
+        '''
+        Impôt net avant réductions suite à l'abrogation de la "réfaction foyers modestes"
+        '''
+        ir_plaf_qf = foyer_fiscal('ir_plaf_qf', period)
+        decote_gain_fiscal = foyer_fiscal('decote_gain_fiscal', period)
+
+        return around(ir_plaf_qf - decote_gain_fiscal)
+
 
 class iaidrdi(Variable):
     value_type = float
@@ -1769,8 +1813,8 @@ class assiette_service(Variable):
 
         return foyer_fiscal.sum(ebic_imps_i)
 
-    # P = _P.impot_revenu.calcul_revenus_imposables.rpns.micro.microentreprise
-    # assert (ebic_imps <= P.servi.max)
+    # parameters = _P.impot_revenu.calcul_revenus_imposables.rpns.micro.microentreprise
+    # assert (ebic_imps <= parameters.servi.max)
 
 
 class assiette_proflib(Variable):
@@ -1789,7 +1833,7 @@ class assiette_proflib(Variable):
         # http://vosdroits.service-public.fr/professionnels-entreprises/F23267.xhtml
         return foyer_fiscal.sum(ebnc_impo_i)
 
-    # assert (ebnc_impo <= P.microentreprise.regime_micro_bnc.services.plafond)
+    # assert (ebnc_impo <= parameters.microentreprise.regime_micro_bic.services.plafond)
 
 
 class microsocial(Variable):
@@ -1828,8 +1872,8 @@ class microentreprise(Variable):
         ebic_impv = foyer_fiscal.sum(ebic_impv_i)
         return (
             max_(0, ebnc_impo - max_(micro.microentreprise.montant_minimum, micro.microentreprise.regime_micro_bnc.taux * ebnc_impo))
-            + max_(0, ebic_imps - max_(micro.microentreprise.montant_minimum, micro.microentreprise.regime_micro_bnc.services.taux * ebic_imps))
-            + max_(0, ebic_impv - max_(micro.microentreprise.montant_minimum, micro.microentreprise.regime_micro_bnc.marchandises.taux * ebic_impv))
+            + max_(0, ebic_imps - max_(micro.microentreprise.montant_minimum, micro.microentreprise.regime_micro_bic.services.taux * ebic_imps))
+            + max_(0, ebic_impv - max_(micro.microentreprise.montant_minimum, micro.microentreprise.regime_micro_bic.marchandises.taux * ebic_impv))
             )
 
 
@@ -1846,9 +1890,9 @@ class tax_rvcm_forfaitaire(Variable):
         Taxation des revenus des valeurs et capitaux mobiliers
         '''
         f2fa = foyer_fiscal('f2fa', period)
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
+        parameters_rvcm = parameters(period).impot_revenu.calcul_revenus_imposables.rvcm
 
-        return f2fa * P.taux_forfaitaire
+        return f2fa * parameters_rvcm.taux_forfaitaire
 
 
 class taxation_plus_values_hors_bareme(Variable):
@@ -1879,7 +1923,7 @@ class taxation_plus_values_hors_bareme(Variable):
         rpns_pvce = foyer_fiscal.sum(rpns_pvce_i)
 
         return round_(
-            pv.plus_values.pvce * rpns_pvce
+            pv.pvce.taux * rpns_pvce
             + pv.pv_cession_valeurs_mobilieres_pv_professionnelles.taux * max_(0, f3vg - f3vh)
             + pv.actions_gratuites.taux2 * glo_taxation_ir_forfaitaire_taux2
             + pv.pv_cession_valeurs_mobilieres_pv_professionnelles.taux * f3vl
@@ -1888,8 +1932,8 @@ class taxation_plus_values_hors_bareme(Variable):
             + pv.plus_values.taux_pv_entrep * f3sa_2012
             + pv.actions_gratuites.taux3 * glo_taxation_ir_forfaitaire_taux3
             + pv.actions_gratuites.taux4 * glo_taxation_ir_forfaitaire_taux4
-            + pv.bspce.plus_3ans.pre_2018 * f3sj
-            + pv.bspce.moins_3ans * f3sk
+            + pv.bspce.taux_plus_3_ans_pre_2018 * f3sj
+            + pv.bspce.taux_moins_3_ans * f3sk
             )
 
     def formula_2013_01_01(foyer_fiscal, period, parameters):
@@ -1909,14 +1953,14 @@ class taxation_plus_values_hors_bareme(Variable):
         pv = parameters(period).impot_revenu.calcul_impot_revenu.pv
 
         return round_(
-            pv.plus_values.pvce * rpns_pvce
+            pv.pvce.taux * rpns_pvce
             + pv.pea.taux_avant_2_ans * f3vm
             + pv.pea.taux_posterieur * f3vt
             + pv.actions_gratuites.taux2 * glo_taxation_ir_forfaitaire_taux2
             + pv.actions_gratuites.taux3 * glo_taxation_ir_forfaitaire_taux3
             + pv.actions_gratuites.taux4 * glo_taxation_ir_forfaitaire_taux4
-            + pv.bspce.plus_3ans.pre_2018 * f3sj
-            + pv.bspce.moins_3ans * f3sk
+            + pv.bspce.taux_plus_3_ans_pre_2018 * f3sj
+            + pv.bspce.taux_moins_3_ans * f3sk
             )
 
     def formula_2016_01_01(foyer_fiscal, period, parameters):
@@ -1938,16 +1982,16 @@ class taxation_plus_values_hors_bareme(Variable):
         pv = parameters(period).impot_revenu.calcul_impot_revenu.pv
 
         return round_(
-            pv.plus_values.pvce * rpns_pvce
+            pv.pvce.taux * rpns_pvce
             + pv.pea.taux_avant_2_ans * f3vm
             + pv.pea.taux_posterieur * f3vt
             + pv.actions_gratuites.taux2 * glo_taxation_ir_forfaitaire_taux2
             + pv.actions_gratuites.taux3 * glo_taxation_ir_forfaitaire_taux3
             + pv.actions_gratuites.taux4 * glo_taxation_ir_forfaitaire_taux4
-            + pv.bspce.plus_3ans.pre_2018 * f3sj
-            + pv.bspce.moins_3ans * f3sk
-            + pv.plus_values.taux_plus_values_report * f3wi
-            + pv.plus_values.taux_plus_values_report_conditionnel * f3wj
+            + pv.bspce.taux_plus_3_ans_pre_2018 * f3sj
+            + pv.bspce.taux_moins_3_ans * f3sk
+            + pv.report_impot_expire.taux_cas_general * f3wi
+            + pv.report_impot_expire.taux_conditionnel * f3wj
             )
 
     def formula_2017_01_01(foyer_fiscal, period, parameters):
@@ -1971,17 +2015,17 @@ class taxation_plus_values_hors_bareme(Variable):
         pv = parameters(period).impot_revenu.calcul_impot_revenu.pv
 
         return round_(
-            pv.plus_values.pvce * rpns_pvce
+            pv.pvce.taux * rpns_pvce
             + pv.actions_gratuites.taux2 * glo_taxation_ir_forfaitaire_taux2
             + pv.actions_gratuites.taux3 * glo_taxation_ir_forfaitaire_taux3
             + pv.actions_gratuites.taux4 * glo_taxation_ir_forfaitaire_taux4
-            + pv.bspce.plus_3ans.pre_2018 * f3sj
-            + pv.bspce.moins_3ans * f3sk
+            + pv.bspce.taux_plus_3_ans_pre_2018 * f3sj
+            + pv.bspce.taux_moins_3_ans * f3sk
             + pv.pea.taux_avant_2_ans * f3vm
             + pv.pea.taux_posterieur * f3vt
-            + pv.plus_values.taux_plus_values_report * f3wi
-            + pv.plus_values.taux_plus_values_report_conditionnel * f3wj
-            + pv.plus_values.taux_plus_values_entc * f3pi
+            + pv.report_impot_expire.taux_cas_general * f3wi
+            + pv.report_impot_expire.taux_conditionnel * f3wj
+            + pv.etnc.taux * f3pi
             )
 
     def formula_2019_01_01(foyer_fiscal, period, parameters):
@@ -2004,19 +2048,19 @@ class taxation_plus_values_hors_bareme(Variable):
         rpns_pvce = foyer_fiscal.sum(rpns_pvce_i)
         rpns_info = foyer_fiscal.sum(rpns_info_i)
         pv = parameters(period).impot_revenu.calcul_impot_revenu.pv
-        P = parameters(period).impot_revenu.calcul_revenus_imposables.rpns
+        parameters_rpns = parameters(period).impot_revenu.calcul_revenus_imposables.rpns
 
         return round_(
-            pv.plus_values.pvce * rpns_pvce
+            pv.pvce.taux * rpns_pvce
             + pv.actions_gratuites.taux2 * glo_taxation_ir_forfaitaire_taux2
             + pv.actions_gratuites.taux3 * glo_taxation_ir_forfaitaire_taux3
             + pv.actions_gratuites.taux4 * glo_taxation_ir_forfaitaire_taux4
-            + P.taux10 * rpns_info
-            + pv.bspce.plus_3ans.pre_2018 * f3sj
-            + pv.bspce.moins_3ans * f3sk
-            + pv.plus_values.taux_plus_values_report * f3wi
-            + pv.plus_values.taux_plus_values_report_conditionnel * f3wj
-            + pv.plus_values.taux_plus_values_entc * f3pi
+            + parameters_rpns.taux10 * rpns_info
+            + pv.bspce.taux_plus_3_ans_pre_2018 * f3sj
+            + pv.bspce.taux_moins_3_ans * f3sk
+            + pv.report_impot_expire.taux_cas_general * f3wi
+            + pv.report_impot_expire.taux_conditionnel * f3wj
+            + pv.etnc.taux * f3pi
             )
 
 
@@ -2264,10 +2308,10 @@ class contribution_exceptionnelle_hauts_revenus(Variable):
         # TODO: Gérer le II.-1 du lissage interannuel ? (problème de non recours)
 
 
-class impot_revenu_restant_a_payer(Variable):
+class impot_revenu_avant_seuils_mise_recouvrement(Variable):
     value_type = float
     entity = FoyerFiscal
-    label = 'Impôt sur le revenu des personnes physiques restant à payer, après prise en compte des éventuels acomptes'
+    label = 'Impôt sur le revenu des personnes physiques restant à payer, après prise en compte des éventuels acomptes et avant application des seuils de mise en recouvrement'
     reference = 'https://www.legifrance.gouv.fr/codes/article_lc/LEGIARTI000041464766'
     definition_period = YEAR
 
@@ -2276,40 +2320,63 @@ class impot_revenu_restant_a_payer(Variable):
         Montant après seuil de recouvrement (hors ppe)
         NB : ce montant l'impôt correspond à une notion administrative :
         dans certains cas, il existe des prélèvements à la source faisant
-        office d'acomptes d'impôt sur le revenu (cf. variable acomptes_ir). Ces acomptes sont comptabilisés
-        dans la feuille d'impôt comme des crédits d'impôt, mais correspondent économiquement à des montants d'impôt dus.
-        Le prélèvement forfaitaire libératoire a déjà été payé il ne doit donc pas être compté dans l'impôt restant à payer.
-        En revanche, il compte dans le calcul du seuil de recouvrement.
+        office d'acomptes d'impôt sur le revenu (cf. variable acomptes_ir). Ces acomptes sont comptabilisés dans la feuille d'impôt comme des crédits d'impôt, mais correspondent économiquement à des montants d'impôt dus.
         '''
         iai = foyer_fiscal('iai', period)
         credits_impot = foyer_fiscal('credits_impot', period)
         acomptes_ir = foyer_fiscal('acomptes_ir', period)
-        contribution_exceptionnelle_hauts_revenus = foyer_fiscal('contribution_exceptionnelle_hauts_revenus', period)
-        prelevement_forfaitaire_unique_ir = foyer_fiscal('prelevement_forfaitaire_unique_ir', period)
-        prelevement_forfaitaire_liberatoire = foyer_fiscal('prelevement_forfaitaire_liberatoire', period)
-        P = parameters(period).impot_revenu.calcul_impot_revenu.recouvrement
+        cehr = foyer_fiscal('contribution_exceptionnelle_hauts_revenus', period)
+        pfu = foyer_fiscal('prelevement_forfaitaire_unique_ir', period)
 
-        pre_result = iai - credits_impot - acomptes_ir + contribution_exceptionnelle_hauts_revenus - prelevement_forfaitaire_unique_ir - prelevement_forfaitaire_liberatoire
-        result = iai - credits_impot - acomptes_ir + contribution_exceptionnelle_hauts_revenus - prelevement_forfaitaire_unique_ir
-        impots_totaux_avant_imputations = iai + contribution_exceptionnelle_hauts_revenus - prelevement_forfaitaire_unique_ir - prelevement_forfaitaire_liberatoire
+        return -(iai + cehr + pfu - credits_impot - acomptes_ir)
 
-        return (
-            (impots_totaux_avant_imputations > P.seuil) * (
-                (pre_result < P.min)
-                * (result > 0)
-                * result
-                * 0
-                + ((pre_result <= 0) + (pre_result >= P.min))
-                * (- result)
-                )
-            + (impots_totaux_avant_imputations <= P.seuil) * (
-                (pre_result < 0)
-                * (-result)
-                + (pre_result >= 0)
-                * 0
-                * result
-                )
-            )
+
+class correction_ir_seuils_recouvrement(Variable):
+    value_type = float
+    entity = FoyerFiscal
+    label = 'Montant de la correction de l´impôt sur le revenu par application des seuils de mise en recouvrement'
+    reference = 'https://bofip.impots.gouv.fr/bofip/2496-PGP.html/identifiant%3DBOI-IR-LIQ-20-20-40-20180704#Franchise_pour_les_impositi_14'
+    definition_period = YEAR
+
+    def formula(foyer_fiscal, period, parameters):
+        iai = foyer_fiscal('iai', period)
+        credits_impot = foyer_fiscal('credits_impot', period)
+        acomptes_ir = foyer_fiscal('acomptes_ir', period)
+        cehr = foyer_fiscal('contribution_exceptionnelle_hauts_revenus', period)
+        pfu = foyer_fiscal('prelevement_forfaitaire_unique_ir', period)
+        pfl = foyer_fiscal('prelevement_forfaitaire_liberatoire', period)
+        impots_avant_seuils_mise_recouvrement = foyer_fiscal('impot_revenu_avant_seuils_mise_recouvrement', period)
+        '''
+        Le prélèvement forfaitaire libératoire a déjà été payé il ne doit donc pas être compté dans l'impôt restant à payer.
+        En revanche, il compte dans le calcul du seuil de recouvrement.
+        '''
+        parameters_recouvrement = parameters(period).impot_revenu.calcul_impot_revenu.recouvrement
+        seuil_avant_imputations = parameters_recouvrement.min_avant_credits_impots
+        seuil_apres_imputations = parameters_recouvrement.min_apres_credits_impots
+
+        impots_totaux_avant_imputations = iai + cehr + pfu - pfl
+        impots_totaux_apres_imputations = iai + cehr + pfu - pfl - credits_impot - acomptes_ir
+
+        condition_1 = (impots_totaux_avant_imputations <= seuil_avant_imputations) * (impots_totaux_apres_imputations > 0)
+        condition_2 = (impots_totaux_avant_imputations > seuil_avant_imputations) * (impots_totaux_apres_imputations > 0) * (impots_totaux_apres_imputations <= seuil_apres_imputations)
+
+        condition_correction_seuils_mise_recouvrement = (condition_1 | condition_2)
+
+        return condition_correction_seuils_mise_recouvrement * impots_avant_seuils_mise_recouvrement
+
+
+class impot_revenu_restant_a_payer(Variable):
+    value_type = float
+    entity = FoyerFiscal
+    label = 'Impôt sur le revenu des personnes physiques restant à payer, après prise en compte des éventuels acomptes et des seuils de mise en recouvrement'
+    reference = 'https://www.legifrance.gouv.fr/codes/article_lc/LEGIARTI000041464766'
+    definition_period = YEAR
+
+    def formula(foyer_fiscal, period, parameters):
+        impots_avant_seuils_recouvrement = foyer_fiscal('impot_revenu_avant_seuils_mise_recouvrement', period)
+        correction_seuils_recouvrement = foyer_fiscal('correction_ir_seuils_recouvrement', period)
+
+        return (impots_avant_seuils_recouvrement - correction_seuils_recouvrement)
 
 
 class foyer_impose(Variable):
@@ -2500,6 +2567,7 @@ class rpns_pvce(Variable):
     entity = Individu
     label = 'Plus values de cession nettes des moins-values - Revenu des professions non salariées'
     definition_period = YEAR
+    # Somme des cases "plus-values nettes à long terme" de la déclaration sur les revenus des professions non-salariées
 
     def formula(individu, period, parameters):
         '''
@@ -2801,15 +2869,15 @@ class defacc(Variable):
         cga = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.cga_taux2
         micro = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.micro
 
-        def abat_rpns(rev, P):
-            return max_(0, rev - min_(rev, max_(P.taux * min_(P.plafond, rev), micro.microentreprise.montant_minimum)))
+        def abat_rpns(rev, parameters):
+            return max_(0, rev - min_(rev, max_(parameters.taux * min_(parameters.plafond, rev), micro.microentreprise.montant_minimum)))
 
         nacc_impn = foyer_fiscal.sum(nacc_impn_i)
         macc_pvct = foyer_fiscal.sum(macc_pvct_i)
         macc_impv = foyer_fiscal.sum(macc_impv_i)
         macc_imps = foyer_fiscal.sum(macc_imps_i)
         aacc_impn = foyer_fiscal.sum(aacc_impn_i)
-        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bnc.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bnc.services)
+        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bic.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bic.services)
         return (
             min_(f5rn + f5ro + f5rp + f5rq + f5rr + f5rw, aacc_impn + macc_pvct + macc_timp + (1 + cga) * nacc_impn)
             )
@@ -2827,14 +2895,14 @@ class defacc(Variable):
         aacc_impn_i = foyer_fiscal.members('aacc_impn', period)
         micro = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.micro
 
-        def abat_rpns(rev, P):
-            return max_(0, rev - min_(rev, max_(P.taux * min_(P.plafond, rev), micro.microentreprise.montant_minimum)))
+        def abat_rpns(rev, parameters):
+            return max_(0, rev - min_(rev, max_(parameters.taux * min_(parameters.plafond, rev), micro.microentreprise.montant_minimum)))
 
         macc_pvct = foyer_fiscal.sum(macc_pvct_i)
         macc_impv = foyer_fiscal.sum(macc_impv_i)
         macc_imps = foyer_fiscal.sum(macc_imps_i)
         aacc_impn = foyer_fiscal.sum(aacc_impn_i)
-        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bnc.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bnc.services)
+        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bic.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bic.services)
         return (
             min_(f5rn + f5ro + f5rp + f5rq + f5rr + f5rw, aacc_impn + macc_pvct + macc_timp)
             )
@@ -2861,15 +2929,15 @@ class defncn(Variable):
         micro = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.micro
         specialbnc = micro.microentreprise.regime_micro_bnc
 
-        def abat_rpns(rev, P):
-            return max_(0, rev - min_(rev, max_(P.taux * min_(P.plafond, rev), micro.microentreprise.montant_minimum)))
+        def abat_rpns(rev, parameters):
+            return max_(0, rev - min_(rev, max_(parameters.taux * min_(parameters.plafond, rev), micro.microentreprise.montant_minimum)))
         cncn_bene = foyer_fiscal.sum(cncn_bene_i)
         mncn_impo = foyer_fiscal.sum(mncn_impo_i)
         mncn_pvct = foyer_fiscal.sum(mncn_pvct_i)
         cncn_aimp = foyer_fiscal.sum(cncn_aimp_i)
         return min_(
             f5ht + f5it + f5jt + f5kt + f5lt + f5mt,
-            abat_rpns(mncn_impo, specialbnc.services) + mncn_pvct + cncn_aimp + (1 + cga) * cncn_bene
+            abat_rpns(mncn_impo, specialbnc) + mncn_pvct + cncn_aimp + (1 + cga) * cncn_bene
             )  # TODO check !
 
     def formula_2023_01_01(foyer_fiscal, period, parameters):
@@ -2885,14 +2953,14 @@ class defncn(Variable):
         micro = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.micro
         specialbnc = micro.microentreprise.regime_micro_bnc
 
-        def abat_rpns(rev, P):
-            return max_(0, rev - min_(rev, max_(P.taux * min_(P.plafond, rev), micro.microentreprise.montant_minimum)))
+        def abat_rpns(rev, parameters):
+            return max_(0, rev - min_(rev, max_(parameters.taux * min_(parameters.plafond, rev), micro.microentreprise.montant_minimum)))
         mncn_impo = foyer_fiscal.sum(mncn_impo_i)
         mncn_pvct = foyer_fiscal.sum(mncn_pvct_i)
         cncn_aimp = foyer_fiscal.sum(cncn_aimp_i)
         return min_(
             f5ht + f5it + f5jt + f5kt + f5lt + f5mt,
-            abat_rpns(mncn_impo, specialbnc.services) + mncn_pvct + cncn_aimp
+            abat_rpns(mncn_impo, specialbnc) + mncn_pvct + cncn_aimp
             )  # TODO check !
 
 
@@ -3136,11 +3204,11 @@ class aacc_timp(Variable):
                 0,
                 (aacc_impn + (aacc_gits > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_gits * (1 - micro.microentreprise.regime_micro_bnc.marchandises.taux)
+                    aacc_gits * (1 - micro.microentreprise.regime_micro_bic.marchandises.taux)
                     ))
                 + (aacc_imps > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_imps * (1 - micro.microentreprise.regime_micro_bnc.services.taux)
+                    aacc_imps * (1 - micro.microentreprise.regime_micro_bic.services.taux)
                     )
                 + max_(0, - alnp_defs) - aacc_defn
                 )
@@ -3161,15 +3229,15 @@ class aacc_timp(Variable):
                 0,
                 (aacc_impn + (aacc_gits > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_gits * (1 - micro.microentreprise.regime_micro_bnc.marchandises.taux)
+                    aacc_gits * (1 - micro.microentreprise.regime_micro_bic.marchandises.taux)
                     ))
                 + (aacc_imps > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_imps * (1 - micro.microentreprise.regime_micro_bnc.services.taux)
+                    aacc_imps * (1 - micro.microentreprise.regime_micro_bic.services.taux)
                     )
                 + (nacc_meup > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    nacc_meup * (1 - micro.microentreprise.regime_micro_bnc.marchandises.taux)
+                    nacc_meup * (1 - micro.microentreprise.regime_micro_bic.marchandises.taux)
                     )
                 + max_(0, nacc_pres - alnp_defs) - aacc_defn
                 )
@@ -3191,19 +3259,19 @@ class aacc_timp(Variable):
                 0,
                 (aacc_impn + (aacc_gits > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_gits * (1 - micro.microentreprise.regime_micro_bnc.marchandises.taux)
+                    aacc_gits * (1 - micro.microentreprise.regime_micro_bic.marchandises.taux)
                     ))
                 + (aacc_imps > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    aacc_imps * (1 - micro.microentreprise.regime_micro_bnc.services.taux)
+                    aacc_imps * (1 - micro.microentreprise.regime_micro_bic.services.taux)
                     )
                 + (nacc_meup > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    nacc_meup * (1 - micro.microentreprise.regime_micro_bnc.marchandises.taux)
+                    nacc_meup * (1 - micro.microentreprise.regime_micro_bic.marchandises.taux)
                     )
                 + (nacc_meuc > 0) * max_(
                     micro.microentreprise.montant_minimum,
-                    nacc_meuc * (1 - micro.microentreprise.regime_micro_bnc.services.taux)
+                    nacc_meuc * (1 - micro.microentreprise.regime_micro_bic.services.taux)
                     )
                 + max_(0, nacc_pres - alnp_defs) - aacc_defn
                 )
@@ -3431,20 +3499,20 @@ class rpns_imposables(Variable):
         locations_pro = individu('locations_pro', period)
         micro = parameters(period).impot_revenu.calcul_revenus_imposables.rpns.micro
 
-        def abat_rpns(rev, P):
-            return max_(0, rev - min_(rev, max_(P.taux * min_(P.plafond, rev), micro.microentreprise.montant_minimum)))
+        def abat_rpns(rev, parameters):
+            return max_(0, rev - min_(rev, max_(parameters.taux * min_(parameters.plafond, rev), micro.microentreprise.montant_minimum)))
 
         # # B revenus industriels et commerciaux professionnels
         # regime micro entreprise
-        mbic_timp = abat_rpns(mbic_impv, micro.microentreprise.regime_micro_bnc.marchandises) + abat_rpns(mbic_imps, micro.microentreprise.regime_micro_bnc.services)
+        mbic_timp = abat_rpns(mbic_impv, micro.microentreprise.regime_micro_bic.marchandises) + abat_rpns(mbic_imps, micro.microentreprise.regime_micro_bic.services)
         # regime micro entreprise
-        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bnc.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bnc.services)
+        macc_timp = abat_rpns(macc_impv, micro.microentreprise.regime_micro_bic.marchandises) + abat_rpns(macc_imps, micro.microentreprise.regime_micro_bic.services)
         # # D revenus non commerciaux professionnels
         # regime déclaratif special ou micro-bnc
-        mbnc_timp = abat_rpns(mbnc_impo, micro.microentreprise.regime_micro_bnc.services)
+        mbnc_timp = abat_rpns(mbnc_impo, micro.microentreprise.regime_micro_bnc)
         # # E revenus non commerciaux non professionnels
         # regime déclaratif special ou micro-bnc
-        mncn_timp = abat_rpns(mncn_impo, micro.microentreprise.regime_micro_bnc.services)
+        mncn_timp = abat_rpns(mncn_impo, micro.microentreprise.regime_micro_bnc)
         macc_mvct = individu.foyer_fiscal('macc_mvct', period) * individu.has_role(FoyerFiscal.DECLARANT_PRINCIPAL)
         mncn_mvct = individu.foyer_fiscal('mncn_mvct', period) * individu.has_role(FoyerFiscal.DECLARANT_PRINCIPAL)
 
@@ -3498,7 +3566,7 @@ class abat_spe(Variable):
         abattements = parameters(period).impot_revenu.calcul_revenus_imposables.abat_rni
 
         # Abattement pour personnes agées de + de 65 ans ou invalide
-        abattement_age_ou_invalidite = abattements.personne_agee_ou_invalide
+        abattement_age_ou_invalidite = abattements.contribuable_age_invalide
 
         # Abattement pour rattachement d'enfants mari·é·s
         abattement_enfant_marie = abattements.enfant_marie
@@ -3512,25 +3580,13 @@ class abat_spe(Variable):
         # Vecteur de montants d'abattement pour personnes âges ou invalides
         as_inv = (
             foyers_eligibles
-            * (
-                (
-                    + abattement_age_ou_invalidite.montant_1
-                    * (revenu_net_global <= abattement_age_ou_invalidite.plafond_1)
-                    )
-                + (
-                    + abattement_age_ou_invalidite.montant_2
-                    * (
-                        (revenu_net_global > abattement_age_ou_invalidite.plafond_1)
-                        & (revenu_net_global <= abattement_age_ou_invalidite.plafond_2)
-                        )
-                    )
-                )
+            * abattement_age_ou_invalidite.calc(revenu_net_global)
             )
 
         # Vecteur de montants d'abattement pour enfants à charge
         as_enf = (
             nombre_enfants
-            * abattement_enfant_marie.montant
+            * abattement_enfant_marie
             )
 
         # Le montant total d'abattement ne peut pas être supérieur au revenu net global
@@ -3634,29 +3690,29 @@ class nbptr(Variable):
 
         # parts des enfants en résidence alternée quand il n'y a que des enfants en résidence alternée
         enf_a = (no_pac & has_alt) * (
-            quotient_familial.enf1 * min_(nbH, 1)
-            + quotient_familial.enf2 * max_(min_(nbH - 1, 1), 0)
-            + quotient_familial.enf3_et_sup * max_(nbH - 2, 0)
+            quotient_familial.cas_general.enf1 * min_(nbH, 1)
+            + quotient_familial.cas_general.enf2 * max_(min_(nbH - 1, 1), 0)
+            + quotient_familial.cas_general.enf3_et_sup * max_(nbH - 2, 0)
             ) * 0.5
         # parts des enfants en résidence alternée quand il y a aussi des enfants entièrement à charge
         enf_b = (has_pac & has_alt) * (
             (nb_pac == 1) * (
-                quotient_familial.enf2 * min_(nbH, 1)
-                + quotient_familial.enf3_et_sup * max_(nbH - 1, 0)
+                quotient_familial.cas_general.enf2 * min_(nbH, 1)
+                + quotient_familial.cas_general.enf3_et_sup * max_(nbH - 1, 0)
                 ) * 0.5
-            + (nb_pac > 1) * (quotient_familial.enf3_et_sup * nbH * 0.5)
+            + (nb_pac > 1) * (quotient_familial.cas_general.enf3_et_sup * nbH * 0.5)
             )
         # parts des enfants entièrement à charge
         enf_c = (
-            quotient_familial.enf1 * min_(nb_pac, 1)
-            + quotient_familial.enf2 * max_(min_(nb_pac - 1, 1), 0)
-            + quotient_familial.enf3_et_sup * max_(nb_pac - 2, 0)
+            quotient_familial.cas_general.enf1 * min_(nb_pac, 1)
+            + quotient_familial.cas_general.enf2 * max_(min_(nb_pac - 1, 1), 0)
+            + quotient_familial.cas_general.enf3_et_sup * max_(nb_pac - 2, 0)
             )
 
         enf = enf_a + enf_b + enf_c
 
         # # note 2 : nombre de parts liées aux invalides (enfant + adulte)
-        n2 = quotient_familial.inv1 * (nbG + nbI / 2) + quotient_familial.inv2 * nbR
+        n2 = quotient_familial.couple_ou_pers_a_charge.inv1 * (nbG + nbI / 2) + quotient_familial.couple_ou_pers_a_charge.inv2 * nbR
 
         # # note 3 : Pas de personne à charge
         # - invalide
@@ -3670,23 +3726,23 @@ class nbptr(Variable):
         n3 = max_(n31, n32)
         # # note 4 Invalidité de la personne ou du conjoint pour les mariés ou
         # # jeunes veuf(ve)s
-        n4 = max_(quotient_familial.not41 * (1 * caseP + 1 * caseF), quotient_familial.not42 * (caseW | caseS))
+        n4 = max_(quotient_familial.couple_ou_pers_a_charge.not41 * (1 * caseP + 1 * caseF), quotient_familial.couple_ou_pers_a_charge.not42 * (caseW | caseS))
 
         # # note 5
         #  - enfant autre et parent isolé
-        n5 = quotient_familial.isol * caseT * (((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2))) + 1 * has_pac)
+        n5 = quotient_familial.couple_ou_pers_a_charge.isol * caseT * (((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2))) + 1 * has_pac)
 
         # # note 6 invalide avec personne à charge
-        n6 = quotient_familial.not6 * (caseP & (has_pac | has_alt))
+        n6 = quotient_familial.couple_ou_pers_a_charge.not6 * (caseP & (has_pac | has_alt))
 
         # # note 7 Parent isolé
-        n7 = quotient_familial.isol * caseT * ((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2)) + 1 * has_pac)
+        n7 = quotient_familial.couple_ou_pers_a_charge.isol * caseT * ((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2)) + 1 * has_pac)
 
         # # Régime des mariés ou pacsés
-        nb_parts_famille = 1 + quotient_familial.conj + enf + n2 + n4
+        nb_parts_famille = 1 + quotient_familial.cas_general.conj + enf + n2 + n4
 
         # # veufs  hors jeune_veuf
-        nb_parts_veuf = 1 + quotient_familial.veuf * (has_pac | has_alt) + enf + n2 + n3 + n5 + n6
+        nb_parts_veuf = 1 + quotient_familial.cas_general.veuf * (has_pac | has_alt) + enf + n2 + n3 + n5 + n6
 
         # # celib div
         nb_parts_celib = 1 + enf + n2 + n3 + n6 + n7
@@ -3699,18 +3755,18 @@ class nbptr(Variable):
 
         note 1 enfants et résidence alternée (formulaire 2041 GV page 10)
 
-        quotient_familial.conj : nb part associées au conjoint d'un couple marié ou pacsé
-        quotient_familial.enf1 : nb part premier enfant
-        quotient_familial.enf2 : nb part deuxième enfant
-        quotient_familial.enf3_et_sup : nb part enfants de rang 3 ou plus
-        quotient_familial.inv1 : nb part supp enfants invalides (I, G)
-        quotient_familial.inv2 : nb part supp adultes invalides (R)
+        quotient_familial.cas_general.conj : nb part associées au conjoint d'un couple marié ou pacsé
+        quotient_familial.cas_general.enf1 : nb part premier enfant
+        quotient_familial.cas_general.enf2 : nb part deuxième enfant
+        quotient_familial.cas_general.enf3_et_sup : nb part enfants de rang 3 ou plus
+        quotient_familial.couple_ou_pers_a_charge.inv1 : nb part supp enfants invalides (I, G)
+        quotient_familial.couple_ou_pers_a_charge.inv2 : nb part supp adultes invalides (R)
         quotient_familial.sans_pers_a_charge.not31 : nb part supp note 3 : cases W ou G pour veuf, celib ou div
         quotient_familial.sans_pers_a_charge.not32 : nb part supp note 3 : personne seule ayant élevé des enfants
-        quotient_familial.not41 : nb part supp adultes invalides (vous et/ou conjoint) note 4
-        quotient_familial.not42 : nb part supp adultes anciens combattants (vous et/ou conjoint) note 4
-        quotient_familial.not6 : nb part supp note 6
-        quotient_familial.isol : demi-part parent isolé (T)
+        quotient_familial.couple_ou_pers_a_charge.not41 : nb part supp adultes invalides (vous et/ou conjoint) note 4
+        quotient_familial.couple_ou_pers_a_charge.not42 : nb part supp adultes anciens combattants (vous et/ou conjoint) note 4
+        quotient_familial.couple_ou_pers_a_charge.not6 : nb part supp note 6
+        quotient_familial.couple_ou_pers_a_charge.isol : demi-part parent isolé (T)
         quotient_familial.edcd : enfant issu du mariage avec conjoint décédé;
         '''
         nb_pac = foyer_fiscal('nb_pac', period)
@@ -3746,29 +3802,29 @@ class nbptr(Variable):
 
         # parts des enfants en résidence alternée quand il n'y a que des enfants en résidence alternée
         enf_a = (no_pac & has_alt) * (
-            quotient_familial.enf1 * min_(nbH, 1)
-            + quotient_familial.enf2 * max_(min_(nbH - 1, 1), 0)
-            + quotient_familial.enf3_et_sup * max_(nbH - 2, 0)
+            quotient_familial.cas_general.enf1 * min_(nbH, 1)
+            + quotient_familial.cas_general.enf2 * max_(min_(nbH - 1, 1), 0)
+            + quotient_familial.cas_general.enf3_et_sup * max_(nbH - 2, 0)
             ) * 0.5
         # parts des enfants en résidence alternée quand il y a aussi des enfants entièrement à charge
         enf_b = (has_pac & has_alt) * (
             (nb_pac == 1) * (
-                quotient_familial.enf2 * min_(nbH, 1)
-                + quotient_familial.enf3_et_sup * max_(nbH - 1, 0)
+                quotient_familial.cas_general.enf2 * min_(nbH, 1)
+                + quotient_familial.cas_general.enf3_et_sup * max_(nbH - 1, 0)
                 ) * 0.5
-            + (nb_pac > 1) * (quotient_familial.enf3_et_sup * nbH * 0.5)
+            + (nb_pac > 1) * (quotient_familial.cas_general.enf3_et_sup * nbH * 0.5)
             )
         # parts des enfants entièrement à charge
         enf_c = (
-            quotient_familial.enf1 * min_(nb_pac, 1)
-            + quotient_familial.enf2 * max_(min_(nb_pac - 1, 1), 0)
-            + quotient_familial.enf3_et_sup * max_(nb_pac - 2, 0)
+            quotient_familial.cas_general.enf1 * min_(nb_pac, 1)
+            + quotient_familial.cas_general.enf2 * max_(min_(nb_pac - 1, 1), 0)
+            + quotient_familial.cas_general.enf3_et_sup * max_(nb_pac - 2, 0)
             )
 
         enf = enf_a + enf_b + enf_c
 
         # # note 2 : nombre de parts liées aux invalides (enfant + adulte)
-        n2 = quotient_familial.inv1 * (nbG + nbI / 2) + quotient_familial.inv2 * nbR
+        n2 = quotient_familial.couple_ou_pers_a_charge.inv1 * (nbG + nbI / 2) + quotient_familial.couple_ou_pers_a_charge.inv2 * nbR
 
         # # note 3 : Pas de personne à charge
         # - invalide
@@ -3782,23 +3838,23 @@ class nbptr(Variable):
         n3 = max_(n31, n32)
         # # note 4 Invalidité de la personne ou du conjoint pour les mariés ou
         # # jeunes veuf(ve)s
-        n4 = max_(quotient_familial.not41 * (1 * caseP + 1 * caseF), quotient_familial.not42 * (caseW | caseS))
+        n4 = max_(quotient_familial.couple_ou_pers_a_charge.not41 * (1 * caseP + 1 * caseF), quotient_familial.couple_ou_pers_a_charge.not42 * (caseW | caseS))
 
         # # note 5
         #  - enfant du conjoint décédé
-        n51 = quotient_familial.cdcd * (caseL & ((nbF + nbJ) > 0))
+        n51 = quotient_familial.couple_ou_pers_a_charge.cdcd * (caseL & ((nbF + nbJ) > 0))
         #  - enfant autre et parent isolé
-        n52 = quotient_familial.isol * caseT * (((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2))) + 1 * has_pac)
+        n52 = quotient_familial.couple_ou_pers_a_charge.isol * caseT * (((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2))) + 1 * has_pac)
         n5 = max_(n51, n52)
 
         # # note 6 invalide avec personne à charge
-        n6 = quotient_familial.not6 * (caseP & (has_pac | has_alt))
+        n6 = quotient_familial.couple_ou_pers_a_charge.not6 * (caseP & (has_pac | has_alt))
 
         # # note 7 Parent isolé
-        n7 = quotient_familial.isol * caseT * ((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2)) + 1 * has_pac)
+        n7 = quotient_familial.couple_ou_pers_a_charge.isol * caseT * ((no_pac & has_alt) * ((nbH == 1) * 0.5 + (nbH >= 2)) + 1 * has_pac)
 
         # # Régime des mariés ou pacsés
-        nb_parts_famille = 1 + quotient_familial.conj + enf + n2 + n4
+        nb_parts_famille = 1 + quotient_familial.cas_general.conj + enf + n2 + n4
 
         # # veufs  hors jeune_veuf
         nb_parts_veuf = 1 + enf + n2 + n3 + n5 + n6