màj

test_bem.xcodeproj/project.pbxproj

@@ -12,7 +12,7 @@
 		0A85289A27D3EDE500D1C30F /* code_louis_function.cpp in Sources */ = {isa = PBXBuildFile; fileRef = 0A85289827D3EDE500D1C30F /* code_louis_function.cpp */; };
 		0A85289D27D6D39C00D1C30F /* BEM_livre.cpp in Sources */ = {isa = PBXBuildFile; fileRef = 0A85289B27D6D39C00D1C30F /* BEM_livre.cpp */; };
 		0A8528A227D6D69F00D1C30F /* code_louis.cpp in Sources */ = {isa = PBXBuildFile; fileRef = 0A8528A027D6D69F00D1C30F /* code_louis.cpp */; };
-		0AAB3F3427DD204300FF9C2C /* main.cpp in Sources */ = {isa = PBXBuildFile; fileRef = 0AAB3F3327DD204300FF9C2C /* main.cpp */; };
+		0AAB3F3827DEB40100FF9C2C /* main.cpp in Sources */ = {isa = PBXBuildFile; fileRef = 0AAB3F3727DEB40100FF9C2C /* main.cpp */; };
 /* End PBXBuildFile section */
 
 /* Begin PBXCopyFilesBuildPhase section */
@@ -41,7 +41,7 @@
 		0A85289C27D6D39C00D1C30F /* BEM_livre.h */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType = sourcecode.c.h; path = BEM_livre.h; sourceTree = "<group>"; };
 		0A8528A027D6D69F00D1C30F /* code_louis.cpp */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType = sourcecode.cpp.cpp; path = code_louis.cpp; sourceTree = "<group>"; };
 		0A8528A127D6D69F00D1C30F /* code_louis.hpp */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType = sourcecode.cpp.h; path = code_louis.hpp; sourceTree = "<group>"; };
-		0AAB3F3327DD204300FF9C2C /* main.cpp */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType = sourcecode.cpp.cpp; path = main.cpp; sourceTree = "<group>"; };
+		0AAB3F3727DEB40100FF9C2C /* main.cpp */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType = sourcecode.cpp.cpp; path = main.cpp; sourceTree = "<group>"; };
 /* End PBXFileReference section */
 
 /* Begin PBXFrameworksBuildPhase section */
@@ -70,7 +70,7 @@
 				0A85287C27D242F200D1C30F /* my_geo.geo */,
 				0A85288C27D395CB00D1C30F /* code_kevin.cpp */,
 				0A85288D27D395CB00D1C30F /* code_kevin.hpp */,
-				0AAB3F3327DD204300FF9C2C /* main.cpp */,
+				0AAB3F3727DEB40100FF9C2C /* main.cpp */,
 				0A85289B27D6D39C00D1C30F /* BEM_livre.cpp */,
 				0A85289C27D6D39C00D1C30F /* BEM_livre.h */,
 				0A85289627D3ECBC00D1C30F /* test.cpp */,
@@ -142,7 +142,7 @@
 			files = (
 				0A8528A227D6D69F00D1C30F /* code_louis.cpp in Sources */,
 				0A85288E27D395CB00D1C30F /* code_kevin.cpp in Sources */,
-				0AAB3F3427DD204300FF9C2C /* main.cpp in Sources */,
+				0AAB3F3827DEB40100FF9C2C /* main.cpp in Sources */,
 				0A85289727D3ECBC00D1C30F /* test.cpp in Sources */,
 				0A85289D27D6D39C00D1C30F /* BEM_livre.cpp in Sources */,
 				0A85289A27D3EDE500D1C30F /* code_louis_function.cpp in Sources */,

test_bem/code_kevin.hpp

@@ -16,5 +16,11 @@ void show_mat(MatrixXd &A);
 
 void show_vec(VectorXd &A);
 
+void range_angle(double &angle);
+
+int code_py(int argc, char **argv);
+
+
+
 
 

test_bem/main.cpp

@@ -18,5 +18,7 @@ int main(int argc, char **argv)
 
     code_kevin(argc,argv);
 
+    //code_py(argc,argv);
+
     return 0;
 }

test_bem/my_geo.geo

@@ -12,12 +12,12 @@ Line(2) = {2, 3};
 Line(3) = {3, 4};
 Line(4) = {4, 1};
 Curve Loop(1) = {4, 1, 2, 3};
-//Plane Surface(1) = {1};
+Plane Surface(1) = {1};
 Transfinite Curve {3, 1} = nx+1 Using Progression 1;
 Transfinite Curve {4, 2} = ny+1 Using Progression 1;
-//Transfinite Surface {1};
+Transfinite Surface {1};
 
-//Recombine Surface {1}; // quads instead of triangles
+Recombine Surface {1}; // quads instead of triangles
 
 //à gauche c'est le tag qu'on donne, il doit être différent pour chaque curve ou surface et à droite c'est l'élément définit correspondant 4 courbes et 1 surface
 Physical Curve("bottom_edge", 1) = {1};
@@ -25,7 +25,7 @@ Physical Curve("right_edge", 2) = {2};
 Physical Curve("top_edge", 3) = {3};
 Physical Curve("left_edge", 4) = {4};
 //Physical Curve("boundary", 5) = {1,2,3,4};
-//Physical Surface("domain", 6) = {1};
+Physical Surface("domain", 6) = {1};
 
 // additional parameters given to the solver
 
@@ -37,19 +37,15 @@ Physical Curve("left_edge", 4) = {4};
 //SetNumber("Volumic Forces/domain/by",0.); //set to -9.81 for gravity
 
 //BC top edge
-SetNumber("Boundary Conditions/top_edge/neumann", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
-//SetNumber("Boundary Conditions/top_edge/neumann_y", 0.); //set to some non-zero value to induce vertical deflection
+SetNumber("Boundary Conditions/top_edge/neumann", 0.);
 
 //BC right edge
-SetNumber("Boundary Conditions/right_edge/dirichlet", 200); // for simple tension conditions
-//SetNumber("Boundary Conditions/right_edge/neumann_y", 0.);
+SetNumber("Boundary Conditions/right_edge/dirichlet", 200); 
 
 //BC bottom edge
 SetNumber("Boundary Conditions/bottom_edge/neumann", 0.);
-//SetNumber("Boundary Conditions/bottom_edge/neumann_y", 0.);
 
 // BC left edge
-SetNumber("Boundary Conditions/left_edge/dirichlet", 100); // HERE YOU DO NOT HAVE TO IMPOSE BOTH ux and uy simultaneously ! (permet aussi de simuler appuis à roulettes)
-//SetNumber("Boundary Conditions/left_edge/dirichlet_y", 0.);
+SetNumber("Boundary Conditions/left_edge/dirichlet", 100);
 //+
 Show "*";