+++ /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
-
-*/
-#include "config.h"
-
-
-#ifdef HAVE_STDLIB_H
-#include <stdlib.h>
-#endif
-
-#include "definitions.h"
-#include "exit_codes.h"
-#include "Erdata.h"
-
-/* We need to set up the permittivity and cell_type arrays. A
-complication arrises if a diaectric is found, that is not vacuum.
-If it is the only dielectric found, we calaculate C by saying C=C*Er. If
-there are more than one dielectric, we need to calculate it twice. */
-
-extern struct pixels Er_in_bitmap[MAX_DIFFERENT_PERMITTIVITIES];
-extern struct pixels Er_on_command_line[MAX_DIFFERENT_PERMITTIVITIES];
-extern int width, height;
-
-extern double **Vij, **Er;
-extern unsigned char *image_data;
-extern unsigned char **cell_type;
-extern int non_vacuum_found;
-extern int coupler;
-extern char *inputfile_name;
-
-void setup_arrays(struct transmission_line_properties *data)
-{
- struct pixels pixels_found;
- int conductor_found;
- int conductors = 0;
- int w,h, offset=-3, colour_mixture, i, z;
- unsigned char red, green, blue;
- int dielectric_found;
- int new_colour_in_image;
- data->dielectrics_in_bitmap=0;
- pixels_found.red=0;
- pixels_found.green=0;
- pixels_found.blue=0;
- pixels_found.white=0;
- pixels_found.other_colour=0;
-
- for(h=0;h<height;h++)
- {
- for (w=0; w<width;++w)
- {
- cell_type[w][height-1-h]=0;
- dielectric_found=FALSE;
- conductor_found=FALSE;
- offset+=3;
- if((w==0) && (offset%4!=0) && (h!=0))
- offset++;
- if((w==0) && (offset%4!=0) && (h!=0))
- offset++;
- if((w==0) && (offset%4!=0) && (h!=0))
- offset++;
- blue=image_data[offset];
- green=image_data[offset+1];
- red=image_data[offset+2];
- colour_mixture=256*256*red+256*green+blue;
- if( colour_mixture == 0xff0000 ) /* +1V red */
- {
- cell_type[w][height-1-h]=CONDUCTOR_PLUS_ONE_V;
- Vij[w][height-1-h]=1.0;
- Er[w][height-1-h]=METAL_ER;
- conductor_found=TRUE;
- pixels_found.red++;
- }
- else if( colour_mixture == 0x00ff00 ) /* 0v green */
- {
- cell_type[w][height-1-h]=CONDUCTOR_ZERO_V;
- Vij[w][height-1-h]=0.0;
- conductor_found=TRUE;
- pixels_found.green++;
- Er[w][height-1-h]=METAL_ER;
- }
- else if( colour_mixture == 0x0000ff ) /* -1V blue */
- {
- cell_type[w][height-1-h]=CONDUCTOR_MINUS_ONE_V;
- Vij[w][height-1-h]=-1.0;
- conductor_found=TRUE;
- pixels_found.blue++;
- Er[w][height-1-h]=METAL_ER;
- coupler=TRUE;
- data->couplerQ=TRUE;
- }
- else /* A dielectric */
- {
- if(colour_mixture == 0xffffff ) /* White */
- pixels_found.white++; /* Vacuum */
- else
- pixels_found.other_colour++; /* Some other dielectric */
- cell_type[w][height-1-h]=DIELECTRIC;
- Vij[w][height-1-h]=0.0; /* XXXXXXXXXXx */
- for(z=0;z<NUMBER_OF_DIELECTRICS_DEFINED;++z)
- {
- /* Check to see if the colour is one of the 10 known
- about, without any need to define on the command line
- */
- if (colour_mixture == colours[z])
- {
- Er[w][height-1-h]=Ers[z];
- dielectric_found=TRUE;
- if(z != 0)
- {
- non_vacuum_found=TRUE;
- data->found_this_dielectric=Ers[z];
- }
- }
- }
- for(i=0;i<data->dielectrics_on_command_line;++i)
- {
- if (Er_on_command_line[i].other_colour == colour_mixture)
- {
- Er[w][height-1-h]=Er_on_command_line[i].epsilon;
- dielectric_found=TRUE;
- data->found_this_dielectric=Er_on_command_line[i].epsilon;
- non_vacuum_found=TRUE;
- }
- }
- }
- if((dielectric_found == FALSE) && (conductor_found==FALSE))
- {
- fprintf(stderr,"Error#7: The colour r=0x%x g=0x%x b=0x%x (0x%02x%02x%02x) exists at pixel %d,%d.\n",red,green,blue,red,green,blue,w,h);
- fprintf(stderr,"However, atlc does not know how to interpret this colour. This is not a\n");
- fprintf(stderr,"conductor (pure red, green or blue), nor is it one of the %d dielectrics that\n",NUMBER_OF_DIELECTRICS_DEFINED);
- fprintf(stderr,"are predefined in Erdata.h, nor is a corresponding dielectric constant defined\n");
- fprintf(stderr,"on the command line with the -d option. Sometimes this occurs when a\n");
- fprintf(stderr,"graphics package is used to make the bitmap, but it performs some form of\n");
- fprintf(stderr,"anti-aliasing or smooting. If this is the case, redraw the image turning such\n");
- fprintf(stderr,"options off. If this is not the case then re-run atlc adding the -d option\n");
- fprintf(stderr,"to define the relative permittivity of the dielectric\n\n");
- fprintf(stderr,"e.g. atlc -d %02x%02x%02x=1.9 %s \n\n",red,green,blue, inputfile_name);
- fprintf(stderr,"if this colour has a permittivity of 1.9. If there are multiple undefined\n");
- fprintf(stderr,"dielectrics, then there will need to be multiple copies of the -d option given.\n");
- exit_with_msg_and_exit_code("",UNDEFINED_COLOUR_IN_BITMAP);
- }
- /* We need to keep a record of the number of dielectrics in the image,
- and determine if they are defined on the command line, or if they are
- defined on in the header file. */
- if (dielectric_found == TRUE)
- {
- new_colour_in_image=TRUE;
- for (i=0; i< data->dielectrics_in_bitmap; ++i)
- {
- if (Er_in_bitmap[i].other_colour == colour_mixture) /* a known colour */
- {
- new_colour_in_image=FALSE;
- if(colour_mixture != 0xffffff)
- {
- non_vacuum_found=TRUE;
- }
- }
- }
- if(new_colour_in_image==TRUE)
- {
- (data->dielectrics_in_bitmap)++;
- Er_in_bitmap[i].other_colour=colour_mixture;
- Er_in_bitmap[i].red=red;
- Er_in_bitmap[i].green=green;
- Er_in_bitmap[i].blue=blue;
- Er_in_bitmap[i].epsilon=Er[w][height-1-h];
- }
- } /* end of if dielctric found */
- } /* end of for w */
- } /*end of for h */
- /* The following prints a lot of data that is generally not wanted
- but is when finding statistics of performance etc. */
- if(data->verbose_level >=3)
- {
- printf("Red (+1 V conductor) pixels found = %8d \n",pixels_found.red);
- printf("Green (0 V conductor) pixels found = %8d \n",pixels_found.green);
- printf("Blue (-1 V conductor) pixels found = %8d \n",pixels_found.blue);
- printf("White (vacuum dielectric) pixels found = %8d \n",pixels_found.white);
- printf("Others (not vacuum dielectic) pixels found = %8d \n",pixels_found.other_colour);
- printf("Width = %8d \n",width);
- printf("Height = %8d \n",height);
- printf("Pixels = %8d \n",width*height);
- printf("Number of Dielectrics found = %8d \n", data->dielectrics_in_bitmap);
- data->non_metallic_pixels=width*height-pixels_found.red-pixels_found.green -pixels_found.blue;
- printf("Number of non-metallic pixels = %8d \n", data->non_metallic_pixels);
- printf("filename = %30s \n",inputfile_name);
- if (pixels_found.red > 0 )
- conductors+=1;
- if (pixels_found.green > 0 )
- conductors+=1;
- if (pixels_found.blue > 0 )
- conductors+=1;
- printf("Number of Conductors = %d \n", conductors);
- }
- /* The following should not be necessary, but may be as a test */
- /* I'd like to Miguel Berg for noticcing a servere bug, where the
- indeces of w and h were transposed, leading to crashes on Windoze
- XP */
- for(h=0;h<height;h++)
- {
- for (w=0; w<width;++w)
- {
- if((Vij[w][h] > 1.0) || (Vij[w][h]<-1.0))
- {
- fprintf(stderr,"Sorry, something is wrong Vij[%d][%d]=%f in %s %d\n",w,h,Vij[w][h], __FILE__,__LINE__);
- exit_with_msg_and_exit_code("Exiting ....",VOLTAGE_OUT_OF_RANGE);
- }
- }
- }
- /* Check two conductors and not next to each other, creating a short */
- check_for_shorts();
-} /* end of function */