git.gag.com Git - debian/atlc/blob - src/non_gui/create_bmp_for_circ_in_circ.c

   1 /* atlc - arbitrary transmission line calculator, for the analysis of
   2 transmission lines are directional couplers.
   3
   4 Copyright (C) 2002. Dr. David Kirkby, PhD (G8WRB).
   5
   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.
  10
  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.
  15
  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,
  19 USA.
  20
  21 Dr. David Kirkby, e-mail drkirkby at ntlworld.com
  22
  23 */
  24 #include "config.h"
  25
  26 #ifdef HAVE_STDLIB_H
  27 #include <stdlib.h>
  28 #endif
  29
  30 #ifdef HAVE_MATH_H
  31 #include <math.h>
  32 #endif
  33
  34 #ifdef HAVE_STRING_H
  35 #include <string.h>
  36 #endif
  37
  38 #ifdef HAVE_STRINGS_H
  39 #include <strings.h>
  40 #endif
  41
  42 #include "exit_codes.h"
  43 #include "definitions.h"
  44
  45 #ifdef WINDOWS
  46 #pragma hrdstop
  47 #include <condefs.h>
  48 #endif
  49
  50 double aa, bb, cc, dd, hh, ww, DD, WW, HH, Er1, Er2;
  51 int a, b, c, d, h, w, D, W, H;
  52
  53 int verbose=FALSE;
  54 int lowest=50, highest=100;
  55
  56 /* The following are only needed so that I can link in both
  57 convert_create_bmp_for_circ_in_circ_dimensions_to_integers.c
  58 and convert_create_bmp_for_rect_in_rect_dimensions_to_integers.c
  59 */
  60
  61
  62 int main(int argc, char **argv) /* Read parameters from command line here   */
  63 {
  64    int  bmp_size=DEFAULT_BMP_SIZE;
  65    int q;
  66    FILE *image_data_fp;
  67    char *filename;
  68    struct transmission_line_properties not_used;
  69    double Zo, x;
  70    filename=string(0,1010);
  71
  72    /* The following line just prevents a warning from the compiler (SGI's to
  73    be precise). There is not reason to set W=1 */
  74    not_used.W=1;
  75
  76    while((q=get_options(argc,argv,"Cb:v")) != -1)
  77    switch (q)
  78    {
  79       case 'C':
  80         print_copyright((char *) "2002");
  81         exit_with_msg_and_exit_code("",OKAY);
  82       break;
  83       case 'b':
  84       bmp_size=atol(my_optarg);
  85       break;
  86       case 'v':
  87       verbose=TRUE;
  88       break;
  89       case '?':
  90       break;
  91    }
  92    if(argc-my_optind == 5)
  93    {
  94       DD=atof(argv[my_optind]);
  95       WW=HH=DD;
  96       dd=atof(argv[my_optind+1]);
  97       hh=atof(argv[my_optind+2]);
  98       Er1=atof(argv[my_optind+3]);
  99       if( dd > DD)
 100         exit_with_msg_and_exit_code("Error: The Inner conductor is larger than the outer conductor!!!", 1);
 101       else if (DD == dd)
 102         exit_with_msg_and_exit_code("Error: The inner and outer conductors are of the same size!!!", 1);
 103       else if (dd/2 + hh >= DD/2)
 104         exit_with_msg_and_exit_code("Error: The inner and outer conductors will touch!!!", 1);
 105       filename=strncpy(filename, argv[my_optind+4],1000);
 106       if( (image_data_fp=fopen(filename,"wb")) == NULL)
 107         exit_with_msg_and_exit_code("Can't open file in create_bmp_for_circ_in_circ.c", CANT_OPEN_FOR_WRITING);
 108       check_parameters_of_create_bmp_for_circ_in_circ();
 109       convert_create_bmp_for_circ_in_circ_dimensions_to_integers(bmp_size);
 110       write_bitmap(image_data_fp, not_used);
 111       x=(double) (d*d+D*D-4*h*h)/(2*D*d);
 112       Zo=59.9585*log(x+sqrt(x*x-1))/sqrt(Er1);
 113       Zo=log(x+sqrt(x*x-1))/(1000*sqrt(EPSILON_0)*sqrt(Er1)*sqrt(10*M_PI));
 114
 115       if(verbose == TRUE)
 116       {
 117         printf("DD=%f dd=%f hh=%f x=%f\n",DD, dd, hh, x);
 118         printf("D=%d d=%d h=%d x=%f\n",D, d, h, x);
 119         printf("Zo is theoretically %f Ohms\n",Zo);
 120       }
 121    }
 122    else
 123       usage_create_bmp_for_circ_in_circ();
 124    free_string(filename,0,1010);
 125    return(OKAY);
 126 }