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

df1b2fn4.cpp (revision 1110)
39 39
      {
40 40
        *zd++ += *xd++;
41 41
      }
42
      
42

  
43 43
      // WRITE WHATEVER ON TAPE
44 44
      //df1b2tape->set_tapeinfo_header(&x,&z,this,xd);
45 45
      // save stuff for first reverse pass
46 46
      // need &x, &z, this,
47
    
47

  
48 48
      // !!! nov 1 02
49 49
      // !!! nov 11 02
50 50
      if (!df1b2_gradlist::no_derivatives)
......
61 61
        (*ncopies)++;
62 62
        u=x.u;
63 63
        u_dot=x.u_dot;
64
        u_bar=x.u_bar;         
64
        u_bar=x.u_bar;
65 65
        u_dot_bar=x.u_dot_bar;
66 66
        u_tilde=x.u_tilde;
67 67
        u_dot_tilde=x.u_dot_tilde;
......
92 92
 * Description not yet available.
93 93
 * \param
94 94
 */
95
int df1b2_gradlist::write_pass1_eq(const df1b2variable * _px, 
95
int df1b2_gradlist::write_pass1_eq(const df1b2variable * _px,
96 96
  df1b2variable * pz)
97 97
{
98 98
  ADUNCONST(df1b2variable*,px)
......
119 119

  
120 120
  memcpy(list,(df1b2_header*)(px),sizeof(df1b2_header));
121 121
  memcpy(list,(df1b2_header*)(pz),sizeof(df1b2_header));
122
  
122

  
123 123
  nlist.bptr->numbytes=adptr_diff(list.bptr,tmpptr);
124 124
  nlist.bptr->pf=(ADrfptr)(&ad_read_pass1_eq);
125 125
  ++nlist;
......
144 144
    read_pass1_eq_3();
145 145
    break;
146 146
  default:
147
    cerr << "illegal value for df1b2variable::pass = " 
147
    cerr << "illegal value for df1b2variable::pass = "
148 148
         << df1b2variable::passnumber << endl;
149 149
    exit(1);
150 150
  }
......
153 153
/**
154 154
 * Description not yet available.
155 155
 * \param
156
 */ 
156
 */
157 157
void read_pass1_eq_1(void)
158 158
{
159 159
  // We are going backword for bptr and forward for bptr2
160 160
  // the current entry+2 in bptr is the size of the record i.e
161 161
  // points to the next record
162 162
  int nvar=df1b2variable::nvar;
163
  fixed_smartlist & nlist=f1b2gradlist->nlist; 
164
  test_smartlist& list=f1b2gradlist->list; 
163
  fixed_smartlist & nlist=f1b2gradlist->nlist;
164
  test_smartlist& list=f1b2gradlist->list;
165 165
   // nlist-=sizeof(int);
166 166
  // get record size
167 167
  int num_bytes=nlist.bptr->numbytes;
......
180 180
  list.bptr+=sizeof(df1b2_header);
181 181
  list.restoreposition(); // save pointer to beginning of record;
182 182
  int i;
183
  
183

  
184 184
  // Do first reverse paSS calculations
185 185
  // ****************************************************************
186 186
  // turn this off if no third derivatives are calculated
......
188 188
  // {
189 189
  // save for second reverse pass
190 190
  // save identifier 1
191
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2; 
192
     test_smartlist& list2=f1b2gradlist->list2; 
191
     fixed_smartlist2& nlist2=f1b2gradlist->nlist2;
192
     test_smartlist& list2=f1b2gradlist->list2;
193 193

  
194 194

  
195 195
     int total_bytes=2*nvar*sizeof(double);
......
214 214
  // ****************************************************************
215 215
#if defined(PRINT_DERS)
216 216
 print_derivatives(" assign ", 1 ,1 ,0, 0,1);
217
 print_derivatives(pz,"z"); 
218
 print_derivatives(px,"x"); 
217
 print_derivatives(pz,"z");
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff