Sacado Package Browser (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
tradoptest_27.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 = 2.; yd = 3.; f = 1.; dx = 0.; dy = 0.;
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
rc
static int rc
Definition: tradoptest_27.cpp:42
cC
const ADT_RAD ConstADvar< double > cC
Definition: tradoptest_01.cpp:38
tol
const double tol
Definition: tradoptest_01.cpp:68
ADT_RAD
#define ADT_RAD
Definition: tradoptest_27.cpp:14
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