From 11a62146a5cad64f026565683c48a1fd12207b64 Mon Sep 17 00:00:00 2001
From: Tristan Carrier Baudouin <tristan.carrier@uclouvain.be>
Date: Tue, 11 Dec 2012 17:59:17 +0000
Subject: [PATCH] periodic Voronoi
---
Mesh/periodical.cpp | 13 ++++++++++++-
1 file changed, 12 insertions(+), 1 deletion(-)
diff --git a/Mesh/periodical.cpp b/Mesh/periodical.cpp
index cfa384f8dd..d214718fcd 100644
--- a/Mesh/periodical.cpp
+++ b/Mesh/periodical.cpp
@@ -129,6 +129,7 @@ void voroMetal3D::execute(std::vector<SPoint3>& vertices,double h)
int val;
int face_number;
int last;
+ int mem;
double x,y,z;
double x1,y1,z1;
double x2,y2,z2;
@@ -205,6 +206,7 @@ void voroMetal3D::execute(std::vector<SPoint3>& vertices,double h)
std::ofstream file("cells.pos");
file << "View \"test\" {\n";
+
std::ofstream file2("cells.geo");
file2 << "c=" << h << ";\n";
@@ -303,6 +305,9 @@ void voroMetal3D::execute(std::vector<SPoint3>& vertices,double h)
for(j=0;j<obj.line_loops2.size();j++){
print_geo_face(get_counter(),obj.line_loops2[j],file2);
obj.faces2.push_back(get_counter());
+ mem = get_counter();
+ increase_counter();
+ print_geo_physical_face(get_counter(),mem,file2);
increase_counter();
}
@@ -358,6 +363,12 @@ void voroMetal3D::print_geo_face(int index1,int index2,std::ofstream& file){
<< "};\n";
}
+void voroMetal3D::print_geo_physical_face(int index1,int index2,std::ofstream& file){
+ file << "Physical Surface(" << index1 << ")={"
+ << index2
+ << "};\n";
+}
+
void voroMetal3D::print_geo_volume(int index1,int index2,std::ofstream& file){
file << "Volume(" << index1 << ")={"
<< index2
@@ -516,7 +527,7 @@ void voroMetal3D::correspondance(double e){
print_segment(p1,p2,file);
- file2 << faces[i]->tag() << " " << faces[j]->tag() << " " << p2.x()-p1.x() << " " << p2.y()-p1.y() << " " << p2.z()-p1.z() << "\n";
+ file2 << faces[i]->physicals[0] << " " << faces[j]->physicals[0] << " " << p2.x()-p1.x() << " " << p2.y()-p1.y() << " " << p2.z()-p1.z() << "\n";
count++;
}
--
GitLab