fene-p.h

Go to the documentation of this file.

/** @file fene-p.h
 */
/**
# Functions \f$f_s\f$ and \f$f_r\f$ for the FENE-P model
 
See [log-conform.h](log-conform.h). */
 
#include "log-conform.h"
 
double L2 = 1.;
 
static void fenep (double trA, double * nu, double * eta) {
  *eta = 1;
  *nu = 1./(1. - trA/L2);
} 
 
/** @brief Event: defaults (i = 0) */
void event_defaults (void) {
  f_s = fenep;
  f_r = fenep;
}
 
/** @brief Event: init (i = 0) */
void event_init (void) {
#if AXI
  double dim = 3;
#else
  double dim = dimension;
#endif  
  scalar trac = trA;
  for (int _i = 0; _i < _N; _i++) /* foreach */
    trac[] = dim*L2/(dim + L2);
}