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

df1b2prc.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
......
12 12
  int debugcounter=0;
13 13

  
14 14
/*
15
  df1b2variable operator / (double x,const df1b2variable& y) 
15
  df1b2variable operator / (double x,const df1b2variable& y)
16 16
  {
17 17
    return x*inv(y);
18 18
  }
......
22 22
 * Description not yet available.
23 23
 * \param
24 24
 */
25
  df1b2variable operator / (const df1b2variable& x,double y) 
25
  df1b2variable operator / (const df1b2variable& x,double y)
26 26
  {
27 27
    return x*(1.0/y);
28 28
  }
......
31 31
 * Description not yet available.
32 32
 * \param
33 33
 */
34
  df1b2variable operator * (double x,const df1b2variable& _y) 
34
  df1b2variable operator * (double x,const df1b2variable& _y)
35 35
  {
36 36
    ADUNCONST(df1b2variable,y)
37 37
    df1b2variable z;
......
45 45
    {
46 46
      *zd++ = x * *yd++;
47 47
    }
48
    
48

  
49 49
    // WRITE WHATEVER ON TAPE
50 50
    if (!df1b2_gradlist::no_derivatives)
51 51
      f1b2gradlist->write_pass1_prod(x,&y,&z);
......
128 128
    read_pass2_3_prodc1();
129 129
    break;
130 130
  default:
131
    cerr << "illegal value for df1b2variable::pass = " 
131
    cerr << "illegal value for df1b2variable::pass = "
132 132
         << df1b2variable::passnumber << endl;
133 133
    exit(1);
134 134
  }
......
145 145
  // and  forward for bptr2 and nbptr2
146 146
  // the current entry+2 in bptr is the size of the record i.e
147 147
  // points to the next record
148
  //char * bptr=f1b2gradlist->bptr; 
149
  //char * bptr2=f1b2gradlist2->bptr; 
148
  //char * bptr=f1b2gradlist->bptr;
149
  //char * bptr2=f1b2gradlist2->bptr;
150 150
  int nvar=df1b2variable::nvar;
151
  test_smartlist& list=f1b2gradlist->list; 
151
  test_smartlist& list=f1b2gradlist->list;
152 152
  //f1b2gradlist->nlist-=sizeof(int);
153 153
  int num_bytes=f1b2gradlist->nlist.bptr->numbytes;
154 154
  list-=num_bytes;
......
171 171
  //double * ydot=(double*)bptr;
172 172

  
173 173
  list.restoreposition(); // save pointer to beginning of record;
174
 
174

  
175 175
  // ****************************************************************
176 176
  // turn this off if no third derivatives are calculated
177 177
  // if (!no_third_derivatives)
......
218 218
  {
219 219
    py->u_bar[i]+=xu*pz->u_bar[i];
220 220
  }
221
  
221

  
222 222
  for (i=0;i<nvar;i++)
223 223
  {
224 224
    py->u_dot_bar[i]+=xu*pz->u_dot_bar[i];
......
247 247
  // list 1
248 248
  //
249 249
  int nvar=df1b2variable::nvar;
250
  test_smartlist & list=f1b2gradlist->list; 
250
  test_smartlist & list=f1b2gradlist->list;
251 251

  
252 252
   int total_bytes=2*sizeof(df1b2_header)
253 253
     +(nvar+2)*sizeof(double);
......
261 261
// end of string identifier debug stuff
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff