octave_packages/secs2d-0.0.8/Utilities/Ustructmesh_right.m

   1 function [p,e,t]=Ustructmesh_right(x,y,region,sides)
   2
   3 % [p,e,t]=Ustructmesh(x,y,region,sides)
   4
   5 % This file is part of
   6 %
   7 %            SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
   8 %         -------------------------------------------------------------------
   9 %            Copyright (C) 2004-2006  Carlo de Falco
  10 %
  11 %
  12 %
  13 %  SECS2D is free software; you can redistribute it and/or modify
  14 %  it under the terms of the GNU General Public License as published by
  15 %  the Free Software Foundation; either version 2 of the License, or
  16 %  (at your option) any later version.
  17 %
  18 %  SECS2D is distributed in the hope that it will be useful,
  19 %  but WITHOUT ANY WARRANTY; without even the implied warranty of
  20 %  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  21 %  GNU General Public License for more details.
  22 %
  23 %  You should have received a copy of the GNU General Public License
  24 %  along with SECS2D; If not, see <http://www.gnu.org/licenses/>.
  25
  26 x = sort(x);
  27 y = sort(y);
  28
  29 nx = length(x);
  30 ny = length(y);
  31 [XX,YY] = meshgrid(x,y);
  32 p = [XX(:),YY(:)]';
  33 iiv (ny,nx)=0;
  34 iiv(:)=1:nx*ny;
  35 iiv(end,:)=[];
  36 iiv(:,end)=[];
  37 iiv=iiv(:)';
  38 t = [[iiv;iiv+ny;iiv+ny+1],[iiv;iiv+ny+1;iiv+1] ];
  39 t (4,:)=region;
  40
  41 l1 = 1+ny*([1:nx]-1);
  42 l4 = 1:ny;
  43 l2 = ny*(nx-1)+1:nx*ny;
  44 l3 = ny + l1 -1;
  45
  46 e = [ l1([1:end-1]) l2([1:end-1]) l3([1:end-1]) l4([1:end-1])
  47        l1([2:end]) l2([2:end]) l3([2:end]) l4([2:end])
  48        [l1([1:end-1]) l2([1:end-1]) l3([1:end-1]) l4([1:end-1])]*0
  49        [l1([1:end-1]) l2([1:end-1]) l3([1:end-1]) l4([1:end-1])]*0
  50        l1([1:end-1])*0+sides(1) l2([1:end-1])*0+sides(2) l3([1:end-1])*0+sides(3) l4([1:end-1])*0+sides(4)
  51            [l1([1:end-1]) l2([1:end-1]) l3([1:end-1]) l4([1:end-1])]*0
  52            [l1([1:end-1]) l2([1:end-1]) l3([1:end-1]) l4([1:end-1])]*0+region
  53            ];
  54