pde.h

Go to the documentation of this file.

#ifndef _PDE_H_
#define _PDE_H_

#include <mc/mc_base.h>

/*
 * ***************************************************************************
 * Class PDE: Parameters and datatypes
 * ***************************************************************************
 */

struct sPDE {

   /*
    * ************************************************************************
    * User-supplied problem information -- ALL must be supplied by the user.
    * ************************************************************************
    */

    int dim;         
    int dimII;            

    int vec;              
    int sym[MAXV][MAXV];  

    double est[MAXV];             
    int bmap[MAXV][VMAX_BDTYPE];  

    void (*initAssemble)(struct sPDE *thee, int ip[], double rp[]);
    void (*initElement)(struct sPDE *thee, int elementType,
        int chart, double tvx[][3], void *data);
    void (*initFace)(struct sPDE *thee, int faceType,
        int chart, double tnvec[]);
    void (*initPoint)(struct sPDE *thee, int pointType,
        int chart, double txq[],
        double U[], double dU[][3]);

    void (*Fu)(struct sPDE *thee, int key, double F[]);
    double (*Ju)(struct sPDE *thee, int key);
    double (*Fu_v)(struct sPDE *thee, int key,
        double V[], double dV[][3]);
    double (*DFu_wv)(struct sPDE *thee, int key,
        double W[], double dW[][3],
        double V[], double dV[][3]);
    double (*p_wv)(struct sPDE *thee,
        int key, double W[], double V[]);

    void (*delta)(struct sPDE *thee, int type,
        int chart, double txq[], void *data, double F[]);
    void (*u_D)(struct sPDE *thee, int type,
        int chart, double txq[], double F[]);
    void (*u_T)(struct sPDE *thee, int type,
                int chart, double txq[], double F[], double dF[][3]);

    void (*bisectEdge)(int dim, int dimII,
        int edgeType, int chart[], double vx[][3]);
    void (*mapBoundary)(int dim, int dimII, int vertexType,
        int chart, double vx[3]);
    int (*markSimplex)(int dim, int dimII,
        int simplexType, int faceType[4], int vertexType[4],
        int chart[], double vx[][3], void *data);
    void (*oneChart)(int dim, int dimII,
        int objType, int chart[], double vx[][3], int dimV);

    int (*simplexBasisInit)(int key, int dim, int comp,
        int *ndof, int dof[]);
    void (*simplexBasisForm)(int key, int dim, int comp,
        int pdkey, double xq[], double basis[]);

    void *user;

};

typedef struct sPDE PDE;

/*
 * ***************************************************************************
 * Class PDE: Inlineable methods (pde.c)
 * ***************************************************************************
 */

#if !defined(VINLINE_PDE)

    void PDE_setDim(PDE *thee, int d);

    void PDE_setDimII(PDE *thee, int d);

    int PDE_dim(PDE *thee);

    int PDE_dimII(PDE *thee);

    int PDE_vec(PDE *thee);
#else /* if defined(VINLINE_PDE) */
#   define PDE_setDim(thee,d)         ( (thee)->dim = (d) )

#   define PDE_setDimII(thee,d)       ( (thee)->dimII = (d) )

#   define PDE_dim(thee)              ( (thee)->dim )

#   define PDE_dimII(thee)            ( (thee)->dimII )

#   define PDE_vec(thee)              ( (thee)->vec )
#endif /* if !defined(VINLINE_PDE) */

/*
 * ***************************************************************************
 * Class PDE: Non-inlineable methods (pde.c)
 * ***************************************************************************
 */

PDE* PDE_ctor_default(void);

void PDE_dtor_default(PDE **thee);

int PDE_sym(PDE *thee, int i, int j);

#endif /* _PDE_H_ */

---