Revision 1109 trunk/src/df1b2-separable/df32fun1.cpp

df32fun1.cpp (revision 1109)
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
......
18 18
 * Description not yet available.
19 19
 * \param
20 20
 */
21
 int df1b2_gradlist::write_pass1(const df1b2variable * _px, 
21
 int df1b2_gradlist::write_pass1(const df1b2variable * _px,
22 22
   const df1b2variable * _py,df1b2variable * pz,double df_x,
23 23
   double df_y,
24 24
   double df_xx,
......
104 104
    read_pass2_3_dvdv();
105 105
    break;
106 106
  default:
107
    cerr << "illegal value for df1b2variable::pass = " 
107
    cerr << "illegal value for df1b2variable::pass = "
108 108
         << df1b2variable::passnumber << endl;
109 109
    exit(1);
110 110
  }
......
120 120
  // and  forward for bptr2 and nbptr2
121 121
  // the current entry+2 in bptr is the size of the record i.e
122 122
  // points to the next record
123
  //char * bptr=f1b2gradlist->bptr; 
124
  //char * bptr2=f1b2gradlist2->bptr; 
123
  //char * bptr=f1b2gradlist->bptr;
124
  //char * bptr2=f1b2gradlist2->bptr;
125 125
  int nvar=df1b2variable::nvar;
126
  test_smartlist& list=f1b2gradlist->list; 
126
  test_smartlist& list=f1b2gradlist->list;
127 127
  //f1b2gradlist->nlist-=sizeof(int);
128 128
  int num_bytes=f1b2gradlist->nlist.bptr->numbytes;
129 129
  list-=num_bytes;
......
190 190
  double * ydot=(double*)bptr;
191 191

  
192 192
  list.restoreposition(); // save pointer to beginning of record;
193
 
193

  
194 194
  // ****************************************************************
195 195
  // turn this off if no third derivatives are calculated
196 196
  // if (!no_third_derivatives)
......
226 226
  (df2),(d2f11),(d2f12),(d2f22),
227 227
  (d3f111),(d3f112),(d3f122),
228 228
  (d3f222),1);
229
 print_derivatives(pz,"z"); 
230
 print_derivatives(px,"x"); 
231
 print_derivatives(py,"y"); 
229
 print_derivatives(pz,"z");
230
 print_derivatives(px,"x");
231
 print_derivatives(py,"y");
232 232
#endif
233 233
#if defined(__DERCHECK__)
234 234
  if (derchecker)
......
245 245
        case 1:
246 246
          if (!derchecker->dotflag)
247 247
            px->u_bar[derchecker->index-1]+=derchecker->delta;
248
          else 
248
          else
249 249
            px->u_dot_bar[derchecker->index-1]+=derchecker->delta;
250 250
          break;
251 251
        case 2:
......
261 261
            pz->u_dot_bar[derchecker->index-1]+=derchecker->delta;
262 262
          break;
263 263
        default:
264
          cerr << "Invalid index value for dercheck_index was " 
264
          cerr << "Invalid index value for dercheck_index was "
265 265
               << derchecker->index << endl;
266 266
          break;
267 267
        }
......
272 272
        case 1:
273 273
          if (!derchecker->dotflag)
274 274
            px->u_bar[derchecker->index-1]-=derchecker->delta;
275
          else 
275
          else
276 276
            px->u_dot_bar[derchecker->index-1]-=derchecker->delta;
277 277
          break;
278 278
        case 2:
......
288 288
            pz->u_dot_bar[derchecker->index-1]-=derchecker->delta;
289 289
          break;
290 290
        default:
291
          cerr << "Invalid index value for dercheck_index was " 
291
          cerr << "Invalid index value for dercheck_index was "
292 292
               << derchecker->index << endl;
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff