root/trunk/rheoplast/rheoplast.h

/***************************************
  $Header$

  This #includes all of the the typedefs for the various field variable
  structures, and their function prototypes, and has the main AppCtx and
  +latex+{\tt DPRINTF()}.
  +html+ <tt>DPRINTF()</tt>.
  ***************************************/


#ifndef RHEOPLAST_H
#define RHEOPLAST_H    /*+ To stop multiple inclusions. +*/

#include <petsc.h>  /* For PetscScalar and PetscPrintf() */
#include <stdlib.h> /* For (s)random (seeded by main()) */
#include <illuminator.h>
#include <petscts.h>       /* For the (semi-)implicit parts */
#include <petscblaslapack.h> /* For the BLAS parts */
#include <stdlib.h>          /* For (s)random */
#include <time.h>            /* For CLK_TCK in tsmonitor */
#include <sys/times.h>       /* For times() in tsmonitor */

/* This includes timestep.h and headers for each of the modules. */
#include "timestep.h"
#include "cahnhill.h"
#include "chternary.h"
#include "vectorphase.h"
#include "heatcond.h"
#include "vortflow.h"
#include "pressflow.h"
#include "shearstrain.h"
#include "electra.h"

/* Redefine DPRINTF to use data->logfile instead of timestep.h's logfile. */
#undef DPRINTF
#ifdef DEBUG
#define DPRINTF(fmt, args...) { \
  ierr = PetscPrintf (PETSC_COMM_WORLD, "%s: " fmt, __FUNCT__, args); \
  CHKERRQ (ierr); \
  if (data->logfile) { \
    ierr = PetscFPrintf (PETSC_COMM_WORLD, data->logfile, "%s: " fmt, __FUNCT__, args); \
    CHKERRQ (ierr); }}
#else
#define DPRINTF(fmt, args...) { \
  if (data->logfile) { \
    ierr = PetscFPrintf (PETSC_COMM_WORLD, data->logfile, "%s: " fmt, __FUNCT__, args); \
    CHKERRQ (ierr); }}
#endif

typedef struct {
  DA theda;
  Vec global, local, localfunc; /* localfunc and parameters are for TS only */
  PetscScalar *parameters, xwid, ywid, zwid;                  /* Domain size */
  PetscScalar deltax_m2, deltay_m2, deltaz_m2, deltat; /* 1/dx^2, 1/dy^2, dt */
  int expsteps, impsteps, current_timestep, monsteps, vars, tempvars, intvars,
    bcflags, *symmtypes;
  char **label, *save_basename;
  FILE *logfile;
  field_plot_type *plot_types;
  PetscTruth contours, jacobian, threedee, load_data;
  EqStyle *timestyle;
  PetscViewer theviewer;
  Mat J;
  Vec intvec;
  CommStyle style;
  /* Flags and structures for each of the modules */
  PetscTruth cahnhill, chternary, vectorphase, heatcond, vortflow, pressflow,
    shearstrain, electra;
  chparm thecahnhill;
  chtparm thechternary;
  vectorphasers thephasers;
  heatparm theheater;
  vortparm thevortex;
  pressparm thepressure;
  strainparm thestrain;
  echemparm thepotential;
} AppCtx;

#define APPCTX_DEFINED

/* This is kind of convoluted, we need vectorphasers etc. typedefs for the
   AppCtx typedef, and in turn need the AppCtx typedef for the function
   prototypes!  So those files are split into three parts; effectively we use
   the first part (typedefs) in the first include above, and the second part
   (prototypes) in the second include here. */
#include "cahnhill.h"
#include "chternary.h"
#include "vectorphase.h"
#include "heatcond.h"
#include "vortflow.h"
#include "pressflow.h"
#include "shearstrain.h"
#include "electra.h"

#endif /* RHEOPLAST_H */