fpdoublewell.h File Reference

Go to the source code of this file.

Data Structures
struct	FPDoubleWell

Macros
#define	_FP_DOUBLE_WELL_   "fp-double-well"
#define	_MODEL_NDIMS_   1
#define	_MODEL_NVARS_   1
#define	drift(x)   (4.0*(x)*(1.0-(x)*(x)))

Functions
int	FPDoubleWellInitialize (void *, void *)
int	FPDoubleWellCleanup (void *)

---

Data Structure Documentation

FPDoubleWell

struct FPDoubleWell

Definition at line 26 of file fpdoublewell.h.

Data Fields
double	q
double	pdf_integral

Macro Definition Documentation

_FP_DOUBLE_WELL_

#define _FP_DOUBLE_WELL_   "fp-double-well"

Definition at line 18 of file fpdoublewell.h.

_MODEL_NDIMS_

#define _MODEL_NDIMS_   1

Definition at line 23 of file fpdoublewell.h.

_MODEL_NVARS_

#define _MODEL_NVARS_   1

Definition at line 24 of file fpdoublewell.h.

drift

#define drift

(

x

)

(4.0*(x)*(1.0-(x)*(x)))

Definition at line 31 of file fpdoublewell.h.

Function Documentation

FPDoubleWellInitialize()

int FPDoubleWellInitialize	(	void *	,
		void *
	)

Definition at line 19 of file FPDoubleWellInitialize.c.

{
  HyPar         *solver  = (HyPar*)         s;
  MPIVariables  *mpi     = (MPIVariables*)  m; 
  FPDoubleWell  *physics = (FPDoubleWell*)  solver->physics;
  int           ferr     = 0;
  _DECLARE_IERR_;

  static int count = 0;

  if (solver->nvars != _MODEL_NVARS_) {
    fprintf(stderr,"Error in FPDoubleWellInitializeO(): nvars has to be %d.\n",_MODEL_NVARS_);
    return(1);
  }
  if (solver->ndims != _MODEL_NDIMS_) {
    fprintf(stderr,"Error in FPDoubleWellInitializeO(): ndims has to be %d.\n",_MODEL_NDIMS_);
    return(1);
  }

  /* reading physical model specific inputs - all processes */
  if (!mpi->rank) {
    FILE *in;
    if (!count) printf("Reading physical model inputs from file \"physics.inp\".\n");
    in = fopen("physics.inp","r");
    if (!in) {
      fprintf(stderr,"Error: File \"physics.inp\" not found.\n");
      return(1);
    } else {
      char word[_MAX_STRING_SIZE_];
      ferr = fscanf(in,"%s",word); if (ferr != 1) return(1);
      if (!strcmp(word, "begin")){
          while (strcmp(word, "end")){
              ferr = fscanf(in,"%s",word); if (ferr != 1) return(1);
          if (!strcmp(word, "q")) {
            /* read diffusion coefficient */
            ferr = fscanf(in,"%lf",&physics->q);
            if (ferr != 1) return(1);
          } else if (strcmp(word,"end")) {
            char useless[_MAX_STRING_SIZE_];
            ferr = fscanf(in,"%s",useless); if (ferr != 1) return(ferr);
            printf("Warning: keyword %s in file \"physics.inp\" with value %s not ",word,useless);
            printf("recognized or extraneous. Ignoring.\n");
          }
        }
        } else {
          fprintf(stderr,"Error: Illegal format in file \"physics.inp\".\n");
        return(1);
        }
    }
    fclose(in);
  }

#ifndef serial
  IERR MPIBroadcast_double(&physics->q,1,0,&mpi->world); CHECKERR(ierr);
#endif

  if (!strcmp(solver->SplitHyperbolicFlux,"yes")) {
    if (!mpi->rank) {
      fprintf(stderr,"Error in FPDoubleWellInitialize: This physical model does not have a splitting ");
      fprintf(stderr,"of the hyperbolic term defined.\n");
    }
    return(1);
  }

  /* initializing physical model-specific functions */
  solver->ComputeCFL         = FPDoubleWellComputeCFL;
  solver->ComputeDiffNumber  = FPDoubleWellComputeDiffNumber;
  solver->FFunction          = FPDoubleWellAdvection;
  solver->GFunction          = FPDoubleWellDiffusion;
  solver->Upwind             = FPDoubleWellUpwind;
  solver->PostStep           = FPDoubleWellPostStep;
  solver->PrintStep          = FPDoubleWellPrintStep;

  /* Calculate and print the PDF integral of the initial solution */
  IERR FPDoubleWellPostStep(solver->u,solver,mpi,0.0,0);CHECKERR(ierr);
  IERR FPDoubleWellPrintStep(solver,mpi,0.0);           CHECKERR(ierr);
  
  count++;
  return(0);
}

FPDoubleWellCleanup()

int FPDoubleWellCleanup

(

void *

)

Definition at line 4 of file FPDoubleWellCleanup.c.

{
  return(0);
}