FPDoubleWellComputeCFL.c File Reference

#include <basic.h>
#include <mathfunctions.h>
#include <physicalmodels/fpdoublewell.h>
#include <mpivars.h>
#include <hypar.h>

Go to the source code of this file.

Functions
double	FPDoubleWellComputeCFL (void *s, void *m, double dt, double t)

Function Documentation

FPDoubleWellComputeCFL()

double FPDoubleWellComputeCFL	(	void *	s,
		void *	m,
		double	dt,
		double	t
	)

Definition at line 7 of file FPDoubleWellComputeCFL.c.

{
  HyPar         *solver = (HyPar*)        s;
  int           d, i, v;

  int     ndims  = solver->ndims;
  int     nvars  = solver->nvars;
  int     ghosts = solver->ghosts;
  int     *dim   = solver->dim_local;
  double  *dxinv = solver->dxinv;

  int     offset  = 0;
  double  max_cfl = 0;
  for (d = 0; d < ndims; d++) {
    for (i = 0; i < dim[d]; i++) {
      for (v = 0; v < nvars; v++) {
        double x; _GetCoordinate_(0,i,dim,ghosts,solver->x,x);
        double local_cfl =  absolute(drift(x)) * dt 
                          * dxinv[offset+ghosts+i];
        if (local_cfl > max_cfl) max_cfl = local_cfl;
      }
    }
    offset += dim[d] + 2*ghosts;
  }

  return(max_cfl);
}