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

df1b2f27.cpp (revision 1110)
26 26
 * Description not yet available.
27 27
 * \param
28 28
 */
29
int df1b2_gradlist::write_pass1x(const df1b2variable * _px, 
29
int df1b2_gradlist::write_pass1x(const df1b2variable * _px,
30 30
  df1b2variable * pz, df1b2function1 * pf)
31 31
{
32 32
  ADUNCONST(df1b2variable*,px)
......
87 87
    read_pass1_3x();
88 88
    break;
89 89
  default:
90
    cerr << "illegal value for df1b2variable::pass = " 
90
    cerr << "illegal value for df1b2variable::pass = "
91 91
         << df1b2variable::passnumber << endl;
92 92
    exit(1);
93 93
  }
......
103 103
  // the current entry+2 in bptr is the size of the record i.e
104 104
  // points to the next record
105 105
  int nvar=df1b2variable::nvar;
106
  fixed_smartlist & nlist=f1b2gradlist->nlist; 
107
  test_smartlist& list=f1b2gradlist->list; 
106
  fixed_smartlist & nlist=f1b2gradlist->nlist;
107
  test_smartlist& list=f1b2gradlist->list;
108 108
   // nlist-=sizeof(int);
109 109
  // get record size
110 110
  int num_bytes=nlist.bptr->numbytes;
......
134 134
  xdot=(double*)list.bptr;
135 135
  list.restoreposition(); // save pointer to beginning of record;
136 136
  int i;
137
  
137

  
138 138
  // Do first reverse paSS calculations
139 139
  // ****************************************************************
140 140
  // turn this off if no third derivatives are calculated
......
142 142
  // {
143 143
  // save for second reverse pass
144 144
  // save identifier 1
145
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2; 
146
     test_smartlist& list2=f1b2gradlist->list2; 
145
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2;
146
     test_smartlist& list2=f1b2gradlist->list2;
147 147

  
148 148

  
149 149
  int total_bytes=2*nvar*sizeof(double);
......
156 156

  
157 157
  list2.check_buffer_size(total_bytes);
158 158
  void * tmpptr2=list2.bptr;
159
  
159

  
160 160
#if defined(SAFE_ALL)
161 161
  memcpy(list2,ids,slen);
162 162
#endif
......
171 171
#if defined(PRINT_DERS)
172 172
 print_derivatives(pf->funname,(pf->df)(xu),(pf->df)(xu),(pf->d2f)(xu),
173 173
  (pf->d3f)(xu),1);
174
 print_derivatives(pz,"z"); 
175
 print_derivatives(px,"x"); 
174
 print_derivatives(pz,"z");
175
 print_derivatives(px,"x");
176 176
#endif
177
 
177

  
178 178
  double df=(pf->df)(xu);
179 179
  double d2f=(pf->d2f)(xu);
180 180
  //double d3f=(pf->d3f)(xu);
......
194 194
    //px->u_dot_bar[i]+=(pf->df)(xu)*pz->u_dot_bar[i];
195 195
    px->u_dot_bar[i]+=df*pz->u_dot_bar[i];
196 196
  }
197
  
197

  
198 198
  // !!!!!!!!!!!!!!!!!!!!!!
199 199
  for (i=0;i<nvar;i++)
200 200
  {
......
204 204
  {
205 205
    pz->u_dot_bar[i]=0;
206 206
  }
207
  
207

  
208 208
#if defined(PRINT_DERS)
209
 print_derivatives(pz,"z"); 
210
 print_derivatives(px,"x"); 
209
 print_derivatives(pz,"z");
210
 print_derivatives(px,"x");
211 211
#endif
212 212
}
213 213

  
......
223 223
  // list 1
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff