+++ /dev/null
-/* atlc - arbitrary transmission line calculator, for the analysis of
-transmission lines are directional couplers.
-
-Copyright (C) 2002. Dr. David Kirkby, PhD (G8WRB).
-
-This program 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 package_version 2
-of the License, or (at your option) any later package_version.
-
-This program 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 this program; if not, write to the Free Software
-Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
-USA.
-
-Dr. David Kirkby, e-mail drkirkby at ntlworld.com
-
-*/
-
-/* The function setup_arrays sets the type of a pixel to be dielectric,
-or metal. This routines refines that, setting the type of a pixel to be
-some number between 1 and 28 according to where it is located, and what
-is around it. See Green paper to understand what they mean */
-
-#include "config.h"
-#include "definitions.h"
-#include "exit_codes.h"
-
-extern int width, height;
-extern unsigned char **cell_type, **oddity;
-extern double **Er;
-
-#ifdef HAVE_STDLIB_H
-#include <stdlib.h>
-#endif
-
-void set_oddity_value(void)
-{
- int i, j;
- double er, ERa, ERb, ERl, ERr;
- int cl, cr, ca, cb;
- unsigned char c;
-
- /* Its easier to set the endge values first, as it
- reduces the amount of checking needed in the main body.
- There are only 11 cases here - 3 types of metal,
- the four corners, and the four sides */
-
- for(i=0;i<width; ++i) {
- for(j=0;j<height;++j) {
-
- oddity[i][j]=UNDEFINED_ODDITY; /* Stick it to some underfined status */
-
- c=cell_type[i][j]; /* Cell type at point (i,j) */
-
-
- /* The 3 metal cases can be quickly checked and the
- oddity value assigned to a fixed value depending on
- whether it's -1, 0 or +1 V */
-
- if(c <= CONDUCTOR_PLUS_ONE_V ) /* a metal */
- oddity[i][j]=c;
-
-
- /* Now do the 4 courners */
- else if( (i == 0) && (j == height-1) )
- oddity[i][j]=BOTTOM_LEFT_CORNER;
-
- else if ( (i == width-1) && (j == height-1) )
- oddity[i][j]=BOTTOM_RIGHT_CORNER;
-
- else if( (i == 0) && (j == 0) )
- oddity[i][j]=TOP_LEFT_CORNER;
-
- else if( (i == width-1) && (j == 0) )
- oddity[i][j]=TOP_RIGHT_CORNER;
-
- /* Now the four edges */
- else if ( i == 0 )
- oddity[i][j]=ORDINARY_POINT_LEFT_EDGE;
-
- else if ( j == 0 )
- oddity[i][j]=ORDINARY_POINT_TOP_EDGE;
-
- else if ( j == height-1)
- oddity[i][j]=ORDINARY_POINT_BOTTOM_EDGE;
-
- else if ( i == width-1 )
- oddity[i][j]=ORDINARY_POINT_RIGHT_EDGE;
-
-
- else if ( (i == 0 || i == width-1 || j == 0 || j == height-1 ) && (oddity[i][j]==UNDEFINED_ODDITY)) {
- fprintf(stderr,"Internal error: one of the edge points (%d,%d) is still undefined\n",i,j);
- fprintf(stderr, "ZZZZZZZZZZZZZ width=%d height=%d\n", width, height);
- fprintf(stderr,"Error set_oddity_value.c\n");
- exit(INTERNAL_ERROR);
- }
- } /* end of for(j=0;j<height-0;++j) { */
- } /* end of for(i=0;i<width-0; ++i) { */
- /* With the oddity values of all the edges now know, the centre
- values can be attempted */
-
- for(i=1;i<width-1; ++i) {
- for(j=1;j<height-1;++j) {
-
- c=cell_type[i][j]; /* Cell type at point (i,j) */
- cl=cell_type[i-1][j]; /* Cell type to left of point (i,j) */
- cr=cell_type[i+1][j]; /* Cell type to right of point (i,j) */
- ca=cell_type[i][j-1]; /* Cell type above point (i,j) */
- cb=cell_type[i][j+1]; /* Cell type below point (i,j) */
-
- ERa=Er[i][j-1];
- ERb=Er[i][j+1];
- ERl=Er[i-1][j];
- ERr=Er[i+1][j];
- er=Er[i][j];
-
-
- /* If the conductor is at a fixed voltage, it must stay there
- so there is nothing to do with it */
-
- if(c == CONDUCTOR_ZERO_V )
- oddity[i][j] = CONDUCTOR_ZERO_V;
-
- else if(c== CONDUCTOR_PLUS_ONE_V)
- oddity[i][j] = CONDUCTOR_PLUS_ONE_V;
-
- else if(c == CONDUCTOR_MINUS_ONE_V )
- oddity[i][j] = CONDUCTOR_MINUS_ONE_V;
-
- else if ( cr <= CONDUCTOR_PLUS_ONE_V && cb<= CONDUCTOR_PLUS_ONE_V)
- oddity[i][j]= METAL_BELOW_AND_RIGHT;
-
- else if ( cr <= CONDUCTOR_PLUS_ONE_V && ca<= CONDUCTOR_PLUS_ONE_V)
- oddity[i][j]= METAL_ABOVE_AND_RIGHT;
-
- else if ( cl <= CONDUCTOR_PLUS_ONE_V && cb<= CONDUCTOR_PLUS_ONE_V)
- oddity[i][j]= METAL_BELOW_AND_LEFT;
-
- else if ( cl <= CONDUCTOR_PLUS_ONE_V && ca<= CONDUCTOR_PLUS_ONE_V)
- oddity[i][j]= METAL_ABOVE_AND_LEFT;
-
- else if ( ca <= CONDUCTOR_PLUS_ONE_V )
- oddity[i][j]= METAL_ABOVE;
-
- else if ( cb <= CONDUCTOR_PLUS_ONE_V )
- oddity[i][j]= METAL_BELOW;
-
- else if ( cl <= CONDUCTOR_PLUS_ONE_V )
- oddity[i][j]= METAL_LEFT;
-
- else if ( cr <= CONDUCTOR_PLUS_ONE_V )
- oddity[i][j]= METAL_RIGHT;
-
- else if ( ERb != ERa)
- oddity[i][j]= DIFFERENT_DIELECTRIC_VERTICALLY;
-
- else if ( ERl != ERr )
- oddity[i][j]= DIFFERENT_DIELECTRIC_HORIZONTALLY;
-/*
- else if ( ERa != er && ERr != er)
- oddity[i][j]= DIFFERENT_DIELECTRIC_ABOVE_AND_RIGHT;
-
- else if ( ERa != er && ERl != er)
- oddity[i][j]= DIFFERENT_DIELECTRIC_ABOVE_AND_LEFT;
-
- else if ( ERb != er && ERl != er)
- oddity[i][j]= DIFFERENT_DIELECTRIC_BELOW_AND_LEFT;
-
- else if ( ERa != er )
- oddity[i][j]= DIFFERENT_DIELECTRIC_ABOVE;
-
- else if ( ERb != er )
- oddity[i][j]= DIFFERENT_DIELECTRIC_BELOW;
- else if ( ERl != er )
- oddity[i][j]= DIFFERENT_DIELECTRIC_LEFT;
-
- else if ( ERr != er )
- oddity[i][j]= DIFFERENT_DIELECTRIC_RIGHT;
-*/
- else
- oddity[i][j]= ORDINARY_INTERIOR_POINT;
- }/* end of for i=0 to width-1 */
- } /* end of for j= 0 to height-1 */
-}
-