particleswarm.h

Go to the documentation of this file.

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

// include UG header files
#include "gm.h"                             /* for data structure                                                */
#include "shapes.h"                         /* general element shape calculations                                */
#include "ugstruct.h"                       /* for GetStringValue                                                */
#include "misc.h"                           /* for MIN, MAX, PI, ...                                             */
#include "ugdevices.h"                      /* for UserWrite, PrintErrorMessage                                  */ 
#include "commands.h"                       /* for GetCurrentMultigrid                                           */
#include "cmdint.h"                         /* for CreateCommand                                                 */
#include "cmdline.h"                        /* for command line processing                                       */
#include "general.h"                        /* for general desclarations                                         */
#include "np.h"                             /* for NUMPROC implementation                                        */
#include "assemble.h"                       /* for CreateCLASS implementation                                    */

// proprietary includes
#include "whitney.h"                        /* include Whitney base function definitions and procedures          */
#include "tdwaveforms.h"                    /* time domain shapes of excitation: gaussian etc.                   */
#include "aqhdSysParams.h"                  /* get system parameters from standard include file                  */
#include "aqhdReleaseInfo.h"                /* get release information                                           */
#include "aqhdPhysicoMath.h"                /* get mathematical & physical constants                             */
#include "aqhdBaseTypes.h"                  /* get mathematical & physical constants                             */
#include "auxiliary.h"


#ifdef HAVE_RLOG
/* Include the files for rlog. */
#include <rlog/rlog.h>
#include <rlog/rloglocation.h>
#include <rlog/Error.h>
#include <rlog/RLogChannel.h>
#include <rlog/StdioNode.h>
// Include this if you want to log the time.
#include <rlog/RLogTime.h>
#endif // HAVE_RLOG

using namespace rlog;
USING_UG_NAMESPACE
USING_UGDIM_NAMESPACE

#ifdef ModelP
USING_PPIF_NAMESPACE
#endif

#ifndef PARTICLESWARM_H_
#define PARTICLESWARM_H_

#define PARTICLE_SWARM_SOURCE_CLASS_NAME "particle_swarm"


/* declare NUMPROC class and methods for particleswarm current density source */
struct np_particle_swarm_source
{
  NP_BASE base;                            // inherits base class

  INT fl;                                  // from level with respect to multigrid
  INT tl;                                  // to level with respect to multigrid
  
  ELEMENT *theElement;                     // pointer to element containing source location
  VECDATA_DESC* force;                     // pointer to excitation aka force vector
  
  std::vector<std::vector<double> > particlepos;      
  // function adapts the discrete system obtained from the FE approach
  // by setting the respective coefficients of the involved matrices
  DOUBLE (*ParticleSwarmSourceSetup)(struct np_particle_swarm_source *);

  // function calculates the value of the time domain excitation and
  // and adapts the right hand side of the linear system to be solved
  DOUBLE (*ParticleSwarmSourceUpdate)(struct np_particle_swarm_source *, DOUBLE);

};
typedef struct np_particle_swarm_source NP_PARTICLE_SWARM_SOURCE;


INT InitParticleSwarmSource(void);

#endif /*PARTICLESWARM_H_*/

---