Revision 1107 trunk/src/df1b2separable/gammdev.cpp
gammdev.cpp (revision 1107)  

2  2 
* $Id$ 
3  3 
* 
4  4 
* Author: David Fournier 
5 
* Copyright (c) 20082012 Regents of the University of California


5 
* Copyright (c) 20082012 Regents of the University of California 

6  6 
*/ 
7  7 
/** 
8  8 
* \file 
...  ...  
57  57 
//init_df3_two_variable aa(2.0); 
58  58 
*xx.get_u_x()=1.0; 
59  59 
*aa.get_u_y()=1.0; 
60 


60  
61  61 
df3_two_variable z=cumd_gamma(xx,aa); 
62  62  
63  63 
// now use the derivatives of z to get the 
...  ...  
66  66 
//double ca=value(a); 
67  67  
68  68 
double F_x=1.0/(*z.get_u_x()); 
69 


69  
70  70 
double F_y=F_x* *z.get_u_y(); 
71  71 
double F_xx=F_x* *z.get_u_xx()/square(*z.get_u_x()); 
72  72  
73  73 
double F_xy=(F_xx * *z.get_u_x() * *z.get_u_y() + F_x * *z.get_u_xy())/ 
74  74 
*z.get_u_x(); 
75  75  
76 
double F_yy=(F_xx * square(*z.get_u_y())


76 
double F_yy=(F_xx * square(*z.get_u_y()) 

77  77 
+2.0* F_xy* *z.get_u_y() 
78  78 
+F_x * *z.get_u_yy()); 
79  79  
...  ...  
84  84 
double F_xxy=(F_xxx * square(*z.get_u_x())* *z.get_u_y() 
85  85 
+ 2.0*F_xx* *z.get_u_x()* *z.get_u_xy() 
86  86 
+ F_xx* *z.get_u_xx() * *z.get_u_y() 
87 
+ F_xy * *z.get_u_xx()


87 
+ F_xy * *z.get_u_xx() 

88  88 
+ F_x * *z.get_u_xxy())/ 
89  89 
square(*z.get_u_x()); 
90  90  
91  91 
double F_xyy=(F_xxx* *z.get_u_x() *square(*z.get_u_y()) 
92  92 
+2.0* F_xx * *z.get_u_xy() * *z.get_u_y() 
93  93 
+2.0*F_xxy * *z.get_u_x() * *z.get_u_y() 
94 
+ 2.0*F_xy * *z.get_u_xy()


94 
+ 2.0*F_xy * *z.get_u_xy() 

95  95 
+ F_xx * *z.get_u_x() * *z.get_u_yy() 
96  96 
+ F_x * *z.get_u_xyy())/ 
97  97 
*z.get_u_x(); 
...  ...  
114  114 
} 
115  115 
if (!df1b2_gradlist::no_derivatives) 
116  116 
{ 
117 
f1b2gradlist>write_pass1(&_y,&_a,&zz,


117 
f1b2gradlist>write_pass1(&_y,&_a,&zz, 

118  118 
F_x,F_y, 
119  119 
F_xx,F_xy,F_yy, 
120  120 
F_xxx,F_xxy,F_xyy,F_yyy); 
Also available in: Unified diff