rcce/lu/verify.c

8 #define xcrref(m) xcrref[m-1]
15         double xcrref[5],xceref[5],xciref, xcrdif[5], xcedif[5], xcidif,
28            xcrref(m) = 1.0;
41          xcrref(1) = 1.6196343210976702e-02;
42          xcrref(2) = 2.1976745164821318e-03;
43          xcrref(3) = 1.5179927653399185e-03;
44          xcrref(4) = 1.5029584435994323e-03;
45          xcrref(5) = 3.4264073155896461e-02;
72            xcrref(1) =   0.1236511638192e+02;
73            xcrref(2) =   0.1317228477799e+01;
74            xcrref(3) =   0.2550120713095e+01;
75            xcrref(4) =   0.2326187750252e+01;
76            xcrref(5) =   0.2826799444189e+02;
102          xcrref(1) = 7.7902107606689367e+02;
103          xcrref(2) = 6.3402765259692870e+01;
104          xcrref(3) = 1.9499249727292479e+02;
105          xcrref(4) = 1.7845301160418537e+02;
106          xcrref(5) = 1.8384760349464247e+03;
134          xcrref(1) = 3.5532672969982736e+03;
135          xcrref(2) = 2.6214750795310692e+02;
136          xcrref(3) = 8.8333721850952190e+02;
137          xcrref(4) = 7.7812774739425265e+02;
138          xcrref(5) = 7.3087969592545314e+03;
166          xcrref(1) = 1.03766980323537846e+04;
167          xcrref(2) = 8.92212458801008552e+02;
168          xcrref(3) = 2.56238814582660871e+03;
169          xcrref(4) = 2.19194343857831427e+03;
170          xcrref(5) = 1.78078057261061185e+04;
198          xcrref(1) = 0.4868417937025e+05;
199          xcrref(2) = 0.4696371050071e+04;
200          xcrref(3) = 0.1218114549776e+05;
201          xcrref(4) = 0.1033801493461e+05;
202          xcrref(5) = 0.7142398413817e+05;
230          xcrref(1) = 0.2099641687874e+06;
231          xcrref(2) = 0.2130403143165e+05;
232          xcrref(3) = 0.5319228789371e+05;
233          xcrref(4) = 0.4509761639833e+05;
234          xcrref(5) = 0.2932360006590e+06;
264            xcrdif(m) = fabs((xcr(m)-xcrref(m))/xcrref(m));
293                       m,xcr(m),xcrref(m),xcrdif(m));
297                       m,xcr(m),xcrref(m),xcrdif(m));