Sacado Package Browser (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
tradoptest_24.cpp
Go to the documentation of this file.
1 
2 
3 /* Try to test all combinations of types and operations */
4 
5 
6 
7 #define ADT_RAD Sacado::Rad::
8 
9 
10 
11 #include "Sacado_trad.hpp"
12 
13 #include <cstdio>
14 
15 using std::printf;
16 
17 
18 
19 typedef ADT_RAD IndepADvar<double> AI;
20 
21 typedef ADT_RAD ADvar<double> A;
22 
23 typedef ADT_RAD ConstADvar<double> C;
24 
25 typedef ADT_RAD ADvari<double> Ai;
26 
27 typedef const ADT_RAD IndepADvar<double> cAI;
28 
29 typedef const ADT_RAD ADvar<double> cA;
30 
31 typedef const ADT_RAD ConstADvar<double> cC;
32 
33 typedef const ADT_RAD ADvari<double> cAi;
34 
35 static int rc;
36 
37 
38 
39 /* This is to be run through an awk program that changes lines */
40 
41 /* with "BINTEST" or "UNOPTEST" at the beginning of the line into */
42 
43 /* a the desired C++ (which we can then inspect). */
44 
45 
46 
47  void
48 
49 botch(const char *what, double wanted, double got)
50 
51 {
52 
53  printf("%s: expected %g, got %g, diff = %.2g\n", what, wanted, got, wanted-got);
54 
55  rc = 1;
56 
57  }
58 
59 
60 
61  const double tol = 5e-16;
62 
63 
64 
65  int
66 
67 differ(double a, double b)
68 
69 {
70 
71  double d = a - b;
72 
73  if (d < 0.)
74 
75  d = -d;
76 
77  if (a < 0.)
78 
79  a = -a;
80 
81  if (b < 0.)
82 
83  b = -b;
84 
85  if (a < b)
86 
87  a = b;
88 
89  if (a > 0.)
90 
91  d /= a;
92 
93  return d > tol;
94 
95  }
96 
97 
98 
99 #ifndef RAD_EQ_ALIAS
100 
101 #define Plus_dx 1.
102 
103 #else
104 
105 #ifdef RAD_AUTO_AD_Const
106 
107 #define Plus_dx 1.
108 
109 #else
110 
111 #define Plus_dx 0.
112 
113 #endif
114 
115 #endif
116 
117 
118 
119  int
120 
121 main(void)
122 
123 {
124 
125  AI xAI, yAI;
126 
127  A fA, xA, yA;
128 
129  C xC, yC;
130 
131  double dx, dy, f, xd, yd;
132 
133  long xL, yL;
134 
135  int xi, yi;
136 
137 
138 
139  rc = 0;
140 
141 
142  /**** Test of - ****/
143 
144  xd = 7.; yd = 4; f = 3.; dx = 1.; dy = -1.;
145  xAI = xd;
146  yAI = yd;
147  fA = xAI - yAI;
148  A::Gradcomp();
149  if (differ(fA.val(), f)) botch("fA = xAI - yAI", f, fA.val());
150  else if (differ(xAI.adj(), dx)) botch("d xAI - yAI/dx", dx, xAI.adj());
151  else if (differ(yAI.adj(), dy)) botch("d xAI - yAI/dy", dy, yAI.adj());
152  {
153  A::aval_reset();
154  cAI xcAI(xd);
155  yAI = yd;
156  fA = xcAI - yAI;
157  A::Gradcomp();
158  if (differ(fA.val(), f)) botch("fA = xcAI - yAI", f, fA.val());
159  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yAI/dx", dx, xcAI.adj());
160  else if (differ(yAI.adj(), dy)) botch("d xcAI - yAI/dy", dy, yAI.adj());
161  }
162  {
163  A::aval_reset();
164  xAI = xd;
165  cAI ycAI(yd);
166  fA = xAI - ycAI;
167  A::Gradcomp();
168  if (differ(fA.val(), f)) botch("fA = xAI - ycAI", f, fA.val());
169  else if (differ(xAI.adj(), dx)) botch("d xAI - ycAI/dx", dx, xAI.adj());
170  else if (differ(ycAI.adj(), dy)) botch("d xAI - ycAI/dy", dy, ycAI.adj());
171  }
172  {
173  A::aval_reset();
174  cAI xcAI(xd);
175  cAI ycAI(yd);
176  fA = xcAI - ycAI;
177  A::Gradcomp();
178  if (differ(fA.val(), f)) botch("fA = xcAI - ycAI", f, fA.val());
179  else if (differ(xcAI.adj(), dx)) botch("d xcAI - ycAI/dx", dx, xcAI.adj());
180  else if (differ(ycAI.adj(), dy)) botch("d xcAI - ycAI/dy", dy, ycAI.adj());
181  }
182  xAI = xd;
183  yA = yd;
184  fA = xAI - yA;
185  A::Gradcomp();
186  if (differ(fA.val(), f)) botch("fA = xAI - yA", f, fA.val());
187  else if (differ(xAI.adj(), dx)) botch("d xAI - yA/dx", dx, xAI.adj());
188  else if (differ(yA.adj(), dy)) botch("d xAI - yA/dy", dy, yA.adj());
189  {
190  A::aval_reset();
191  cAI xcAI(xd);
192  yA = yd;
193  fA = xcAI - yA;
194  A::Gradcomp();
195  if (differ(fA.val(), f)) botch("fA = xcAI - yA", f, fA.val());
196  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yA/dx", dx, xcAI.adj());
197  else if (differ(yA.adj(), dy)) botch("d xcAI - yA/dy", dy, yA.adj());
198  }
199  {
200  A::aval_reset();
201  xAI = xd;
202  cA ycA(yd);
203  fA = xAI - ycA;
204  A::Gradcomp();
205  if (differ(fA.val(), f)) botch("fA = xAI - ycA", f, fA.val());
206  else if (differ(xAI.adj(), dx)) botch("d xAI - ycA/dx", dx, xAI.adj());
207  else if (differ(ycA.adj(), dy)) botch("d xAI - ycA/dy", dy, ycA.adj());
208  }
209  {
210  A::aval_reset();
211  cAI xcAI(xd);
212  cA ycA(yd);
213  fA = xcAI - ycA;
214  A::Gradcomp();
215  if (differ(fA.val(), f)) botch("fA = xcAI - ycA", f, fA.val());
216  else if (differ(xcAI.adj(), dx)) botch("d xcAI - ycA/dx", dx, xcAI.adj());
217  else if (differ(ycA.adj(), dy)) botch("d xcAI - ycA/dy", dy, ycA.adj());
218  }
219  xAI = xd;
220  yC = yd;
221  fA = xAI - yC;
222  A::Gradcomp();
223  if (differ(fA.val(), f)) botch("fA = xAI - yC", f, fA.val());
224  else if (differ(xAI.adj(), dx)) botch("d xAI - yC/dx", dx, xAI.adj());
225  else if (differ(yC.adj(), dy)) botch("d xAI - yC/dy", dy, yC.adj());
226  {
227  A::aval_reset();
228  cAI xcAI(xd);
229  yC = yd;
230  fA = xcAI - yC;
231  A::Gradcomp();
232  if (differ(fA.val(), f)) botch("fA = xcAI - yC", f, fA.val());
233  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yC/dx", dx, xcAI.adj());
234  else if (differ(yC.adj(), dy)) botch("d xcAI - yC/dy", dy, yC.adj());
235  }
236  {
237  A::aval_reset();
238  xAI = xd;
239  cC ycC(yd);
240  fA = xAI - ycC;
241  A::Gradcomp();
242  if (differ(fA.val(), f)) botch("fA = xAI - ycC", f, fA.val());
243  else if (differ(xAI.adj(), dx)) botch("d xAI - ycC/dx", dx, xAI.adj());
244  else if (differ(ycC.adj(), dy)) botch("d xAI - ycC/dy", dy, ycC.adj());
245  }
246  {
247  A::aval_reset();
248  cAI xcAI(xd);
249  cC ycC(yd);
250  fA = xcAI - ycC;
251  A::Gradcomp();
252  if (differ(fA.val(), f)) botch("fA = xcAI - ycC", f, fA.val());
253  else if (differ(xcAI.adj(), dx)) botch("d xcAI - ycC/dx", dx, xcAI.adj());
254  else if (differ(ycC.adj(), dy)) botch("d xcAI - ycC/dy", dy, ycC.adj());
255  }
256  {
257  xAI = xd;
258  Ai yAi(yd);
259  fA = xAI - yAi;
260  A::Gradcomp();
261  if (differ(fA.val(), f)) botch("fA = xAI - yAi", f, fA.val());
262  else if (differ(xAI.adj(), dx)) botch("d xAI - yAi/dx", dx, xAI.adj());
263  else if (differ(yAi.aval, dy)) botch("d xAI - yAi/dy", dy, yAi.aval);
264  }
265  {
266  A::aval_reset();
267  cAI xcAI(xd);
268  Ai yAi(yd);
269  fA = xcAI - yAi;
270  A::Gradcomp();
271  if (differ(fA.val(), f)) botch("fA = xcAI - yAi", f, fA.val());
272  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yAi/dx", dx, xcAI.adj());
273  else if (differ(yAi.aval, dy)) botch("d xcAI - yAi/dy", dy, yAi.aval);
274  }
275  {
276  A::aval_reset();
277  xAI = xd;
278  cAi ycAi(yd);
279  fA = xAI - ycAi;
280  A::Gradcomp();
281  if (differ(fA.val(), f)) botch("fA = xAI - ycAi", f, fA.val());
282  else if (differ(xAI.adj(), dx)) botch("d xAI - ycAi/dx", dx, xAI.adj());
283  else if (differ(ycAi.aval, dy)) botch("d xAI - ycAi/dy", dy, ycAi.aval);
284  }
285  {
286  A::aval_reset();
287  cAI xcAI(xd);
288  cAi ycAi(yd);
289  fA = xcAI - ycAi;
290  A::Gradcomp();
291  if (differ(fA.val(), f)) botch("fA = xcAI - ycAi", f, fA.val());
292  else if (differ(xcAI.adj(), dx)) botch("d xcAI - ycAi/dx", dx, xcAI.adj());
293  else if (differ(ycAi.aval, dy)) botch("d xcAI - ycAi/dy", dy, ycAi.aval);
294  }
295  xAI = xd;
296  fA = xAI - yd;
297  A::Gradcomp();
298  if (differ(fA.val(), f)) botch("fA = xAI - yd", f, fA.val());
299  else if (differ(xAI.adj(), dx)) botch("d xAI - yd/dx", dx, xAI.adj());
300  {
301  A::aval_reset();
302  cAI xcAI(xd);
303  fA = xcAI - yd;
304  A::Gradcomp();
305  if (differ(fA.val(), f)) botch("fA = xcAI - yd", f, fA.val());
306  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yd/dx", dx, xcAI.adj());
307  }
308  xAI = xd;
309  yL = (long)yd;
310  fA = xAI - yL;
311  A::Gradcomp();
312  if (differ(fA.val(), f)) botch("fA = xAI - yL", f, fA.val());
313  else if (differ(xAI.adj(), dx)) botch("d xAI - yL/dx", dx, xAI.adj());
314  {
315  A::aval_reset();
316  cAI xcAI(xd);
317  yL = (long)yd;
318  fA = xcAI - yL;
319  A::Gradcomp();
320  if (differ(fA.val(), f)) botch("fA = xcAI - yL", f, fA.val());
321  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yL/dx", dx, xcAI.adj());
322  }
323  xAI = xd;
324  yi = (int)yd;
325  fA = xAI - yi;
326  A::Gradcomp();
327  if (differ(fA.val(), f)) botch("fA = xAI - yi", f, fA.val());
328  else if (differ(xAI.adj(), dx)) botch("d xAI - yi/dx", dx, xAI.adj());
329  {
330  A::aval_reset();
331  cAI xcAI(xd);
332  yi = (int)yd;
333  fA = xcAI - yi;
334  A::Gradcomp();
335  if (differ(fA.val(), f)) botch("fA = xcAI - yi", f, fA.val());
336  else if (differ(xcAI.adj(), dx)) botch("d xcAI - yi/dx", dx, xcAI.adj());
337  }
338  xA = xd;
339  yAI = yd;
340  fA = xA - yAI;
341  A::Gradcomp();
342  if (differ(fA.val(), f)) botch("fA = xA - yAI", f, fA.val());
343  else if (differ(xA.adj(), dx)) botch("d xA - yAI/dx", dx, xA.adj());
344  else if (differ(yAI.adj(), dy)) botch("d xA - yAI/dy", dy, yAI.adj());
345  {
346  A::aval_reset();
347  cA xcA(xd);
348  yAI = yd;
349  fA = xcA - yAI;
350  A::Gradcomp();
351  if (differ(fA.val(), f)) botch("fA = xcA - yAI", f, fA.val());
352  else if (differ(xcA.adj(), dx)) botch("d xcA - yAI/dx", dx, xcA.adj());
353  else if (differ(yAI.adj(), dy)) botch("d xcA - yAI/dy", dy, yAI.adj());
354  }
355  {
356  A::aval_reset();
357  xA = xd;
358  cAI ycAI(yd);
359  fA = xA - ycAI;
360  A::Gradcomp();
361  if (differ(fA.val(), f)) botch("fA = xA - ycAI", f, fA.val());
362  else if (differ(xA.adj(), dx)) botch("d xA - ycAI/dx", dx, xA.adj());
363  else if (differ(ycAI.adj(), dy)) botch("d xA - ycAI/dy", dy, ycAI.adj());
364  }
365  {
366  A::aval_reset();
367  cA xcA(xd);
368  cAI ycAI(yd);
369  fA = xcA - ycAI;
370  A::Gradcomp();
371  if (differ(fA.val(), f)) botch("fA = xcA - ycAI", f, fA.val());
372  else if (differ(xcA.adj(), dx)) botch("d xcA - ycAI/dx", dx, xcA.adj());
373  else if (differ(ycAI.adj(), dy)) botch("d xcA - ycAI/dy", dy, ycAI.adj());
374  }
375  xA = xd;
376  yA = yd;
377  fA = xA - yA;
378  A::Gradcomp();
379  if (differ(fA.val(), f)) botch("fA = xA - yA", f, fA.val());
380  else if (differ(xA.adj(), dx)) botch("d xA - yA/dx", dx, xA.adj());
381  else if (differ(yA.adj(), dy)) botch("d xA - yA/dy", dy, yA.adj());
382  {
383  A::aval_reset();
384  cA xcA(xd);
385  yA = yd;
386  fA = xcA - yA;
387  A::Gradcomp();
388  if (differ(fA.val(), f)) botch("fA = xcA - yA", f, fA.val());
389  else if (differ(xcA.adj(), dx)) botch("d xcA - yA/dx", dx, xcA.adj());
390  else if (differ(yA.adj(), dy)) botch("d xcA - yA/dy", dy, yA.adj());
391  }
392  {
393  A::aval_reset();
394  xA = xd;
395  cA ycA(yd);
396  fA = xA - ycA;
397  A::Gradcomp();
398  if (differ(fA.val(), f)) botch("fA = xA - ycA", f, fA.val());
399  else if (differ(xA.adj(), dx)) botch("d xA - ycA/dx", dx, xA.adj());
400  else if (differ(ycA.adj(), dy)) botch("d xA - ycA/dy", dy, ycA.adj());
401  }
402  {
403  A::aval_reset();
404  cA xcA(xd);
405  cA ycA(yd);
406  fA = xcA - ycA;
407  A::Gradcomp();
408  if (differ(fA.val(), f)) botch("fA = xcA - ycA", f, fA.val());
409  else if (differ(xcA.adj(), dx)) botch("d xcA - ycA/dx", dx, xcA.adj());
410  else if (differ(ycA.adj(), dy)) botch("d xcA - ycA/dy", dy, ycA.adj());
411  }
412  xA = xd;
413  yC = yd;
414  fA = xA - yC;
415  A::Gradcomp();
416  if (differ(fA.val(), f)) botch("fA = xA - yC", f, fA.val());
417  else if (differ(xA.adj(), dx)) botch("d xA - yC/dx", dx, xA.adj());
418  else if (differ(yC.adj(), dy)) botch("d xA - yC/dy", dy, yC.adj());
419  {
420  A::aval_reset();
421  cA xcA(xd);
422  yC = yd;
423  fA = xcA - yC;
424  A::Gradcomp();
425  if (differ(fA.val(), f)) botch("fA = xcA - yC", f, fA.val());
426  else if (differ(xcA.adj(), dx)) botch("d xcA - yC/dx", dx, xcA.adj());
427  else if (differ(yC.adj(), dy)) botch("d xcA - yC/dy", dy, yC.adj());
428  }
429  {
430  A::aval_reset();
431  xA = xd;
432  cC ycC(yd);
433  fA = xA - ycC;
434  A::Gradcomp();
435  if (differ(fA.val(), f)) botch("fA = xA - ycC", f, fA.val());
436  else if (differ(xA.adj(), dx)) botch("d xA - ycC/dx", dx, xA.adj());
437  else if (differ(ycC.adj(), dy)) botch("d xA - ycC/dy", dy, ycC.adj());
438  }
439  {
440  A::aval_reset();
441  cA xcA(xd);
442  cC ycC(yd);
443  fA = xcA - ycC;
444  A::Gradcomp();
445  if (differ(fA.val(), f)) botch("fA = xcA - ycC", f, fA.val());
446  else if (differ(xcA.adj(), dx)) botch("d xcA - ycC/dx", dx, xcA.adj());
447  else if (differ(ycC.adj(), dy)) botch("d xcA - ycC/dy", dy, ycC.adj());
448  }
449  {
450  xA = xd;
451  Ai yAi(yd);
452  fA = xA - yAi;
453  A::Gradcomp();
454  if (differ(fA.val(), f)) botch("fA = xA - yAi", f, fA.val());
455  else if (differ(xA.adj(), dx)) botch("d xA - yAi/dx", dx, xA.adj());
456  else if (differ(yAi.aval, dy)) botch("d xA - yAi/dy", dy, yAi.aval);
457  }
458  {
459  A::aval_reset();
460  cA xcA(xd);
461  Ai yAi(yd);
462  fA = xcA - yAi;
463  A::Gradcomp();
464  if (differ(fA.val(), f)) botch("fA = xcA - yAi", f, fA.val());
465  else if (differ(xcA.adj(), dx)) botch("d xcA - yAi/dx", dx, xcA.adj());
466  else if (differ(yAi.aval, dy)) botch("d xcA - yAi/dy", dy, yAi.aval);
467  }
468  {
469  A::aval_reset();
470  xA = xd;
471  cAi ycAi(yd);
472  fA = xA - ycAi;
473  A::Gradcomp();
474  if (differ(fA.val(), f)) botch("fA = xA - ycAi", f, fA.val());
475  else if (differ(xA.adj(), dx)) botch("d xA - ycAi/dx", dx, xA.adj());
476  else if (differ(ycAi.aval, dy)) botch("d xA - ycAi/dy", dy, ycAi.aval);
477  }
478  {
479  A::aval_reset();
480  cA xcA(xd);
481  cAi ycAi(yd);
482  fA = xcA - ycAi;
483  A::Gradcomp();
484  if (differ(fA.val(), f)) botch("fA = xcA - ycAi", f, fA.val());
485  else if (differ(xcA.adj(), dx)) botch("d xcA - ycAi/dx", dx, xcA.adj());
486  else if (differ(ycAi.aval, dy)) botch("d xcA - ycAi/dy", dy, ycAi.aval);
487  }
488  xA = xd;
489  fA = xA - yd;
490  A::Gradcomp();
491  if (differ(fA.val(), f)) botch("fA = xA - yd", f, fA.val());
492  else if (differ(xA.adj(), dx)) botch("d xA - yd/dx", dx, xA.adj());
493  {
494  A::aval_reset();
495  cA xcA(xd);
496  fA = xcA - yd;
497  A::Gradcomp();
498  if (differ(fA.val(), f)) botch("fA = xcA - yd", f, fA.val());
499  else if (differ(xcA.adj(), dx)) botch("d xcA - yd/dx", dx, xcA.adj());
500  }
501  xA = xd;
502  yL = (long)yd;
503  fA = xA - yL;
504  A::Gradcomp();
505  if (differ(fA.val(), f)) botch("fA = xA - yL", f, fA.val());
506  else if (differ(xA.adj(), dx)) botch("d xA - yL/dx", dx, xA.adj());
507  {
508  A::aval_reset();
509  cA xcA(xd);
510  yL = (long)yd;
511  fA = xcA - yL;
512  A::Gradcomp();
513  if (differ(fA.val(), f)) botch("fA = xcA - yL", f, fA.val());
514  else if (differ(xcA.adj(), dx)) botch("d xcA - yL/dx", dx, xcA.adj());
515  }
516  xA = xd;
517  yi = (int)yd;
518  fA = xA - yi;
519  A::Gradcomp();
520  if (differ(fA.val(), f)) botch("fA = xA - yi", f, fA.val());
521  else if (differ(xA.adj(), dx)) botch("d xA - yi/dx", dx, xA.adj());
522  {
523  A::aval_reset();
524  cA xcA(xd);
525  yi = (int)yd;
526  fA = xcA - yi;
527  A::Gradcomp();
528  if (differ(fA.val(), f)) botch("fA = xcA - yi", f, fA.val());
529  else if (differ(xcA.adj(), dx)) botch("d xcA - yi/dx", dx, xcA.adj());
530  }
531  xC = xd;
532  yAI = yd;
533  fA = xC - yAI;
534  A::Gradcomp();
535  if (differ(fA.val(), f)) botch("fA = xC - yAI", f, fA.val());
536  else if (differ(xC.adj(), dx)) botch("d xC - yAI/dx", dx, xC.adj());
537  else if (differ(yAI.adj(), dy)) botch("d xC - yAI/dy", dy, yAI.adj());
538  {
539  A::aval_reset();
540  cC xcC(xd);
541  yAI = yd;
542  fA = xcC - yAI;
543  A::Gradcomp();
544  if (differ(fA.val(), f)) botch("fA = xcC - yAI", f, fA.val());
545  else if (differ(xcC.adj(), dx)) botch("d xcC - yAI/dx", dx, xcC.adj());
546  else if (differ(yAI.adj(), dy)) botch("d xcC - yAI/dy", dy, yAI.adj());
547  }
548  {
549  A::aval_reset();
550  xC = xd;
551  cAI ycAI(yd);
552  fA = xC - ycAI;
553  A::Gradcomp();
554  if (differ(fA.val(), f)) botch("fA = xC - ycAI", f, fA.val());
555  else if (differ(xC.adj(), dx)) botch("d xC - ycAI/dx", dx, xC.adj());
556  else if (differ(ycAI.adj(), dy)) botch("d xC - ycAI/dy", dy, ycAI.adj());
557  }
558  {
559  A::aval_reset();
560  cC xcC(xd);
561  cAI ycAI(yd);
562  fA = xcC - ycAI;
563  A::Gradcomp();
564  if (differ(fA.val(), f)) botch("fA = xcC - ycAI", f, fA.val());
565  else if (differ(xcC.adj(), dx)) botch("d xcC - ycAI/dx", dx, xcC.adj());
566  else if (differ(ycAI.adj(), dy)) botch("d xcC - ycAI/dy", dy, ycAI.adj());
567  }
568  xC = xd;
569  yA = yd;
570  fA = xC - yA;
571  A::Gradcomp();
572  if (differ(fA.val(), f)) botch("fA = xC - yA", f, fA.val());
573  else if (differ(xC.adj(), dx)) botch("d xC - yA/dx", dx, xC.adj());
574  else if (differ(yA.adj(), dy)) botch("d xC - yA/dy", dy, yA.adj());
575  {
576  A::aval_reset();
577  cC xcC(xd);
578  yA = yd;
579  fA = xcC - yA;
580  A::Gradcomp();
581  if (differ(fA.val(), f)) botch("fA = xcC - yA", f, fA.val());
582  else if (differ(xcC.adj(), dx)) botch("d xcC - yA/dx", dx, xcC.adj());
583  else if (differ(yA.adj(), dy)) botch("d xcC - yA/dy", dy, yA.adj());
584  }
585  {
586  A::aval_reset();
587  xC = xd;
588  cA ycA(yd);
589  fA = xC - ycA;
590  A::Gradcomp();
591  if (differ(fA.val(), f)) botch("fA = xC - ycA", f, fA.val());
592  else if (differ(xC.adj(), dx)) botch("d xC - ycA/dx", dx, xC.adj());
593  else if (differ(ycA.adj(), dy)) botch("d xC - ycA/dy", dy, ycA.adj());
594  }
595  {
596  A::aval_reset();
597  cC xcC(xd);
598  cA ycA(yd);
599  fA = xcC - ycA;
600  A::Gradcomp();
601  if (differ(fA.val(), f)) botch("fA = xcC - ycA", f, fA.val());
602  else if (differ(xcC.adj(), dx)) botch("d xcC - ycA/dx", dx, xcC.adj());
603  else if (differ(ycA.adj(), dy)) botch("d xcC - ycA/dy", dy, ycA.adj());
604  }
605  xC = xd;
606  yC = yd;
607  fA = xC - yC;
608  A::Gradcomp();
609  if (differ(fA.val(), f)) botch("fA = xC - yC", f, fA.val());
610  else if (differ(xC.adj(), dx)) botch("d xC - yC/dx", dx, xC.adj());
611  else if (differ(yC.adj(), dy)) botch("d xC - yC/dy", dy, yC.adj());
612  {
613  A::aval_reset();
614  cC xcC(xd);
615  yC = yd;
616  fA = xcC - yC;
617  A::Gradcomp();
618  if (differ(fA.val(), f)) botch("fA = xcC - yC", f, fA.val());
619  else if (differ(xcC.adj(), dx)) botch("d xcC - yC/dx", dx, xcC.adj());
620  else if (differ(yC.adj(), dy)) botch("d xcC - yC/dy", dy, yC.adj());
621  }
622  {
623  A::aval_reset();
624  xC = xd;
625  cC ycC(yd);
626  fA = xC - ycC;
627  A::Gradcomp();
628  if (differ(fA.val(), f)) botch("fA = xC - ycC", f, fA.val());
629  else if (differ(xC.adj(), dx)) botch("d xC - ycC/dx", dx, xC.adj());
630  else if (differ(ycC.adj(), dy)) botch("d xC - ycC/dy", dy, ycC.adj());
631  }
632  {
633  A::aval_reset();
634  cC xcC(xd);
635  cC ycC(yd);
636  fA = xcC - ycC;
637  A::Gradcomp();
638  if (differ(fA.val(), f)) botch("fA = xcC - ycC", f, fA.val());
639  else if (differ(xcC.adj(), dx)) botch("d xcC - ycC/dx", dx, xcC.adj());
640  else if (differ(ycC.adj(), dy)) botch("d xcC - ycC/dy", dy, ycC.adj());
641  }
642  {
643  xC = xd;
644  Ai yAi(yd);
645  fA = xC - yAi;
646  A::Gradcomp();
647  if (differ(fA.val(), f)) botch("fA = xC - yAi", f, fA.val());
648  else if (differ(xC.adj(), dx)) botch("d xC - yAi/dx", dx, xC.adj());
649  else if (differ(yAi.aval, dy)) botch("d xC - yAi/dy", dy, yAi.aval);
650  }
651  {
652  A::aval_reset();
653  cC xcC(xd);
654  Ai yAi(yd);
655  fA = xcC - yAi;
656  A::Gradcomp();
657  if (differ(fA.val(), f)) botch("fA = xcC - yAi", f, fA.val());
658  else if (differ(xcC.adj(), dx)) botch("d xcC - yAi/dx", dx, xcC.adj());
659  else if (differ(yAi.aval, dy)) botch("d xcC - yAi/dy", dy, yAi.aval);
660  }
661  {
662  A::aval_reset();
663  xC = xd;
664  cAi ycAi(yd);
665  fA = xC - ycAi;
666  A::Gradcomp();
667  if (differ(fA.val(), f)) botch("fA = xC - ycAi", f, fA.val());
668  else if (differ(xC.adj(), dx)) botch("d xC - ycAi/dx", dx, xC.adj());
669  else if (differ(ycAi.aval, dy)) botch("d xC - ycAi/dy", dy, ycAi.aval);
670  }
671  {
672  A::aval_reset();
673  cC xcC(xd);
674  cAi ycAi(yd);
675  fA = xcC - ycAi;
676  A::Gradcomp();
677  if (differ(fA.val(), f)) botch("fA = xcC - ycAi", f, fA.val());
678  else if (differ(xcC.adj(), dx)) botch("d xcC - ycAi/dx", dx, xcC.adj());
679  else if (differ(ycAi.aval, dy)) botch("d xcC - ycAi/dy", dy, ycAi.aval);
680  }
681  xC = xd;
682  fA = xC - yd;
683  A::Gradcomp();
684  if (differ(fA.val(), f)) botch("fA = xC - yd", f, fA.val());
685  else if (differ(xC.adj(), dx)) botch("d xC - yd/dx", dx, xC.adj());
686  {
687  A::aval_reset();
688  cC xcC(xd);
689  fA = xcC - yd;
690  A::Gradcomp();
691  if (differ(fA.val(), f)) botch("fA = xcC - yd", f, fA.val());
692  else if (differ(xcC.adj(), dx)) botch("d xcC - yd/dx", dx, xcC.adj());
693  }
694  xC = xd;
695  yL = (long)yd;
696  fA = xC - yL;
697  A::Gradcomp();
698  if (differ(fA.val(), f)) botch("fA = xC - yL", f, fA.val());
699  else if (differ(xC.adj(), dx)) botch("d xC - yL/dx", dx, xC.adj());
700  {
701  A::aval_reset();
702  cC xcC(xd);
703  yL = (long)yd;
704  fA = xcC - yL;
705  A::Gradcomp();
706  if (differ(fA.val(), f)) botch("fA = xcC - yL", f, fA.val());
707  else if (differ(xcC.adj(), dx)) botch("d xcC - yL/dx", dx, xcC.adj());
708  }
709  xC = xd;
710  yi = (int)yd;
711  fA = xC - yi;
712  A::Gradcomp();
713  if (differ(fA.val(), f)) botch("fA = xC - yi", f, fA.val());
714  else if (differ(xC.adj(), dx)) botch("d xC - yi/dx", dx, xC.adj());
715  {
716  A::aval_reset();
717  cC xcC(xd);
718  yi = (int)yd;
719  fA = xcC - yi;
720  A::Gradcomp();
721  if (differ(fA.val(), f)) botch("fA = xcC - yi", f, fA.val());
722  else if (differ(xcC.adj(), dx)) botch("d xcC - yi/dx", dx, xcC.adj());
723  }
724  {
725  Ai xAi(xd);
726  yAI = yd;
727  fA = xAi - yAI;
728  A::Gradcomp();
729  if (differ(fA.val(), f)) botch("fA = xAi - yAI", f, fA.val());
730  else if (differ(xAi.aval, dx)) botch("d xAi - yAI/dx", dx, xAi.aval);
731  else if (differ(yAI.adj(), dy)) botch("d xAi - yAI/dy", dy, yAI.adj());
732  }
733  {
734  A::aval_reset();
735  cAi xcAi(xd);
736  yAI = yd;
737  fA = xcAi - yAI;
738  A::Gradcomp();
739  if (differ(fA.val(), f)) botch("fA = xcAi - yAI", f, fA.val());
740  else if (differ(xcAi.aval, dx)) botch("d xcAi - yAI/dx", dx, xcAi.aval);
741  else if (differ(yAI.adj(), dy)) botch("d xcAi - yAI/dy", dy, yAI.adj());
742  }
743  {
744  A::aval_reset();
745  Ai xAi(xd);
746  cAI ycAI(yd);
747  fA = xAi - ycAI;
748  A::Gradcomp();
749  if (differ(fA.val(), f)) botch("fA = xAi - ycAI", f, fA.val());
750  else if (differ(xAi.aval, dx)) botch("d xAi - ycAI/dx", dx, xAi.aval);
751  else if (differ(ycAI.adj(), dy)) botch("d xAi - ycAI/dy", dy, ycAI.adj());
752  }
753  {
754  Ai xAi(xd);
755  yA = yd;
756  fA = xAi - yA;
757  A::Gradcomp();
758  if (differ(fA.val(), f)) botch("fA = xAi - yA", f, fA.val());
759  else if (differ(xAi.aval, dx)) botch("d xAi - yA/dx", dx, xAi.aval);
760  else if (differ(yA.adj(), dy)) botch("d xAi - yA/dy", dy, yA.adj());
761  }
762  {
763  A::aval_reset();
764  cAi xcAi(xd);
765  yA = yd;
766  fA = xcAi - yA;
767  A::Gradcomp();
768  if (differ(fA.val(), f)) botch("fA = xcAi - yA", f, fA.val());
769  else if (differ(xcAi.aval, dx)) botch("d xcAi - yA/dx", dx, xcAi.aval);
770  else if (differ(yA.adj(), dy)) botch("d xcAi - yA/dy", dy, yA.adj());
771  }
772  {
773  A::aval_reset();
774  Ai xAi(xd);
775  cA ycA(yd);
776  fA = xAi - ycA;
777  A::Gradcomp();
778  if (differ(fA.val(), f)) botch("fA = xAi - ycA", f, fA.val());
779  else if (differ(xAi.aval, dx)) botch("d xAi - ycA/dx", dx, xAi.aval);
780  else if (differ(ycA.adj(), dy)) botch("d xAi - ycA/dy", dy, ycA.adj());
781  }
782  {
783  Ai xAi(xd);
784  yC = yd;
785  fA = xAi - yC;
786  A::Gradcomp();
787  if (differ(fA.val(), f)) botch("fA = xAi - yC", f, fA.val());
788  else if (differ(xAi.aval, dx)) botch("d xAi - yC/dx", dx, xAi.aval);
789  else if (differ(yC.adj(), dy)) botch("d xAi - yC/dy", dy, yC.adj());
790  }
791  {
792  A::aval_reset();
793  cAi xcAi(xd);
794  yC = yd;
795  fA = xcAi - yC;
796  A::Gradcomp();
797  if (differ(fA.val(), f)) botch("fA = xcAi - yC", f, fA.val());
798  else if (differ(xcAi.aval, dx)) botch("d xcAi - yC/dx", dx, xcAi.aval);
799  else if (differ(yC.adj(), dy)) botch("d xcAi - yC/dy", dy, yC.adj());
800  }
801  {
802  A::aval_reset();
803  Ai xAi(xd);
804  cC ycC(yd);
805  fA = xAi - ycC;
806  A::Gradcomp();
807  if (differ(fA.val(), f)) botch("fA = xAi - ycC", f, fA.val());
808  else if (differ(xAi.aval, dx)) botch("d xAi - ycC/dx", dx, xAi.aval);
809  else if (differ(ycC.adj(), dy)) botch("d xAi - ycC/dy", dy, ycC.adj());
810  }
811  {
812  Ai xAi(xd);
813  Ai yAi(yd);
814  fA = xAi - yAi;
815  A::Gradcomp();
816  if (differ(fA.val(), f)) botch("fA = xAi - yAi", f, fA.val());
817  else if (differ(xAi.aval, dx)) botch("d xAi - yAi/dx", dx, xAi.aval);
818  else if (differ(yAi.aval, dy)) botch("d xAi - yAi/dy", dy, yAi.aval);
819  }
820  {
821  A::aval_reset();
822  cAi xcAi(xd);
823  Ai yAi(yd);
824  fA = xcAi - yAi;
825  A::Gradcomp();
826  if (differ(fA.val(), f)) botch("fA = xcAi - yAi", f, fA.val());
827  else if (differ(xcAi.aval, dx)) botch("d xcAi - yAi/dx", dx, xcAi.aval);
828  else if (differ(yAi.aval, dy)) botch("d xcAi - yAi/dy", dy, yAi.aval);
829  }
830  {
831  A::aval_reset();
832  Ai xAi(xd);
833  cAi ycAi(yd);
834  fA = xAi - ycAi;
835  A::Gradcomp();
836  if (differ(fA.val(), f)) botch("fA = xAi - ycAi", f, fA.val());
837  else if (differ(xAi.aval, dx)) botch("d xAi - ycAi/dx", dx, xAi.aval);
838  else if (differ(ycAi.aval, dy)) botch("d xAi - ycAi/dy", dy, ycAi.aval);
839  }
840  {
841  Ai xAi(xd);
842  fA = xAi - yd;
843  A::Gradcomp();
844  if (differ(fA.val(), f)) botch("fA = xAi - yd", f, fA.val());
845  else if (differ(xAi.aval, dx)) botch("d xAi - yd/dx", dx, xAi.aval);
846  }
847  {
848  A::aval_reset();
849  cAi xcAi(xd);
850  fA = xcAi - yd;
851  A::Gradcomp();
852  if (differ(fA.val(), f)) botch("fA = xcAi - yd", f, fA.val());
853  else if (differ(xcAi.aval, dx)) botch("d xcAi - yd/dx", dx, xcAi.aval);
854  }
855  {
856  Ai xAi(xd);
857  yL = (long)yd;
858  fA = xAi - yL;
859  A::Gradcomp();
860  if (differ(fA.val(), f)) botch("fA = xAi - yL", f, fA.val());
861  else if (differ(xAi.aval, dx)) botch("d xAi - yL/dx", dx, xAi.aval);
862  }
863  {
864  A::aval_reset();
865  cAi xcAi(xd);
866  yL = (long)yd;
867  fA = xcAi - yL;
868  A::Gradcomp();
869  if (differ(fA.val(), f)) botch("fA = xcAi - yL", f, fA.val());
870  else if (differ(xcAi.aval, dx)) botch("d xcAi - yL/dx", dx, xcAi.aval);
871  }
872  {
873  Ai xAi(xd);
874  yi = (int)yd;
875  fA = xAi - yi;
876  A::Gradcomp();
877  if (differ(fA.val(), f)) botch("fA = xAi - yi", f, fA.val());
878  else if (differ(xAi.aval, dx)) botch("d xAi - yi/dx", dx, xAi.aval);
879  }
880  {
881  A::aval_reset();
882  cAi xcAi(xd);
883  yi = (int)yd;
884  fA = xcAi - yi;
885  A::Gradcomp();
886  if (differ(fA.val(), f)) botch("fA = xcAi - yi", f, fA.val());
887  else if (differ(xcAi.aval, dx)) botch("d xcAi - yi/dx", dx, xcAi.aval);
888  }
889  yAI = yd;
890  fA = xd - yAI;
891  A::Gradcomp();
892  if (differ(fA.val(), f)) botch("fA = xd - yAI", f, fA.val());
893  else if (differ(yAI.adj(), dy)) botch("d xd - yAI/dy", dy, yAI.adj());
894  {
895  A::aval_reset();
896  cAI ycAI(yd);
897  fA = xd - ycAI;
898  A::Gradcomp();
899  if (differ(fA.val(), f)) botch("fA = xd - ycAI", f, fA.val());
900  else if (differ(ycAI.adj(), dy)) botch("d xd - ycAI/dy", dy, ycAI.adj());
901  }
902  yA = yd;
903  fA = xd - yA;
904  A::Gradcomp();
905  if (differ(fA.val(), f)) botch("fA = xd - yA", f, fA.val());
906  else if (differ(yA.adj(), dy)) botch("d xd - yA/dy", dy, yA.adj());
907  {
908  A::aval_reset();
909  cA ycA(yd);
910  fA = xd - ycA;
911  A::Gradcomp();
912  if (differ(fA.val(), f)) botch("fA = xd - ycA", f, fA.val());
913  else if (differ(ycA.adj(), dy)) botch("d xd - ycA/dy", dy, ycA.adj());
914  }
915  yC = yd;
916  fA = xd - yC;
917  A::Gradcomp();
918  if (differ(fA.val(), f)) botch("fA = xd - yC", f, fA.val());
919  else if (differ(yC.adj(), dy)) botch("d xd - yC/dy", dy, yC.adj());
920  {
921  A::aval_reset();
922  cC ycC(yd);
923  fA = xd - ycC;
924  A::Gradcomp();
925  if (differ(fA.val(), f)) botch("fA = xd - ycC", f, fA.val());
926  else if (differ(ycC.adj(), dy)) botch("d xd - ycC/dy", dy, ycC.adj());
927  }
928  {
929  Ai yAi(yd);
930  fA = xd - yAi;
931  A::Gradcomp();
932  if (differ(fA.val(), f)) botch("fA = xd - yAi", f, fA.val());
933  else if (differ(yAi.aval, dy)) botch("d xd - yAi/dy", dy, yAi.aval);
934  }
935  {
936  A::aval_reset();
937  cAi ycAi(yd);
938  fA = xd - ycAi;
939  A::Gradcomp();
940  if (differ(fA.val(), f)) botch("fA = xd - ycAi", f, fA.val());
941  else if (differ(ycAi.aval, dy)) botch("d xd - ycAi/dy", dy, ycAi.aval);
942  }
943  xL = (long)xd;
944  yAI = yd;
945  fA = xL - yAI;
946  A::Gradcomp();
947  if (differ(fA.val(), f)) botch("fA = xL - yAI", f, fA.val());
948  else if (differ(yAI.adj(), dy)) botch("d xL - yAI/dy", dy, yAI.adj());
949  {
950  A::aval_reset();
951  xL = (long)xd;
952  cAI ycAI(yd);
953  fA = xL - ycAI;
954  A::Gradcomp();
955  if (differ(fA.val(), f)) botch("fA = xL - ycAI", f, fA.val());
956  else if (differ(ycAI.adj(), dy)) botch("d xL - ycAI/dy", dy, ycAI.adj());
957  }
958  xL = (long)xd;
959  yA = yd;
960  fA = xL - yA;
961  A::Gradcomp();
962  if (differ(fA.val(), f)) botch("fA = xL - yA", f, fA.val());
963  else if (differ(yA.adj(), dy)) botch("d xL - yA/dy", dy, yA.adj());
964  {
965  A::aval_reset();
966  xL = (long)xd;
967  cA ycA(yd);
968  fA = xL - ycA;
969  A::Gradcomp();
970  if (differ(fA.val(), f)) botch("fA = xL - ycA", f, fA.val());
971  else if (differ(ycA.adj(), dy)) botch("d xL - ycA/dy", dy, ycA.adj());
972  }
973  xL = (long)xd;
974  yC = yd;
975  fA = xL - yC;
976  A::Gradcomp();
977  if (differ(fA.val(), f)) botch("fA = xL - yC", f, fA.val());
978  else if (differ(yC.adj(), dy)) botch("d xL - yC/dy", dy, yC.adj());
979  {
980  A::aval_reset();
981  xL = (long)xd;
982  cC ycC(yd);
983  fA = xL - ycC;
984  A::Gradcomp();
985  if (differ(fA.val(), f)) botch("fA = xL - ycC", f, fA.val());
986  else if (differ(ycC.adj(), dy)) botch("d xL - ycC/dy", dy, ycC.adj());
987  }
988  {
989  xL = (long)xd;
990  Ai yAi(yd);
991  fA = xL - yAi;
992  A::Gradcomp();
993  if (differ(fA.val(), f)) botch("fA = xL - yAi", f, fA.val());
994  else if (differ(yAi.aval, dy)) botch("d xL - yAi/dy", dy, yAi.aval);
995  }
996  {
997  A::aval_reset();
998  xL = (long)xd;
999  cAi ycAi(yd);
1000  fA = xL - ycAi;
1001  A::Gradcomp();
1002  if (differ(fA.val(), f)) botch("fA = xL - ycAi", f, fA.val());
1003  else if (differ(ycAi.aval, dy)) botch("d xL - ycAi/dy", dy, ycAi.aval);
1004  }
1005  xi = (int)xd;
1006  yAI = yd;
1007  fA = xi - yAI;
1008  A::Gradcomp();
1009  if (differ(fA.val(), f)) botch("fA = xi - yAI", f, fA.val());
1010  else if (differ(yAI.adj(), dy)) botch("d xi - yAI/dy", dy, yAI.adj());
1011  {
1012  A::aval_reset();
1013  xi = (int)xd;
1014  cAI ycAI(yd);
1015  fA = xi - ycAI;
1016  A::Gradcomp();
1017  if (differ(fA.val(), f)) botch("fA = xi - ycAI", f, fA.val());
1018  else if (differ(ycAI.adj(), dy)) botch("d xi - ycAI/dy", dy, ycAI.adj());
1019  }
1020  xi = (int)xd;
1021  yA = yd;
1022  fA = xi - yA;
1023  A::Gradcomp();
1024  if (differ(fA.val(), f)) botch("fA = xi - yA", f, fA.val());
1025  else if (differ(yA.adj(), dy)) botch("d xi - yA/dy", dy, yA.adj());
1026  {
1027  A::aval_reset();
1028  xi = (int)xd;
1029  cA ycA(yd);
1030  fA = xi - ycA;
1031  A::Gradcomp();
1032  if (differ(fA.val(), f)) botch("fA = xi - ycA", f, fA.val());
1033  else if (differ(ycA.adj(), dy)) botch("d xi - ycA/dy", dy, ycA.adj());
1034  }
1035  xi = (int)xd;
1036  yC = yd;
1037  fA = xi - yC;
1038  A::Gradcomp();
1039  if (differ(fA.val(), f)) botch("fA = xi - yC", f, fA.val());
1040  else if (differ(yC.adj(), dy)) botch("d xi - yC/dy", dy, yC.adj());
1041  {
1042  A::aval_reset();
1043  xi = (int)xd;
1044  cC ycC(yd);
1045  fA = xi - ycC;
1046  A::Gradcomp();
1047  if (differ(fA.val(), f)) botch("fA = xi - ycC", f, fA.val());
1048  else if (differ(ycC.adj(), dy)) botch("d xi - ycC/dy", dy, ycC.adj());
1049  }
1050  {
1051  xi = (int)xd;
1052  Ai yAi(yd);
1053  fA = xi - yAi;
1054  A::Gradcomp();
1055  if (differ(fA.val(), f)) botch("fA = xi - yAi", f, fA.val());
1056  else if (differ(yAi.aval, dy)) botch("d xi - yAi/dy", dy, yAi.aval);
1057  }
1058  {
1059  A::aval_reset();
1060  xi = (int)xd;
1061  cAi ycAi(yd);
1062  fA = xi - ycAi;
1063  A::Gradcomp();
1064  if (differ(fA.val(), f)) botch("fA = xi - ycAi", f, fA.val());
1065  else if (differ(ycAi.aval, dy)) botch("d xi - ycAi/dy", dy, ycAi.aval);
1066  }
1067 
1068 
1069  if (!rc) // chatter for cppunit test, which cannot tolerate silence
1070 
1071  printf("OK\n");
1072 
1073  return rc;
1074 
1075  }
Ai
ADT_RAD ADvari< double > Ai
Definition: tradoptest_01.cpp:25
botch
int botch(ExpectedAnswer *e, const char *partial, double got, double wanted)
Definition: tradvec_example.cpp:74
f
void f()
dx
expr expr dx(i)
Definition: Sacado_ELRFad_Ops.hpp:205
AI
ADT_RAD IndepADvar< double > AI
Definition: tradoptest_01.cpp:19
cAI
const ADT_RAD IndepADvar< double > cAI
Definition: tradoptest_01.cpp:27
Ai
#define Ai
Definition: Sacado_rad.hpp:575
a
a
Definition: Sacado_CacheFad_Ops.hpp:426
C
#define C(x)
Definition: Sacado_trad.hpp:2627
A
#define A
Definition: Sacado_rad.hpp:572
main
int main()
Definition: ad_example.cpp:191
Sacado_trad.hpp
cAi
const ADT_RAD ADvari< double > cAi
Definition: tradoptest_01.cpp:33
differ
static int differ(double a, double b)
Definition: hesopcheck.cpp:147
ADT_RAD
#define ADT_RAD
Definition: tradoptest_24.cpp:7
cC
const ADT_RAD ConstADvar< double > cC
Definition: tradoptest_01.cpp:31
tol
const double tol
Definition: tradoptest_01.cpp:61
rc
static int rc
Definition: tradoptest_24.cpp:35
AI
#define AI
Definition: Sacado_rad.hpp:576
cA
const ADT_RAD ADvar< double > cA
Definition: tradoptest_01.cpp:29
Generated by   doxygen 1.8.5