MoM example.

docs/source/comparisons_numerical.rst

@@ -7,20 +7,20 @@ This section presents some comparisons of models using different numerical model
 FDTD/MoM
 ========
 
-The Finite-Difference Time-Domain (FDTD) technique using gprMax is compared with the Method of Moments (MoM) using the MATLAB antenna toolbox (http://uk.mathworks.com/products/antenna/).
+The Finite-Difference Time-Domain (FDTD) method from gprMax is compared with the Method of Moments (MoM) from the MATLAB antenna toolbox (http://uk.mathworks.com/products/antenna/).
 
 Bowtie antenna in free space
 ----------------------------
 
-:download:`hertzian_dipole_fs.in <../../tests/numerical/vs_MoM_MATLAB/antenna_bowtie_fs/antenna_bowtie_fs.in>`
+:download:`antenna_bowtie_fs.in <../../tests/numerical/vs_MoM_MATLAB/antenna_bowtie_fs/antenna_bowtie_fs.in>`
 
-This example considers the input impedance of a planar bowtie antenna in free space. The length and height of the bowtie are 100mm, giving a flare angle of 90$^\circ$.
+This example considers the input impedance of a planar bowtie antenna in free space. The length and height of the bowtie are 100mm, giving a flare angle of :math:`90^\circ`.
 
 .. literalinclude:: ../../tests/numerical/vs_MoM_MATLAB/antenna_bowtie_fs/antenna_bowtie_fs.in
     :language: none
     :linenos:
 
-The MoM was created in MATLAB using the ``bowtieTriangular`` class:
+For the MoM, the bowtie antenna was created in MATLAB using the ``bowtieTriangular`` class:
 
 .. code-block:: matlab
 

tests/numerical/vs_MoM_MATLAB/antenna_bowtie_fs/antenna_bowtie_fs.in

@@ -0,0 +1,36 @@
+#python:
+
+from gprMax.input_cmd_funcs import *
+
+title = 'antenna_bowtie_fs'
+print('#title: {}'.format(title))
+
+domain = domain(0.200, 0.120, 0.120)
+dxdydz = dx_dy_dz(0.001, 0.001, 0.001)
+time_window = time_window(30e-9)
+bowtie_dims = (0.050, 0.100) # Length, height
+tx_pos = (domain[0]/2, domain[1]/2, 0.020)
+
+# Source excitation and type
+print('#waveform: gaussian 1 1.5e9 mypulse')
+print('#transmission_line: x {:g} {:g} {:g} 50 mypulse'.format(tx_pos[0], tx_pos[1], tx_pos[2]))
+
+# Output point
+print('#rx: {:g} {:g} {:g}'.format(tx_pos[0], tx_pos[1], tx_pos[2] + 0.020))
+
+# Bowtie - upper x half
+triangle(tx_pos[0] + dxdydz[0], tx_pos[1] - dxdydz[1], tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] - bowtie_dims[1]/2 - dxdydz[1], tx_pos[2], tx_pos[0] + bowtie_dims[0] + dxdydz[0], tx_pos[1] + bowtie_dims[1]/2 - dxdydz[1], tx_pos[2], 0, 'pec')
+
+# Bowtie - lower x half
+triangle(tx_pos[0], tx_pos[1], tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] - bowtie_dims[1]/2, tx_pos[2], tx_pos[0] - bowtie_dims[0], tx_pos[1] + bowtie_dims[1]/2, tx_pos[2], 0, 'pec')
+
+# Uncomment to temporarily check location of source
+#edge(tx_pos[0], tx_pos[1], tx_pos[2], tx_pos[0] + dxdydz[0], tx_pos[1], tx_pos[2], 'pec')
+
+# Geometry view of entire domain
+geometry_view(0, 0, 0, domain[0], domain[1], domain[2], dxdydz[0], dxdydz[1], dxdydz[2], title)
+
+# Detailed geometry view of PCB and bowties
+geometry_view(tx_pos[0] - bowtie_dims[0] - 2*dxdydz[0], tx_pos[1] - bowtie_dims[1]/2  - 2*dxdydz[1], tx_pos[2] - 2*dxdydz[2], tx_pos[0] + bowtie_dims[0] + 2*dxdydz[0], tx_pos[1] + bowtie_dims[1]/2  + 2*dxdydz[1], tx_pos[2] + 2*dxdydz[2], dxdydz[0], dxdydz[1], dxdydz[2], title + '_pcb', type='f')
+
+#end_python: