octave_packages/secs2d-0.0.8/Utilities/Umsh2pdetool.m

   1 function msh=Umsh2pdetool(filename);
   2
   3 ##
   4 ##
   5 ##  loadgmshmesh(filename);
   6 ##
   7 ##
   8
   9 awk_command = "BEGIN {  filename = ARGV[1] ; gsub(/\\./,""_"",filename) }\n\
  10 \n\
  11 /\\$NOD/,/\\$ENDNOD/ { \n\
  12     if ( FNR>2 ) \n\
  13     {  \n\
  14         if($0 ~ /^[^\\$]/  )  \n\
  15         {\n\
  16             print ""p ( "" $1 "" ,:) = ["" $2 "" ""  $3""];""  > filename ""_p.m"" \n\
  17         }\n\
  18     } \n\
  19 } \n\
  20 \n\
  21 /\\$ELM/,/\\$ENDNELM/ { \n\
  22     if (  $1 ~ /\\$ELM/   )\n\
  23     {\n\
  24         gsub(/\\$ELM/,""t=["")\n\
  25         print > filename ""_t.m""\n\
  26         gsub(/t=\\[/,""e=["")\n\
  27         print > filename ""_e.m""\n\
  28 \n\
  29     } else if  ($1 ~ /\\$ENDELM/ ){\n\
  30                 gsub(/\\$ENDELM/,""];"")\n\
  31                 print > filename ""_t.m""\n\
  32                 print > filename ""_e.m""\n\
  33     }\n\
  34     else if ( $2 == ""2"" )\n\
  35     {\n\
  36         print ( $6 "" "" $7 "" "" $8 "" "" $4)  > filename ""_t.m"" \n\
  37     }\n\
  38     else if ( $2 == ""1"" )\n\
  39     {\n\
  40         print ( $6 "" "" $7 "" 0 0 "" $4 "" 0 0"")  > filename ""_e.m"" \n\
  41     }\n\
  42     else if ( $2 == ""9"" )\n\
  43     {\n\
  44      print ( $6 "" "" $7 "" "" $8 "" "" $9 "" "" $10 "" "" $11 "" "" \
  45             $4)  > filename ""_t.m"" \n\
  46     }\n\
  47     else if ( $2 == ""8"" )\n\
  48     {\n\
  49      print ( $6 "" "" $7 "" "" $8 "" 0 "" $4)  > filename ""_e.m"" \n\
  50     }\n\
  51 }\n\
  52 \n\
  53 { }"
  54
  55 system(["awk '" awk_command "' " filename ".msh"]);
  56 eval([ filename "_msh_p"]);
  57 eval([ filename "_msh_e"]);
  58 eval([ filename "_msh_t"]);
  59
  60
  61 msh=struct("p",p',"t",t',"e",e');