X-Git-Url: https://git.gag.com/?p=debian%2Fatlc;a=blobdiff_plain;f=src%2Fwrite_fields_for_two_conductor_lines.c;fp=src%2Fwrite_fields_for_two_conductor_lines.c;h=c41b4538da88e9ce5839c3b8e2f91935cc6e991f;hp=0000000000000000000000000000000000000000;hb=6b66971c81ba66fc706c04decfa5a69a19e57c4d;hpb=a51f9bc3e458756855df5ab51a7e7ec015a9aebd

diff --git a/src/write_fields_for_two_conductor_lines.c b/src/write_fields_for_two_conductor_lines.c
new file mode 100644
index 0000000..c41b453
--- /dev/null
+++ b/src/write_fields_for_two_conductor_lines.c
@@ -0,0 +1,267 @@
+/* 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"
+
+#include "definitions.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
+
+#ifdef HAVE_STDIO_H
+#include <stdio.h>
+#endif
+
+#include "exit_codes.h"
+
+extern double **Vij;
+extern double **Er;
+extern unsigned char *bmp_buff;
+extern int width, height, errno;
+#ifdef WRITE_ODDITY_DATA
+extern unsigned char **oddity;
+#endif
+
+/* Write the following files, assuming an input of example.bmp 
+where example.bmp is a 2 conductor transmission lines. For
+3 conductor transmission lines (couplers) the function
+
+
+example.E.bmp   Grayscale Bitmap of |E-field|, normallised to 1,
+but corrected for Gamma
+example.Ex.bmp  Colour Bitmap of x-directed E-field, normallised to 1, 
+but corrected for Gamma
+example.Ey.bmp  Colour Bitmap of y-directed E-field, normallised to 1, 
+but corrected for Gamma
+example.V.bmp   Colour Bitmap of Voltage field, normallised to 1, but 
+corrected for Gamma
+eexample.U.bmp  Grayscale bitmap, with just the energy (U=CV^2).
+
+example.Ex.bin  binary file, with just the x-directed E-field 
+(in volts/pixel) as doubles 
+example.Ey.bin  binary file, with just the y-directed E-field 
+(in volts/pixel) as doubles 
+example.E.bin   binary file, with just the E-field {sqrt(Ex^2+Ey^2)} 
+(in volts/pixel) as doubles 
+example.V.bin   binary file, with just the Voltage as doubles 
+eexample.U.bin  binary file, with just the energy (U=CV^2).
+
+*/
+
+extern double image_fiddle_factor;
+
+
+void write_fields_for_two_conductor_lines(char * filename, struct transmission_line_properties data, size_t size)
+{
+  FILE *Ex_bin_fp=NULL, *Ey_bin_fp=NULL;
+  FILE *E_bin_fp=NULL, *V_bin_fp, *U_bin_fp=NULL;
+  FILE *Ex_bmp_fp=NULL, *Ey_bmp_fp=NULL;
+  FILE *E_bmp_fp=NULL, *V_bmp_fp=NULL, *U_bmp_fp=NULL;
+#ifdef WRITE_ODDITY_DATA
+  FILE *oddity_bmp_fp=NULL;
+  double odd;
+#endif
+  unsigned char r, g, b;
+  FILE *permittivity_bin_fp=NULL, *permittivity_bmp_fp=NULL;
+  struct max_values maximum_values;
+  int offset=-3, w, h;
+  double V, E, Ex, Ey, U;
+
+  unsigned char *image_data_Ex=NULL; 
+  unsigned char *image_data_Ey=NULL;
+  unsigned char *image_data_E=NULL;
+  unsigned char *image_data_U=NULL; 
+  unsigned char *image_data_V=NULL;
+  unsigned char *image_data_Er=NULL;
+#ifdef WRITE_ODDITY_DATA
+  unsigned char *image_data_oddity=NULL;
+#endif 
+
+
+  if(data.write_binary_field_imagesQ==TRUE)
+  {
+    Ex_bin_fp=get_file_pointer_with_right_filename(filename,".Ex.bin");
+    Ey_bin_fp=get_file_pointer_with_right_filename(filename,".Ey.bin");
+    E_bin_fp=get_file_pointer_with_right_filename(filename,".E.bin");
+    V_bin_fp=get_file_pointer_with_right_filename(filename,".V.bin");
+    U_bin_fp=get_file_pointer_with_right_filename(filename,".U.bin");
+    permittivity_bin_fp=get_file_pointer_with_right_filename(filename,".Er.bin");
+    for(h=height-1;h>=0;h--)
+    {
+      for(w=0;w<width;++w)
+      {
+	V=Vij[w][h];
+        Ex=find_Ex(w,h);
+        Ey=find_Ey(w,h);
+        E=find_E(w,h); 
+        U=find_energy_per_metre(w,h);
+        if( fwrite((void *) &Ex,sizeof(double), 1, Ex_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#1: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+        if( fwrite((void *) &Ey,sizeof(double), 1, Ey_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#2: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+        if( fwrite((void *) &E,sizeof(double), 1, E_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#3: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+        if( fwrite((void *) &V,sizeof(double), 1, V_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#4: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+        if( fwrite((void *) &U,sizeof(double), 1, U_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#5: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+        if( fwrite((void *) &Er[w][h],sizeof(double), 1, permittivity_bin_fp) != 1)
+          exit_with_msg_and_exit_code("Error#6: Failed to write binary file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+      }
+    }
+  } /* end of writing binary files for 2 conductor lines */
+
+  if(data.write_bitmap_field_imagesQ==TRUE)
+  {
+    find_maximum_values(&(maximum_values), ZERO_ELEMENTS_FIRST); /* sets stucture maximum_values */
+
+    /* Allocate ram to store the bitmaps before they are written to disk */
+    image_data_Ex=ustring(0L,(long) size);
+    image_data_Ey=ustring(0L,(long) size);
+    image_data_E=ustring(0L,(long) size);
+    image_data_V=ustring(0L,(long) size);
+    image_data_Er=ustring(0L,(long) size);
+    image_data_U=ustring(0L,(long) size);
+#ifdef WRITE_ODDITY_DATA
+    image_data_oddity=ustring(0L,(long) size);
+#endif
+
+    Ex_bmp_fp=get_file_pointer_with_right_filename(filename,".Ex.bmp");
+    Ey_bmp_fp=get_file_pointer_with_right_filename(filename,".Ey.bmp");
+    E_bmp_fp=get_file_pointer_with_right_filename(filename,".E.bmp");
+    V_bmp_fp=get_file_pointer_with_right_filename(filename,".V.bmp");
+    U_bmp_fp=get_file_pointer_with_right_filename(filename,".U.bmp");
+    permittivity_bmp_fp=get_file_pointer_with_right_filename(filename,".Er.bmp");
+#ifdef WRITE_ODDITY_DATA
+    oddity_bmp_fp=get_file_pointer_with_right_filename(filename,".oddity.bmp");
+#endif
+
+    if( fwrite(bmp_buff,0x36,1,Ex_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer Ex_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fwrite(bmp_buff,0x36,1,Ey_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer Ey_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fwrite(bmp_buff,0x36,1,E_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer E_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fwrite(bmp_buff,0x36,1,V_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer V_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fwrite(bmp_buff,0x36,1,U_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer U_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fwrite(bmp_buff,0x36,1,permittivity_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer permittivity_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+
+#ifdef WRITE_ODDITY_DATA
+    if( fwrite(bmp_buff,0x36,1,oddity_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("failed to write file with the pointer oddity_bmp_fp in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+#endif 
+
+    offset=-3;
+    for(h=height-1;h>=0;h--)
+    {
+      for(w=0;w<width;++w)
+      {
+	offset+=3;
+	if((w==0) && (offset%4!=0) && (h!=0)) 
+	  offset+=(4-offset%4);
+        Ex=find_Ex(w,h);
+        Ey=find_Ey(w,h);
+        E=find_E(w,h); 
+        U=find_energy_per_metre(w,h);
+        calculate_colour_data(Ex, maximum_values.Ex_or_Ey_max, w, h, offset,image_data_Ex, COLOUR,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+        calculate_colour_data(Ey, maximum_values.Ex_or_Ey_max, w, h, offset,image_data_Ey, COLOUR,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+        calculate_colour_data(E, maximum_values.E_max, w, h, offset,image_data_E, MONOCHROME,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+        calculate_colour_data(U, maximum_values.U_max, w, h, offset,image_data_U, MONOCHROME,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+        calculate_colour_data(Vij[w][h], maximum_values.V_max, w, h, offset,image_data_V, COLOUR,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+        calculate_colour_data(Er[w][h], MAX_ER, w, h, offset,image_data_Er, MIXED,&r,&g,&b,IMAGE_FIDDLE_FACTOR);
+#ifdef WRITE_ODDITY_DATA
+	odd=(double) oddity[w][h];
+        calculate_colour_data(odd, 255 , w, h, offset,image_data_oddity, MONOCHROME,&r,&g,&b,1.0);
+#endif
+
+      }
+    }
+
+    /* if( fwrite((void *) &(image_data_Ex[0]),size, 1, Ex_bmp_fp) != 1) */
+    if( fwrite((void *) image_data_Ex,size, 1, Ex_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#7: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(Ex_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#8: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+    /* if( fwrite((void *) &(image_data_Ey[0]),size, 1, Ey_bmp_fp) != 1) */
+    if( fwrite((void *) image_data_Ey,size, 1, Ey_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#9: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(Ey_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#10: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+    if( fwrite((void *) &(image_data_E[0]),size, 1, E_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#11: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(E_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#12: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+    if( fwrite((void *) &(image_data_U[0]),size, 1, U_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#13: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(U_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#14: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+    if( fwrite((void *) &(image_data_V[0]),size, 1, V_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#15: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(V_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#16: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+    if( fwrite((void *) &(image_data_Er[0]),size, 1, permittivity_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#17: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(permittivity_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#18: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+
+#ifdef WRITE_ODDITY_DATA
+    if( fwrite((void *) &(image_data_oddity[0]),size, 1, oddity_bmp_fp) != 1)
+      exit_with_msg_and_exit_code("Error#19: Failed to write bitmap file in write_fields_for_two_conductor_lines.c",WRITE_FAILURE);
+    if( fclose(oddity_bmp_fp) != 0)
+      exit_with_msg_and_exit_code("Error#20: Unable to close file in write_fields_for_two_conductor_lines.c",CANT_CLOSE_FILE);
+#endif
+
+    free_ustring(image_data_Ex,0L,(long)size);
+    free_ustring(image_data_Ey,0L,(long)size);
+    free_ustring(image_data_E,0L,(long)size);
+    free_ustring(image_data_V,0L,(long)size);
+    free_ustring(image_data_U,0L,(long)size);
+    free_ustring(image_data_Er,0L,(long)size);
+#ifdef WRITE_ODDITY_DATA
+    free_ustring(image_data_oddity,0L,(long)size);
+#endif
+  }
+}