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
43 #include "definitions.h"
44 #include "exit_codes.h"
46 /* names, colours and Ers are all arrays of 10. It would be better they were
47 in a structure as they are all linked closely, but they are not and I
48 can't be bothered to change it */
50 extern double WW, HH, aa, bb, cc, dd, ww, hh, Er1, Er2;
51 extern int W, H, a, b, c, d, w, h;
56 void convert_create_bmp_for_rect_in_rect_dimensions_to_integers(int accuracy_level)
58 /* A number of paramters are set to -1, as otherwise the compiler
59 complaings they may be used unitialised. It is wise to check they
60 are not -1 before assining a variable to their contents,
61 but I can't be bothered. */
62 int best_W=-1, best_H=-1, best_a=-1, best_b=-1, best_c=-1;
63 int best_d=-1, best_w=-1, best_h=-1;
66 double max_gridsize, min_gridsize;
67 double error, error_min=VERY_LARGE, gridsize=-1, best_grid_size=-1;
69 /* scale grid size, according to a command line option */
70 /*lowest*=(int) 0.5+pow(2,accuracy_level); */
71 /*highest*=(int) 0.5+pow(2,accuracy_level); */
72 lowest*=1<<accuracy_level;
73 highest*=1<<accuracy_level;
75 max_gridsize=sqrt(WW*HH/(double )lowest); /* minimum dimension in m */
76 min_gridsize=sqrt(WW*HH/(double )highest);
89 /* Try various combinations for H,a,b,c, etc */
94 W=(int)(WW/gridsize + 0.5);
95 H=(int)(HH/gridsize + 0.5);
96 a=(int)(aa/gridsize + 0.5);
97 b=(int)(bb/gridsize + 0.5);
98 c=(int)(cc/gridsize + 0.5);
99 d=(int)(dd/gridsize + 0.5);
100 w=(int)(ww/gridsize + 0.5);
101 h=(int)(hh/gridsize + 0.5);
103 error+=pow((WW-W*gridsize)/WW,2.0); /* relative error in W */
104 error+=pow((HH-H*gridsize)/HH,2.0);
106 error+=pow((aa-a*gridsize)/aa,2.0);
108 error+=pow((bb-b*gridsize)/bb,2.0);
110 error+=pow((cc-c*gridsize)/cc,2.0);
112 error+=pow((dd-d*gridsize)/dd,2.0);
114 error+=pow((ww-w*gridsize)/ww,2.0);
116 error+=pow((hh-h*gridsize)/hh,2.0);
117 if(error < (error_min-TINY))
128 best_grid_size=gridsize;
141 fprintf(stderr,"error_min=%.16f\n",error_min);
142 fprintf(stderr,"User requested: WW=%f HH=%f a=%f b=%f c=%f\nd=%f w=%f h=%fEr1=%f Er2=%f\n\n",WW,HH,aa,bb,cc,dd,ww,hh,Er1,Er2);
143 fprintf(stderr,"Internally the program is using the following grid:\n");
144 fprintf(stderr,"W=%d H=%d a=%d b=%d c=%d d=%d w=%d h=%d. \nThe grid size is %f mm, inches or whatever\n\n", W, H, a, b, c, d, w, h,best_grid_size);
147 fprintf(stderr,"This means we are simulating a transmission line with these dimensions:\n");
148 fprintf(stderr,"W=%f H=%f a=%f b=%f c=%f d=%f\nw=%f h=%f (mm, inches or whatever)\n",W*best_grid_size,H*best_grid_size,a*best_grid_size,b*best_grid_size,c*best_grid_size,d*best_grid_size,w*best_grid_size,h*best_grid_size);
149 fprintf(stderr,"\nThese are slightly different to what you indicated on the command line,\n");
150 fprintf(stderr,"but they are the best approximation possible, given the grid size\n");
155 fprintf(stderr,"Error #12. The gap between the two conductors (b) is too small. Either increase the bitmap size (-b option), or make b at least %f\n",gridsize);
156 exit_with_msg_and_exit_code("",GAP_BETWEEN_CONDUCTORS_TOO_SMALL);
160 fprintf(stderr,"Error #13. The thickness of the inner conductor (c) is too small. Either increase the bitmap size (-b option), or make c at least %f\n",gridsize);
161 exit_with_msg_and_exit_code("",GAP_BETWEEN_CONDUCTORS_TOO_SMALL);
165 fprintf(stderr,"Error #14. The thickness of the dielectric of Er2 (h) is too small. Either increase the bitmap size (-b option), or make h at least %f\n",gridsize);
166 exit_with_msg_and_exit_code("",GAP_BETWEEN_CONDUCTORS_TOO_SMALL);
168 check_error(WW,W,best_grid_size,"W");
169 check_error(HH,H,best_grid_size,"H");
170 check_error(aa,a,best_grid_size,"a");
171 check_error(bb,b,best_grid_size,"b");
172 check_error(cc,c,best_grid_size,"c");
173 check_error(dd,d,best_grid_size,"d");
174 check_error(ww,w,best_grid_size,"w");
175 check_error(hh,h,best_grid_size,"h");