--- /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 gmail.com
+
+*/
+
+#include "config.h"
+
+#ifdef HAVE_STRING_H
+#include <string.h>
+#endif
+
+#ifdef HAVE_STRINGS_H
+#include <strings.h>
+#endif
+
+#ifdef HAVE_STDLIB_H
+#include <stdlib.h>
+#endif
+#include "definitions.h"
+
+/* names, colours and Ers are all arrays of 10. It would be better they were
+in a structure as they are all linked closely, but they are not and I
+can't be bothered to change it */
+
+extern double DD, dd, hh, Er1;
+extern int D, d, h, W, H;
+extern int lowest;
+extern int highest;
+extern int verbose;
+
+void convert_create_bmp_for_circ_in_circ_dimensions_to_integers(int accuracy_level)
+{
+ int best_D=-1, best_d=-1, best_h=-1;
+ int min, max;
+ int i;
+ double max_gridsize, min_gridsize;
+ double error, error_min=VERY_LARGE, gridsize, best_grid_size=-1;
+
+ /* scale grid size, according to a command line option */
+ lowest*=(int) 0.5+sqrt(pow(2.0,(double)accuracy_level));
+ highest*=(int) 0.5+sqrt(pow(2.0,(double)accuracy_level));
+
+ max_gridsize=DD/lowest; /* minimum dimension in m */
+ min_gridsize=DD/highest;
+ min=DD/max_gridsize;
+ max=DD/min_gridsize;
+ for(i=min;i<=max;++i)
+ {
+ /* Try various combinations for D,d,o */
+ gridsize=DD/i;
+ D=(int)(DD/gridsize + 0.5);
+ d=(int)(dd/gridsize + 0.5);
+ h=(int)(hh/gridsize + 0.5);
+ error=0.0;
+ error+=pow((DD-D*gridsize)/DD,2.0);
+ error+=pow((dd-d*gridsize)/dd,2.0);
+ if(h>0)
+ error+=pow((hh-h*gridsize)/hh,2.0);
+ if(error < (error_min-TINY))
+ {
+ error_min=error;
+ best_D=D;
+ best_d=d;
+ best_h=h;
+ best_grid_size=gridsize;
+ }
+ }
+ D=W=H=best_D;
+ d=best_d;
+ h=best_h;
+ if(verbose==TRUE)
+ {
+ fprintf(stderr,"error_min=%.16f\n",error_min);
+ fprintf(stderr,"User requested: DD=%f dd=%f hh=%f Er1=%f\n\n",DD,dd,hh,Er1);
+ fprintf(stderr,"Internally the program is using the following grid:\n");
+ fprintf(stderr,"D=%d d=%d h=%d \nThe grid size is %f mm, inches or whatever\n\n", D,d,h,best_grid_size);
+ if(error_min > TINY)
+ {
+ fprintf(stderr,"This means we are simulating a transmission line with these dimensions:\n");
+ fprintf(stderr,"D=%f d=%f h=%f (mm, inches or whatever)\n",D*best_grid_size,d*best_grid_size,h*best_grid_size);
+ fprintf(stderr,"\nThese are slightly different to what you indicated on the command line,\n");
+ fprintf(stderr,"but they are the best approximation possible, given the grid size\n");
+ }
+ }
+ check_error(DD,D,best_grid_size,"D");
+ check_error(dd,d,best_grid_size,"d");
+ if(h > 0.0)
+ check_error(hh,h,best_grid_size,"h");
+}
+