diff --git a/srcs/FEM/solverFEM.cpp b/srcs/FEM/solverFEM.cpp
index 2449aa9d805320566ed2f217d85e529db82683a8..088dc6aa5bec43af8f5c496b7e69c91e4378d368 100644
--- a/srcs/FEM/solverFEM.cpp
+++ b/srcs/FEM/solverFEM.cpp
@@ -250,7 +250,7 @@ void solverFEMnonLinear(std::map<int, double> &electrostaticPressure,
                 relativeForces = sqrt(relativeForces/(2*FEM_nodeTags.size())/currentMaxNodalForce); // norm made relative
             else //to avoid division by zero
                 relativeForces = 1e17; 
-            std::cout << "step: " << step << ", relativeForces: " << relativeForces << "\n";
+            //std::cout << "step: " << step << ", relativeForces: " << relativeForces << "\n";
 
             if(previousRelativeForces < relativeForces) //the current step induces an increase of nodal difference
             {
diff --git a/srcs/MEMS.geo b/srcs/MEMS.geo
index a918f61c3068a787a1b17244848dfd096e19f79e..1caef0d5d166b5ad1b61ec8fcf77a4c51d5462cb 100644
--- a/srcs/MEMS.geo
+++ b/srcs/MEMS.geo
@@ -9,7 +9,7 @@ Lx_diaph = 2*scale;
 nx = 20;
 ny = 50;
 
-// monter à 500+ steps en FEM non linéaire car converge très très lentement
+// THIS .geo FILE is mentioned once IN THE REPORT, in section 5.3.
 
 Point(1) = {0, 0, 0, 2};
 Point(2) = {Lx_foot, 0, 0, 2};
@@ -90,7 +90,7 @@ SetNumber("Materials/FEM_domain/Young", 210e3);
 SetNumber("Materials/FEM_domain/Poisson", 0.3);
 SetNumber("Materials/FEM_domain/rho",7800);
 
-SetNumber("Boundary Conditions/top_edge/tx", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
+SetNumber("Boundary Conditions/top_edge/tx", 0.);
 SetNumber("Boundary Conditions/top_edge/ty", 0.);
 SetNumber("Volumic Forces/FEM_domain/bx",0.);
 SetNumber("Volumic Forces/FEM_domain/by",0.);
diff --git a/srcs/beamActuation.geo b/srcs/beamActuation.geo
index a37ee6d35de807819cba90a988a91c0491262946..fbec4f3386b33cc611be8ea64f4faf12c9e2671e 100644
--- a/srcs/beamActuation.geo
+++ b/srcs/beamActuation.geo
@@ -1,5 +1,7 @@
 scale = 1e-6;
 
+// THIS .geo FILE IS NOT USED IN THE REPORT
+
 bt = 1; // Beam tickness
 bh = 1; // Beam height (position)
 
diff --git a/srcs/beamActuationGeneral.geo b/srcs/beamActuationGeneral.geo
deleted file mode 100644
index 3ad6cf8a80d55447284b1467115cedc983abd0e0..0000000000000000000000000000000000000000
--- a/srcs/beamActuationGeneral.geo
+++ /dev/null
@@ -1,229 +0,0 @@
-scale = 1e-6;
-
-bt = 0.2; // Beam tickness
-bh = 1; // Beam height (position)
-
-eL = 3; // Electrode base length
-el = 1.5; // Electrode tip length
-et = 1; // Electrode tickness
-es = 1.5; // Electrode spacing
-
-nbElectrodes = 10;
-
-W = nbElectrodes*(es + eL) + 2*es;    // Width of the space
-H = W;  // Height of the space
-
-nEs = 12;
-nEL = 15;
-nBt = 8;
-nEt = 3;
-nBh = 3;
-
-bt = bt * scale;
-bh = bh * scale;
-eL = eL * scale;
-el = el * scale;
-et = et * scale;
-es = es * scale;
-W = W * scale;
-H = H * scale;
-
-e = (eL - el)/2;    // Electrode constant
-
-/*
- *  Points
- */
-
-For i In {0:nbElectrodes-1}
-    Point(6*i + 1) = {(es+eL)*i, bh + bt, 0, 1};
-    Point(6*i + 2) = {(es+eL)*i + es, bh + bt, 0, 1};
-    Point(6*i + 3) = {(es+eL)*i + es+e, bh + bt + et, 0, 1};
-    Point(6*i + 4) = {(es+eL)*i + es+e+el, bh + bt + et, 0, 1};
-    Point(6*i + 5) = {(es+eL)*i + es, bh, 0, 1};
-    Point(6*i + 6) = {(es+eL)*i, bh, 0, 1};
-EndFor
-Point(6*nbElectrodes + 1) = {nbElectrodes*(es+eL), bh + bt, 0, 1};
-Point(6*nbElectrodes + 2) = {nbElectrodes*(es+eL) + es, bh + bt, 0, 1};
-Point(6*nbElectrodes + 3) = {nbElectrodes*(es+eL) + es, bh, 0, 1};
-Point(6*nbElectrodes + 6) = {nbElectrodes*(es+eL), bh, 0, 1};
-
-Point(6*nbElectrodes + 4) = {0, 0, 0, 1};
-Point(6*nbElectrodes + 5) = {W, 0, 0, 1};
-Point(6*nbElectrodes + 7) = {W, H, 0, 1};
-Point(6*nbElectrodes + 8) = {0, H, 0, 1};
-
-For i In {0:nbElectrodes-1}
-    Point(6*nbElectrodes + 8 + i + 1) = {i*(es+eL) + es + e + el/2, bh + bt + et - ((el*e)/(2*et)), 0, 1};
-EndFor
-
-/*
- *  Lines
- */
-Line(9*nbElectrodes + 1) = {6*nbElectrodes + 4, 6*nbElectrodes + 5};
-Line(9*nbElectrodes + 2) = {6*nbElectrodes + 5, 6*nbElectrodes + 7};
-Line(9*nbElectrodes + 3) = {6*nbElectrodes + 7, 6*nbElectrodes + 8};
-Line(9*nbElectrodes + 4) = {6*nbElectrodes + 8, 1};
-Line(9*nbElectrodes + 5) = {6, 6*nbElectrodes + 4};
-
-Line(9*nbElectrodes + 9) = {6*nbElectrodes + 1, 6*nbElectrodes + 2};
-Line(9*nbElectrodes + 6) = {6*nbElectrodes + 2, 6*nbElectrodes + 3};
-Line(9*nbElectrodes + 7) = {6*nbElectrodes + 3, 6*nbElectrodes + 6};
-Line(9*nbElectrodes + 8) = {6*nbElectrodes + 6, 6*nbElectrodes + 1};
-For i In {0:nbElectrodes-1}
-    Line(9*i + 1) = {6*i + 1, 6*i + 2};
-    Line(9*i + 2) = {6*i + 2, 6*i + 3};
-    // Line(9*i + 3) = {6*i + 3, 6*i + 4};
-    Circle(9*i + 3) = {6*i + 3, 6*nbElectrodes + 8 + i + 1, 6*i + 4};
-    Line(9*i + 4) = {6*i + 4, 6*(i+1) + 1};
-    Line(9*i + 5) = {6*(i+1) + 1, 6*i + 2};
-    Line(9*i + 6) = {6*i + 5, 6*i + 2};
-    Line(9*i + 7) = {6*i + 5, 6*i + 6};
-    Line(9*i + 8) = {6*i + 6, 6*i + 1};
-    Line(9*i + 9) = {6*(i+1) + 6, 6*i + 5};
-EndFor
-
-/*
- *  Transfinite Curve
- */
-For i In {0:nbElectrodes-1}
-    Transfinite Curve{9*i + 1, 9*i + 7} = nEs + 1 Using Progression 1;
-    Transfinite Curve{9*i + 2, 9*i + 4} = nEt + 1 Using Progression 1;
-    Transfinite Curve{9*i + 3, 9*i + 5, 9*i + 9} = nEL + 1 Using Progression 1;
-    Transfinite Curve{9*i + 8, 9*i + 6} = nBt + 1 Using Progression 1;
-EndFor
-Transfinite Curve{9*nbElectrodes + 9, 9*nbElectrodes + 7} = nEs + 1 Using Progression 1;
-Transfinite Curve{9*nbElectrodes + 8, 9*nbElectrodes + 6} = nBt + 1 Using Progression 1;
-
-Transfinite Curve{9*nbElectrodes + 1, 9*nbElectrodes + 3} = nbElectrodes*(nEs + nEL) + 2*nEs + 1 Using Progression 1;
-Transfinite Curve{9*nbElectrodes + 2} = nbElectrodes*(nEs + nEL) + 2*nEs + 1 Using Progression 1;
-Transfinite Curve{9*nbElectrodes + 4} = nbElectrodes*(nEs + nEL) + 2*nEs - nBt - nBh + 1 Using Progression 1;
-Transfinite Curve{9*nbElectrodes + 5} = nBh + 1 Using Progression 1;
-
-/*
- *  Curve Loop
- */
-
-For i In {0:nbElectrodes-1}
-    Curve Loop(3*i + 1) = {9*i + 1, -(9*i + 6), 9*i + 7, 9*i + 8};
-    Curve Loop(3*i + 2) = {9*i + 3, 9*i + 4, 9*i + 5, 9*i + 2};
-    Curve Loop(3*i + 3) = {-(9*i + 5), -(9*(i+1) + 8), 9*i + 9, 9*i + 6};
-EndFor
-Curve Loop(3*nbElectrodes + 1) = {9*nbElectrodes + 9, 9*nbElectrodes + 6, 9*nbElectrodes + 7, 9*nbElectrodes + 8};
-// How to fill a vector in a loop ?
-
-bemBoundary = {1:6*nbElectrodes + 8};
-bemBoundaryWithoutE = {1:4*nbElectrodes + 8 + 2};
-bemFemBoundary = {1:6*nbElectrodes + 3};
-leftElectrodesFull = {1:nbElectrodes};
-rightElectrodesFull = {1:nbElectrodes};
-For i In {0:nbElectrodes-1}
-    bemBoundary[4*i+1 - 1] = 9*i + 1;
-    bemBoundary[4*i+2 - 1] = 9*i + 2;
-    bemBoundary[4*i+3 - 1] = 9*i + 3;
-    bemBoundary[4*i+4 - 1] = 9*i + 4;
-
-    bemBoundaryWithoutE[2*i+1 - 1] = 9*i + 1;
-    bemBoundaryWithoutE[2*i+2 - 1] = 9*i + 3;
-
-    bemFemBoundary[4*i+1 - 1] = 9*i + 1;
-    bemFemBoundary[4*i+2 - 1] = 9*i + 2;
-    bemFemBoundary[4*i+3 - 1] = 9*i + 3;
-    bemFemBoundary[4*i+4 - 1] = 9*i + 4;
-
-    leftElectrodesFull[i] = 9*i + 2;
-    rightElectrodesFull[i] = 9*i + 4;
-EndFor
-bemBoundary[4*nbElectrodes + 1 - 1] = 9*nbElectrodes + 9;
-bemBoundary[4*nbElectrodes + 2 - 1] = 9*nbElectrodes + 6;
-bemBoundary[4*nbElectrodes + 3 - 1] = 9*nbElectrodes + 7;
-
-bemBoundaryWithoutE[2*nbElectrodes + 1 - 1] = 9*nbElectrodes + 9;
-bemBoundaryWithoutE[2*nbElectrodes + 2 - 1] = 9*nbElectrodes + 6;
-bemBoundaryWithoutE[2*nbElectrodes + 3 - 1] = 9*nbElectrodes + 7;
-
-bemFemBoundary[4*nbElectrodes + 1 - 1] = 9*nbElectrodes + 9;
-bemFemBoundary[4*nbElectrodes + 2 - 1] = 9*nbElectrodes + 6;
-bemFemBoundary[4*nbElectrodes + 3 - 1] = 9*nbElectrodes + 7;
-For i In {0:nbElectrodes-1}
-    bemBoundary[4*nbElectrodes + 3 + 2*i+1 - 1] = 9*i + 9;
-    bemBoundary[4*nbElectrodes + 3 + 2*i+2 - 1] = 9*i + 7;
-
-    bemBoundaryWithoutE[2*nbElectrodes + 3 + 2*i+1 - 1] = 9*i + 9;
-    bemBoundaryWithoutE[2*nbElectrodes + 3 + 2*i+2 - 1] = 9*i + 7;
-
-    bemFemBoundary[4*nbElectrodes + 3 + 2*i+1 - 1] = 9*i + 9;
-    bemFemBoundary[4*nbElectrodes + 3 + 2*i+2 - 1] = 9*i + 7;
-EndFor
-bemBoundary[6*nbElectrodes + 4 - 1] = 9*nbElectrodes + 5;
-bemBoundary[6*nbElectrodes + 5 - 1] = 9*nbElectrodes + 1;
-bemBoundary[6*nbElectrodes + 6 - 1] = 9*nbElectrodes + 2;
-bemBoundary[6*nbElectrodes + 7 - 1] = 9*nbElectrodes + 3;
-bemBoundary[6*nbElectrodes + 8 - 1] = 9*nbElectrodes + 4;
-
-bemBoundaryWithoutE[4*nbElectrodes + 4 - 1] = 9*nbElectrodes + 5;
-bemBoundaryWithoutE[4*nbElectrodes + 5 - 1] = 9*nbElectrodes + 1;
-bemBoundaryWithoutE[4*nbElectrodes + 6 - 1] = 9*nbElectrodes + 2;
-bemBoundaryWithoutE[4*nbElectrodes + 7 - 1] = 9*nbElectrodes + 3;
-bemBoundaryWithoutE[4*nbElectrodes + 8 - 1] = 9*nbElectrodes + 4;
-
-leftElectrodes = {1:nbElectrodes-1};
-rightElectrodes = {1:nbElectrodes-1};
-For i In {0:nbElectrodes-2}
-    leftElectrodes[i] = leftElectrodesFull[i+1];
-    rightElectrodes[i] = rightElectrodesFull[i];
-EndFor
-
-bemBoundaryWithoutE[4*nbElectrodes + 9 - 1] = 2;
-bemBoundaryWithoutE[4*nbElectrodes + 10 - 1] = 9*nbElectrodes + 4;
-
-Curve Loop(3*nbElectrodes + 2) = bemBoundary[];
-
-/*
- *  Plane Surface
- */
-For i In {0:nbElectrodes-1}
-    Plane Surface(3*i + 1) = {3*i + 1};
-    Plane Surface(3*i + 2) = {3*i + 2};
-    Plane Surface(3*i + 3) = {3*i + 3};
-EndFor
-Plane Surface(3*nbElectrodes + 1) = {3*nbElectrodes + 1};
-Plane Surface(3*nbElectrodes + 2) = {3*nbElectrodes + 2};
-
-/*
- *  Surface operations
- */
-Recombine Surface{1:3*nbElectrodes + 1};
-Transfinite Surface{1:3*nbElectrodes + 1};
-
-/*
- *  Physical Curve
- */
-Physical Curve("clamp", 1) = {8};
-Physical Curve("BEM_boundary_without_e", 2) = bemBoundaryWithoutE[];
-Physical Curve("BEM_FEM_boundary", 3) = bemFemBoundary[];
-Physical Curve("left_electrodes", 4) = leftElectrodes[];
-Physical Curve("right_electrodes", 5) = rightElectrodes[];
-
-Physical Surface("FEM_domain") = {1:3*nbElectrodes + 1};
-Physical Surface("BEM_domain_1") = {3*nbElectrodes + 2};
-
-/*
- *  Parameters
- */
- SetNumber("Boundary Conditions/clamp/ux", 0.0);
- SetNumber("Boundary Conditions/clamp/uy", 0.0);
- 
-//  SetNumber("Materials/FEM_domain/Young", 3e7);
-SetNumber("Materials/FEM_domain/Young", 1e7);
- SetNumber("Materials/FEM_domain/Poisson", 0.3);
- SetNumber("Materials/FEM_domain/rho",7800);
- 
- SetNumber("Volumic Forces/FEM_domain/bx", 0.0);
- SetNumber("Volumic Forces/FEM_domain/by", 0.0);
- 
- SetNumber("Boundary Conditions/BEM_boundary_without_e/BEM_domain_1/neumann", 0);
- SetNumber("Boundary Conditions/left_electrodes/BEM_domain_1/dirichlet", 8);
- SetNumber("Boundary Conditions/right_electrodes/BEM_domain_1/dirichlet", 0);
- SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
- 
- SetNumber("Non_linear_solver", 1);
\ No newline at end of file
diff --git a/srcs/clampedMicroBeam.geo b/srcs/clampedMicroBeam.geo
index 5e94f6d3570ef5bd1d8605fdd2833d07b4ee2f59..8b19c2193394b977c8f093c4bd56cf441d51a9cc 100644
--- a/srcs/clampedMicroBeam.geo
+++ b/srcs/clampedMicroBeam.geo
@@ -10,13 +10,13 @@ n = 6; // FEM MESH DENSITY
 // additional parameters given to the solver
 SetNumber("Boundary Conditions/left_edge/ux", 0.); 
 SetNumber("Boundary Conditions/left_edge/uy", 0);
-SetNumber("Materials/domain/Young", 150e9);
-SetNumber("Materials/domain/Poisson", 0.27);
-SetNumber("Materials/domain/rho",2300);
+SetNumber("Materials/FEM_domain/Young", 150e9);
+SetNumber("Materials/FEM_domain/Poisson", 0.27);
+SetNumber("Materials/FEM_domain/rho",2300);
 SetNumber("Volumic Forces/FEM_domain/bx",0.);
-SetNumber("Volumic Forces/FEM_domain/by",0.); //set to -9.81 for gravity
+SetNumber("Volumic Forces/FEM_domain/by",0.);
 
-phi_top = 112; // à modifier valou
+phi_top = 111.2; 
 SetNumber("Boundary Conditions/mass/BEM_domain_1/dirichlet", 0);
 SetNumber("Boundary Conditions/BEM_FEM_boundary/BEM_domain_1/dirichlet", phi_top);
 SetNumber("Boundary Conditions/rest_of_outside/BEM_domain_1/neumann", 0);
diff --git a/srcs/coupling_validation.geo b/srcs/coupling_validation.geo
index e5510d1e2ee35de5eebe1b186d0f0be9839fefe1..78440018e8001ebf4a75ebb44cf8abcf605323eb 100644
--- a/srcs/coupling_validation.geo
+++ b/srcs/coupling_validation.geo
@@ -44,19 +44,19 @@ Physical Surface("FEM_domain", 2) = {1};
 Physical Point("fixed_node", 3) = {1};
 
 // additional parameters given to the solver
-SetNumber("Boundary Conditions/left_edge/ux", 0.); // HERE YOU DO NOT HAVE TO IMPOSE BOTH ux and uy simultaneously ! (permet aussi de simuler appuis à roulettes)
+SetNumber("Boundary Conditions/left_edge/ux", 0.); 
 //SetNumber("Boundary Conditions/left_edge/uy", 0.);
 SetNumber("Boundary Conditions/fixed_node/uy", 0.);
 SetNumber("Materials/FEM_domain/Young", 210e3);
 SetNumber("Materials/FEM_domain/Poisson", 0.3);
-SetNumber("Materials/FEM_domain/rho",7800); //volumic mass of acier
+SetNumber("Materials/FEM_domain/rho",7800); 
 
 Physical Curve("BEM_FEM_boundary", 4) = {2};
 Physical Surface("BEM_domain_1", 5) = {2};
 Physical Curve("right_BEM", 6) = {6};
 Physical Curve("homogeneous_field", 7) = {5, 7};
 
-phi = 30;
+phi = 90;
 SetNumber("Boundary Conditions/right_BEM/BEM_domain_1/dirichlet", 0);
 SetNumber("Boundary Conditions/BEM_FEM_boundary/BEM_domain_1/dirichlet", phi);
 SetNumber("Boundary Conditions/homogeneous_field/BEM_domain_1/neumann", 0);
diff --git a/srcs/hybrid_geo.geo b/srcs/hybrid_geo.geo
index 53ba5b7907d95934a50e77a778638bc84538074a..ef858c7410c7b7db6318730c3a5832abe20d0f95 100644
--- a/srcs/hybrid_geo.geo
+++ b/srcs/hybrid_geo.geo
@@ -4,9 +4,11 @@ Ly_poutre = 2*scale;
 h_poutre = 2*scale;
 h_tot = 10*scale;
 width = 10*scale;
-nx = 20; // prend beaucoup de temps àpd de 200x40
+nx = 20;
 ny = 8;
 
+// THIS .geo FILE IS NOT USED IN THE REPORT
+
 Point(1) = {0, h_poutre, 0, 2};
 Point(2) = {Lx_poutre, h_poutre, 0, 2};
 Point(3) = {Lx_poutre, h_poutre + Ly_poutre, 0, 2};
@@ -59,17 +61,17 @@ Physical Point("fixed_node", 9) = {1};
 Physical Surface("BEM_domain_1", 10) = {2};
 
 // additional parameters given to the solver
-SetNumber("Boundary Conditions/left_edge/ux", 0.); // HERE YOU DO NOT HAVE TO IMPOSE BOTH ux and uy simultaneously ! (permet aussi de simuler appuis à roulettes)
+SetNumber("Boundary Conditions/left_edge/ux", 0.);
 SetNumber("Boundary Conditions/left_edge/uy", 0);
 SetNumber("Materials/domain/Young", 210e3);
 SetNumber("Materials/domain/Poisson", 0.3);
-SetNumber("Materials/domain/rho",7800); //volumic mass of acier
-SetNumber("Boundary Conditions/top_edge/tx", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
-SetNumber("Boundary Conditions/top_edge/ty", 0.); //set to some non-zero value to induce vertical deflection
-//SetNumber("Boundary Conditions/right_edge/tx", 21e3); // for simple tension conditions
+SetNumber("Materials/domain/rho",7800);
+SetNumber("Boundary Conditions/top_edge/tx", 0.); 
+SetNumber("Boundary Conditions/top_edge/ty", 0.); 
+//SetNumber("Boundary Conditions/right_edge/tx", 21e3);
 SetNumber("Boundary Conditions/right_edge/ty", 0.);
 SetNumber("Volumic Forces/FEM_domain/bx",0.);
-SetNumber("Volumic Forces/FEM_domain/by",0.); //set to -9.81 for gravity
+SetNumber("Volumic Forces/FEM_domain/by",0.);
 //SetNumber("Boundary Conditions/fixed_node/uy",2.);
 
 // BEM geometry
@@ -81,8 +83,5 @@ SetNumber("Boundary Conditions/mass/BEM_domain_1/dirichlet", 0);
 SetNumber("Boundary Conditions/BEM_FEM_boundary/BEM_domain_1/dirichlet", phi_top);
 SetNumber("Boundary Conditions/rest_of_outside/BEM_domain_1/neumann", 0);
 SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
-//Physical Curve("BEM_boundary", 14) = {1,3,5,6,7,8,9,10,2};
 
-Physical Curve("Electrode", 15) = {1};
-//SetNumber("Boundary Conditions/bottom_edge/tx", 0);
-//SetNumber("Boundary Conditions/bottom_edge/ty", -25000);
\ No newline at end of file
+Physical Curve("Electrode", 15) = {1};
\ No newline at end of file
diff --git a/srcs/longitudinalCombDevice.geo b/srcs/longitudinalCombDevice.geo
index 098922a6252baea88986052ff3e46fad2b5af584..8398239152c737c0473d347a3969cd5b81471d3b 100644
--- a/srcs/longitudinalCombDevice.geo
+++ b/srcs/longitudinalCombDevice.geo
@@ -1,13 +1,15 @@
 scale = 2e-6;
 
 // USE WITH MINIMUM 2 FINS, else use longitudinal_comb.geo
-N_fins = 3; // number of fins on one side of the comb // MAX 12 si on change ps la longueur
+N_fins = 4; // number of fins on one side of the comb // MAX 12 si on change ps la longueur
 
 // WARNING: when using more fins the pull-in voltage decreases
 
-n = 1; // FEM elements density
+n = 6; // FEM elements density
 nBEM = 1; // BEM elements density
 
+a = 0; // acceleration verticale
+
 // mechanical properties and boundary conditions
 SetNumber("Boundary Conditions/left/ux", 0.); // encastrement
 SetNumber("Boundary Conditions/left/uy", 0);
@@ -17,10 +19,12 @@ SetNumber("Materials/FEM_domain/Young", 150e9);
 SetNumber("Materials/FEM_domain/Poisson", 0.27);
 SetNumber("Materials/FEM_domain/rho",2300); 
 SetNumber("Volumic Forces/FEM_domain/bx",0);
-SetNumber("Volumic Forces/FEM_domain/by",00); // acceleration of accelerometer
+SetNumber("Volumic Forces/FEM_domain/by",a); // acceleration of accelerometer
+
+SetNumber("Non_linear_solver", 1);
 
 // BEM properties in bottom domain
-phi_1 = 50;
+phi_1 = 20;
 SetNumber("Boundary Conditions/BEM_FEM_boundary_1/BEM_domain_1/dirichlet", 0);
 SetNumber("Boundary Conditions/electrode_1/BEM_domain_1/dirichlet", phi_1);
 SetNumber("Boundary Conditions/outside_1/BEM_domain_1/neumann", 0);
@@ -42,7 +46,7 @@ h_fin_elec = 2.4*scale; // length of the fins of the electrode, can be longer th
 
 //l_bord = 10*scale;
 //l_tot = 34.4*scale;
-l_tot = 40*scale;
+l_tot = 60*scale;
 l_fin = 0.8*scale;
 l_space = 0.8*scale;
 t_electrode = 1*scale; // width of one electrode
diff --git a/srcs/longitudinal_comb.geo b/srcs/longitudinal_comb.geo
deleted file mode 100644
index 6c89eec081168abdcadd2fbd9403a00c051875fd..0000000000000000000000000000000000000000
--- a/srcs/longitudinal_comb.geo
+++ /dev/null
@@ -1,147 +0,0 @@
-scale = 1e-5;
-
-h_tot = 20*scale;
-h_base = 1*scale;
-h_pin = 2.5*scale;
-h_space = 1.2*scale;
-
-l_tot = 20*scale;
-l_pin = 0.8*scale;
-l_space = 0.5*scale;
-
-t_electrode = 1*scale; // width of fixed electrodes
-// FEM domain
-Point(1) = {-l_tot/2, h_base/2, 0, 0.2*scale};
-Point(2) = {-l_pin/2, h_base/2, 0, 0.2*scale};
-Point(3) = {-l_pin/2, h_base/2 + h_pin, 0, 0.2*scale};
-Point(4) = {l_pin/2, h_base/2 + h_pin, 0, 0.2*scale};
-Point(5) = {l_pin/2, h_base/2, 0, 0.2*scale};
-Point(6) = {l_tot/2, h_base/2, 0, 0.2*scale};
-Point(7) = {l_tot/2, -h_base/2, 0, 0.2*scale};
-Point(8) = {l_pin/2, -h_base/2, 0, 0.2*scale};
-Point(9) = {l_pin/2, -(h_base/2 + h_pin), 0, 0.2*scale};
-Point(10) = {-l_pin/2, -(h_base/2 + h_pin), 0, 0.2*scale};
-Point(11) = {-l_pin/2, -h_base/2, 0, 0.2*scale};
-Point(12) = {-l_tot/2, -h_base/2, 0, 0.2*scale};
-
-// FEM domain
-Line(1) = {1, 2};
-Line(2) = {2, 3};
-Line(3) = {3, 4};
-Line(4) = {4, 5};
-Line(5) = {5, 6};
-Line(6) = {6, 7};
-Line(7) = {7, 8};
-Line(8) = {8, 9};
-Line(9) = {9, 10};
-Line(10) = {10, 11};
-Line(11) = {11, 12};
-Line(12) = {12, 1};
-
-Curve Loop(1) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
-Plane Surface(1) = {1};
-
-Recombine Surface {1}; // quads instead of triangles
-
-Physical Curve("left", 1) = {12};
-Physical Curve("right", 2) = {6};
-
-Physical Surface("FEM_domain", 3) = {1};
-Physical Curve("BEM_FEM_boundary", 4) = {1, 2, 3, 4, 5, 7, 8, 9, 10, 11};
-
-// mechanical properties and boundary conditions
-SetNumber("Boundary Conditions/left/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/left/uy", 0);
-SetNumber("Boundary Conditions/right/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/right/uy", 0);
-SetNumber("Materials/domain/Young", 210e3); // A DETERMINER PRECISEMENT
-SetNumber("Materials/domain/Poisson", 0.3);
-SetNumber("Materials/domain/rho",7800); //volumic mass of acier
-SetNumber("Volumic Forces/FEM_domain/bx",0); // acceleration of accelerometer
-SetNumber("Volumic Forces/FEM_domain/by",-200.); 
-
-// BEM domain 1
-Point(13) = {-(l_pin/2 + l_space + t_electrode), -(h_base/2 + h_pin + h_space + t_electrode), 0, 0.2*scale}; // first fixed electrode
-Point(14) = {-(l_pin/2 + l_space + t_electrode), -(h_base/2 + h_space), 0, 0.2*scale};
-Point(15) = {-(l_pin/2 + l_space), -(h_base/2 + h_space), 0, 0.2*scale};
-Point(16) = {-(l_pin/2 + l_space), -(h_base/2 + h_pin + h_space), 0, 0.2*scale};
-Point(17) = {l_pin/2 + l_space, -(h_base/2 + h_pin + h_space), 0, 0.2*scale};
-Point(18) = {l_pin/2 + l_space, -(h_base/2 + h_space), 0, 0.2*scale};
-Point(19) = {l_pin/2 + l_space + t_electrode, -(h_base/2 + h_space), 0, 0.2*scale};
-Point(20) = {l_pin/2 + l_space + t_electrode, -(h_base/2 + h_pin + h_space + t_electrode), 0, 0.2*scale};
-
-Line(13) = {13, 14};
-Line(14) = {14, 15};
-Line(15) = {15, 16};
-Line(16) = {16, 17};
-Line(17) = {17, 18};
-Line(18) = {18, 19};
-Line(19) = {19, 20};
-Line(20) = {20, 13};
-
-Point(29) = {-l_tot/2, -h_tot/2, 0, 0.2*scale};
-Point(30) = {l_tot/2, -h_tot/2, 0, 0.2*scale};
-
-Line(29) = {12, 29};
-Line(30) = {29, 30};
-Line(31) = {30, 7};
-
-Curve Loop(2) = {29, 30, 31, 7, 8, 9, 10, 11};
-Curve Loop(3) = {13, 14, 15, 16, 17, 18, 19, 20};
-Plane Surface(2) = {2, 3};
-
-Physical Surface("BEM_domain_1", 5) = {2};
-Physical Curve("BEM_FEM_boundary_1", 6) = {7, 8, 9, 10, 11};
-Physical Curve("electrode_1", 7) = {13, 14, 15, 16, 17, 18, 19, 20};
-Physical Curve("outside_1", 8) = {29, 30, 31};
-
-phi_1 = 0.1;
-SetNumber("Boundary Conditions/BEM_FEM_boundary_1/BEM_domain_1/dirichlet", 0);
-SetNumber("Boundary Conditions/electrode_1/BEM_domain_1/dirichlet", phi_1);
-SetNumber("Boundary Conditions/outside_1/BEM_domain_1/neumann", 0);
-SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
-
-// BEM domain 2
-Point(21) = {-(l_pin/2 + l_space + t_electrode), h_base/2 + h_space, 0, 0.2*scale}; // second fixed electrode
-Point(22) = {-(l_pin/2 + l_space + t_electrode), h_base/2 + h_pin + h_space + t_electrode, 0, 0.2*scale}; 
-Point(23) = {l_pin/2 + l_space + t_electrode, h_base/2 + h_pin + h_space + t_electrode, 0, 0.2*scale};
-Point(24) = {l_pin/2 + l_space + t_electrode, h_base/2 + h_space, 0, 0.2*scale};
-Point(25) = {l_pin/2 + l_space, h_base/2 + h_space, 0, 0.2*scale};
-Point(26) = {l_pin/2 + l_space, h_base/2 + h_pin + h_space, 0, 0.2*scale};
-Point(27) = {-(l_pin/2 + l_space), h_base/2 + h_pin + h_space, 0, 0.2*scale};
-Point(28) = {-(l_pin/2 + l_space), h_base/2 + h_space, 0, 0.2*scale};
-
-
-Line(21) = {21, 22};
-Line(22) = {22, 23};
-Line(23) = {23, 24};
-Line(24) = {24, 25};
-Line(25) = {25, 26};
-Line(26) = {26, 27};
-Line(27) = {27, 28};
-Line(28) = {28, 21};
-
-Point(31) = {l_tot/2, h_tot/2, 0, 0.2*scale};
-Point(32) = {-l_tot/2, h_tot/2, 0, 0.2*scale};
-
-Line(32) = {6, 31};
-Line(33) = {31, 32};
-Line(34) = {32, 1};
-
-Curve Loop(4) = {1, 2, 3, 4, 5, 32, 33, 34};
-Curve Loop(5) = {21, 22, 23, 24, 25, 26, 27, 28};
-Plane Surface(3) = {4, 5};
-
-Physical Surface("BEM_domain_2", 9) = {3};
-Physical Curve("BEM_FEM_boundary_2", 10) = {1, 2, 3, 4, 5};
-Physical Curve("electrode_2", 11) = {21, 22, 23, 24, 25, 26, 27, 28};
-Physical Curve("outside_2", 12) = {32, 33, 34};
-
-phi_2 = 25;
-SetNumber("Boundary Conditions/BEM_FEM_boundary_2/BEM_domain_2/dirichlet", 0);
-SetNumber("Boundary Conditions/electrode_2/BEM_domain_2/dirichlet", phi_2);
-SetNumber("Boundary Conditions/outside_2/BEM_domain_2/neumann", 0);
-SetNumber("Materials/BEM_domain_2/Epsilon", 8.8541878128e-12); // dielectric permittivity
-
-
-Mesh.Algorithm = 8;
\ No newline at end of file
diff --git a/srcs/main.cpp b/srcs/main.cpp
index 9441044f46fb72de62c7539748e282f8a59afcb3..817334806b0ad1d63d3f7b4069e6c25a9b3f25f5 100644
--- a/srcs/main.cpp
+++ b/srcs/main.cpp
@@ -33,6 +33,9 @@ int main(int argc, char **argv)
     if(nonLinearSolver)
     {
 
+        std::cout << "-----------------------------------------\n";
+        std::cout << "NON-LINEAR TWO-WAY COUPLED ITERATIVE SOLVER\n";
+
         bool postProcessing = 0; // pass it as an argument to the solver, only compute post processing at last iteration
 
         std::map<int, double> electrostaticPressure;
@@ -114,6 +117,9 @@ int main(int argc, char **argv)
     else
     {
 
+        std::cout << "-----------------------------------------\n";
+        std::cout << "LINEAR ONE-WAY COUPLED SOLVER\n";
+
         std::map<int, double> electrostaticPressure;
 
         std::map<int,std::pair<double,double>> boundaryDisplacementMap;
diff --git a/srcs/squeeze.geo b/srcs/squeeze.geo
index ad1a4a9bd13f24d997c7c2e461eb74453df44709..fc23efd5c4281e22c98503bab48d0d68d5ad7be6 100644
--- a/srcs/squeeze.geo
+++ b/srcs/squeeze.geo
@@ -1,25 +1,36 @@
 scale = 4e-6;
 
-h_tot = 5*scale;
-w_tot = 10*scale;
-h_foot = 4*scale;
-h_diaph = 4.4*scale;
-h_end_diaph = 4.6*scale;
-Lx_diaph = 1*scale;
-
-nx = 10; // prend beaucoup de temps àpd de 200x40
-ny = 4;
-
-Point(1) = {0, 0, 0, 2};
-Point(2) = {w_tot, 0, 0, 2};
-Point(3) = {w_tot, h_foot, 0, 2};
+h_tot0 = 5;
+w_tot0 = 10;
+h_foot0 = 4;
+h_diaph0 = 4.4;
+h_end_diaph0 = 4.6;
+Lx_diaph0 = 1;
+
+unit_l = 0.2; // serves as a reference
+
+h_tot = h_tot0*scale;
+w_tot = w_tot0*scale;
+h_foot = h_foot0*scale;
+h_diaph = h_diaph0*scale;
+h_end_diaph = h_end_diaph0*scale;
+Lx_diaph = Lx_diaph0*scale;
+
+n = 8; // FEM elements density
+nBEM = 0.25; // BEM elements density
+
+// THIS .geo FILE IS NOT USED IN THE REPORT
+
+Point(1) = {0, 0, 0, 0.8*nBEM*scale};
+Point(2) = {w_tot, 0, 0, 0.8*nBEM*scale};
+Point(3) = {w_tot, h_foot, 0, 0.4*nBEM*scale};
 Point(4) = {w_tot, h_tot, 0, 2};
 Point(5) = {w_tot-Lx_diaph, h_end_diaph, 0, 2};
 Point(6) = {Lx_diaph, h_end_diaph, 0, 2};
 Point(7) = {0, h_tot, 0, 2};
-Point(8) = {0, h_foot, 0, 2};
-Point(9) = {Lx_diaph, h_diaph, 0, 2};
-Point(10) = {w_tot-Lx_diaph, h_diaph, 0, 2};
+Point(8) = {0, h_foot, 0, 0.4*nBEM*scale};
+Point(9) = {Lx_diaph, h_diaph, 0, 0.1*nBEM*scale};
+Point(10) = {w_tot-Lx_diaph, h_diaph, 0, 0.1*nBEM*scale};
 
 Line(1) = {9, 8};
 Line(2) = {10, 9};
@@ -34,22 +45,31 @@ Line(9) = {8, 1};
 Line(10) = {1, 2};
 Line(11) = {2 ,3};
 
-Curve Loop(1) = {1, 2, 3, 4, 5, 6, 7, 8};
+Line(13) = {9, 6};
+Line(14) = {5, 10};
+
+Curve Loop(1) = {-13, 1, 8, 7};
 Plane Surface(1) = {1};
 
-Curve Loop(2) = {9, 10, 11, 3, 2, 1};
-Plane Surface(2) = {2}; 
+Curve Loop(2) = {2, 13, 6, 14};
+Plane Surface(2) = {2};
+
+Curve Loop(3) = {-14, 5, 4, 3};
+Plane Surface(3) = {3};
+
+Curve Loop(4) = {9, 10, 11, 3, 2, 1};
+Plane Surface(4) = {4}; 
 
-Transfinite Curve {10} = 10*nx+1 Using Progression 1;
-Transfinite Curve {1, 2, 3} = 3*nx+1 Using Progression 1;
-Transfinite Curve {7, 6, 5} = 3*nx+1 Using Progression 1;
-Transfinite Curve {8, 4} = 2*ny+1 Using Progression 1;
-Transfinite Curve {7, 1, 5, 3} = 1*ny+1 Using Progression 1;
-Transfinite Curve {9, 11}  = 4*ny+1 Using Progression 1;
+Transfinite Curve {13, 14} = (h_end_diaph0 - h_diaph0)/unit_l*n+1 Using Progression 1;
+Transfinite Curve {2, 6} = (w_tot0 - 2*Lx_diaph0)/unit_l*n+1 Using Progression 1;
+Transfinite Curve {4, 8} = (h_tot0-h_foot0)/unit_l*n+1 Using Progression 1;
+Transfinite Curve {1, 3, 5, 7} = (h_tot0-h_foot0)/unit_l*n+1 Using Progression 1;
 
-//Transfinite Surface {1};
+Transfinite Surface {2};
 
 Recombine Surface {1};
+Recombine Surface {2};
+Recombine Surface {3};
 
 Mesh.ElementOrder = 1;
 
@@ -62,8 +82,8 @@ Physical Curve ("right_up_foot", 10) = {5};
 Physical Curve ("up_diaph", 11) = {6};
 Physical Curve ("low_diaph", 12) = {2};
 
-Physical Surface("FEM_domain", 13) = {1};
-Physical Surface("BEM_domain_1", 14) = {2};
+Physical Surface("FEM_domain", 13) = {1,2,3};
+Physical Surface("BEM_domain_1", 14) = {4};
 
 
 SetNumber("Boundary Conditions/left_edge_foot/ux", 0.); 
@@ -79,9 +99,9 @@ SetNumber("Boundary Conditions/right_low_foot/uy", 0);
 SetNumber("Boundary Conditions/right_up_foot/ux", 0.);
 SetNumber("Boundary Conditions/right_up_foot/uy", 0);
 
-SetNumber("Materials/FEM_domain/Young", 210e3);
-SetNumber("Materials/FEM_domain/Poisson", 0.3);
-SetNumber("Materials/FEM_domain/rho",7800);
+SetNumber("Materials/FEM_domain/Young", 150e9);
+SetNumber("Materials/FEM_domain/Poisson", 0.27);
+SetNumber("Materials/FEM_domain/rho",2300);
 
 SetNumber("Volumic Forces/FEM_domain/bx",0.);
 SetNumber("Volumic Forces/FEM_domain/by",0.);
diff --git a/srcs/transverse_comb.geo b/srcs/transverse_comb.geo
deleted file mode 100644
index 21122d65ca96ae9a77466cb6eb8c97cf2f302e65..0000000000000000000000000000000000000000
--- a/srcs/transverse_comb.geo
+++ /dev/null
@@ -1,157 +0,0 @@
-scale = 1e-5;
-
-h_tot = 20*scale;
-h_base = 3*scale;
-h_pin = 2.5*scale;
-h_space = 1*scale;
-
-l_tot = 20*scale;
-l_base = 1*scale;
-l_pin = 0.8*scale;
-l_space = 1.2*scale;
-
-t_electrode = 1*scale; // width of fixed electrodes
-// FEM domain
-Point(1) = {-l_tot/2, -h_tot/2, 0, 0.2*scale};
-Point(2) = {-l_tot/2, h_tot/2, 0, 0.2*scale};
-Point(3) = {-(l_tot/2 - l_base), h_tot/2, 0, 0.2*scale};
-Point(4) = {-(l_tot/2 - l_base), h_base/2, 0, 0.2*scale};
-Point(5) = {-l_pin/2, h_base/2, 0, 0.2*scale};
-Point(6) = {-l_pin/2, h_base/2 + h_pin, 0, 0.2*scale};
-Point(7) = {l_pin/2, h_base/2 + h_pin, 0, 0.2*scale};
-Point(8) = {l_pin/2, h_base/2, 0, 0.2*scale};
-Point(9) = {(l_tot/2 - l_base), h_base/2, 0, 0.2*scale};
-Point(10) = {(l_tot/2 - l_base), h_tot/2, 0, 0.2*scale};
-Point(11) = {l_tot/2, h_tot/2, 0, 0.2*scale};
-Point(12) = {l_tot/2, -h_tot/2, 0, 0.2*scale};
-Point(13) = {(l_tot/2 - l_base), -h_tot/2, 0, 0.2*scale};
-Point(14) = {(l_tot/2 - l_base), -h_base/2, 0, 0.2*scale};
-Point(15) = {l_pin/2, -h_base/2, 0, 0.2*scale};
-Point(16) = {l_pin/2, -(h_base/2 + h_pin), 0, 0.2*scale};
-Point(17) = {-l_pin/2, -(h_base/2 + h_pin), 0, 0.2*scale};
-Point(18) = {-l_pin/2, -h_base/2, 0, 0.2*scale};
-Point(19) = {-(l_tot/2 - l_base), -h_base/2, 0, 0.2*scale};
-Point(20) = {-(l_tot/2 - l_base), -h_tot/2, 0, 0.2*scale};
-
-// FEM domain
-Line(1) = {1, 2};
-Line(2) = {2, 3};
-Line(3) = {3, 4};
-Line(4) = {4, 5};
-Line(5) = {5, 6};
-Line(6) = {6, 7};
-Line(7) = {7, 8};
-Line(8) = {8, 9};
-Line(9) = {9, 10};
-Line(10) = {10, 11};
-Line(11) = {11, 12};
-Line(12) = {12, 13};
-Line(13) = {13, 14};
-Line(14) = {14, 15};
-Line(15) = {15, 16};
-Line(16) = {16, 17};
-Line(17) = {17, 18};
-Line(18) = {18, 19};
-Line(19) = {19, 20};
-Line(20) = {20, 1};
-
-Curve Loop(1) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20};
-Plane Surface(1) = {1};
-
-Recombine Surface {1}; // quads instead of triangles
-
-Physical Curve("bottom_left", 1) = {20};
-Physical Curve("top_left", 2) = {2};
-Physical Curve("top_right", 3) = {10};
-Physical Curve("bottom_right", 4) = {12};
-
-Physical Surface("FEM_domain", 5) = {1};
-Physical Curve("BEM_FEM_boundary", 6) = {3, 4, 5, 6, 7, 8, 9, 13, 14, 15, 16, 17, 18, 19};
-
-// mechanical properties and boundary conditions
-SetNumber("Boundary Conditions/bottom_left/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/bottom_left/uy", 0);
-SetNumber("Boundary Conditions/top_left/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/top_left/uy", 0);
-SetNumber("Boundary Conditions/top_right/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/top_right/uy", 0);
-SetNumber("Boundary Conditions/bottom_right/ux", 0.); // encastrement
-SetNumber("Boundary Conditions/bottom_right/uy", 0);
-SetNumber("Materials/domain/Young", 210e3); // A DETERMINER PRECISEMENT
-SetNumber("Materials/domain/Poisson", 0.3);
-SetNumber("Materials/domain/rho",7800); //volumic mass of acier
-SetNumber("Volumic Forces/FEM_domain/bx",50); // acceleration of accelerometer
-SetNumber("Volumic Forces/FEM_domain/by",0.); 
-
-// BEM domain 1
-Point(21) = {-(l_pin/2 + l_space + t_electrode), -(h_base/2 + h_pin + h_space + t_electrode), 0, 0.2*scale}; // first fixed electrode
-Point(22) = {-(l_pin/2 + l_space + t_electrode), -(h_base/2 + h_space), 0, 0.2*scale};
-Point(23) = {-(l_pin/2 + l_space), -(h_base/2 + h_space), 0, 0.2*scale};
-Point(24) = {-(l_pin/2 + l_space), -(h_base/2 + h_pin + h_space), 0, 0.2*scale};
-Point(25) = {l_pin/2 + l_space, -(h_base/2 + h_pin + h_space), 0, 0.2*scale};
-Point(26) = {l_pin/2 + l_space, -(h_base/2 + h_space), 0, 0.2*scale};
-Point(27) = {l_pin/2 + l_space + t_electrode, -(h_base/2 + h_space), 0, 0.2*scale};
-Point(28) = {l_pin/2 + l_space + t_electrode, -(h_base/2 + h_pin + h_space + t_electrode), 0, 0.2*scale};
-
-Line(21) = {21, 22};
-Line(22) = {22, 23};
-Line(23) = {23, 24};
-Line(24) = {24, 25};
-Line(25) = {25, 26};
-Line(26) = {26, 27};
-Line(27) = {27, 28};
-Line(28) = {28, 21};
-
-Line(37) = {20, 13};
-
-Curve Loop(2) = {37, 13, 14, 15, 16, 17, 18, 19};
-Curve Loop(3) = {21, 22, 23, 24, 25, 26, 27, 28};
-Plane Surface(2) = {2, 3};
-
-Physical Surface("BEM_domain_1", 7) = {2};
-Physical Curve("BEM_FEM_boundary_1", 8) = {13, 14, 15, 16, 17, 18, 19};
-Physical Curve("electrode_1", 9) = {21, 22, 23, 24, 25, 26, 27, 28};
-Physical Curve("outside_1", 10) = {37};
-
-phi_1 = 30;
-SetNumber("Boundary Conditions/BEM_FEM_boundary_1/BEM_domain_1/dirichlet", 0);
-SetNumber("Boundary Conditions/electrode_1/BEM_domain_1/dirichlet", phi_1);
-SetNumber("Boundary Conditions/outside_1/BEM_domain_1/neumann", 0);
-SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
-
-// BEM domain 2
-Point(29) = {-(l_pin/2 + l_space + t_electrode), h_base/2 + h_space, 0, 0.2*scale}; // second fixed electrode
-Point(30) = {-(l_pin/2 + l_space + t_electrode), h_base/2 + h_pin + h_space + t_electrode, 0, 0.2*scale}; 
-Point(31) = {l_pin/2 + l_space + t_electrode, h_base/2 + h_pin + h_space + t_electrode, 0, 0.2*scale};
-Point(32) = {l_pin/2 + l_space + t_electrode, h_base/2 + h_space, 0, 0.2*scale};
-Point(33) = {l_pin/2 + l_space, h_base/2 + h_space, 0, 0.2*scale};
-Point(34) = {l_pin/2 + l_space, h_base/2 + h_pin + h_space, 0, 0.2*scale};
-Point(35) = {-(l_pin/2 + l_space), h_base/2 + h_pin + h_space, 0, 0.2*scale};
-Point(36) = {-(l_pin/2 + l_space), h_base/2 + h_space, 0, 0.2*scale};
-
-
-Line(29) = {29, 30};
-Line(30) = {30, 31};
-Line(31) = {31, 32};
-Line(32) = {32, 33};
-Line(33) = {33, 34};
-Line(34) = {34, 35};
-Line(35) = {35, 36};
-Line(36) = {36, 29};
-
-Line(38) = {10, 3};
-
-Curve Loop(4) = {38, 3, 4, 5, 6, 7, 8, 9};
-Curve Loop(5) = {29, 30, 31, 32, 33, 34, 35, 36};
-Plane Surface(3) = {4, 5};
-
-Physical Surface("BEM_domain_2", 11) = {3};
-Physical Curve("BEM_FEM_boundary_2", 12) = {3, 4, 5, 6, 7, 8, 9};
-Physical Curve("electrode_2", 13) = {29, 30, 31, 32, 33, 34, 35, 36};
-Physical Curve("outside_2", 14) = {38};
-
-phi_2 = 30;
-SetNumber("Boundary Conditions/BEM_FEM_boundary_2/BEM_domain_2/dirichlet", 0);
-SetNumber("Boundary Conditions/electrode_2/BEM_domain_2/dirichlet", phi_2);
-SetNumber("Boundary Conditions/outside_2/BEM_domain_2/neumann", 0);
-SetNumber("Materials/BEM_domain_2/Epsilon", 8.8541878128e-12); // dielectric permittivity
\ No newline at end of file
diff --git a/srcs/two_BEM_test.geo b/srcs/two_BEM_test.geo
index 89926c27ade94d3510c858b0f78b71bc8ac2cc14..449a3efd30ec9a30eccfceae47ea7e4505dbf785 100644
--- a/srcs/two_BEM_test.geo
+++ b/srcs/two_BEM_test.geo
@@ -7,6 +7,8 @@ h_FEM = 4.5*scale;
 nx = 40;
 ny = 4;
 
+// THIS .geo FILE IS NOT USED IN THE REPORT
+
 Point(1) = {0, 0, 0, 2};
 Point(2) = {w_tot, 0, 0, 2};
 Point(3) = {w_tot, h_FEM, 0, 2};