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