7#define ADT_RAD Sacado::Rad::
49botch(
const char *what,
double wanted,
double got)
53 printf(
"%s: expected %g, got %g, diff = %.2g\n", what, wanted, got, wanted-got);
61 const double tol = 5e-16;
105#ifdef RAD_AUTO_AD_Const
131 double dx, dy, f, xd, yd;
144 xd = 7.; yd = 4; f = 3.;
dx = 1.; dy = -1.;
147 fA = operator-(xAI,yAI);
149 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yAI)", f, fA.val());
150 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yAI)/dx",
dx, xAI.adj());
151 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xAI,yAI)/dy", dy, yAI.adj());
156 fA = operator-(xcAI,yAI);
158 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yAI)", f, fA.val());
159 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yAI)/dx",
dx, xcAI.adj());
160 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xcAI,yAI)/dy", dy, yAI.adj());
166 fA = operator-(xAI,ycAI);
168 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,ycAI)", f, fA.val());
169 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,ycAI)/dx",
dx, xAI.adj());
170 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xAI,ycAI)/dy", dy, ycAI.adj());
176 fA = operator-(xcAI,ycAI);
178 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,ycAI)", f, fA.val());
179 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,ycAI)/dx",
dx, xcAI.adj());
180 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xcAI,ycAI)/dy", dy, ycAI.adj());
184 fA = operator-(xAI,yA);
186 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yA)", f, fA.val());
187 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yA)/dx",
dx, xAI.adj());
188 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xAI,yA)/dy", dy, yA.adj());
193 fA = operator-(xcAI,yA);
195 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yA)", f, fA.val());
196 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yA)/dx",
dx, xcAI.adj());
197 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xcAI,yA)/dy", dy, yA.adj());
203 fA = operator-(xAI,ycA);
205 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,ycA)", f, fA.val());
206 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,ycA)/dx",
dx, xAI.adj());
207 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xAI,ycA)/dy", dy, ycA.adj());
213 fA = operator-(xcAI,ycA);
215 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,ycA)", f, fA.val());
216 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,ycA)/dx",
dx, xcAI.adj());
217 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xcAI,ycA)/dy", dy, ycA.adj());
221 fA = operator-(xAI,yC);
223 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yC)", f, fA.val());
224 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yC)/dx",
dx, xAI.adj());
225 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xAI,yC)/dy", dy, yC.adj());
230 fA = operator-(xcAI,yC);
232 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yC)", f, fA.val());
233 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yC)/dx",
dx, xcAI.adj());
234 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xcAI,yC)/dy", dy, yC.adj());
240 fA = operator-(xAI,ycC);
242 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,ycC)", f, fA.val());
243 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,ycC)/dx",
dx, xAI.adj());
244 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xAI,ycC)/dy", dy, ycC.adj());
250 fA = operator-(xcAI,ycC);
252 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,ycC)", f, fA.val());
253 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,ycC)/dx",
dx, xcAI.adj());
254 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xcAI,ycC)/dy", dy, ycC.adj());
259 fA = operator-(xAI,yAi);
261 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yAi)", f, fA.val());
262 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yAi)/dx",
dx, xAI.adj());
263 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xAI,yAi)/dy", dy, yAi.aval);
269 fA = operator-(xcAI,yAi);
271 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yAi)", f, fA.val());
272 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yAi)/dx",
dx, xcAI.adj());
273 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xcAI,yAi)/dy", dy, yAi.aval);
279 fA = operator-(xAI,ycAi);
281 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,ycAi)", f, fA.val());
282 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,ycAi)/dx",
dx, xAI.adj());
283 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xAI,ycAi)/dy", dy, ycAi.aval);
289 fA = operator-(xcAI,ycAi);
291 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,ycAi)", f, fA.val());
292 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,ycAi)/dx",
dx, xcAI.adj());
293 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xcAI,ycAi)/dy", dy, ycAi.aval);
296 fA = operator-(xAI,yd);
298 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yd)", f, fA.val());
299 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yd)/dx",
dx, xAI.adj());
303 fA = operator-(xcAI,yd);
305 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yd)", f, fA.val());
306 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yd)/dx",
dx, xcAI.adj());
310 fA = operator-(xAI,yL);
312 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yL)", f, fA.val());
313 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yL)/dx",
dx, xAI.adj());
318 fA = operator-(xcAI,yL);
320 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yL)", f, fA.val());
321 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yL)/dx",
dx, xcAI.adj());
325 fA = operator-(xAI,yi);
327 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAI,yi)", f, fA.val());
328 else if (
differ(xAI.adj(),
dx))
botch(
"d operator-(xAI,yi)/dx",
dx, xAI.adj());
333 fA = operator-(xcAI,yi);
335 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAI,yi)", f, fA.val());
336 else if (
differ(xcAI.adj(),
dx))
botch(
"d operator-(xcAI,yi)/dx",
dx, xcAI.adj());
340 fA = operator-(xA,yAI);
342 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yAI)", f, fA.val());
343 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yAI)/dx",
dx, xA.adj());
344 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xA,yAI)/dy", dy, yAI.adj());
349 fA = operator-(xcA,yAI);
351 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yAI)", f, fA.val());
352 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yAI)/dx",
dx, xcA.adj());
353 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xcA,yAI)/dy", dy, yAI.adj());
359 fA = operator-(xA,ycAI);
361 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,ycAI)", f, fA.val());
362 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,ycAI)/dx",
dx, xA.adj());
363 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xA,ycAI)/dy", dy, ycAI.adj());
369 fA = operator-(xcA,ycAI);
371 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,ycAI)", f, fA.val());
372 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,ycAI)/dx",
dx, xcA.adj());
373 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xcA,ycAI)/dy", dy, ycAI.adj());
377 fA = operator-(xA,yA);
379 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yA)", f, fA.val());
380 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yA)/dx",
dx, xA.adj());
381 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xA,yA)/dy", dy, yA.adj());
386 fA = operator-(xcA,yA);
388 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yA)", f, fA.val());
389 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yA)/dx",
dx, xcA.adj());
390 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xcA,yA)/dy", dy, yA.adj());
396 fA = operator-(xA,ycA);
398 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,ycA)", f, fA.val());
399 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,ycA)/dx",
dx, xA.adj());
400 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xA,ycA)/dy", dy, ycA.adj());
406 fA = operator-(xcA,ycA);
408 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,ycA)", f, fA.val());
409 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,ycA)/dx",
dx, xcA.adj());
410 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xcA,ycA)/dy", dy, ycA.adj());
414 fA = operator-(xA,yC);
416 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yC)", f, fA.val());
417 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yC)/dx",
dx, xA.adj());
418 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xA,yC)/dy", dy, yC.adj());
423 fA = operator-(xcA,yC);
425 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yC)", f, fA.val());
426 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yC)/dx",
dx, xcA.adj());
427 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xcA,yC)/dy", dy, yC.adj());
433 fA = operator-(xA,ycC);
435 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,ycC)", f, fA.val());
436 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,ycC)/dx",
dx, xA.adj());
437 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xA,ycC)/dy", dy, ycC.adj());
443 fA = operator-(xcA,ycC);
445 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,ycC)", f, fA.val());
446 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,ycC)/dx",
dx, xcA.adj());
447 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xcA,ycC)/dy", dy, ycC.adj());
452 fA = operator-(xA,yAi);
454 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yAi)", f, fA.val());
455 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yAi)/dx",
dx, xA.adj());
456 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xA,yAi)/dy", dy, yAi.aval);
462 fA = operator-(xcA,yAi);
464 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yAi)", f, fA.val());
465 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yAi)/dx",
dx, xcA.adj());
466 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xcA,yAi)/dy", dy, yAi.aval);
472 fA = operator-(xA,ycAi);
474 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,ycAi)", f, fA.val());
475 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,ycAi)/dx",
dx, xA.adj());
476 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xA,ycAi)/dy", dy, ycAi.aval);
482 fA = operator-(xcA,ycAi);
484 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,ycAi)", f, fA.val());
485 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,ycAi)/dx",
dx, xcA.adj());
486 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xcA,ycAi)/dy", dy, ycAi.aval);
489 fA = operator-(xA,yd);
491 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yd)", f, fA.val());
492 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yd)/dx",
dx, xA.adj());
496 fA = operator-(xcA,yd);
498 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yd)", f, fA.val());
499 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yd)/dx",
dx, xcA.adj());
503 fA = operator-(xA,yL);
505 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yL)", f, fA.val());
506 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yL)/dx",
dx, xA.adj());
511 fA = operator-(xcA,yL);
513 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yL)", f, fA.val());
514 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yL)/dx",
dx, xcA.adj());
518 fA = operator-(xA,yi);
520 if (
differ(fA.val(), f))
botch(
"fA = operator-(xA,yi)", f, fA.val());
521 else if (
differ(xA.adj(),
dx))
botch(
"d operator-(xA,yi)/dx",
dx, xA.adj());
526 fA = operator-(xcA,yi);
528 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcA,yi)", f, fA.val());
529 else if (
differ(xcA.adj(),
dx))
botch(
"d operator-(xcA,yi)/dx",
dx, xcA.adj());
533 fA = operator-(xC,yAI);
535 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yAI)", f, fA.val());
536 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yAI)/dx",
dx, xC.adj());
537 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xC,yAI)/dy", dy, yAI.adj());
542 fA = operator-(xcC,yAI);
544 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yAI)", f, fA.val());
545 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yAI)/dx",
dx, xcC.adj());
546 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xcC,yAI)/dy", dy, yAI.adj());
552 fA = operator-(xC,ycAI);
554 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,ycAI)", f, fA.val());
555 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,ycAI)/dx",
dx, xC.adj());
556 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xC,ycAI)/dy", dy, ycAI.adj());
562 fA = operator-(xcC,ycAI);
564 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,ycAI)", f, fA.val());
565 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,ycAI)/dx",
dx, xcC.adj());
566 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xcC,ycAI)/dy", dy, ycAI.adj());
570 fA = operator-(xC,yA);
572 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yA)", f, fA.val());
573 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yA)/dx",
dx, xC.adj());
574 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xC,yA)/dy", dy, yA.adj());
579 fA = operator-(xcC,yA);
581 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yA)", f, fA.val());
582 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yA)/dx",
dx, xcC.adj());
583 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xcC,yA)/dy", dy, yA.adj());
589 fA = operator-(xC,ycA);
591 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,ycA)", f, fA.val());
592 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,ycA)/dx",
dx, xC.adj());
593 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xC,ycA)/dy", dy, ycA.adj());
599 fA = operator-(xcC,ycA);
601 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,ycA)", f, fA.val());
602 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,ycA)/dx",
dx, xcC.adj());
603 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xcC,ycA)/dy", dy, ycA.adj());
607 fA = operator-(xC,yC);
609 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yC)", f, fA.val());
610 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yC)/dx",
dx, xC.adj());
611 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xC,yC)/dy", dy, yC.adj());
616 fA = operator-(xcC,yC);
618 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yC)", f, fA.val());
619 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yC)/dx",
dx, xcC.adj());
620 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xcC,yC)/dy", dy, yC.adj());
626 fA = operator-(xC,ycC);
628 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,ycC)", f, fA.val());
629 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,ycC)/dx",
dx, xC.adj());
630 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xC,ycC)/dy", dy, ycC.adj());
636 fA = operator-(xcC,ycC);
638 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,ycC)", f, fA.val());
639 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,ycC)/dx",
dx, xcC.adj());
640 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xcC,ycC)/dy", dy, ycC.adj());
645 fA = operator-(xC,yAi);
647 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yAi)", f, fA.val());
648 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yAi)/dx",
dx, xC.adj());
649 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xC,yAi)/dy", dy, yAi.aval);
655 fA = operator-(xcC,yAi);
657 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yAi)", f, fA.val());
658 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yAi)/dx",
dx, xcC.adj());
659 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xcC,yAi)/dy", dy, yAi.aval);
665 fA = operator-(xC,ycAi);
667 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,ycAi)", f, fA.val());
668 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,ycAi)/dx",
dx, xC.adj());
669 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xC,ycAi)/dy", dy, ycAi.aval);
675 fA = operator-(xcC,ycAi);
677 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,ycAi)", f, fA.val());
678 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,ycAi)/dx",
dx, xcC.adj());
679 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xcC,ycAi)/dy", dy, ycAi.aval);
682 fA = operator-(xC,yd);
684 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yd)", f, fA.val());
685 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yd)/dx",
dx, xC.adj());
689 fA = operator-(xcC,yd);
691 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yd)", f, fA.val());
692 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yd)/dx",
dx, xcC.adj());
696 fA = operator-(xC,yL);
698 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yL)", f, fA.val());
699 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yL)/dx",
dx, xC.adj());
704 fA = operator-(xcC,yL);
706 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yL)", f, fA.val());
707 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yL)/dx",
dx, xcC.adj());
711 fA = operator-(xC,yi);
713 if (
differ(fA.val(), f))
botch(
"fA = operator-(xC,yi)", f, fA.val());
714 else if (
differ(xC.adj(),
dx))
botch(
"d operator-(xC,yi)/dx",
dx, xC.adj());
719 fA = operator-(xcC,yi);
721 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcC,yi)", f, fA.val());
722 else if (
differ(xcC.adj(),
dx))
botch(
"d operator-(xcC,yi)/dx",
dx, xcC.adj());
727 fA = operator-(xAi,yAI);
729 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yAI)", f, fA.val());
730 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yAI)/dx",
dx, xAi.aval);
731 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xAi,yAI)/dy", dy, yAI.adj());
737 fA = operator-(xcAi,yAI);
739 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yAI)", f, fA.val());
740 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yAI)/dx",
dx, xcAi.aval);
741 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xcAi,yAI)/dy", dy, yAI.adj());
747 fA = operator-(xAi,ycAI);
749 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,ycAI)", f, fA.val());
750 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,ycAI)/dx",
dx, xAi.aval);
751 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xAi,ycAI)/dy", dy, ycAI.adj());
756 fA = operator-(xAi,yA);
758 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yA)", f, fA.val());
759 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yA)/dx",
dx, xAi.aval);
760 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xAi,yA)/dy", dy, yA.adj());
766 fA = operator-(xcAi,yA);
768 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yA)", f, fA.val());
769 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yA)/dx",
dx, xcAi.aval);
770 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xcAi,yA)/dy", dy, yA.adj());
776 fA = operator-(xAi,ycA);
778 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,ycA)", f, fA.val());
779 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,ycA)/dx",
dx, xAi.aval);
780 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xAi,ycA)/dy", dy, ycA.adj());
785 fA = operator-(xAi,yC);
787 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yC)", f, fA.val());
788 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yC)/dx",
dx, xAi.aval);
789 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xAi,yC)/dy", dy, yC.adj());
795 fA = operator-(xcAi,yC);
797 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yC)", f, fA.val());
798 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yC)/dx",
dx, xcAi.aval);
799 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xcAi,yC)/dy", dy, yC.adj());
805 fA = operator-(xAi,ycC);
807 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,ycC)", f, fA.val());
808 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,ycC)/dx",
dx, xAi.aval);
809 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xAi,ycC)/dy", dy, ycC.adj());
814 fA = operator-(xAi,yAi);
816 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yAi)", f, fA.val());
817 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yAi)/dx",
dx, xAi.aval);
818 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xAi,yAi)/dy", dy, yAi.aval);
824 fA = operator-(xcAi,yAi);
826 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yAi)", f, fA.val());
827 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yAi)/dx",
dx, xcAi.aval);
828 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xcAi,yAi)/dy", dy, yAi.aval);
834 fA = operator-(xAi,ycAi);
836 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,ycAi)", f, fA.val());
837 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,ycAi)/dx",
dx, xAi.aval);
838 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xAi,ycAi)/dy", dy, ycAi.aval);
842 fA = operator-(xAi,yd);
844 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yd)", f, fA.val());
845 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yd)/dx",
dx, xAi.aval);
850 fA = operator-(xcAi,yd);
852 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yd)", f, fA.val());
853 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yd)/dx",
dx, xcAi.aval);
858 fA = operator-(xAi,yL);
860 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yL)", f, fA.val());
861 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yL)/dx",
dx, xAi.aval);
867 fA = operator-(xcAi,yL);
869 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yL)", f, fA.val());
870 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yL)/dx",
dx, xcAi.aval);
875 fA = operator-(xAi,yi);
877 if (
differ(fA.val(), f))
botch(
"fA = operator-(xAi,yi)", f, fA.val());
878 else if (
differ(xAi.aval,
dx))
botch(
"d operator-(xAi,yi)/dx",
dx, xAi.aval);
884 fA = operator-(xcAi,yi);
886 if (
differ(fA.val(), f))
botch(
"fA = operator-(xcAi,yi)", f, fA.val());
887 else if (
differ(xcAi.aval,
dx))
botch(
"d operator-(xcAi,yi)/dx",
dx, xcAi.aval);
890 fA = operator-(xd,yAI);
892 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,yAI)", f, fA.val());
893 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xd,yAI)/dy", dy, yAI.adj());
897 fA = operator-(xd,ycAI);
899 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,ycAI)", f, fA.val());
900 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xd,ycAI)/dy", dy, ycAI.adj());
903 fA = operator-(xd,yA);
905 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,yA)", f, fA.val());
906 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xd,yA)/dy", dy, yA.adj());
910 fA = operator-(xd,ycA);
912 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,ycA)", f, fA.val());
913 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xd,ycA)/dy", dy, ycA.adj());
916 fA = operator-(xd,yC);
918 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,yC)", f, fA.val());
919 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xd,yC)/dy", dy, yC.adj());
923 fA = operator-(xd,ycC);
925 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,ycC)", f, fA.val());
926 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xd,ycC)/dy", dy, ycC.adj());
930 fA = operator-(xd,yAi);
932 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,yAi)", f, fA.val());
933 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xd,yAi)/dy", dy, yAi.aval);
938 fA = operator-(xd,ycAi);
940 if (
differ(fA.val(), f))
botch(
"fA = operator-(xd,ycAi)", f, fA.val());
941 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xd,ycAi)/dy", dy, ycAi.aval);
945 fA = operator-(xL,yAI);
947 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,yAI)", f, fA.val());
948 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xL,yAI)/dy", dy, yAI.adj());
953 fA = operator-(xL,ycAI);
955 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,ycAI)", f, fA.val());
956 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xL,ycAI)/dy", dy, ycAI.adj());
960 fA = operator-(xL,yA);
962 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,yA)", f, fA.val());
963 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xL,yA)/dy", dy, yA.adj());
968 fA = operator-(xL,ycA);
970 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,ycA)", f, fA.val());
971 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xL,ycA)/dy", dy, ycA.adj());
975 fA = operator-(xL,yC);
977 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,yC)", f, fA.val());
978 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xL,yC)/dy", dy, yC.adj());
983 fA = operator-(xL,ycC);
985 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,ycC)", f, fA.val());
986 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xL,ycC)/dy", dy, ycC.adj());
991 fA = operator-(xL,yAi);
993 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,yAi)", f, fA.val());
994 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xL,yAi)/dy", dy, yAi.aval);
1000 fA = operator-(xL,ycAi);
1002 if (
differ(fA.val(), f))
botch(
"fA = operator-(xL,ycAi)", f, fA.val());
1003 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xL,ycAi)/dy", dy, ycAi.aval);
1007 fA = operator-(xi,yAI);
1009 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,yAI)", f, fA.val());
1010 else if (
differ(yAI.adj(), dy))
botch(
"d operator-(xi,yAI)/dy", dy, yAI.adj());
1015 fA = operator-(xi,ycAI);
1017 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,ycAI)", f, fA.val());
1018 else if (
differ(ycAI.adj(), dy))
botch(
"d operator-(xi,ycAI)/dy", dy, ycAI.adj());
1022 fA = operator-(xi,yA);
1024 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,yA)", f, fA.val());
1025 else if (
differ(yA.adj(), dy))
botch(
"d operator-(xi,yA)/dy", dy, yA.adj());
1030 fA = operator-(xi,ycA);
1032 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,ycA)", f, fA.val());
1033 else if (
differ(ycA.adj(), dy))
botch(
"d operator-(xi,ycA)/dy", dy, ycA.adj());
1037 fA = operator-(xi,yC);
1039 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,yC)", f, fA.val());
1040 else if (
differ(yC.adj(), dy))
botch(
"d operator-(xi,yC)/dy", dy, yC.adj());
1045 fA = operator-(xi,ycC);
1047 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,ycC)", f, fA.val());
1048 else if (
differ(ycC.adj(), dy))
botch(
"d operator-(xi,ycC)/dy", dy, ycC.adj());
1053 fA = operator-(xi,yAi);
1055 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,yAi)", f, fA.val());
1056 else if (
differ(yAi.aval, dy))
botch(
"d operator-(xi,yAi)/dy", dy, yAi.aval);
1062 fA = operator-(xi,ycAi);
1064 if (
differ(fA.val(), f))
botch(
"fA = operator-(xi,ycAi)", f, fA.val());
1065 else if (
differ(ycAi.aval, dy))
botch(
"d operator-(xi,ycAi)/dy", dy, ycAi.aval);
ADT_RAD IndepADvar< double > AI
const ADT_RAD ConstADvar< double > cC
const ADT_RAD IndepADvar< double > cAI
const ADT_RAD ADvari< double > cAi
const ADT_RAD ADvar< double > cA
ADT_RAD ADvari< double > Ai
ADT_RAD IndepADvar< double > AI
const ADT_RAD ConstADvar< double > cC
const ADT_RAD IndepADvar< double > cAI
ADT_RAD ADvar< double > A
const ADT_RAD ADvari< double > cAi
ADT_RAD ConstADvar< double > C
void botch(const char *what, double wanted, double got)
const ADT_RAD ADvar< double > cA
ADT_RAD ADvari< double > Ai
int differ(double a, double b)