1 /* atlc - arbitrary transmission line calculator, for the analysis of
2 transmission lines are directional couplers.
4 Copyright (C) 2002. Dr. David Kirkby, PhD (G8WRB).
6 This program is free software; you can redistribute it and/or
7 modify it under the terms of the GNU General Public License
8 as published by the Free Software Foundation; either package_version 2
9 of the License, or (at your option) any later package_version.
11 This program is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
16 You should have received a copy of the GNU General Public License
17 along with this program; if not, write to the Free Software
18 Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
21 Dr. David Kirkby, e-mail drkirkby at ntlworld.com
38 #include "definitions.h"
40 /* names, colours and Ers are all arrays of 10. It would be better they were
41 in a structure as they are all linked closely, but they are not and I
42 can't be bothered to change it */
44 extern double DD, ww, hh, xx, yy, Er1;
45 extern int D, w, h, x, y;
47 extern int highest, verbose;
49 void convert_create_bmp_for_rect_in_circ_dimensions_to_integers(int accuracy_level)
51 int best_D=-1, best_w=-1, best_h=-1, best_x=-1, best_y=-1;
54 double error, error_min=VERY_LARGE, gridsize, best_grid_size=-1;
56 /* scale grid size, according to a command line option */
57 lowest*=(int) 0.5+pow(2.0,(double) accuracy_level);
58 highest*=(int) 0.5+pow(2.0,(double) accuracy_level);
62 /* Try various combinations for D, w, h, x, & y etc */
64 D=(int)(DD/gridsize + 0.5);
65 w=(int)(ww/gridsize + 0.5);
66 h=(int)(hh/gridsize + 0.5);
67 x=(int)(xx/gridsize + 0.5);
68 y=(int)(yy/gridsize + 0.5);
70 error+=pow((DD-D*gridsize)/DD,2.0); /* relative error in D */
71 error+=pow((ww-w*gridsize)/ww,2.0);
72 error+=pow((hh-h*gridsize)/hh,2.0);
74 error+=pow((xx-x*gridsize)/xx,2.0);
76 error+=pow((yy-y*gridsize)/yy,2.0);
77 if(error < error_min-TINY)
85 best_grid_size=gridsize;
95 fprintf(stderr,"error_min=%.16f\n",error_min);
96 fprintf(stderr,"User requested: D=%f w=%f h=%f x=%f y=%f Er1=%f\n\n",DD,ww,hh,xx,yy,Er1);
97 fprintf(stderr,"Internally the program is using the following grid:\n");
98 fprintf(stderr,"D=%d w=%d h=%d x=%d y=%d \n", D, w, h, x, y);
99 fprintf(stderr,"which equates to dimensions of D=%f w=%f h=%f x=%f y=%f\n", D*gridsize, w*gridsize, h*gridsize, x*gridsize,
103 fprintf(stderr,"This means we are simulating a transmission line with these dimensions:\n");
104 fprintf(stderr,"W=%f w=%f h=%f x=%f y=%f (mm, inches or whatever)\n",D*best_grid_size,w*best_grid_size,h*best_grid_size,x*best_grid_size,y*best_grid_size);
105 fprintf(stderr,"\nThese are slightly different to what you indicated on the command line,\n");
106 fprintf(stderr,"but they are the best approximation possible, given the grid size\n");
109 check_error(DD,D,best_grid_size,"D");
110 check_error(ww,w,best_grid_size,"w");
111 check_error(hh,h,best_grid_size,"h");
112 check_error(xx,x,best_grid_size,"x");
113 check_error(yy,y,best_grid_size,"y");