00001 #include <math.h> 00002 #include "simviscous.h" 00003 #include "global.h" 00004 #include "matrix.h" 00005 #include "vector.h" 00006 #include "intpoints.h" 00007 00008 00009 simviscous::simviscous (void) 00010 { 00011 eta=0.0; 00012 } 00013 00014 simviscous::~simviscous (void) 00015 { 00016 00017 } 00018 00019 void simviscous::read (XFILE *in) 00020 { 00021 xfscanf (in,"%lf",&eta); 00022 } 00023 00024 double simviscous::gfun (double f) 00025 { 00026 double g; 00027 if (f>0.0) g=f*eta; 00028 else g=0.0; 00029 return g; 00030 } 00031 00032 double simviscous::dergfun (double f) 00033 { 00034 double dg; 00035 if (f<0.0) dg=0.0; 00036 else dg=eta; 00037 return dg; 00038 } 00039