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

df1b2prb.cpp (revision 1110)
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
......
14 14
 * Description not yet available.
15 15
 * \param
16 16
 */
17
  df1b2variable operator * (const df1b2variable& _x, double y) 
17
  df1b2variable operator * (const df1b2variable& _x, double y)
18 18
  {
19 19
    ADUNCONST(df1b2variable,x)
20 20
    df1b2variable z;
......
28 28
    {
29 29
      *zd++ = y * *xd++;
30 30
    }
31
    
31

  
32 32
    // WRITE WHATEVER ON TAPE
33 33
    if (!df1b2_gradlist::no_derivatives)
34 34
      f1b2gradlist->write_pass1_prod(&x,y,&z);
35 35
    return z;
36 36
  }
37 37

  
38

  
39 38
void ad_read_pass2_prodc2(void);
40 39

  
41 40
/**
......
111 110
    read_pass2_3_prodc2();
112 111
    break;
113 112
  default:
114
    cerr << "illegal value for df1b2variable::pass = " 
113
    cerr << "illegal value for df1b2variable::pass = "
115 114
         << df1b2variable::passnumber << endl;
116 115
    exit(1);
117 116
  }
......
127 126
  // and  forward for bptr2 and nbptr2
128 127
  // the current entry+2 in bptr is the size of the record i.e
129 128
  // points to the next record
130
  //char * bptr=f1b2gradlist->bptr; 
131
  //char * bptr2=f1b2gradlist2->bptr; 
129
  //char * bptr=f1b2gradlist->bptr;
130
  //char * bptr2=f1b2gradlist2->bptr;
132 131
  int nvar=df1b2variable::nvar;
133
  test_smartlist& list=f1b2gradlist->list; 
132
  test_smartlist& list=f1b2gradlist->list;
134 133
  //f1b2gradlist->nlist-=sizeof(int);
135 134
  int num_bytes=f1b2gradlist->nlist.bptr->numbytes;
136 135
  list-=num_bytes;
......
153 152
  //double * xdot=(double*)bptr;
154 153

  
155 154
  list.restoreposition(); // save pointer to beginning of record;
156
 
155

  
157 156
  // ****************************************************************
158 157
  // turn this off if no third derivatives are calculated
159 158
  // if (!no_third_derivatives)
......
188 187
  {
189 188
    px->u_bar[i]+=yu*pz->u_bar[i];
190 189
  }
191
  
190

  
192 191
  for (i=0;i<nvar;i++)
193 192
  {
194 193
    px->u_dot_bar[i]+=yu*pz->u_dot_bar[i];
......
217 216
  // list 1
218 217
  //
219 218
  int nvar=df1b2variable::nvar;
220
  test_smartlist & list=f1b2gradlist->list; 
219
  test_smartlist & list=f1b2gradlist->list;
221 220

  
222 221
   int total_bytes=2*sizeof(df1b2_header)
223 222
     +(nvar+2)*sizeof(double);
......
231 230
// end of string identifier debug stuff
232 231

  
233 232
  list.saveposition(); // save pointer to beginning of record;
234
  fixed_smartlist & nlist=f1b2gradlist->nlist; 
233
  fixed_smartlist & nlist=f1b2gradlist->nlist;
235 234
   // nlist-=sizeof(int);
236 235
  // get record size
237 236
  int num_bytes=nlist.bptr->numbytes;
238 237
  //
239 238
  // list 2
240 239
  //
241
  test_smartlist & list2=f1b2gradlist->list2; 
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff