cg_rc_prb_output.txt

14 January 2013 02:34:44 PM

CG_RC_PRB:
  C++ version
  Test the CG_RC library.

TEST01
  Use CG_RC on the 1, -2, 1 matrix.

  Number of iterations was 9
  Estimated error is 4.21596e-16
  Loo error is 6.66134e-16

     I      X(I)         X_EXACT(I)        B(I)

     0     7.80626e-17               0        0.309017
     1        0.309017        0.309017      -0.0302487
     2        0.587785        0.587785      -0.0575365
     3        0.809017        0.809017      -0.0791922
     4        0.951057        0.951057       -0.093096
     5               1               1       -0.097887
     6        0.951057        0.951057       -0.093096
     7        0.809017        0.809017      -0.0791922
     8        0.587785        0.587785      -0.0575365
     9        0.309017        0.309017      -0.0302487
    10    -4.79828e-17     1.22465e-16     3.88578e-16
    11       -0.309017       -0.309017       0.0302487
    12       -0.587785       -0.587785       0.0575365
    13       -0.809017       -0.809017       0.0791922
    14       -0.951057       -0.951057        0.093096
    15              -1              -1        0.097887
    16       -0.951057       -0.951057        0.093096
    17       -0.809017       -0.809017       0.0791922
    18       -0.587785       -0.587785       0.0575365
    19       -0.309017       -0.309017       0.0302487
    20    -4.10262e-16    -2.44929e-16       -0.309017

TEST02
  Use CG_RC to solve a linear system
  involving the Wathen matrix.

  NX = 5
  NY = 4
  N  = 21

  Number of iterations was 18
  Estimated error is 0.00730674
  Loo error is 4.22526e-05

     I      X(I)         X_EXACT(I)        B(I)

     0        0.218432        0.218418        -19.3113
     1        0.956316        0.956318         88.7317
     2        0.829512        0.829509        -36.3017
     3        0.561714        0.561695         90.0773
     4        0.415273        0.415307        -43.7437
     5        0.066156       0.0661187         54.0011
     6         0.25756        0.257578        -23.2089
     7        0.109949        0.109957         32.2168
     8       0.0438096        0.043829        -33.6528
     9        0.633973        0.633966         65.2136
    10         0.06171       0.0617272        -13.8091
    11        0.449528        0.449539          70.363
    12        0.401292        0.401306         154.809
    13        0.754653        0.754673         168.305
    14        0.797282        0.797287         118.754
    15      0.00182286      0.00183837         72.6417
    16        0.897471        0.897504         74.6033
    17        0.350743        0.350752        -37.7129
    18       0.0945624       0.0945448          107.05
    19       0.0136232       0.0136169        -109.428
    20        0.859112        0.859097         207.905
    21        0.840875        0.840847         -64.218
    22        0.123101        0.123104         99.2914
    23      0.00751016      0.00751236        -67.2724
    24        0.260308        0.260303         101.661
    25        0.912477        0.912484        -42.2354
    26        0.113677        0.113664         99.6324
    27        0.351621        0.351629        -30.9739
    28        0.822884        0.822887         78.3529
    29         0.26713        0.267132         144.803
    30        0.692067        0.692066         166.674
    31         0.56167        0.561662          111.98
    32         0.86121        0.861216         146.417
    33        0.453772        0.453794         49.2886
    34        0.911983        0.911977        -10.5401
    35        0.597914        0.597917         114.666
    36         0.18896        0.188955        -79.7521
    37        0.761504        0.761492         176.378
    38        0.396996        0.396988        -81.2853
    39        0.185316        0.185314         108.356
    40        0.574347        0.574366        -72.6881
    41        0.367041        0.367027         156.156
    42         0.61722        0.617205        -72.2069
    43        0.361536        0.361529         124.845
    44        0.212947         0.21293        -41.1752
    45        0.714457        0.714471         46.6628
    46        0.117699        0.117707         83.0738
    47        0.299328        0.299329         131.474
    48        0.824989        0.825003         189.243
    49        0.824647         0.82466         191.901
    50       0.0618704       0.0618618         62.9693
    51        0.710738        0.710781        -1.73137
    52       0.0882712       0.0882833         31.9862
    53         0.77798        0.777994        -19.4569
    54        0.745317        0.745303         131.412
    55        0.308673        0.308675        -61.8462
    56        0.899362        0.899373         151.142
    57        0.763542        0.763537         -50.183
    58        0.761737        0.761731         163.776
    59        0.406984         0.40697        -88.3359
    60        0.938776        0.938749         202.851
    61        0.562085        0.562088        -47.2525
    62       0.0178467         0.01782        -2.62116
    63        0.501099        0.501103         66.5544
    64       0.0418957       0.0419093         68.9155
    65        0.368859        0.368851           96.42
    66        0.271712        0.271724          138.43
    67        0.858573        0.858573         114.962
    68       0.0290067       0.0290366        -1.53808
    69       0.0174254       0.0174423        -1.44348
    70        0.152396        0.152384        -11.5233
    71        0.114317        0.114319         32.6938
    72        0.353917        0.353907        -24.8832
    73        0.119309        0.119308         32.1821
    74        0.206644        0.206653        -28.0064
    75        0.212946        0.212924          38.958
    76        0.612979        0.612948        -40.6021
    77        0.809533        0.809519         107.951
    78        0.587116         0.58709        -27.9974

CG_RC_PRB:
  Normal end of execution.

14 January 2013 02:34:44 PM