src/graphic/vtk_plot.cc

   1 /*****************************************************************************/
   2 /*              8888888         88888888         88888888                    */
   3 /*            8                8                8                            */
   4 /*           8                 8                8                            */
   5 /*           8                  88888888         88888888                    */
   6 /*           8      8888                8                8                   */
   7 /*            8       8                 8                8                   */
   8 /*              888888         888888888        888888888                    */
   9 /*                                                                           */
  10 /*       A Two-Dimensional General Purpose Semiconductor Simulator.          */
  11 /*                                                                           */
  12 /*  GSS 0.4x                                                                 */
  13 /*  Last update: Oct 30, 2006                                                */
  14 /*                                                                           */
  15 /*  Gong Ding                                                                */
  16 /*  gdiso@ustc.edu                                                           */
  17 /*  NINT, No.69 P.O.Box, Xi'an City, China                                   */
  18 /*                                                                           */
  19 /*****************************************************************************/
  20 #define HAVE_VTK
  21 #ifdef HAVE_VTK
  22
  23 #include <unistd.h>
  24 #include "bsolver.h"
  25 #include "log.h"
  26
  27 int BSolver::vtk_display_data(PlotDefine & plot_define)
  28 {
  29     // This is the child process.
  30     char filename[32];
  31     char *vtk_wish = new char [strlen(getenv("GSS_DIR"))+32];
  32     char *vtk_script = new char [strlen(getenv("GSS_DIR"))+32];
  33     sprintf(vtk_wish,"%s/bin/vtk",getenv("GSS_DIR"));
  34     sprintf(vtk_script,"%s/bin/vtk_plot.tcl",getenv("GSS_DIR"));
  35
  36     if ( plot_define.Variable == Potential)
  37     {
  38       sprintf(filename,"data_%s.vtk","Potential");
  39       vtk_output_file(filename);
  40       if(execl(vtk_wish, "vtk",vtk_script,filename,"Potential",0)==-1)
  41       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  42     }
  43     else if (plot_define.Variable == DopingNa)
  44     {
  45       sprintf(filename,"data_%s.vtk","Na");
  46       vtk_output_file(filename);
  47       if(execl(vtk_wish, "vtk",vtk_script,filename,"Na",0)==-1)
  48       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  49     }
  50     else if (plot_define.Variable == DopingNd)
  51     {
  52       sprintf(filename,"data_%s.vtk","Nd");
  53       vtk_output_file(filename);
  54       if(execl(vtk_wish, "vtk",vtk_script,filename,"Nd",0)==-1)
  55       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  56     }
  57     else if (plot_define.Variable == ElecDensity)
  58     {
  59       sprintf(filename,"data_%s.vtk","ElecDensity");
  60       vtk_output_file(filename);
  61       if(execl(vtk_wish, "vtk",vtk_script,filename,"ElecDensity",0)==-1)
  62       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  63     }
  64     else if (plot_define.Variable == HoleDensity)
  65     {
  66       sprintf(filename,"data_%s.vtk","HoleDensity");
  67       vtk_output_file(filename);
  68       if(execl(vtk_wish, "vtk",vtk_script,filename,"HoleDensity",0)==-1)
  69       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  70     }
  71     else if (plot_define.Variable == Temperature)
  72     {
  73       sprintf(filename,"data_%s.vtk","Temperature");
  74       vtk_output_file(filename);
  75       if(execl(vtk_wish, "vtk",vtk_script,filename,"Temperature",0)==-1)
  76       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  77     }
  78     else if (plot_define.Variable == OpticalEz)
  79     {
  80       sprintf(filename,"data_%s.vtk","OpticalEz");
  81       vtk_output_file(filename);
  82       if(execl(vtk_wish, "vtk",vtk_script,filename,"OpticalEz",0)==-1)
  83       { gss_log.string_buf()<<"VTK plot error!";   gss_log.record();}
  84     }
  85     delete [] vtk_wish;
  86     delete [] vtk_script;
  87     return 0;
  88 }
  89
  90
  91
  92 int BSolver::vtk_output_file(char *filename)
  93 {
  94   FILE *fp;
  95   fp = fopen(filename,"w");
  96   if(fp==NULL)
  97   {
  98     sprintf(log_buf,"Error, can't open file for vtk data output\n");
  99     GSS_LOG();
 100     return 1;
 101   }
 102
 103   fprintf(fp,"# vtk DataFile Version 2.0\n");
 104   fprintf(fp,"GSS data file\n");
 105   fprintf(fp,"ASCII\n");
 106   fprintf(fp,"DATASET UNSTRUCTURED_GRID\n\n");
 107   fprintf(fp,"POINTS %d double\n",gnode.size());
 108   for(int i=0;i<gnode.size();i++)
 109     fprintf(fp,"%e %e %lf\n",gnode[i].x/scale_unit.s_micron,gnode[i].y/scale_unit.s_micron,0.0);
 110   fprintf(fp,"\nCELLS %d %d\n",gtri.size(),4*gtri.size());
 111   for(int i=0;i<gtri.size();i++)
 112     fprintf(fp,"%d %d %d %d\n",3,gtri[i].g_node[0],gtri[i].g_node[1],gtri[i].g_node[2]);
 113   fprintf(fp,"CELL_TYPES %d\n",gtri.size());
 114   for(int i=0;i<gtri.size();i++)
 115     fprintf(fp,"%d\n",5);
 116   fprintf(fp,"\n");
 117
 118   fprintf(fp,"\nPOINT_DATA %d\n",gnode.size());
 119
 120   fprintf(fp,"SCALARS Potential double 1\n");
 121   fprintf(fp,"LOOKUP_TABLE default\n");
 122   for(int i=0;i<gnode.size();i++)
 123   {
 124     int zone_index = gnode[i].zone_index;
 125     int local_index = gnode[i].local_index;
 126     if(zonedata[zone_index]->material_type == Semiconductor)
 127     {
 128       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 129       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].P/scale_unit.s_volt));
 130     }
 131     else if(zonedata[zone_index]->material_type == Insulator)
 132     {
 133       ISZone *pzonedata = dynamic_cast< ISZone * >(zonedata[zone_index]);
 134       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].P/scale_unit.s_volt));
 135     }
 136     else if(zonedata[zone_index]->material_type == Conductor)
 137     {
 138       ElZone *pzonedata = dynamic_cast< ElZone * >(zonedata[zone_index]);
 139       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].P/scale_unit.s_volt));
 140     }
 141     else
 142       fprintf(fp,"%e\n",0.0);
 143   }
 144
 145   fprintf(fp,"SCALARS ElecDensity double 1\n");
 146   fprintf(fp,"LOOKUP_TABLE default\n");
 147   for(int i=0;i<gnode.size();i++)
 148   {
 149     int zone_index = gnode[i].zone_index;
 150     int local_index = gnode[i].local_index;
 151     if(zonedata[zone_index]->material_type == Semiconductor)
 152     {
 153       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 154       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].n*pow(scale_unit.s_centimeter,3)));
 155     }
 156     else
 157       fprintf(fp,"%e\n",0.0);
 158   }
 159
 160   fprintf(fp,"SCALARS HoleDensity double 1\n");
 161   fprintf(fp,"LOOKUP_TABLE default\n");
 162   for(int i=0;i<gnode.size();i++)
 163   {
 164     int zone_index = gnode[i].zone_index;
 165     int local_index = gnode[i].local_index;
 166     if(zonedata[zone_index]->material_type == Semiconductor)
 167     {
 168       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 169       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].p*pow(scale_unit.s_centimeter,3)));
 170     }
 171     else
 172       fprintf(fp,"%e\n",0.0);
 173   }
 174
 175
 176   fprintf(fp,"SCALARS Nd double 1\n");
 177   fprintf(fp,"LOOKUP_TABLE default\n");
 178   for(int i=0;i<gnode.size();i++)
 179   {
 180     int zone_index = gnode[i].zone_index;
 181     int local_index = gnode[i].local_index;
 182     if(zonedata[zone_index]->material_type == Semiconductor)
 183     {
 184       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 185       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].Nd*pow(scale_unit.s_centimeter,3)));
 186     }
 187     else
 188       fprintf(fp,"%e\n",0.0);
 189   }
 190
 191   fprintf(fp,"SCALARS Na double 1\n");
 192   fprintf(fp,"LOOKUP_TABLE default\n");
 193   for(int i=0;i<gnode.size();i++)
 194   {
 195     int zone_index = gnode[i].zone_index;
 196     int local_index = gnode[i].local_index;
 197     if(zonedata[zone_index]->material_type == Semiconductor)
 198     {
 199       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 200       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].Na*pow(scale_unit.s_centimeter,3)));
 201     }
 202     else
 203       fprintf(fp,"%e\n",0.0);
 204   }
 205
 206
 207   fprintf(fp,"SCALARS Temperature double 1\n");
 208   fprintf(fp,"LOOKUP_TABLE default\n");
 209   for(int i=0;i<gnode.size();i++)
 210   {
 211     int zone_index = gnode[i].zone_index;
 212     int local_index = gnode[i].local_index;
 213     if(zonedata[zone_index]->material_type == Semiconductor)
 214     {
 215       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 216       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].T/scale_unit.s_kelvin));
 217     }
 218     else if(zonedata[zone_index]->material_type == Insulator)
 219     {
 220       ISZone *pzonedata = dynamic_cast< ISZone * >(zonedata[zone_index]);
 221       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].T/scale_unit.s_kelvin));
 222     }
 223     else if(zonedata[zone_index]->material_type == Conductor)
 224     {
 225       ElZone *pzonedata = dynamic_cast< ElZone * >(zonedata[zone_index]);
 226       fprintf(fp,"%e\n",double(pzonedata->fs[local_index].T/scale_unit.s_kelvin));
 227     }
 228     else
 229       fprintf(fp,"%e\n",0.0);
 230   }
 231
 232   fprintf(fp,"SCALARS opticalEz double 1\n");
 233   fprintf(fp,"LOOKUP_TABLE default\n");
 234   for(int i=0;i<gnode.size();i++)
 235   {
 236     int zone_index = gnode[i].zone_index;
 237     int local_index = gnode[i].local_index;
 238     if(zonedata[zone_index]->material_type == Semiconductor)
 239     {
 240       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 241       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpEz/scale_unit.s_volt*scale_unit.s_centimeter));
 242     }
 243     else if(zonedata[zone_index]->material_type == Insulator)
 244     {
 245      ISZone *pzonedata = dynamic_cast< ISZone * >(zonedata[zone_index]);
 246      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpEz/scale_unit.s_volt*scale_unit.s_centimeter));
 247     }
 248     else if(zonedata[zone_index]->material_type == Conductor)
 249     {
 250      ElZone *pzonedata = dynamic_cast< ElZone * >(zonedata[zone_index]);
 251      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpEz/scale_unit.s_volt*scale_unit.s_centimeter));
 252     }
 253     else if(zonedata[zone_index]->material_type == Vacuum)
 254     {
 255      VacuumZone *pzonedata = dynamic_cast< VacuumZone * >(zonedata[zone_index]);
 256      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpEz/scale_unit.s_volt*scale_unit.s_centimeter));
 257     }
 258     else if(zonedata[zone_index]->material_type == PML)
 259     {
 260       PMLZone *pzonedata = dynamic_cast< PMLZone * >(zonedata[zone_index]);
 261       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpEz/scale_unit.s_volt*scale_unit.s_centimeter));
 262     }
 263     else
 264       fprintf(fp,"%e\n",0.0);
 265   }
 266
 267   fprintf(fp,"SCALARS OptG double 1\n");
 268   fprintf(fp,"LOOKUP_TABLE default\n");
 269   for(int i=0;i<gnode.size();i++)
 270   {
 271     int zone_index = gnode[i].zone_index;
 272     int local_index = gnode[i].local_index;
 273     if(zonedata[zone_index]->material_type == Semiconductor)
 274     {
 275       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 276       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OptG*pow(scale_unit.s_centimeter,3)*scale_unit.s_second));
 277     }
 278     else
 279       fprintf(fp,"%e\n",0.0);
 280   }
 281
 282   fprintf(fp,"SCALARS opticalHz double 1\n");
 283   fprintf(fp,"LOOKUP_TABLE default\n");
 284   for(int i=0;i<gnode.size();i++)
 285   {
 286     int zone_index = gnode[i].zone_index;
 287     int local_index = gnode[i].local_index;
 288     if(zonedata[zone_index]->material_type == Semiconductor)
 289     {
 290       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 291       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpHz/scale_unit.s_A*scale_unit.s_centimeter));
 292     }
 293     else if(zonedata[zone_index]->material_type == Insulator)
 294     {
 295      ISZone *pzonedata = dynamic_cast< ISZone * >(zonedata[zone_index]);
 296      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpHz/scale_unit.s_A*scale_unit.s_centimeter));
 297     }
 298     else if(zonedata[zone_index]->material_type == Conductor)
 299     {
 300      ElZone *pzonedata = dynamic_cast< ElZone * >(zonedata[zone_index]);
 301      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpHz/scale_unit.s_A*scale_unit.s_centimeter));
 302     }
 303     else if(zonedata[zone_index]->material_type == Vacuum)
 304     {
 305      VacuumZone *pzonedata = dynamic_cast< VacuumZone * >(zonedata[zone_index]);
 306      fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpHz/scale_unit.s_A*scale_unit.s_centimeter));
 307     }
 308     else if(zonedata[zone_index]->material_type == PML)
 309     {
 310       PMLZone *pzonedata = dynamic_cast< PMLZone * >(zonedata[zone_index]);
 311       fprintf(fp,"%e\n",double(pzonedata->aux[local_index].OpHz/scale_unit.s_A*scale_unit.s_centimeter));
 312     }
 313   }
 314
 315   fprintf(fp,"VECTORS EField double \n");
 316   for(int i=0;i<gnode.size();i++)
 317   {
 318     int zone_index = gnode[i].zone_index;
 319     int local_index = gnode[i].local_index;
 320     if(zonedata[zone_index]->material_type == Semiconductor)
 321     {
 322       SMCZone *pzonedata = dynamic_cast< SMCZone * >(zonedata[zone_index]);
 323       fprintf(fp,"%e %e 0\n",double(pzonedata->aux[local_index].Ex), double(pzonedata->aux[local_index].Ey));
 324     }
 325     else if(zonedata[zone_index]->material_type == Insulator)
 326     {
 327       ISZone *pzonedata = dynamic_cast< ISZone * >(zonedata[zone_index]);
 328       fprintf(fp,"%e %e 0\n",double(pzonedata->aux[local_index].Ex), double(pzonedata->aux[local_index].Ey));
 329     }
 330     else
 331       fprintf(fp,"0 0 0\n");
 332   }
 333
 334
 335
 336   fclose(fp);
 337   return 0;
 338 }
 339
 340 #endif