c_ggrid.h

/*
  Copyright (C) 2004-2015  Timothy C.A. Molteno
  
  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 version 2 of the License, or
  (at your option) any later 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
*/
#pragma once
 
#include "math_util.h"
#include "common.h"
#include "misc.h"
#include "c_evlcom.h"
 
using namespace std;
 
 
class c_ggrid
{
public:
  static int m_nxa[3], m_nya[3];
  static nec_float m_dxa[3], m_dya[3];
  static nec_float m_xsa[3], m_ysa[3];
  
  nec_complex m_epscf;
  complex_array m_ar1, m_ar2, m_ar3;
  
  c_evlcom m_evlcom;
  
  void initialize()  {
    m_ar1.resize(11*10*4);
    m_ar2.resize(17*5*4);
    m_ar3.resize(9*8*4);
 
    // initialize interpolation cache
    m_ip_ixs = -10;
    m_ip_iys = -10;
    m_ip_igrs = -10;
    m_ip_ixeg = 0;
    m_ip_iyeg = 0;
    m_ip_dx = 1.0;
    m_ip_dy = 1.0;
    m_ip_xs = 0.0;
    m_ip_ys = 0.0;
  }
 
  void interpolate( nec_float x, nec_float y, 
    nec_complex *f1, nec_complex *f2,
    nec_complex *f3, nec_complex *f4 );
 
  void sommerfeld( nec_float epr, nec_float sig, nec_float wavelength );
 
private:
  // interpolation cache (per-instance, thread-safe)
  int m_ip_ix, m_ip_iy, m_ip_ixs, m_ip_iys, m_ip_igrs, m_ip_ixeg, m_ip_iyeg;
  int m_ip_nxm2, m_ip_nym2, m_ip_nxms, m_ip_nyms, m_ip_nd, m_ip_ndp;
  nec_float m_ip_dx, m_ip_dy, m_ip_xs, m_ip_ys, m_ip_xz, m_ip_yz;
  nec_complex m_ip_a[4][4], m_ip_b[4][4], m_ip_c[4][4], m_ip_d[4][4];
};
 
class c_ground_wave
{
public:
  /* common  /gwav/ */
  nec_float r1, r2, zmh, zph;
  nec_complex u, u2, xx1, xx2;
  
  void set_u(nec_complex in_u)  {
    u = in_u;
    u2 = u * u;
  }
};
 
void gwave( nec_complex& erv, nec_complex& ezv,
  nec_complex& erh, nec_complex& ezh, nec_complex& eph,
  c_ground_wave& ground_wave);