Revision 1107 trunk/src/df1b2-separable/gammdev.cpp

gammdev.cpp (revision 1107)
2 2
 * $Id$
3 3
 *
4 4
 * Author: David Fournier
5
 * Copyright (c) 2008-2012 Regents of the University of California 
5
 * Copyright (c) 2008-2012 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