Merge commit 'upstream/4.6.1'

[debian/atlc] / src / swap_conductor_voltages.c

diff --git a/src/swap_conductor_voltages.c b/src/swap_conductor_voltages.c
new file mode 100644 (file)
index 0000000..0cec911
--- /dev/null
+++ b/src/swap_conductor_voltages.c
@@ -0,0 +1,61 @@
+/* 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 
+
+*/
+
+/* When calculationg the properties of couplers it is necessary to swap
+around the voltages on the conductors. First we consider the voltages as
+present, which will be positive and negative. Next the negative pixels
+must be changed to positive, to calculate the even mode impedance. This
+routine does the swapping around, based on the value of the parameter
+oddity, which will be set to CONDUCTOR_PLUS_ONE_V or
+
+CONDUCTOR_MINUS_ONE_V in definitions.h */
+
+#include "config.h"
+
+#include "definitions.h"
+
+extern int width, height;
+unsigned char **oddity;
+extern double **Vij;
+
+/* This swaps any conductor orignally at -1 V to +1 V. Any
+dielectric is set to zero volts. This is used in computation
+of couplers. */
+
+void swap_conductor_voltages()
+{
+  int i, j;
+  int pixel;
+  for(i=0;i<width;++i)
+  {
+    for(j=0;j<height;++j)
+    {
+      pixel=oddity[i][j];
+      if(pixel==CONDUCTOR_MINUS_ONE_V){
+        Vij[i][j]=+1.0;
+       oddity[i][j]=CONDUCTOR_PLUS_ONE_V;
+      }
+    }
+  }
+}