Revision 1488
trunk/contrib/statslib/logistic.cpp (revision 1488)  

1 
#include "statsLib.h" 

1 
#include "statsLib.h"


2  2  
3  3 
/** 
4  4 
* @file logistic.cpp 
trunk/src/df1b2separable/df3fun.cpp (revision 1488)  

378  378 
+ ::exp(*x.get_u())* *x.get_udot2(); 
379  379  
380  380 
*z.get_udot3() = ::exp(*x.get_u()) * cube(*x.get_udot()) 
381 
+ 3.0 * ::exp(*x.get_u()) * *x.get_udot() * *x.get_udot2();


381 
+ 3.0 * ::exp(*x.get_u()) * *x.get_udot() * *x.get_udot2() 

382  382 
+ ::exp(*x.get_u()) * *x.get_udot3(); 
383  383 
return z; 
384  384 
} 
...  ...  
399  399 
+ xp * *x.get_udot2(); 
400  400  
401  401 
*z.get_udot3() = xp3 * cube(*x.get_udot()) 
402 
+ 3.0 * xp2 * *x.get_udot() * *x.get_udot2();


402 
+ 3.0 * xp2 * *x.get_udot() * *x.get_udot2() 

403  403 
+ xp * *x.get_udot3(); 
404  404 
return z; 
405  405 
} 
...  ...  
779  779 
if (value(x)>=0) 
780  780 
{ 
781  781 
df3_one_variable u1=p*x; 
782 
df3_one_variable u2=1.+p*x;


783 
df3_one_variable u=1./(1.+p*x);


782 
df3_one_variable u2=1.+u1;


783 
df3_one_variable u=1./u2;


784  784 
df3_one_variable y= ((((b5*u+b4)*u+b3)*u+b2)*u+b1)*u; 
785  785 
df3_one_variable tmp1=0.3989422804*exp(.5*x*x); 
786  786 
z=1.0+tmp1*y; 
...  ...  
811  811 
if (value(x)>=0) 
812  812 
{ 
813  813 
df3_one_variable u1=p*x; 
814 
df3_one_variable u2=1.+p*x;


815 
df3_one_variable u=1./(1.+p*x);


814 
df3_one_variable u2=1.+u1;


815 
df3_one_variable u=1./u2;


816  816 
df3_one_variable y= ((((b5*u+b4)*u+b3)*u+b2)*u+b1)*u; 
817  817 
df3_one_variable tmp1=0.3989422804*exp(.5*x*x); 
818  818 
z=1.0+tmp1*y; 
Also available in: Unified diff