X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fsecs2d-0.0.8%2FUtilities%2FUmsh2pdetool.m;fp=octave_packages%2Fsecs2d-0.0.8%2FUtilities%2FUmsh2pdetool.m;h=63c2ba47331cf6f27dd00d6c73973c572fe66a85;hp=0000000000000000000000000000000000000000;hb=c880e8788dfc484bf23ce13fa2787f2c6bca4863;hpb=1705066eceaaea976f010f669ce8e972f3734b05

diff --git a/octave_packages/secs2d-0.0.8/Utilities/Umsh2pdetool.m b/octave_packages/secs2d-0.0.8/Utilities/Umsh2pdetool.m
new file mode 100644
index 0000000..63c2ba4
--- /dev/null
+++ b/octave_packages/secs2d-0.0.8/Utilities/Umsh2pdetool.m
@@ -0,0 +1,61 @@
+function msh=Umsh2pdetool(filename);
+
+##
+##
+##  loadgmshmesh(filename);
+##
+##
+
+awk_command = "BEGIN {  filename = ARGV[1] ; gsub(/\\./,""_"",filename) }\n\
+\n\
+/\\$NOD/,/\\$ENDNOD/ { \n\
+    if ( FNR>2 ) \n\
+    {  \n\
+	if($0 ~ /^[^\\$]/  )  \n\
+	{\n\
+	    print ""p ( "" $1 "" ,:) = ["" $2 "" ""  $3""];""  > filename ""_p.m"" \n\
+	}\n\
+    } \n\
+} \n\
+\n\
+/\\$ELM/,/\\$ENDNELM/ { \n\
+    if (  $1 ~ /\\$ELM/   )\n\
+    {\n\
+	gsub(/\\$ELM/,""t=["")\n\
+	print > filename ""_t.m""\n\
+	gsub(/t=\\[/,""e=["")\n\
+	print > filename ""_e.m""\n\
+\n\
+    } else if  ($1 ~ /\\$ENDELM/ ){\n\
+		gsub(/\\$ENDELM/,""];"")\n\
+		print > filename ""_t.m""\n\
+		print > filename ""_e.m""\n\
+    }\n\
+    else if ( $2 == ""2"" )\n\
+    {\n\
+	print ( $6 "" "" $7 "" "" $8 "" "" $4)  > filename ""_t.m"" \n\
+    }\n\
+    else if ( $2 == ""1"" )\n\
+    {\n\
+	print ( $6 "" "" $7 "" 0 0 "" $4 "" 0 0"")  > filename ""_e.m"" \n\
+    }\n\
+    else if ( $2 == ""9"" )\n\
+    {\n\
+     print ( $6 "" "" $7 "" "" $8 "" "" $9 "" "" $10 "" "" $11 "" "" \
+	    $4)  > filename ""_t.m"" \n\
+    }\n\
+    else if ( $2 == ""8"" )\n\
+    {\n\
+     print ( $6 "" "" $7 "" "" $8 "" 0 "" $4)  > filename ""_e.m"" \n\
+    }\n\
+}\n\
+\n\
+{ }"
+
+system(["awk '" awk_command "' " filename ".msh"]);
+eval([ filename "_msh_p"]);
+eval([ filename "_msh_e"]);
+eval([ filename "_msh_t"]);
+
+
+msh=struct("p",p',"t",t',"e",e');
\ No newline at end of file