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

df1b2f12.cpp (revision 1110)
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,df1b2function2 * pf)
23 23
 {
24 24
   ADUNCONST(df1b2variable*,px)
......
84 84
    read_pass2_3();
85 85
    break;
86 86
  default:
87
    cerr << "illegal value for df1b2variable::pass = " 
87
    cerr << "illegal value for df1b2variable::pass = "
88 88
         << df1b2variable::passnumber << endl;
89 89
    exit(1);
90 90
  }
......
100 100
  // and  forward for bptr2 and nbptr2
101 101
  // the current entry+2 in bptr is the size of the record i.e
102 102
  // points to the next record
103
  //char * bptr=f1b2gradlist->bptr; 
104
  //char * bptr2=f1b2gradlist2->bptr; 
103
  //char * bptr=f1b2gradlist->bptr;
104
  //char * bptr2=f1b2gradlist2->bptr;
105 105
  int nvar=df1b2variable::nvar;
106
  test_smartlist& list=f1b2gradlist->list; 
106
  test_smartlist& list=f1b2gradlist->list;
107 107
  //f1b2gradlist->nlist-=sizeof(int);
108 108
  int num_bytes=f1b2gradlist->nlist.bptr->numbytes;
109 109
  list-=num_bytes;
......
134 134
  double * ydot=(double*)bptr;
135 135

  
136 136
  list.restoreposition(); // save pointer to beginning of record;
137
 
137

  
138 138
  // ****************************************************************
139 139
  // turn this off if no third derivatives are calculated
140 140
  // if (!no_third_derivatives)
......
170 170
  (pf->df2)(xu,yu),(pf->d2f11)(xu,yu),(pf->d2f12)(xu,yu),(pf->d2f22)(xu,yu),
171 171
  (pf->d3f111)(xu,yu),(pf->d3f112)(xu,yu),(pf->d3f122)(xu,yu),
172 172
  (pf->d3f222)(xu,yu),1);
173
 print_derivatives(pz,"z"); 
174
 print_derivatives(px,"x"); 
175
 print_derivatives(py,"y"); 
173
 print_derivatives(pz,"z");
174
 print_derivatives(px,"x");
175
 print_derivatives(py,"y");
176 176
#endif
177 177
#if defined(__DERCHECK__)
178 178
  if (derchecker)
......
189 189
        case 1:
190 190
          if (!derchecker->dotflag)
191 191
            px->u_bar[derchecker->index-1]+=derchecker->delta;
192
          else 
192
          else
193 193
            px->u_dot_bar[derchecker->index-1]+=derchecker->delta;
194 194
          break;
195 195
        case 2:
......
205 205
            pz->u_dot_bar[derchecker->index-1]+=derchecker->delta;
206 206
          break;
207 207
        default:
208
          cerr << "Invalid index value for dercheck_index was " 
208
          cerr << "Invalid index value for dercheck_index was "
209 209
               << derchecker->index << endl;
210 210
          break;
211 211
        }
......
216 216
        case 1:
217 217
          if (!derchecker->dotflag)
218 218
            px->u_bar[derchecker->index-1]-=derchecker->delta;
219
          else 
219
          else
220 220
            px->u_dot_bar[derchecker->index-1]-=derchecker->delta;
221 221
          break;
222 222
        case 2:
......
232 232
            pz->u_dot_bar[derchecker->index-1]-=derchecker->delta;
233 233
          break;
234 234
        default:
235
          cerr << "Invalid index value for dercheck_index was " 
235
          cerr << "Invalid index value for dercheck_index was "
236 236
               << derchecker->index << endl;
237 237
          break;
238 238
        }
239 239
        break;
240 240
      }
241
    }  
241
    }
242 242
  }
243 243
#endif
244 244

  
... This diff was truncated because it exceeds the maximum size that can be displayed.

Also available in: Unified diff