diff --git a/src/Computation.cxx b/src/Computation.cxx
index 7896092..dc05e57 100644
--- a/src/Computation.cxx
+++ b/src/Computation.cxx
@@ -375,36 +375,6 @@ std::vector<double> Computation::computeAugmentedPerturbationSaturatedAndDissipa
     }
 }
 
-//TODO test if it works in the correct way
-std::vector<double> Computation::computeAugmentedPerturbationEnhanced1(double timeStep, bool saturation, const std::vector<double>& saturationsVector){
-    if(saturation){
-        if (saturationsVector.size() ) {
-            if (saturationsVector.size() == InputAugmentedArma.n_elem) {
-                arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-                for(uint i = 0;i<outputArma.n_elem;i++){
-                    outputArma[i] = hyperbolicTangentScaled(outputArma[i], saturationsVector[i]);
-                }
-                outputAugmented = armaColumnToVector(outputArma);
-                return outputAugmented;
-            } else{
-                throw std::invalid_argument("[ERROR] Computation::computeAugmentedPerturbationEnhanced1: saturationVector is not of the same size as output vector. abort");
-            }
-        }
-        else {
-            arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-            for(uint i = 0;i<outputArma.n_elem;i++){
-                outputArma[i] = hyperbolicTangentScaled(outputArma[i], 1);
-            }
-            outputAugmented = armaColumnToVector(outputArma);
-            return outputAugmented;
-        }
-    } else {
-        arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-        outputAugmented = armaColumnToVector(outputArma);
-        return outputAugmented;
-    }
-}
-
 //TODO test if it works in the correct way
 std::vector<double> Computation::computeAugmentedPerturbationEnhanced2(double timeStep, bool saturation, const std::vector<double>& saturationsVector,const std::vector<double>& qVector){
     if (saturation) {
diff --git a/src/Computation.h b/src/Computation.h
index 2660d18..701ba5f 100644
--- a/src/Computation.h
+++ b/src/Computation.h
@@ -69,7 +69,6 @@ class Computation{
         std::vector<double> computeAugmentedPerturbationDissipatedAfterCompute(double timeStep); //taking into account dissipation after every iteration(Dissipation model), dissipation after the computation of the perturbated value
         std::vector<double> computeAugmentedPerturbationDissipatedBeforeCompute(double timeStep); //taking into account dissipation after every iteration (Dissipation model), dissipation before the computation of the perturbated value
         std::vector<double> computeAugmentedPerturbationSaturatedAndDissipatedBeforeCompute(double timeStep,const std::vector<double>& saturationsVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling) and dissipation after every iteration
-        std::vector<double> computeAugmentedPerturbationEnhanced1(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling) and dissipation after every iteration
         std::vector<double> computeAugmentedPerturbationEnhanced2(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>(),const std::vector<double>& qVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling), dissipation and conservation after every iteration
         std::vector<double> computeAugmentedPerturbationEnhanced3(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>(),const std::vector<double>& qVector = std::vector<double>(), std::function<double(double)> propagationScaleFunction = [](double time)-> double{return 1.0;}); //taking into account scaling
         std::vector<double> computeAugmentedPerturbationEnhanced4(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>(),const std::vector<double>& qVector = std::vector<double>()); //all the models
diff --git a/src/ComputationOptimized.cxx b/src/ComputationOptimized.cxx
index 9166acc..a515ef4 100644
--- a/src/ComputationOptimized.cxx
+++ b/src/ComputationOptimized.cxx
@@ -315,35 +315,6 @@ std::vector<double> ComputationOptimized::computeAugmentedPerturbationSaturatedA
     }
 }
 
-//TODO test if it works in the correct way
-std::vector<double> ComputationOptimized::computeAugmentedPerturbationEnhanced1(double timeStep, bool saturation, const std::vector<double>& saturationsVector){
-    if(saturation){
-        if (saturationsVector.size() ) {
-            if (saturationsVector.size() == InputAugmentedArma.n_elem) {
-                arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-                for(uint i = 0;i<outputArma.n_elem;i++){
-                    outputArma[i] = hyperbolicTangentScaled(outputArma[i], saturationsVector[i]);
-                }
-                outputAugmented = armaColumnToVector(outputArma);
-                return outputAugmented;
-            } else{
-                throw std::invalid_argument("[ERROR] ComputationOptimized::computeAugmentedPerturbationEnhanced1: saturationVector is not of the same size as output vector. abort");
-            }
-        }
-        else {
-            arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-            for(uint i = 0;i<outputArma.n_elem;i++){
-                outputArma[i] = hyperbolicTangentScaled(outputArma[i], 1);
-            }
-            outputAugmented = armaColumnToVector(outputArma);
-            return outputAugmented;
-        }
-    } else {
-        arma::Col<double> outputArma =  pseudoInverseAugmentedArma * dissipationModel->dissipate(InputAugmentedArma, timeStep);
-        outputAugmented = armaColumnToVector(outputArma);
-        return outputAugmented;
-    }
-}
 
 //TODO test if it works in the correct way
 std::vector<double> ComputationOptimized::computeAugmentedPerturbationEnhanced2(double timeStep, bool saturation, const std::vector<double>& saturationsVector,const std::vector<double>& qVector){
diff --git a/src/ComputationOptimized.h b/src/ComputationOptimized.h
index 3357a1f..063a3aa 100644
--- a/src/ComputationOptimized.h
+++ b/src/ComputationOptimized.h
@@ -66,7 +66,6 @@ class ComputationOptimized{
         std::vector<double> computeAugmentedPerturbationDissipatedAfterCompute(double timeStep); //taking into account dissipation after every iteration(Dissipation model), dissipation after the ComputationOptimized of the perturbated value
         std::vector<double> computeAugmentedPerturbationDissipatedBeforeCompute(double timeStep); //taking into account dissipation after every iteration (Dissipation model), dissipation before the ComputationOptimized of the perturbated value
         std::vector<double> computeAugmentedPerturbationSaturatedAndDissipatedBeforeCompute(double timeStep,const std::vector<double>& saturationsVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling) and dissipation after every iteration
-        std::vector<double> computeAugmentedPerturbationEnhanced1(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling) and dissipation after every iteration
         std::vector<double> computeAugmentedPerturbationEnhanced2(double timeStep, bool saturation = true, const std::vector<double>& saturationsVector = std::vector<double>(),const std::vector<double>& qVector = std::vector<double>()); //taking into account saturation(hyperbolic tangent and scaling), dissipation and conservation after every iteration
         std::pair<std::string,double> getMapVirtualOutputsToCellInputs(); //TODO
         void updateInput(const std::vector<double>& newInp = std::vector<double>(), bool augmented = false);