Sacado Package Browser (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
tradoptest_25.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 = 6.; yd = 7.; f = 42.; dx = 7.; dy = 6.;
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_25.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_25.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