fene-p.h

📄 View in API Reference (Doxygen) · View on basilisk.fr

#include "log-conform.h" [api]

double L2 = 1.;

static void fenep (double trA, double * nu, double * eta) {
  *eta = 1;
  *nu = 1./(1. - trA/L2);
} 

event defaults (i = 0) {
  f_s = fenep;
  f_r = fenep;
}

event init (i = 0) {
#if AXI
  double dim = 3;
#else
  double dim = dimension;
#endif  
  scalar trac = trA;
  foreach()
    trac[] = dim*L2/(dim + L2);
}