Revision 1110 trunk/src/df1b2-separable/df1b2f11.cpp

df1b2f11.cpp (revision 1110)
32 32
 * Description not yet available.
33 33
 * \param
34 34
 */
35
int df1b2_gradlist::write_pass1(const df1b2variable * _px, 
35
int df1b2_gradlist::write_pass1(const df1b2variable * _px,
36 36
  df1b2variable * pz, df1b2function1 * pf)
37 37
{
38 38
  ADUNCONST(df1b2variable*,px)
......
78 78
/**
79 79
 * Description not yet available.
80 80
 * \param
81
 */ 
81
 */
82 82
void ad_read_pass1(void)
83 83
{
84 84
  switch(df1b2variable::passnumber)
......
93 93
    read_pass1_3();
94 94
    break;
95 95
  default:
96
    cerr << "illegal value for df1b2variable::pass = " 
96
    cerr << "illegal value for df1b2variable::pass = "
97 97
         << df1b2variable::passnumber << endl;
98 98
    exit(1);
99 99
  }
......
102 102
/**
103 103
 * Description not yet available.
104 104
 * \param
105
 */    
105
 */
106 106
void read_pass1_1(void)
107 107
{
108 108
  // We are going backword for bptr and forward for bptr2
109 109
  // the current entry+2 in bptr is the size of the record i.e
110 110
  // points to the next record
111 111
  int nvar=df1b2variable::nvar;
112
  fixed_smartlist & nlist=f1b2gradlist->nlist; 
113
  test_smartlist& list=f1b2gradlist->list; 
112
  fixed_smartlist & nlist=f1b2gradlist->nlist;
113
  test_smartlist& list=f1b2gradlist->list;
114 114
   // nlist-=sizeof(int);
115 115
  // get record size
116 116
  int num_bytes=nlist.bptr->numbytes;
......
140 140
  xdot=(double*)list.bptr;
141 141
  list.restoreposition(); // save pointer to beginning of record;
142 142
  int i;
143
  
143

  
144 144
  // Do first reverse paSS calculations
145 145
  // ****************************************************************
146 146
  // turn this off if no third derivatives are calculated
......
148 148
  // {
149 149
  // save for second reverse pass
150 150
  // save identifier 1
151
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2; 
152
     test_smartlist& list2=f1b2gradlist->list2; 
151
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2;
152
     test_smartlist& list2=f1b2gradlist->list2;
153 153

  
154 154

  
155 155
  int total_bytes=2*nvar*sizeof(double);
......
162 162

  
163 163
  list2.check_buffer_size(total_bytes);
164 164
  void * tmpptr2=list2.bptr;
165
  
165

  
166 166
#if defined(SAFE_ALL)
167 167
  memcpy(list2,ids,slen);
168 168
#endif
......
177 177
#if defined(PRINT_DERS)
178 178
 print_derivatives(pf->funname,(pf->df)(xu),(pf->df)(xu),(pf->d2f)(xu),
179 179
  (pf->d3f)(xu),1);
180
 print_derivatives(pz,"z"); 
181
 print_derivatives(px,"x"); 
180
 print_derivatives(pz,"z");
181
 print_derivatives(px,"x");
182 182
#endif
183
 
183

  
184 184
  double df=(pf->df)(xu);
185 185
  double d2f=(pf->d2f)(xu);
186 186
  //double d3f=(pf->d3f)(xu);
......
200 200
    //px->u_dot_bar[i]+=(pf->df)(xu)*pz->u_dot_bar[i];
201 201
    px->u_dot_bar[i]+=df*pz->u_dot_bar[i];
202 202
  }
203
  
203

  
204 204
  // !!!!!!!!!!!!!!!!!!!!!!
205 205
  for (i=0;i<nvar;i++)
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff