Add a useful package (from Source forge) for octave

[CreaPhase.git] / octave_packages / secs2d-0.0.8 / Utilities / UDXoutput2Ddata.m

diff --git a/octave_packages/secs2d-0.0.8/Utilities/UDXoutput2Ddata.m b/octave_packages/secs2d-0.0.8/Utilities/UDXoutput2Ddata.m
new file mode 100644 (file)
index 0000000..754a92c
--- /dev/null
+++ b/octave_packages/secs2d-0.0.8/Utilities/UDXoutput2Ddata.m
@@ -0,0 +1,115 @@
+function UDXoutput2Ddata(filename,p,t,u,attr_name,attr_rank,attr_shape,endfile)
+
+##
+##   UDXoutput2Ddata(filename,p,t,u,attr_name,attr_rank,attr_shape,endfile)
+##
+##   Outputs data in DX form.
+##   Only one variable can be written to the file
+##   variable must be a scalar, vector or tensor of doubles   
+##
+##   x
+##   attr_name  = name of the variable                   (type string)
+##   attr_rank  = rank of variable data                  (0 for scalar, 1 for vector, etc.)
+##   attr_shape = number of components of variable data  (assumed 1 for scalar)
+##
+
+
+% This file is part of 
+%
+%            SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
+%         -------------------------------------------------------------------
+%            Copyright (C) 2004-2006  Carlo de Falco
+%
+%
+%
+%  SECS2D is free software; you can redistribute it and/or modify
+%  it under the terms of the GNU General Public License as published by
+%  the Free Software Foundation; either version 2 of the License, or
+%  (at your option) any later version.
+%
+%  SECS2D is distributed in the hope that it will be useful,
+%  but WITHOUT ANY WARRANTY; without even the implied warranty of
+%  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+%  GNU General Public License for more details.
+%
+%  You should have received a copy of the GNU General Public License
+%  along with SECS2D; If not, see <http://www.gnu.org/licenses/>.
+
+
+
+
+  p = p';
+  t = t';
+  t = t(:,1:3);
+
+  %eval(['!rm -f ',filename]);
+  
+  fid=fopen (filename,'w');
+  Nnodi = size(p,1);
+  Ntriangoli = size(t,1);
+  Ndati = size(u,1);
+
+  fprintf(fid,'object "pos"\nclass array type float rank 1 shape 2 items %d data follows',Nnodi);
+
+  for i=1:Nnodi
+    fprintf(fid,'\n');
+    fprintf(fid,'    %1.7e',p(i,:));
+  end
+  
+  if (min(min(t))==1)
+    t=t-1;
+  elseif(min(min(t))~=0)
+    disp('WARNING: check triangle structure')
+  end                    
+                               # In DX format nodes are 
+                               # numbered starting from zero,
+                               # instead we want to number
+                               # them starting from 1!
+                               # Here we restore the DX
+                               # format
+
+  fprintf(fid,'\n\nobject "con"\nclass array type int rank 1 shape 3 items %d data follows',Ntriangoli);
+  for i=1:Ntriangoli
+    fprintf(fid,'\n');
+    fprintf(fid,'      %d',t(i,:));
+  end
+  
+  fprintf(fid,'\nattribute "element type" string "triangles"\nattribute "ref" string "positions"\n\n');
+  
+  if ((attr_rank==0) & (min(size(u))==1))
+    fprintf(fid,'object "%s.data"\nclass array type double rank 0 items %d data follows',attr_name,Ndati);
+    fprintf(fid,'\n %1.7e',u);
+
+  else
+    fprintf(fid,'object "%s.data"\nclass array type double rank %d shape %d items %d data follows', ...
+            attr_name,attr_rank,attr_shape,Ndati);
+    for i=1:Ndati
+      fprintf(fid,'\n');
+      fprintf(fid,'    %1.7e',u(i,:));
+    end
+
+end
+
+if Ndati==Nnodi
+  fprintf(fid,['\nattribute "dep" string "positions"\n\n' ...
+               'object "%s" class field\n'...
+               'component "positions" value "pos"\n'...
+               'component "connections" value "con"\n'...
+               'component "data" value "%s.data"\n'],...
+         attr_name,attr_name);
+  elseif Ndati==Ntriangoli
+ fprintf(fid,['\nattribute "dep" string "connections"\n\n' ...
+               'object "%s" class field\n'...
+               'component "positions" value "pos"\n'...
+               'component "connections" value "con"\n'...
+               'component "data" value "%s.data"\n'],...
+         attr_name,attr_name);
+end
+
+  if(endfile)
+    fprintf(fid,'\nend\n');
+  end
+  
+  fclose (fid);
+
+