5 This test allows to analyze the quality of an adjoint operator associated
6 to some given direct operator. If the adjoint operator is approximated and
7 not given, the test measures the quality of the automatic approximation.
9 Using the "ScalarProduct" formula, one observes the residue R which is the
10 difference of two scalar products:
12 R(Alpha) = | < TangentF_X(dX) , Y > - < dX , AdjointF_X(Y) > |
14 which must remain constantly equal to zero to the accuracy of the calculation.
15 One takes dX0 = Normal(0,X) and dX = Alpha*dX0, where F is the calculation
16 operator. If it is given, Y must be in the image of F. If it is not given,
19 (Remark: numbers that are (about) under 2e-16 represent 0 to machine precision)
21 -------------------------------------------------------------
22 i Alpha ||X|| ||Y|| ||dX|| R(Alpha)
23 -------------------------------------------------------------
24 0 1e+00 2.236e+00 1.910e+01 3.536e+00 0.000e+00
25 1 1e-01 2.236e+00 1.910e+01 3.536e-01 0.000e+00
26 2 1e-02 2.236e+00 1.910e+01 3.536e-02 0.000e+00
27 3 1e-03 2.236e+00 1.910e+01 3.536e-03 0.000e+00
28 4 1e-04 2.236e+00 1.910e+01 3.536e-04 0.000e+00
29 5 1e-05 2.236e+00 1.910e+01 3.536e-05 0.000e+00
30 6 1e-06 2.236e+00 1.910e+01 3.536e-06 0.000e+00
31 7 1e-07 2.236e+00 1.910e+01 3.536e-07 0.000e+00
32 8 1e-08 2.236e+00 1.910e+01 3.536e-08 0.000e+00
33 9 1e-09 2.236e+00 1.910e+01 3.536e-09 0.000e+00
34 10 1e-10 2.236e+00 1.910e+01 3.536e-10 0.000e+00
35 11 1e-11 2.236e+00 1.910e+01 3.536e-11 0.000e+00
36 12 1e-12 2.236e+00 1.910e+01 3.536e-12 0.000e+00
37 -------------------------------------------------------------