From 658611b29633f32d37bfe080598976bebbe0f6c1 Mon Sep 17 00:00:00 2001
From: Tristan Carrier Baudouin <tristan.carrier@uclouvain.be>
Date: Thu, 7 Feb 2013 10:38:17 +0000
Subject: [PATCH] Supplementary option to generate a Voronoi geo file
---
Common/CommandLine.cpp | 33 +++++++++++++++++++++++++++++++++
1 file changed, 33 insertions(+)
diff --git a/Common/CommandLine.cpp b/Common/CommandLine.cpp
index 3ec33d7a0a..a48ac2d2d4 100644
--- a/Common/CommandLine.cpp
+++ b/Common/CommandLine.cpp
@@ -17,6 +17,8 @@
#include "GModel.h"
#include "CreateFile.h"
#include "OS.h"
+#include <fstream>
+#include <periodical.h>
#if defined(HAVE_FLTK)
#include <FL/Fl.H>
@@ -77,6 +79,7 @@ std::vector<std::pair<std::string, std::string> > GetUsage()
s.push_back(mp("-hoElasticity float","Poisson ration for elasticity analogy (nu in [-1.0,0.5])"));
s.push_back(mp("-optimize[_netgen]", "Optimize quality of tetrahedral elements"));
s.push_back(mp("-optimize_lloyd", "Optimize 2D meshes using Lloyd algorithm"));
+ s.push_back(mp("-microstructure", "Generate polycrystal Voronoi geometry"));
s.push_back(mp("-clscale float", "Set global mesh element size scaling factor"));
s.push_back(mp("-clmin float", "Set minimum mesh element size"));
s.push_back(mp("-clmax float", "Set maximum mesh element size"));
@@ -468,6 +471,36 @@ void GetOptions(int argc, char *argv[])
CTX::instance()->mesh.optimizeLloyd = atoi(argv[i++]);
else
Msg::Fatal("Missing number of lloyd iterations");
+ }
+ else if(!strcmp(argv[i] + 1, "microstructure")) {
+ i++;
+
+ int j;
+ double number;
+ double temp;
+ std::vector<double> coordinates;
+
+ if(argv[i]){
+ std::ifstream file(argv[i++]);
+ file >> number;
+ coordinates.clear();
+ coordinates.resize(3*number);
+
+ for(j=0;j<number;j++){
+ file >> coordinates[3*j];
+ file >> coordinates[3*j+1];
+ file >> coordinates[3*j+2];
+ file >> temp;
+ }
+
+ voroMetal3D vm1;
+ vm1.execute(coordinates,0.1);
+
+ GModel::current()->load("MicrostructurePolycrystal3D.geo");
+
+ voroMetal3D vm2;
+ vm2.correspondance(0.00001);
+ }
}
else if(!strcmp(argv[i] + 1, "nopopup")) {
CTX::instance()->noPopup = 1;
--
GitLab