cases.h

// rights - copyright by benedikt oswald, 2002-2003.
// rights - universitas ruperto carola heidelbergensis
// project - aqhdes2003
// file name - cases.h
// file type - C include file
// objective - install a specific problem into UG
// creation - 2003 feb 10 ~ 17:25:00 by benedikt oswald - implemented problem description
// modified - 2003 sep 17 ~ 21:33:00 by benedikt oswald - implemented PDE coefficient functions
// required software - UG 3.1 or higher
// rights for UG - cf. paper by Bastian et al., 1997
// file id -
// inheritance - none
// feature - implement simple benchmark problems to test the
// feature - implementation of the whitney discretization
// feature - of Maxwell's equations; additionally implements
// feature - coefficient functions for the material of the PDE to be solved
// feature - PROBLEM 0 : InitRectangularCavityResonanceProblem : calculates cavity resonances
// feature - PROBLEM 1 : InitCylindricalCavityResonanceProblem : calculates cavity resonances

/****************************************************************************/
/*                                                                          */
/* auto include mechanism and other include files                           */
/*                                                                          */
/****************************************************************************/

#ifndef __PROBLEMS__
#define __PROBLEMS__

#ifndef __COMPILER__
#include "compiler.h"
#endif

#ifndef __GM__
#include "gm.h"
#endif

// include standard C header files
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>

// include UG specific header files
#include "compiler.h"
#include "evm.h"         /* for V_COPY,M_COPY */
#include "gm.h"
#include "misc.h"
#include "cmdint.h"
#include "commands.h"    /* for GetCurrentMultigrid (Quickhack) */
#include "ugstruct.h"
#include "initug.h"
#include "np.h"          /* for NUMPROC implementation                                        */
#include "scan.h"        /* for ReadArgvVecDesc */
#include "std_domain.h"
#include "general.h"
#include "udm.h"         /* for VECDATA_DESC */
#include "ugdevices.h"
#include "shapes.h"

//  proprietary includes
#include "casesdim.h"
#include "linmaterials.h"

/****************************************************************************/
/*                                                                          */
/* defines in the following order                                           */
/*                                                                          */
/*        compile time constants defining static data size (i.e. arrays)    */
/*        other constants                                                   */
/*        macros                                                            */
/*                                                                          */
/****************************************************************************/

#define RECTANGULAR_CAVITY_ELECTROMANGETIC_RESONANCES 0
#define CYLINDRICAL_CAVITY_ELECTROMAGNETIC_RESONANCES 1

/****************************************************************************/
/*                                                                          */
/* data structures exported by the corresponding source file                */
/*                                                                          */
/****************************************************************************/

/****************************************************************************/
/*                                                                          */
/* data structures used in this source file (exported data structures are   */
/*        in the corresponding include file!)                               */
/*                                                                          */
/****************************************************************************/

/* declare NUMPROC class and methods for 'NP_DISCRETE_MATERIAL_LAYERS';
   which implements a simple model for the situation of a GPR antenna
   lying on soil which variable electrical material properties in the z
   direction, i.e. in the direction going into the soil; we use the
   parameters given below to describe the distribution of material
   properties in the air-soil configuration:
   
   - electric material properties vary in the z dimension only
   - z axis is directed as in cartesian coordinate system, i.e z=0 at bottom
   - fraction of interval filled with air material properties
   - use NUMPROC NP_LINEAR_MATERIALS, defined in linmaterials.h
*/

struct np_discrete_layers_z
{
    NP_LINEAR_MATERIALS linmat;  // uses NP_LINEARMATERIALS NUMPROC

    INT nl;                      // number of layers in z direction
    DOUBLE hz;                   // length of interval over which material properties vary

    DOUBLE *lthickness;          // tickness of layer in absolute dimensions
    DOUBLE *er;                  // array stores relative permittivity values of discrete layers
    DOUBLE *mr;                  // array stores relative permeability values of discrete layers
    DOUBLE *sigma;               // array stores ohmic conductivity of discrete layers
    DOUBLE *zintrfc;             // maximum z coordinate of a layer 
};
typedef struct np_discrete_layers_z NP_DISCRETE_LAYERS_Z;




/****************************************************************************/
/*                                                                          */
/* definition of exported global variables                                  */
/*                                                                          */
/****************************************************************************/

/****************************************************************************/
/*                                                                          */
/* function declarations                                                    */
/*                                                                          */
/****************************************************************************/

/* INT InitRectangularCavityResonanceProblem(void); */

INT InitProblems(void);

#endif

---