c_ggrid.h

/*
        Copyright (C) 2004-2006  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
*/
#ifndef __c_ggrid__
#define __c_ggrid__

#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);
        }

        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 freq_mhz );
        
};

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);


#endif