Bug#1014514: Armadillo library doesnt work correct anymore after update to Debian 11
Schlager, Andreas
andreas.schlager at airbus.com
Thu Jul 7 12:51:47 BST 2022
Package: libarmadillo10
Version: 1:10.1.2+dfsg-6+b1
After updating to Debian 11, the armadillo library for C++ doesn't work properly anymore (e.g. matrix inversions deliver obviously wrong results), as can be seen in the example code and results attached.
Short note about the problem: Assume you have a square matrix B, then B*inv(B) should yield the unit matrix. In armadillo this worked without problems until the update to Debian 11. Now I get a completely random result for B*inv(B). The larger the dimension of B the more severe the problem is. There might be a numerical issue that didn't exist in older Debian versions.
According to the Armadillo-team, this is not a problem of Armadillo itself but might be the problem of the versions of libopenblas, liblapack or libblas used by Debian 11.
One of these libraries might be corrupt or incompatible with Debian 11.
Operating System: Debian GNU/Linux 11 (bullseye)
Kernel: Linux 5.10.0-14-amd64
Library Version Description
------------------------------------------------------------------------------------------------------------------------------------------------------
libarmadillo10 1:10.1.2+dfsg-6+b1 streamlined C++ linear algebra library
libblas3:amd64 3.9.0-3 Basic Linear Algebra Reference implementations, shared library
liblapack3:amd64 3.9.0-3 Library of linear algebra routines 3 - shared version
libopenblas0-pthread:amd64 0.3.13+ds-3 Optimized BLAS (linear algebra) library (shared lib, pthread)
Regards
Andreas Schlager
The information in this e-mail is confidential. The contents may not be disclosed or used by anyone other than the addressee. Access to this e-mail by anyone else is unauthorised.
If you are not the intended recipient, please notify Airbus immediately and delete this e-mail.
Airbus cannot accept any responsibility for the accuracy or completeness of this e-mail as it has been sent over public networks. If you have any concerns over the content of this message or its Accuracy or Integrity, please contact Airbus immediately.
All outgoing e-mails from Airbus are checked using regularly updated virus scanning software but you should take whatever measures you deem to be appropriate to ensure that this message and any attachments are virus free.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-science-maintainers/attachments/20220707/f65c3db2/attachment-0001.htm>
-------------- next part --------------
//===========================================================================================================================================================================================================================================================
// After update to debian 11
// B*inv(B) does NOT yield the unit matrix. This is an error
//===========================================================================================================================================================================================================================================================
B:
4.4143e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 2.4141e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 6.2101e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 4.2099e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.5046e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.3045e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 4.2407e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.2405e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 4.3532e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 2.3530e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.8471e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.6470e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.7360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.5360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03
5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 6.6301e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.9830e+01 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 8.4259e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 1.9941e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.7262e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.0830e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 6.4565e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.4639e+00 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 6.5690e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.3717e+01 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 5.0686e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.4254e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.9576e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.3144e+03
2.4141e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 9.0354e+02 4.7152e+02 8.2781e+02 4.3081e+02 3.8315e+02 9.8281e+02 6.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 4.2099e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 4.7152e+02 1.0831e+03 1.0903e+03 4.6861e+02 4.4006e+02 1.3395e+03 8.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.3045e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 8.2781e+02 1.0903e+03 2.9667e+03 7.8489e+02 7.1454e+02 3.2970e+03 1.8912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 2.2405e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 4.3081e+02 4.6861e+02 7.8489e+02 5.2406e+02 3.8445e+02 9.2970e+02 6.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 2.3530e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 3.8315e+02 4.4006e+02 7.1454e+02 3.8445e+02 8.9743e+02 1.0317e+03 7.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.6470e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 9.8281e+02 1.3395e+03 3.2970e+03 9.2970e+02 1.0317e+03 5.3091e+03 2.8217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.5360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 6.2067e+02 8.6604e+02 1.8912e+03 6.0393e+02 7.1263e+02 2.8217e+03 2.1981e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03
1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.9830e+01 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.4124e+03 4.7152e+02 8.2781e+02 4.3081e+02 3.8315e+02 9.8281e+02 6.2067e+02
3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 1.9941e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 4.7152e+02 1.5920e+03 1.0903e+03 4.6861e+02 4.4006e+02 1.3395e+03 8.6604e+02
7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.0830e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 8.2781e+02 1.0903e+03 3.4755e+03 7.8489e+02 7.1454e+02 3.2970e+03 1.8912e+03
3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.4639e+00 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 4.3081e+02 4.6861e+02 7.8489e+02 7.4564e+02 3.8445e+02 9.2970e+02 6.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.3717e+01 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 3.8315e+02 4.4006e+02 7.1454e+02 3.8445e+02 1.4063e+03 1.0317e+03 7.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.4254e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 9.8281e+02 1.3395e+03 3.2970e+03 9.2970e+02 1.0317e+03 5.8180e+03 2.8217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.3144e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 6.2067e+02 8.6604e+02 1.8912e+03 6.0393e+02 7.1263e+02 2.8217e+03 2.7069e+03
inv(B):
1.3229e+01 6.1365e+00 -5.9469e+00 -1.0600e+01 -4.7321e+00 3.5751e+00 -1.4661e+00 -7.3821e+00 -3.4458e+00 3.3376e+00 5.8779e+00 2.6101e+00 -2.0024e+00 8.5273e-01 -1.2172e-02 1.0486e-02 -2.4158e-03 3.8827e-03 -3.7162e-03 1.8848e-03 -3.1579e-03 -1.3512e-02 -2.7177e-03 -6.1167e-03 -2.1554e-02 1.1714e-03 3.6646e-03 1.6240e-03
6.1365e+00 4.8781e+00 -4.3877e+00 -4.7961e+00 -2.7973e+00 3.3776e+00 -2.4618e+00 -3.4076e+00 -2.7415e+00 2.4588e+00 2.6680e+00 1.5555e+00 -1.8617e+00 1.3795e+00 3.8998e-03 -3.7796e-03 -2.1500e-03 -1.3552e-03 -4.2272e-04 -7.8284e-04 3.5852e-04 -3.3354e-03 -4.4774e-05 3.6002e-04 -5.9148e-03 3.9876e-04 2.4089e-03 1.6210e-04
-5.9469e+00 -4.3877e+00 6.8499e+00 -1.1783e+00 8.9728e+00 -3.4478e+00 -7.6279e-01 3.3257e+00 2.4527e+00 -3.8004e+00 6.4188e-01 -5.0250e+00 1.9375e+00 4.1745e-01 4.7591e-03 2.5658e-03 -3.0040e-03 -5.4352e-03 5.1286e-03 4.0578e-04 2.0698e-03 -1.5411e-02 -2.2870e-04 4.7385e-03 3.5973e-03 2.0548e-03 -1.2520e-03 -3.8662e-06
-1.0600e+01 -4.7961e+00 -1.1783e+00 2.0821e+01 -9.7552e+00 -1.8556e+00 7.2419e+00 5.9078e+00 2.6154e+00 6.4116e-01 -1.1644e+01 5.3921e+00 1.0147e+00 -3.9795e+00 1.7110e-02 1.6246e-02 -9.5099e-03 6.1883e-03 -8.1867e-03 -6.1449e-04 -8.6453e-05 1.9084e-03 -8.7856e-05 -7.5568e-03 -2.9955e-02 7.6663e-04 1.9201e-03 5.0212e-03
-4.7321e+00 -2.7973e+00 8.9728e+00 -9.7552e+00 1.6717e+01 -2.7346e+00 -5.7230e+00 2.6346e+00 1.5279e+00 -5.0386e+00 5.4287e+00 -9.3211e+00 1.5099e+00 3.2328e+00 4.9808e-03 4.5813e-03 -4.6681e-03 1.0125e-02 -9.5433e-03 -2.4357e-04 -1.7293e-03 1.3301e-02 1.8273e-04 -9.3013e-03 -1.9442e-02 -1.6138e-03 7.6740e-04 8.7713e-05
3.5751e+00 3.3776e+00 -3.4478e+00 -1.8556e+00 -2.7346e+00 2.5620e+00 -1.5558e+00 -2.0095e+00 -1.8844e+00 1.9158e+00 1.0433e+00 1.5103e+00 -1.4358e+00 8.7284e-01 2.5269e-02 6.1179e-03 -3.4263e-03 3.4803e-03 1.6352e-04 1.7624e-03 7.5574e-04 -4.6109e-03 1.7212e-05 -1.6043e-03 1.9502e-03 -1.7683e-04 -5.4544e-03 -1.2426e-04
-1.4661e+00 -2.4618e+00 -7.6279e-01 7.2419e+00 -5.7230e+00 -1.5558e+00 5.0047e+00 8.2085e-01 1.4275e+00 4.3465e-01 -4.0025e+00 3.2403e+00 8.8892e-01 -2.8381e+00 -1.8514e-02 -1.7075e-02 6.2865e-03 -1.0390e-03 3.9425e-03 6.8232e-04 8.2321e-05 -1.3381e-03 2.9790e-04 3.2438e-03 2.6199e-02 -2.6035e-03 -2.1052e-03 -6.7402e-03
-7.3821e+00 -3.4076e+00 3.3257e+00 5.9078e+00 2.6346e+00 -2.0095e+00 8.2085e-01 4.1249e+00 1.9169e+00 -1.8614e+00 -3.2673e+00 -1.4522e+00 1.1167e+00 -4.7629e-01 9.1172e-03 -5.2233e-03 3.8083e-03 -3.6260e-04 2.0864e-03 6.5066e-04 -3.4155e-04 2.4501e-03 6.8927e-05 1.4161e-03 1.0482e-02 -6.5311e-04 -2.0424e-03 -9.0674e-04
-3.4458e+00 -2.7415e+00 2.4527e+00 2.6154e+00 1.5279e+00 -1.8844e+00 1.4275e+00 1.9169e+00 1.5364e+00 -1.3699e+00 -1.4844e+00 -8.6705e-01 1.0359e+00 -7.6747e-01 -4.1559e-03 3.5662e-03 3.0048e-04 1.3526e-03 -3.6208e-04 4.3636e-04 -1.0177e-04 1.3434e-03 2.5172e-05 -7.9789e-04 1.7694e-03 -2.2268e-04 -1.3427e-03 -9.2365e-05
3.3376e+00 2.4588e+00 -3.8004e+00 6.4116e-01 -5.0386e+00 1.9158e+00 4.3465e-01 -1.8614e+00 -1.3699e+00 2.1189e+00 -3.5979e-01 2.8012e+00 -1.0800e+00 -2.2915e-01 -2.3024e-03 -1.6080e-03 2.3161e-03 3.6248e-03 -3.4531e-03 -2.2554e-04 -1.7026e-03 1.2364e-02 1.2753e-04 -3.2358e-03 -3.5271e-03 -1.1443e-03 6.9590e-04 7.7225e-08
5.8779e+00 2.6680e+00 6.4188e-01 -1.1644e+01 5.4287e+00 1.0433e+00 -4.0025e+00 -3.2673e+00 -1.4844e+00 -3.5979e-01 6.5035e+00 -2.9970e+00 -5.6676e-01 2.1976e+00 -1.2764e-02 -1.6680e-02 1.1613e-02 -5.1704e-03 7.5732e-03 8.3710e-04 5.4409e-05 -4.4344e-04 2.2009e-04 6.7144e-03 3.5769e-02 -1.9601e-03 -2.5948e-03 -4.3116e-03
2.6101e+00 1.5555e+00 -5.0250e+00 5.3921e+00 -9.3211e+00 1.5103e+00 3.2403e+00 -1.4522e+00 -8.6705e-01 2.8012e+00 -2.9970e+00 5.2023e+00 -8.3464e-01 -1.8018e+00 -2.8990e-03 -1.9630e-03 2.6014e-03 -7.3632e-03 6.7632e-03 1.3469e-04 9.6387e-04 -7.5882e-03 -2.7304e-04 6.6401e-03 9.3054e-03 2.4378e-03 -4.2846e-04 -5.0937e-05
-2.0024e+00 -1.8617e+00 1.9375e+00 1.0147e+00 1.5099e+00 -1.4358e+00 8.8892e-01 1.1167e+00 1.0359e+00 -1.0800e+00 -5.6676e-01 -8.3464e-01 8.0764e-01 -4.8638e-01 -2.0813e-02 -4.3170e-03 1.9150e-03 -2.4374e-03 -5.8877e-04 -1.4782e-03 -4.2434e-04 2.5721e-03 -9.1203e-06 8.9343e-04 -2.6098e-03 9.7822e-05 4.5728e-03 6.7067e-05
8.5273e-01 1.3795e+00 4.1745e-01 -3.9795e+00 3.2328e+00 8.7284e-01 -2.8381e+00 -4.7629e-01 -7.6747e-01 -2.2915e-01 2.1976e+00 -1.8018e+00 -4.8638e-01 1.5930e+00 1.3908e-02 1.1183e-02 -3.4990e-03 5.6162e-04 -2.1970e-03 -3.7956e-04 -5.0280e-05 7.5070e-04 -1.6787e-04 -1.8101e-03 -1.6110e-02 1.4488e-03 1.1710e-03 5.2895e-03
-1.2172e-02 3.8998e-03 4.7591e-03 1.7110e-02 4.9808e-03 2.5269e-02 -1.8514e-02 9.1172e-03 -4.1559e-03 -2.3024e-03 -1.2764e-02 -2.8990e-03 -2.0813e-02 1.3908e-02 2.3675e-03 -2.2598e-04 -4.1555e-04 1.0719e-03 -1.0728e-03 -9.0756e-08 1.7624e-04 -9.2985e-05 3.1469e-07 -1.0718e-03 4.9973e-06 -8.8867e-08 2.8622e-07 -1.1323e-06
1.0486e-02 -3.7796e-03 2.5658e-03 1.6246e-02 4.5813e-03 6.1179e-03 -1.7075e-02 -5.2233e-03 3.5662e-03 -1.6080e-03 -1.6680e-02 -1.9630e-03 -4.3170e-03 1.1183e-02 -2.2598e-04 2.5313e-03 -4.1386e-04 1.0710e-03 -1.0716e-03 -7.3301e-09 1.7613e-04 -9.2538e-05 3.3559e-07 -1.0707e-03 9.1700e-06 -2.9057e-07 1.4066e-08 -1.5326e-06
-2.4158e-03 -2.1500e-03 -3.0040e-03 -9.5099e-03 -4.6681e-03 -3.4263e-03 6.2865e-03 3.8083e-03 3.0048e-04 2.3161e-03 1.1613e-02 2.6014e-03 1.9150e-03 -3.4990e-03 -4.1555e-04 -4.1386e-04 2.7569e-03 -2.7465e-03 7.6294e-07 -1.9081e-06 -2.2263e-06 -8.4830e-07 -1.0308e-06 -3.2837e-06 1.2909e-05 6.0563e-07 -1.3513e-06 -1.5844e-06
3.8827e-03 -1.3552e-03 -5.4352e-03 6.1883e-03 1.0125e-02 3.4803e-03 -1.0390e-03 -3.6260e-04 1.3526e-03 3.6248e-03 -5.1704e-03 -7.3632e-03 -2.4374e-03 5.6162e-04 1.0719e-03 1.0710e-03 -2.7465e-03 2.6449e-02 -2.7548e-03 -2.7469e-03 -2.7362e-03 3.9487e-06 7.2871e-06 1.0514e-05 -1.0665e-04 4.4140e-06 1.0248e-05 1.8061e-05
-3.7162e-03 -4.2272e-04 5.1286e-03 -8.1867e-03 -9.5433e-03 1.6352e-04 3.9425e-03 2.0864e-03 -3.6208e-04 -3.4531e-03 7.5732e-03 6.7632e-03 -5.8877e-04 -2.1970e-03 -1.0728e-03 -1.0716e-03 7.6294e-07 -2.7548e-03 2.7600e-03 -7.7078e-07 -2.1912e-06 -8.8867e-08 -2.9060e-07 6.0486e-07 5.6031e-06 -1.0640e-06 -5.1822e-07 -1.5585e-06
1.8848e-03 -7.8284e-04 4.0578e-04 -6.1449e-04 -2.4357e-04 1.7624e-03 6.8232e-04 6.5066e-04 4.3636e-04 -2.2554e-04 8.3710e-04 1.3469e-04 -1.4782e-03 -3.7956e-04 -9.0756e-08 -7.3301e-09 -1.9081e-06 -2.7469e-03 -7.7078e-07 2.7554e-03 -2.3240e-06 2.8680e-07 1.4528e-08 -1.3513e-06 1.2343e-05 -5.1831e-07 -4.3770e-06 -1.6559e-06
-3.1579e-03 3.5852e-04 2.0698e-03 -8.6453e-05 -1.7293e-03 7.5574e-04 8.2321e-05 -3.4155e-04 -1.0177e-04 -1.7026e-03 5.4409e-05 9.6387e-04 -4.2434e-04 -5.0280e-05 1.7624e-04 1.7613e-04 -2.2263e-06 -2.7362e-03 -2.1912e-06 -2.3240e-06 2.7548e-03 -1.1327e-06 -1.5331e-06 -1.5849e-06 2.2613e-05 -1.5592e-06 -1.6564e-06 -4.8645e-06
-1.3512e-02 -3.3354e-03 -1.5411e-02 1.9084e-03 1.3301e-02 -4.6109e-03 -1.3381e-03 2.4501e-03 1.3434e-03 1.2364e-02 -4.4344e-04 -7.5882e-03 2.5721e-03 7.5070e-04 -9.2985e-05 -9.2538e-05 -8.4830e-07 3.9487e-06 -8.8867e-08 2.8680e-07 -1.1327e-06 1.5390e-03 2.0032e-07 -6.5445e-07 -1.5358e-03 -9.8188e-08 1.7669e-07 -8.6277e-07
-2.7177e-03 -4.4774e-05 -2.2870e-04 -8.7856e-05 1.8273e-04 1.7212e-05 2.9790e-04 6.8927e-05 2.5172e-05 1.2753e-04 2.2009e-04 -2.7304e-04 -9.1203e-06 -1.6787e-04 3.1469e-07 3.3559e-07 -1.0308e-06 7.2871e-06 -2.9060e-07 1.4528e-08 -1.5331e-06 2.0032e-07 1.5381e-03 -7.8820e-07 -1.5328e-03 -2.4602e-07 -2.2844e-08 -1.1562e-06
-6.1167e-03 3.6002e-04 4.7385e-03 -7.5568e-03 -9.3013e-03 -1.6043e-03 3.2438e-03 1.4161e-03 -7.9789e-04 -3.2358e-03 6.7144e-03 6.6401e-03 8.9343e-04 -1.8101e-03 -1.0718e-03 -1.0707e-03 -3.2837e-06 1.0514e-05 6.0486e-07 -1.3513e-06 -1.5849e-06 -6.5445e-07 -7.8820e-07 1.5375e-03 -1.5300e-03 4.1020e-07 -1.0230e-06 -1.1941e-06
-2.1554e-02 -5.9148e-03 3.5973e-03 -2.9955e-02 -1.9442e-02 1.9502e-03 2.6199e-02 1.0482e-02 1.7694e-03 -3.5271e-03 3.5769e-02 9.3054e-03 -2.6098e-03 -1.6110e-02 4.9973e-06 9.1700e-06 1.2909e-05 -1.0665e-04 5.6031e-06 1.2343e-05 2.2613e-05 -1.5358e-03 -1.5328e-03 -1.5300e-03 1.9133e-02 -1.5354e-03 -1.5304e-03 -1.5229e-03
1.1714e-03 3.9876e-04 2.0548e-03 7.6663e-04 -1.6138e-03 -1.7683e-04 -2.6035e-03 -6.5311e-04 -2.2268e-04 -1.1443e-03 -1.9601e-03 2.4378e-03 9.7822e-05 1.4488e-03 -8.8867e-08 -2.9057e-07 6.0563e-07 4.4140e-06 -1.0640e-06 -5.1831e-07 -1.5592e-06 -9.8188e-08 -2.4602e-07 4.1020e-07 -1.5354e-03 1.5391e-03 -4.1282e-07 -1.1751e-06
3.6646e-03 2.4089e-03 -1.2520e-03 1.9201e-03 7.6740e-04 -5.4544e-03 -2.1052e-03 -2.0424e-03 -1.3427e-03 6.9590e-04 -2.5948e-03 -4.2846e-04 4.5728e-03 1.1710e-03 2.8622e-07 1.4066e-08 -1.3513e-06 1.0248e-05 -5.1822e-07 -4.3770e-06 -1.6564e-06 1.7669e-07 -2.2844e-08 -1.0230e-06 -1.5304e-03 -4.1282e-07 1.5367e-03 -1.2465e-06
1.6240e-03 1.6210e-04 -3.8662e-06 5.0212e-03 8.7713e-05 -1.2426e-04 -6.7402e-03 -9.0674e-04 -9.2365e-05 7.7225e-08 -4.3116e-03 -5.0937e-05 6.7067e-05 5.2895e-03 -1.1323e-06 -1.5326e-06 -1.5844e-06 1.8061e-05 -1.5585e-06 -1.6559e-06 -4.8645e-06 -8.6277e-07 -1.1562e-06 -1.1941e-06 -1.5229e-03 -1.1751e-06 -1.2465e-06 1.5363e-03
B*inv(B):
-8.6239e+00 4.0418e+01 3.7016e+01 -5.5621e+01 -5.1818e+01 -2.3570e+00 6.5160e+01 7.0817e+00 -2.6662e+01 -1.2374e+01 1.9205e+01 1.2250e+01 5.8784e+00 3.1939e+00 4.1642e+00 1.5815e+00 7.9734e-01 4.7704e+00 -8.7043e-01 3.1439e+00 -4.9325e-01 -7.9659e+00 -1.0875e+00 -1.9912e+00 5.7766e-15 2.8004e-16 -1.7911e-15 3.0686e-16
-1.2848e+01 5.2823e+01 4.8160e+01 -7.5292e+01 -6.8726e+01 2.6273e+00 8.6947e+01 7.8110e+00 -2.7286e+01 -1.2997e+01 2.0633e+01 1.6962e+01 9.6514e+00 4.3904e+00 4.5580e+00 1.1475e+00 4.2932e-01 3.3061e+00 -7.5204e-01 3.9645e+00 6.1312e-01 -8.8211e+00 -1.0073e+00 -1.4869e+00 4.9045e-15 2.8621e-16 3.2859e-17 -4.2695e-16
1.1841e+01 1.4646e+02 1.2645e+02 -1.2689e+02 -1.5169e+02 4.3478e+00 1.5269e+02 -1.2555e+01 -5.3975e+01 -3.3236e+01 4.1101e+01 4.1586e+01 2.8565e+01 8.0148e+00 1.2227e+01 1.7378e+00 1.7179e+00 1.4641e+00 -4.3988e-01 1.0379e+01 2.4606e+00 -1.9825e+01 -1.4268e+00 -1.7454e+00 2.3064e-14 -1.2842e-15 -2.0370e-15 -4.3584e-17
-1.5795e+01 3.5490e+01 2.7663e+01 -5.1653e+01 -4.4687e+01 3.3549e+00 6.2308e+01 1.1298e+01 -3.0438e+01 -9.8948e+00 1.0486e+01 1.0257e+01 4.2712e+00 4.1290e+00 3.0513e+00 -8.3572e-02 1.6609e+00 1.7227e+00 -1.8818e-01 3.0170e+00 1.4393e-01 -6.0358e+00 -4.8833e-01 -1.0574e+00 7.1701e-15 -2.9828e-17 -1.5541e-15 -1.2715e-16
-1.4766e+01 3.8065e+01 2.2980e+01 -5.0518e+01 -3.2228e+01 -4.0516e-01 6.4356e+01 6.5565e+00 -2.3016e+01 -1.2103e+01 1.9747e+01 1.1738e+01 8.1216e+00 4.7298e+00 3.6147e+00 9.5523e-01 1.9767e-01 3.0908e+00 -6.7761e-01 2.9683e+00 3.9961e-01 -5.6991e+00 -6.6846e-01 -1.3730e+00 3.9242e-15 -5.9606e-17 2.2345e-16 3.1530e-17
1.8906e+01 2.1507e+02 1.6910e+02 -1.5515e+02 -1.8348e+02 6.1315e+00 1.9848e+02 -2.9806e+01 -6.4654e+01 -5.1866e+01 6.3957e+01 6.8659e+01 4.7613e+01 7.3706e+00 1.5143e+01 3.5595e-01 2.5997e+00 -2.9864e+00 4.9046e-01 1.5046e+01 4.8931e+00 -2.6833e+01 -1.6400e+00 -1.1879e+00 3.1730e-14 -1.5636e-15 -6.3486e-15 -1.0312e-15
1.4824e+01 1.2124e+02 9.0283e+01 -7.9270e+01 -9.1183e+01 4.4628e+00 1.1031e+02 -1.7649e+01 -3.0745e+01 -2.6479e+01 2.9577e+01 4.0475e+01 2.9510e+01 5.5066e+00 9.8163e+00 5.0701e-01 1.3246e+00 -7.0271e-01 1.5603e-01 8.7891e+00 2.5627e+00 -1.5109e+01 -9.2719e-01 -1.1424e+00 1.0587e-14 -4.3039e-16 -1.8259e-15 -2.4099e-15
-6.8748e+00 4.2440e+01 4.0370e+01 -5.7384e+01 -5.4557e+01 -6.0706e+00 6.7333e+01 7.4278e+00 -2.7074e+01 -1.4092e+01 2.1584e+01 1.3947e+01 7.2410e+00 1.8470e+00 4.3214e+00 1.7010e+00 1.1807e+00 4.9411e+00 -7.8320e-01 3.3859e+00 -8.0610e-01 -8.9004e+00 -1.2933e+00 -2.1897e+00 5.7766e-15 5.0208e-16 -2.0132e-15 1.1950e-15
-1.6898e+01 4.9926e+01 4.8988e+01 -8.4812e+01 -7.3566e+01 1.8812e+00 9.6130e+01 1.0709e+01 -2.6390e+01 -1.2811e+01 2.2316e+01 1.7156e+01 9.7358e+00 3.7328e+00 4.2938e+00 1.1543e+00 3.3454e-01 3.3051e+00 -7.5203e-01 3.9645e+00 6.1332e-01 -8.8872e+00 -1.2261e+00 -1.4868e+00 5.3486e-15 5.0825e-16 -1.2994e-15 2.3918e-16
1.6275e+01 1.4854e+02 1.2838e+02 -1.2725e+02 -1.5548e+02 4.0721e+00 1.5309e+02 -1.4221e+01 -5.4341e+01 -3.2795e+01 3.8925e+01 4.0590e+01 2.8233e+01 9.0442e+00 1.2462e+01 1.8976e+00 1.5911e+00 1.7653e+00 -5.2472e-01 1.0380e+01 2.3828e+00 -1.9288e+01 -1.4266e+00 -2.0486e+00 2.1287e-14 -9.5113e-16 -2.6066e-16 8.4459e-16
-1.7418e+01 3.5671e+01 2.5459e+01 -5.3333e+01 -4.4018e+01 4.0744e+00 6.3539e+01 1.3507e+01 -3.4574e+01 -1.0021e+01 7.6154e+00 1.0749e+01 4.6722e+00 3.4426e+00 2.5065e+00 -1.2545e+00 2.7750e+00 -6.0507e-01 4.5777e-01 3.3035e+00 3.5865e-01 -5.7294e+00 -2.4652e-01 -4.8462e-01 4.9496e-15 -2.9828e-17 -1.3321e-15 9.4890e-17
-1.8567e+01 3.7510e+01 1.8954e+01 -5.6476e+01 -3.1637e+01 -2.3874e+00 7.1552e+01 8.8092e+00 -2.5140e+01 -1.2713e+01 2.3558e+01 1.1790e+01 9.1521e+00 4.7207e+00 3.6823e+00 1.1243e+00 1.9743e-01 3.0642e+00 -6.9100e-01 2.9683e+00 3.9999e-01 -5.7238e+00 -6.9275e-01 -1.1661e+00 3.4801e-15 -1.1598e-16 6.6754e-16 2.5357e-16
1.6750e+01 2.1766e+02 1.7160e+02 -1.5811e+02 -1.8552e+02 5.5885e+00 2.0186e+02 -2.9980e+01 -6.6502e+01 -5.3740e+01 6.6328e+01 6.9683e+01 4.8172e+01 7.2244e+00 1.4184e+01 2.2670e-01 2.6009e+00 -2.8408e+00 4.1978e-01 1.4757e+01 4.8930e+00 -2.6832e+01 -1.6400e+00 -1.1877e+00 2.6401e-14 -2.8728e-16 -8.5413e-15 1.6334e-15
1.9920e+01 1.2227e+02 8.9031e+01 -7.1960e+01 -8.4881e+01 5.2467e+00 1.0428e+02 -2.0173e+01 -2.6731e+01 -2.4483e+01 2.5500e+01 4.0964e+01 3.0815e+01 6.3771e+00 1.0313e+01 7.3045e-01 1.3257e+00 -4.9069e-01 1.5644e-01 8.7896e+00 2.3442e+00 -1.5108e+01 -9.2716e-01 -1.1424e+00 8.7060e-15 -5.3369e-16 2.0796e-16 -2.3566e-15
-1.9292e+01 3.5691e+01 3.1529e+01 -4.1802e+01 -4.2283e+01 1.9209e+01 5.0561e+01 1.1226e+01 -2.6678e+01 -9.5488e+00 6.5383e+00 6.4851e+00 -7.7071e+00 1.3493e+01 4.8165e+00 1.4040e+00 -4.0144e-01 6.4885e+00 -1.6927e+00 2.5339e+00 3.9596e-01 -5.6606e+00 -5.6788e-01 -2.0922e+00 8.4411e-15 -5.3027e-17 -1.7911e-15 -5.8132e-16
-4.9628e-02 5.7785e+01 4.8143e+01 -4.2757e+01 -5.4203e+01 9.7582e+00 5.5328e+01 -3.2176e-02 -2.7844e+01 -1.4661e+01 8.5966e+00 1.5411e+01 5.2692e+00 1.1896e+01 5.3316e+00 2.2542e+00 4.3635e-01 5.4575e+00 -1.3543e+00 3.9650e+00 7.1278e-01 -8.7060e+00 -4.5500e-01 -2.0886e+00 5.7927e-15 2.8621e-16 -1.2994e-15 2.3918e-16
3.3661e-01 1.4082e+02 1.1990e+02 -1.2708e+02 -1.4509e+02 7.8334e+00 1.5152e+02 -6.5986e+00 -5.3279e+01 -3.3162e+01 4.9224e+01 4.5766e+01 2.6933e+01 6.0304e+00 1.1673e+01 1.3913e+00 2.9545e+00 1.4719e+00 -4.4084e-01 1.0378e+01 2.6918e+00 -2.1196e+01 -1.4271e+00 -1.1969e+00 1.4182e-14 -9.5113e-16 4.1802e-15 4.0050e-16
-1.1169e+01 3.5087e+01 3.3987e+01 -4.7383e+01 -4.5050e+01 5.1386e+00 5.5922e+01 6.1126e+00 -2.0171e+01 -9.9745e+00 1.7434e+01 8.3030e+00 3.4449e-01 8.2081e+00 4.6132e+00 3.0734e+00 -1.2288e+00 9.3519e+00 -2.0300e+00 2.2951e+00 -3.6103e-01 -6.8269e+00 -1.0975e+00 -2.7145e+00 5.3937e-15 -2.9828e-17 -4.4391e-16 -9.0431e-16
-6.5772e+00 3.9169e+01 3.4739e+01 -3.5373e+01 -3.8497e+01 8.3369e+00 4.5654e+01 2.3434e+00 -1.7842e+01 -1.1645e+01 1.0886e+01 1.4815e+01 2.3452e+00 5.5862e+00 3.2452e+00 3.4712e-01 1.1678e-01 3.9202e+00 1.4250e-01 2.9685e+00 4.3443e-01 -5.6564e+00 -6.0682e-01 -2.1074e+00 4.8124e-15 1.4899e-16 8.8958e-16 -7.9492e-17
2.3734e+01 2.0710e+02 1.6435e+02 -1.4542e+02 -1.7749e+02 1.1001e+01 1.8796e+02 -2.7728e+01 -5.9332e+01 -4.7969e+01 5.6608e+01 6.5625e+01 4.3161e+01 9.3671e+00 1.7742e+01 8.8710e-01 2.5144e+00 -2.5881e+00 4.5414e-01 1.6771e+01 4.9286e+00 -2.6851e+01 -1.6397e+00 -1.4023e+00 2.5512e-14 -1.7965e-15 -7.2043e-15 -1.0312e-15
5.4758e-01 1.1876e+02 9.4958e+01 -9.6766e+01 -1.0643e+02 6.4207e+00 1.2217e+02 -1.0567e+01 -4.1023e+01 -3.2532e+01 3.9057e+01 3.8899e+01 2.3112e+01 5.6994e+00 8.8204e+00 2.1699e-01 1.2380e+00 -4.6720e-01 -6.0557e-02 8.7872e+00 4.1475e+00 -1.5130e+01 -9.2734e-01 -1.3571e+00 1.4446e-14 -7.9116e-16 -1.6006e-15 -1.8110e-15
-2.9048e+01 2.9723e+01 1.5332e+01 -4.9575e+01 -3.4513e+01 6.2518e+00 5.8529e+01 9.5108e+00 -2.3815e+01 1.7888e+00 1.3101e+01 2.0089e+00 3.5735e+00 7.3684e+00 3.4358e+00 9.8364e-01 -5.5192e-01 4.1661e+00 -1.1769e+00 2.2922e+00 6.0893e-01 -3.6741e+00 -3.6221e-01 -1.2916e+00 1.7555e-15 -2.4294e-16 -9.9580e-16 -2.5086e-16
-5.0010e+00 6.1676e+01 4.5092e+01 -4.3923e+01 -5.2758e+01 3.7688e+00 5.6868e+01 -7.6626e-01 -2.9954e+01 -1.3406e+01 1.7796e+01 1.6876e+01 1.0180e+01 5.4514e+00 5.3715e+00 1.0089e+00 7.6317e-01 3.3095e+00 -7.5198e-01 3.9652e+00 6.1241e-01 -8.5877e+00 7.6234e-01 -1.4873e+00 5.7684e-15 -2.3666e-16 -2.8203e-16 2.3658e-16
-1.0083e+01 1.3983e+02 1.2044e+02 -1.3302e+02 -1.4811e+02 4.0893e+00 1.5525e+02 -5.4709e+00 -5.3713e+01 -3.5466e+01 5.6313e+01 5.1595e+01 3.0415e+01 2.1738e+00 1.0586e+01 3.6364e-01 2.1638e+00 4.0650e-01 -1.4178e-01 1.0377e+01 2.7338e+00 -2.1714e+01 -1.4276e+00 3.1970e-01 1.3047e-14 -5.8560e-16 2.5331e-15 -2.6823e-16
-1.2438e+01 3.4883e+01 3.4928e+01 -5.1270e+01 -4.7366e+01 6.4764e-02 5.8724e+01 3.8267e+00 -1.6154e+01 -9.0049e+00 2.5463e+01 8.8229e+00 3.2845e+00 6.1232e+00 4.8512e+00 3.9215e+00 -2.2596e+00 9.9132e+00 -2.4611e+00 2.0088e+00 -6.1112e-01 -7.1148e+00 -1.3395e+00 -3.0728e+00 1.0000e+00 -4.9718e-16 -9.2669e-17 -7.9589e-16
-3.7448e+00 4.0051e+01 3.6661e+01 -3.0560e+01 -3.9158e+01 5.8181e+00 3.9779e+01 -1.0634e+00 -1.5815e+01 -9.5064e+00 7.6828e+00 1.5371e+01 4.9170e+00 4.3347e+00 3.2596e+00 2.4351e-01 2.0007e-01 3.1810e+00 -6.2951e-01 2.9686e+00 3.9849e-01 -5.6136e+00 -5.8277e-01 -2.0998e+00 6.3424e-15 1.0000e+00 -4.8816e-16 -8.2095e-17
2.4123e+01 2.0596e+02 1.5983e+02 -1.4578e+02 -1.7663e+02 2.7792e+00 1.8732e+02 -2.8874e+01 -5.8423e+01 -4.4826e+01 5.6973e+01 6.5353e+01 4.9575e+01 7.4572e+00 1.8393e+01 6.9366e-01 2.5961e+00 -3.4984e+00 7.3915e-01 1.6059e+01 4.8930e+00 -2.6832e+01 -1.6400e+00 -1.1884e+00 3.8145e-14 -3.8327e-15 1.0000e+00 -3.0322e-15
-5.4992e+00 1.1766e+02 9.4196e+01 -1.0606e+02 -1.1371e+02 9.4771e-01 1.2780e+02 -8.4420e+00 -4.5166e+01 -3.3070e+01 4.5306e+01 3.9047e+01 2.5334e+01 5.5287e+00 7.9537e+00 -3.9257e-01 1.3192e+00 -1.4425e+00 1.5304e-01 8.7866e+00 3.3298e+00 -1.5113e+01 -9.2775e-01 -1.1433e+00 1.1161e-14 -4.5653e-16 4.8958e-15 1.0000e+00
//===========================================================================================================================================================================================================================================================
// Reference: Before update
// B*inv(B) yields the unit matrix like it should
//===========================================================================================================================================================================================================================================================
B:
4.4143e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 2.4141e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 6.2101e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 4.2099e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.5046e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.3045e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 4.2407e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.2405e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 4.3532e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 2.3530e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.8471e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.6470e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.7360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.5360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03
5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 6.6301e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.9830e+01 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 8.4259e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 1.9941e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.7262e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.0830e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 6.4565e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.4639e+00 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 6.5690e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.3717e+01 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 5.0686e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.4254e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.9576e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.3144e+03
2.4141e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 9.0354e+02 4.7152e+02 8.2781e+02 4.3081e+02 3.8315e+02 9.8281e+02 6.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02
3.7152e+02 4.2099e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 4.7152e+02 1.0831e+03 1.0903e+03 4.6861e+02 4.4006e+02 1.3395e+03 8.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02
7.2781e+02 9.9030e+02 2.3045e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 8.2781e+02 1.0903e+03 2.9667e+03 7.8489e+02 7.1454e+02 3.2970e+03 1.8912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03
3.3081e+02 3.6861e+02 6.8489e+02 2.2405e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 4.3081e+02 4.6861e+02 7.8489e+02 5.2406e+02 3.8445e+02 9.2970e+02 6.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 2.3530e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 3.8315e+02 4.4006e+02 7.1454e+02 3.8445e+02 8.9743e+02 1.0317e+03 7.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.6470e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 9.8281e+02 1.3395e+03 3.2970e+03 9.2970e+02 1.0317e+03 5.3091e+03 2.8217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.5360e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 6.2067e+02 8.6604e+02 1.8912e+03 6.0393e+02 7.1263e+02 2.8217e+03 2.1981e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03
1.6208e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.9830e+01 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 5.2074e+02 3.7152e+02 7.2781e+02 3.3081e+02 2.8315e+02 8.8281e+02 5.2067e+02 1.4124e+03 4.7152e+02 8.2781e+02 4.3081e+02 3.8315e+02 9.8281e+02 6.2067e+02
3.7152e+02 3.4166e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 1.9941e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 3.7152e+02 7.0031e+02 9.9030e+02 3.6861e+02 3.4006e+02 1.2395e+03 7.6604e+02 4.7152e+02 1.5920e+03 1.0903e+03 4.6861e+02 4.4006e+02 1.3395e+03 8.6604e+02
7.2781e+02 9.9030e+02 2.2252e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.0830e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 7.2781e+02 9.9030e+02 2.5839e+03 6.8489e+02 6.1454e+02 3.1970e+03 1.7912e+03 8.2781e+02 1.0903e+03 3.4755e+03 7.8489e+02 7.1454e+02 3.2970e+03 1.8912e+03
3.3081e+02 3.6861e+02 6.8489e+02 1.4472e+02 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 2.4639e+00 2.8445e+02 8.2970e+02 5.0393e+02 3.3081e+02 3.6861e+02 6.8489e+02 5.0337e+02 2.8445e+02 8.2970e+02 5.0393e+02 4.3081e+02 4.6861e+02 7.8489e+02 7.4564e+02 3.8445e+02 9.2970e+02 6.0393e+02
2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.5597e+02 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 1.3717e+01 9.3167e+02 6.1263e+02 2.8315e+02 3.4006e+02 6.1454e+02 2.8445e+02 5.1463e+02 9.3167e+02 6.1263e+02 3.8315e+02 4.4006e+02 7.1454e+02 3.8445e+02 1.4063e+03 1.0317e+03 7.1263e+02
8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.5677e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.4254e+03 2.7217e+03 8.8281e+02 1.2395e+03 3.1970e+03 8.2970e+02 9.3167e+02 4.9263e+03 2.7217e+03 9.8281e+02 1.3395e+03 3.2970e+03 9.2970e+02 1.0317e+03 5.8180e+03 2.8217e+03
5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.4567e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.3144e+03 5.2067e+02 7.6604e+02 1.7912e+03 5.0393e+02 6.1263e+02 2.7217e+03 1.8153e+03 6.2067e+02 8.6604e+02 1.8912e+03 6.0393e+02 7.1263e+02 2.8217e+03 2.7069e+03
inv(B):
1.3303e+01 6.1467e+00 -5.9733e+00 -1.0521e+01 -4.6845e+00 3.5674e+00 -1.5387e+00 -7.4260e+00 -3.4200e+00 3.3285e+00 5.8313e+00 2.6133e+00 -1.9812e+00 8.6037e-01 -3.1381e-03 4.7923e-03 -7.0594e-05 -2.0879e-02 1.5986e-03 5.0027e-03 2.2158e-03 -2.1363e-03 3.5115e-03 -5.0233e-05 -2.1554e-02 1.1714e-03 3.6646e-03 1.6240e-03
6.1467e+00 4.9201e+00 -4.4002e+00 -4.7735e+00 -2.7847e+00 3.3398e+00 -2.4843e+00 -3.4224e+00 -2.7582e+00 2.4509e+00 2.6527e+00 1.5530e+00 -1.8566e+00 1.3848e+00 3.0187e-03 -6.8869e-03 -8.6937e-04 -4.0744e-03 5.4422e-04 3.2890e-03 2.2060e-04 2.2121e-03 -4.8836e-03 -6.3596e-04 -5.9148e-03 3.9876e-04 2.4089e-03 1.6210e-04
-5.9733e+00 -4.4002e+00 6.8725e+00 -1.1838e+00 8.9635e+00 -3.4444e+00 -7.5598e-01 3.3279e+00 2.4513e+00 -3.8443e+00 6.6686e-01 -4.9928e+00 1.9169e+00 4.2129e-01 -1.3297e-03 -1.0325e-03 -5.9826e-03 6.0715e-03 2.8040e-03 -1.7103e-03 -5.4490e-06 -9.7444e-04 -7.5645e-04 -4.2225e-03 3.5973e-03 2.0548e-03 -1.2520e-03 -3.8662e-06
-1.0521e+01 -4.7735e+00 -1.1838e+00 2.0942e+01 -9.6832e+00 -1.8630e+00 7.1199e+00 5.8656e+00 2.6637e+00 6.6384e-01 -1.1729e+01 5.3991e+00 1.0415e+00 -3.9589e+00 1.0268e-03 2.4700e-03 3.3573e-03 -3.2332e-02 1.2694e-03 2.8424e-03 7.0772e-03 5.8885e-04 1.6467e-03 2.2970e-03 -2.9955e-02 7.6663e-04 1.9201e-03 5.0212e-03
-4.6845e+00 -2.7847e+00 8.9635e+00 -9.6832e+00 1.6767e+01 -2.7328e+00 -5.7832e+00 2.6130e+00 1.5545e+00 -4.9857e+00 5.3678e+00 -9.3563e+00 1.5240e+00 3.2231e+00 1.3635e-03 1.8023e-03 7.1092e-03 -1.8659e-02 -2.4251e-03 1.0459e-03 1.1825e-04 9.9907e-04 1.3208e-03 5.2072e-03 -1.9442e-02 -1.6138e-03 7.6740e-04 8.7713e-05
3.5674e+00 3.3398e+00 -3.4444e+00 -1.8630e+00 -2.7328e+00 2.6030e+00 -1.5542e+00 -1.9866e+00 -1.8596e+00 1.9178e+00 1.0429e+00 1.5230e+00 -1.4671e+00 8.6592e-01 1.5034e-03 1.5458e-03 -1.1950e-03 4.3606e-03 -2.4126e-04 -7.6649e-03 -1.6995e-04 1.1017e-03 1.1327e-03 -8.7478e-04 1.9502e-03 -1.7683e-04 -5.4544e-03 -1.2426e-04
-1.5387e+00 -2.4843e+00 -7.5598e-01 7.1199e+00 -5.7832e+00 -1.5542e+00 5.1248e+00 8.5464e-01 1.3812e+00 4.1826e-01 -3.9256e+00 3.2191e+00 8.6225e-01 -2.8750e+00 -2.3815e-03 -2.7223e-03 -2.3295e-03 3.0359e-02 -3.5527e-03 -2.8737e-03 -9.4195e-03 -1.7452e-03 -1.9948e-03 -1.7068e-03 2.6199e-02 -2.6035e-03 -2.1052e-03 -6.7402e-03
-7.4260e+00 -3.4224e+00 3.3279e+00 5.8656e+00 2.6130e+00 -1.9866e+00 8.5464e-01 4.1510e+00 1.9042e+00 -1.8544e+00 -3.2535e+00 -1.4577e+00 1.1033e+00 -4.7788e-01 3.2869e-03 -2.6708e-03 3.7822e-05 1.0106e-02 -8.9130e-04 -2.7881e-03 -1.2372e-03 2.7290e-03 -1.9570e-03 2.6873e-05 1.0482e-02 -6.5311e-04 -2.0424e-03 -9.0674e-04
-3.4200e+00 -2.7582e+00 2.4513e+00 2.6637e+00 1.5545e+00 -1.8596e+00 1.3812e+00 1.9042e+00 1.5517e+00 -1.3654e+00 -1.4826e+00 -8.6691e-01 1.0337e+00 -7.6990e-01 -1.6821e-03 5.3741e-03 4.8273e-04 7.4452e-04 -3.0389e-04 -1.8333e-03 -1.2570e-04 -1.2326e-03 4.2586e-03 3.5312e-04 1.7694e-03 -2.2268e-04 -1.3427e-03 -9.2365e-05
3.3285e+00 2.4509e+00 -3.8443e+00 6.6384e-01 -4.9857e+00 1.9178e+00 4.1826e-01 -1.8544e+00 -1.3654e+00 2.1558e+00 -3.7630e-01 2.7771e+00 -1.0673e+00 -2.3309e-01 7.3944e-04 5.7348e-04 4.8675e-03 -4.9040e-03 -1.5615e-03 9.5063e-04 2.1073e-07 5.4188e-04 4.2016e-04 3.8882e-03 -3.5271e-03 -1.1443e-03 6.9590e-04 7.7225e-08
5.8313e+00 2.6527e+00 6.6686e-01 -1.1729e+01 5.3678e+00 1.0429e+00 -3.9256e+00 -3.2535e+00 -1.4826e+00 -3.7630e-01 6.5988e+00 -2.9953e+00 -5.8537e-01 2.1804e+00 -2.1037e-03 -2.9015e-03 -3.3899e-03 3.2663e-02 -2.2380e-03 -3.1037e-03 -5.4480e-03 -1.8617e-03 -2.4464e-03 -2.8044e-03 3.5769e-02 -1.9601e-03 -2.5948e-03 -4.3116e-03
2.6133e+00 1.5530e+00 -4.9928e+00 5.3991e+00 -9.3563e+00 1.5230e+00 3.2191e+00 -1.4577e+00 -8.6691e-01 2.7771e+00 -2.9953e+00 5.2265e+00 -8.4933e-01 -1.7940e+00 -7.6021e-04 -1.0052e-03 -3.9618e-03 8.8691e-03 2.8896e-03 -5.8395e-04 -6.8695e-05 -5.5705e-04 -7.3665e-04 -2.9019e-03 9.3054e-03 2.4378e-03 -4.2846e-04 -5.0937e-05
-1.9812e+00 -1.8566e+00 1.9169e+00 1.0415e+00 1.5240e+00 -1.4671e+00 8.6225e-01 1.1033e+00 1.0337e+00 -1.0673e+00 -5.8537e-01 -8.4933e-01 8.3223e-01 -4.8040e-01 -8.3719e-04 -8.6131e-04 6.6338e-04 -3.9507e-03 1.3347e-04 5.8024e-03 9.1742e-05 -6.1350e-04 -6.3110e-04 4.8563e-04 -2.6098e-03 9.7822e-05 4.5728e-03 6.7067e-05
8.6037e-01 1.3848e+00 4.2129e-01 -3.9589e+00 3.2231e+00 8.6592e-01 -2.8750e+00 -4.7788e-01 -7.6990e-01 -2.3309e-01 2.1804e+00 -1.7940e+00 -4.8040e-01 1.6183e+00 1.3252e-03 1.5144e-03 1.2954e-03 -1.8426e-02 1.9771e-03 1.5985e-03 6.7804e-03 9.7107e-04 1.1097e-03 9.4912e-04 -1.6110e-02 1.4488e-03 1.1710e-03 5.2895e-03
-3.1381e-03 3.0187e-03 -1.3297e-03 1.0268e-03 1.3635e-03 1.5034e-03 -2.3815e-03 3.2869e-03 -1.6821e-03 7.3944e-04 -2.1037e-03 -7.6021e-04 -8.3719e-04 1.3252e-03 2.7598e-03 3.7127e-07 -1.2211e-06 -2.7555e-03 -1.8474e-07 3.2789e-07 -1.6092e-06 -1.1625e-06 3.1847e-07 -8.4796e-07 4.9973e-06 -8.8867e-08 2.8622e-07 -1.1323e-06
4.7923e-03 -6.8869e-03 -1.0325e-03 2.4700e-03 1.8023e-03 1.5458e-03 -2.7223e-03 -2.6708e-03 5.3741e-03 5.7348e-04 -2.9015e-03 -1.0052e-03 -8.6131e-04 1.5144e-03 3.7127e-07 2.7582e-03 -1.4701e-06 -2.7509e-03 -4.6001e-07 -4.3488e-08 -2.1556e-06 3.1860e-07 -2.3912e-06 -1.0305e-06 9.1700e-06 -2.9057e-07 1.4066e-08 -1.5326e-06
-7.0594e-05 -8.6937e-04 -5.9826e-03 3.3573e-03 7.1092e-03 -1.1950e-03 -2.3295e-03 3.7822e-05 4.8273e-04 4.8675e-03 -3.3899e-03 -3.9618e-03 6.6338e-04 1.2954e-03 -1.2211e-06 -1.4701e-06 2.7569e-03 -2.7465e-03 7.6294e-07 -1.9081e-06 -2.2263e-06 -8.4830e-07 -1.0308e-06 -3.2837e-06 1.2909e-05 6.0563e-07 -1.3513e-06 -1.5844e-06
-2.0879e-02 -4.0744e-03 6.0715e-03 -3.2332e-02 -1.8659e-02 4.3606e-03 3.0359e-02 1.0106e-02 7.4452e-04 -4.9040e-03 3.2663e-02 8.8691e-03 -3.9507e-03 -1.8426e-02 -2.7555e-03 -2.7509e-03 -2.7465e-03 2.6449e-02 -2.7548e-03 -2.7469e-03 -2.7362e-03 3.9487e-06 7.2871e-06 1.0514e-05 -1.0665e-04 4.4140e-06 1.0248e-05 1.8061e-05
1.5986e-03 5.4422e-04 2.8040e-03 1.2694e-03 -2.4251e-03 -2.4126e-04 -3.5527e-03 -8.9130e-04 -3.0389e-04 -1.5615e-03 -2.2380e-03 2.8896e-03 1.3347e-04 1.9771e-03 -1.8474e-07 -4.6001e-07 7.6294e-07 -2.7548e-03 2.7600e-03 -7.7078e-07 -2.1912e-06 -8.8867e-08 -2.9060e-07 6.0486e-07 5.6031e-06 -1.0640e-06 -5.1822e-07 -1.5585e-06
5.0027e-03 3.2890e-03 -1.7103e-03 2.8424e-03 1.0459e-03 -7.6649e-03 -2.8737e-03 -2.7881e-03 -1.8333e-03 9.5063e-04 -3.1037e-03 -5.8395e-04 5.8024e-03 1.5985e-03 3.2789e-07 -4.3488e-08 -1.9081e-06 -2.7469e-03 -7.7078e-07 2.7554e-03 -2.3240e-06 2.8680e-07 1.4528e-08 -1.3513e-06 1.2343e-05 -5.1831e-07 -4.3770e-06 -1.6559e-06
2.2158e-03 2.2060e-04 -5.4490e-06 7.0772e-03 1.1825e-04 -1.6995e-04 -9.4195e-03 -1.2372e-03 -1.2570e-04 2.1073e-07 -5.4480e-03 -6.8695e-05 9.1742e-05 6.7804e-03 -1.6092e-06 -2.1556e-06 -2.2263e-06 -2.7362e-03 -2.1912e-06 -2.3240e-06 2.7548e-03 -1.1327e-06 -1.5331e-06 -1.5849e-06 2.2613e-05 -1.5592e-06 -1.6564e-06 -4.8645e-06
-2.1363e-03 2.2121e-03 -9.7444e-04 5.8885e-04 9.9907e-04 1.1017e-03 -1.7452e-03 2.7290e-03 -1.2326e-03 5.4188e-04 -1.8617e-03 -5.5705e-04 -6.1350e-04 9.7107e-04 -1.1625e-06 3.1860e-07 -8.4830e-07 3.9487e-06 -8.8867e-08 2.8680e-07 -1.1327e-06 1.5390e-03 2.0032e-07 -6.5445e-07 -1.5358e-03 -9.8188e-08 1.7669e-07 -8.6277e-07
3.5115e-03 -4.8836e-03 -7.5645e-04 1.6467e-03 1.3208e-03 1.1327e-03 -1.9948e-03 -1.9570e-03 4.2586e-03 4.2016e-04 -2.4464e-03 -7.3665e-04 -6.3110e-04 1.1097e-03 3.1847e-07 -2.3912e-06 -1.0308e-06 7.2871e-06 -2.9060e-07 1.4528e-08 -1.5331e-06 2.0032e-07 1.5381e-03 -7.8820e-07 -1.5328e-03 -2.4602e-07 -2.2844e-08 -1.1562e-06
-5.0233e-05 -6.3596e-04 -4.2225e-03 2.2970e-03 5.2072e-03 -8.7478e-04 -1.7068e-03 2.6873e-05 3.5312e-04 3.8882e-03 -2.8044e-03 -2.9019e-03 4.8563e-04 9.4912e-04 -8.4796e-07 -1.0305e-06 -3.2837e-06 1.0514e-05 6.0486e-07 -1.3513e-06 -1.5849e-06 -6.5445e-07 -7.8820e-07 1.5375e-03 -1.5300e-03 4.1020e-07 -1.0230e-06 -1.1941e-06
-2.1554e-02 -5.9148e-03 3.5973e-03 -2.9955e-02 -1.9442e-02 1.9502e-03 2.6199e-02 1.0482e-02 1.7694e-03 -3.5271e-03 3.5769e-02 9.3054e-03 -2.6098e-03 -1.6110e-02 4.9973e-06 9.1700e-06 1.2909e-05 -1.0665e-04 5.6031e-06 1.2343e-05 2.2613e-05 -1.5358e-03 -1.5328e-03 -1.5300e-03 1.9133e-02 -1.5354e-03 -1.5304e-03 -1.5229e-03
1.1714e-03 3.9876e-04 2.0548e-03 7.6663e-04 -1.6138e-03 -1.7683e-04 -2.6035e-03 -6.5311e-04 -2.2268e-04 -1.1443e-03 -1.9601e-03 2.4378e-03 9.7822e-05 1.4488e-03 -8.8867e-08 -2.9057e-07 6.0563e-07 4.4140e-06 -1.0640e-06 -5.1831e-07 -1.5592e-06 -9.8188e-08 -2.4602e-07 4.1020e-07 -1.5354e-03 1.5391e-03 -4.1282e-07 -1.1751e-06
3.6646e-03 2.4089e-03 -1.2520e-03 1.9201e-03 7.6740e-04 -5.4544e-03 -2.1052e-03 -2.0424e-03 -1.3427e-03 6.9590e-04 -2.5948e-03 -4.2846e-04 4.5728e-03 1.1710e-03 2.8622e-07 1.4066e-08 -1.3513e-06 1.0248e-05 -5.1822e-07 -4.3770e-06 -1.6564e-06 1.7669e-07 -2.2844e-08 -1.0230e-06 -1.5304e-03 -4.1282e-07 1.5367e-03 -1.2465e-06
1.6240e-03 1.6210e-04 -3.8662e-06 5.0212e-03 8.7713e-05 -1.2426e-04 -6.7402e-03 -9.0674e-04 -9.2365e-05 7.7225e-08 -4.3116e-03 -5.0937e-05 6.7067e-05 5.2895e-03 -1.1323e-06 -1.5326e-06 -1.5844e-06 1.8061e-05 -1.5585e-06 -1.6559e-06 -4.8645e-06 -8.6277e-07 -1.1562e-06 -1.1941e-06 -1.5229e-03 -1.1751e-06 -1.2465e-06 1.5363e-03
B*inv(B):
1.0000e+00 -7.1042e-12 -7.7739e-12 2.7194e-11 -2.3030e-11 -2.3622e-12 1.6791e-11 1.6855e-12 3.7279e-12 5.5280e-12 -1.3805e-11 1.2505e-11 -8.6503e-13 -8.0886e-12 -9.9400e-15 -7.8110e-15 -1.1100e-14 1.8792e-14 -1.8207e-14 6.6412e-14 -1.1054e-14 -2.4236e-12 -1.2365e-12 8.8498e-13 2.7821e-12 -6.7431e-13 -1.4959e-12 1.1687e-12
-4.5257e-12 1.0000e+00 -6.2867e-12 2.1806e-11 -1.9814e-11 -3.1426e-12 1.3381e-11 1.8667e-12 2.1778e-12 3.2970e-12 -1.0187e-11 1.1729e-11 -8.7694e-13 -5.3948e-12 -7.4806e-15 -5.0470e-15 -8.9468e-15 -3.6776e-16 -1.2099e-14 5.2685e-14 -5.3920e-15 -2.7615e-12 -1.9950e-12 8.4913e-13 3.7064e-12 -7.5460e-13 -1.7282e-12 1.5787e-12
-1.1031e-11 -1.0031e-11 1.0000e+00 4.5382e-11 -4.1160e-11 -6.7466e-12 2.6525e-11 4.3665e-12 4.6908e-12 7.5800e-12 -2.1138e-11 2.1814e-11 -1.9157e-12 -1.1950e-11 -1.5790e-14 -1.0877e-14 -1.5367e-14 1.0332e-14 -2.2665e-14 1.1742e-13 -1.3428e-14 -6.1566e-12 -4.1473e-12 1.1836e-12 8.4515e-12 -1.6522e-12 -3.6782e-12 3.1664e-12
-2.2167e-12 -2.5437e-12 -7.0957e-12 1.0000e+00 -2.1778e-11 -1.1150e-12 1.2101e-11 1.4671e-12 1.1074e-12 5.7059e-12 -1.0292e-11 1.2269e-11 -2.6327e-12 -5.3846e-12 -6.2472e-15 -3.9990e-15 -9.6031e-15 -3.4469e-15 -1.2838e-14 4.8775e-14 3.4734e-15 -2.2713e-12 -1.5696e-12 8.6787e-13 2.7558e-12 -3.0104e-13 -1.6093e-12 1.1786e-12
-3.0340e-12 -2.9504e-12 -5.8059e-12 1.8982e-11 1.0000e+00 -1.4307e-12 9.3445e-12 1.2044e-12 6.3975e-13 3.3456e-12 -8.1881e-12 7.8565e-12 -1.3266e-12 -4.0066e-12 -5.5172e-15 -3.7445e-15 -7.1578e-15 4.2310e-15 -9.2363e-15 4.1462e-14 -5.2957e-15 -2.4404e-12 -1.5724e-12 7.6826e-13 3.1429e-12 -7.9435e-13 -1.3891e-12 1.4597e-12
-1.0055e-11 -9.0164e-12 -1.4599e-11 5.7568e-11 -5.0646e-11 1.0000e+00 3.5321e-11 8.8218e-12 6.3849e-12 7.6999e-12 -2.1991e-11 2.4622e-11 -3.0038e-12 -1.4866e-11 -1.9473e-14 -1.2548e-14 -2.0109e-14 4.0593e-15 -2.5478e-14 1.4512e-13 -1.6254e-14 -8.6573e-12 -6.1884e-12 1.4937e-12 1.2771e-11 -3.6273e-12 -4.8475e-12 4.9774e-12
-8.2125e-12 -4.9051e-12 -8.9749e-12 3.4555e-11 -3.0434e-11 -4.9936e-12 1.0000e+00 3.6371e-12 2.6510e-12 7.8761e-12 -1.4589e-11 1.5353e-11 -1.3104e-12 -8.4555e-12 -1.1127e-14 -6.0208e-15 -1.6476e-14 -5.4436e-15 -1.4043e-14 8.2779e-14 -1.1197e-14 -5.3263e-12 -3.9053e-12 1.1227e-12 7.7258e-12 -2.3144e-12 -2.9013e-12 3.1042e-12
-6.2678e-12 -9.5027e-12 -8.5018e-12 3.3789e-11 -2.5807e-11 -3.2876e-12 2.0525e-11 1.0000e+00 5.5665e-12 5.7359e-12 -1.7676e-11 1.4079e-11 -2.7617e-13 -1.0118e-11 -1.2117e-14 -1.0340e-14 -1.2661e-14 3.4457e-14 -2.2402e-14 7.6406e-14 -1.7423e-14 -2.5625e-12 -1.1477e-12 1.0294e-12 2.7333e-12 -8.1825e-13 -1.5376e-12 1.2160e-12
-4.0017e-12 -3.4176e-12 -6.4315e-12 2.2980e-11 -2.0443e-11 -3.4322e-12 1.3285e-11 9.9454e-13 1.0000e+00 3.9019e-12 -9.7602e-12 1.1719e-11 -1.0581e-12 -5.8069e-12 -6.6598e-15 -6.1450e-15 -7.4454e-15 3.8060e-15 -1.2972e-14 5.4803e-14 -5.9358e-15 -2.7008e-12 -2.0714e-12 8.9231e-13 3.6677e-12 -7.5282e-13 -1.7146e-12 1.6276e-12
-9.9134e-12 -9.3101e-12 -1.2142e-11 4.5842e-11 -3.8311e-11 -6.6613e-12 2.8383e-11 4.5690e-12 4.8631e-12 1.0000e+00 -2.3014e-11 2.2987e-11 -2.5179e-12 -1.2159e-11 -1.5606e-14 -1.0342e-14 -1.5362e-14 1.1276e-14 -2.4345e-14 1.1252e-13 -1.5939e-14 -6.1757e-12 -4.1339e-12 1.1534e-12 8.3822e-12 -1.4377e-12 -3.7119e-12 3.0882e-12
-2.2504e-12 -1.6870e-12 -8.9242e-12 2.3104e-11 -2.5179e-11 -6.1226e-13 1.3221e-11 7.8720e-13 3.6042e-13 6.7045e-12 1.0000e+00 1.4157e-11 -3.2895e-12 -5.6333e-12 -6.5759e-15 -3.8266e-15 -1.0180e-14 -1.0656e-14 -1.3471e-14 4.9540e-14 9.5770e-15 -2.3437e-12 -1.7700e-12 9.9326e-13 2.7505e-12 -1.5760e-13 -1.7756e-12 1.2510e-12
-3.6842e-12 -3.6320e-12 -4.7137e-12 1.8469e-11 -1.5228e-11 -1.3605e-12 9.3907e-12 1.2044e-12 8.2759e-13 3.1800e-12 -7.4429e-12 1.0000e+00 -1.4367e-12 -3.3502e-12 -5.4836e-15 -3.6976e-15 -6.5658e-15 9.1368e-15 -9.7834e-15 4.2265e-14 -6.8136e-15 -2.5760e-12 -1.6171e-12 8.2945e-13 3.1678e-12 -7.6515e-13 -1.4253e-12 1.5964e-12
-8.8329e-12 -1.0580e-11 -1.5227e-11 5.7698e-11 -5.2448e-11 -7.9762e-12 3.5037e-11 7.5855e-12 6.5590e-12 6.7793e-12 -2.3274e-11 2.2759e-11 1.0000e+00 -1.4134e-11 -1.9091e-14 -1.2049e-14 -2.0990e-14 -1.0839e-14 -2.5257e-14 1.4123e-13 -1.3279e-14 -8.6810e-12 -6.2228e-12 1.4760e-12 1.2916e-11 -3.7325e-12 -4.8591e-12 5.0333e-12
-9.4973e-12 -5.1183e-12 -9.1022e-12 3.5091e-11 -3.0639e-11 -5.3080e-12 2.1211e-11 3.6811e-12 2.2679e-12 7.9187e-12 -1.4846e-11 1.4913e-11 -1.0085e-12 1.0000e+00 -1.1224e-14 -6.2775e-15 -1.4602e-14 -1.3264e-14 -1.3966e-14 8.4222e-14 -9.8298e-15 -5.4508e-12 -4.0364e-12 1.1927e-12 7.8710e-12 -2.4320e-12 -2.9226e-12 3.1695e-12
-1.1741e-12 1.0697e-13 -6.3777e-12 1.3633e-11 -1.6451e-11 1.9580e-13 9.0496e-12 -7.3830e-15 3.0712e-13 3.8702e-12 -5.1590e-12 9.0772e-12 -2.6689e-12 -3.2090e-12 1.0000e+00 -9.4586e-16 -5.5808e-15 -7.3951e-15 -7.7955e-15 3.5834e-14 8.0751e-16 -2.0715e-12 -1.4600e-12 5.1972e-13 2.9109e-12 -3.1221e-13 -1.3915e-12 1.0484e-12
-3.7992e-12 -3.5064e-12 -7.3148e-12 2.3960e-11 -2.2260e-11 -2.3548e-12 1.3896e-11 1.2912e-12 2.1760e-12 4.2451e-12 -1.0395e-11 1.2443e-11 -1.2020e-12 -6.3967e-12 -7.5309e-15 1.0000e+00 -8.9212e-15 3.1780e-15 -1.2584e-14 5.5833e-14 -6.7771e-15 -2.9147e-12 -1.8017e-12 7.4288e-13 3.7876e-12 -7.5837e-13 -1.7650e-12 1.4535e-12
-1.1336e-11 -9.7737e-12 -1.4348e-11 4.7411e-11 -4.0520e-11 -7.1161e-12 2.6734e-11 3.1053e-12 4.2183e-12 9.1361e-12 -2.0537e-11 2.3185e-11 -2.0436e-12 -1.1044e-11 -1.4708e-14 -1.0017e-14 1.0000e+00 3.2821e-15 -2.4157e-14 1.0748e-13 -1.5373e-14 -6.1073e-12 -4.1810e-12 1.2644e-12 8.6113e-12 -2.1920e-12 -3.5831e-12 3.3582e-12
-3.8767e-12 -5.1638e-12 -3.8178e-12 1.5878e-11 -1.4871e-11 -2.7519e-12 1.0756e-11 1.6780e-12 3.0107e-12 3.1496e-12 -8.1481e-12 7.5349e-12 -9.5007e-13 -4.6962e-12 -6.0973e-15 -5.1194e-15 -7.3274e-15 1.0000e+00 -1.1670e-14 4.5239e-14 -1.2323e-14 -2.0880e-12 -1.0655e-12 5.5013e-13 2.7505e-12 -6.6270e-13 -1.1923e-12 9.9820e-13
-3.2463e-12 -2.5347e-12 -4.8818e-12 1.7880e-11 -1.5306e-11 -1.4102e-12 1.0773e-11 1.5090e-12 9.4994e-13 3.5517e-12 -8.0860e-12 8.2207e-12 -1.2849e-12 -4.7651e-12 -5.7797e-15 -3.1447e-15 -6.9626e-15 2.9091e-15 1.0000e+00 3.9936e-14 -3.2565e-15 -2.0987e-12 -1.4605e-12 6.1494e-13 3.0643e-12 -8.6699e-13 -1.2963e-12 1.1156e-12
-9.9236e-12 -1.0679e-11 -1.7181e-11 6.0162e-11 -5.0798e-11 -7.7950e-12 3.6366e-11 8.2623e-12 6.7296e-12 7.7879e-12 -2.1652e-11 2.3778e-11 -3.3520e-12 -1.4776e-11 -1.9208e-14 -1.1504e-14 -2.2320e-14 6.4809e-15 -2.5323e-14 1.0000e+00 -1.7165e-14 -8.5964e-12 -6.1011e-12 1.5432e-12 1.2451e-11 -3.3654e-12 -4.8120e-12 4.8326e-12
-7.9736e-12 -5.9497e-12 -7.8256e-12 3.3211e-11 -3.0869e-11 -4.5578e-12 2.0364e-11 3.4868e-12 2.8063e-12 7.2890e-12 -1.3877e-11 1.4847e-11 -1.3961e-12 -8.6278e-12 -1.0377e-14 -6.5252e-15 -1.6551e-14 -2.3245e-15 -1.5064e-14 8.1491e-14 1.0000e+00 -5.0126e-12 -3.5739e-12 9.4679e-13 7.3817e-12 -2.0204e-12 -2.8490e-12 2.9359e-12
-4.6541e-13 1.1668e-12 -5.2408e-12 8.8662e-12 -1.1928e-11 7.5845e-13 5.2482e-12 1.5210e-13 -1.3223e-12 3.2739e-12 -2.6557e-12 7.1431e-12 -3.1267e-12 -1.7488e-12 -1.4173e-15 8.8905e-16 -2.6951e-15 -2.2825e-14 -3.7984e-15 2.6624e-14 6.8404e-15 1.0000e+00 -1.5484e-12 3.7534e-13 2.9623e-12 -1.6886e-13 -1.3489e-12 1.0004e-12
-3.3387e-12 -4.5302e-12 -6.7537e-12 2.3946e-11 -2.3715e-11 -2.5445e-12 1.5582e-11 1.1616e-12 2.8821e-12 3.4129e-12 -1.0394e-11 1.2489e-11 -1.1660e-12 -6.9527e-12 -7.7486e-15 -5.2874e-15 -8.8673e-15 -4.7531e-16 -1.2741e-14 5.6802e-14 -5.8452e-15 -2.9748e-12 1.0000e+00 7.0064e-13 3.8276e-12 -7.6057e-13 -1.7799e-12 1.4049e-12
-1.0296e-11 -9.7246e-12 -1.2250e-11 4.6951e-11 -3.5838e-11 -7.0890e-12 2.7500e-11 3.5045e-12 4.1852e-12 9.3692e-12 -2.0979e-11 2.0675e-11 -2.2935e-12 -1.0708e-11 -1.5243e-14 -1.0186e-14 -1.9057e-14 9.4369e-15 -2.3517e-14 1.1502e-13 -1.3121e-14 -6.0876e-12 -4.1948e-12 1.0000e+00 8.6744e-12 -2.4070e-12 -3.5418e-12 3.4364e-12
-5.0180e-12 -6.1638e-12 -2.0068e-12 1.5116e-11 -1.1361e-11 -3.2968e-12 1.0437e-11 2.2885e-12 3.9633e-12 2.2054e-12 -7.4203e-12 5.6679e-12 -1.8836e-13 -5.0777e-12 -5.8087e-15 -5.0952e-15 -6.5056e-15 1.6195e-14 -1.1183e-14 4.5113e-14 -1.6629e-14 -2.0147e-12 -8.6451e-13 4.2429e-13 1.0000e+00 -8.0606e-13 -1.0255e-12 9.2737e-13
-3.4044e-12 -2.4647e-12 -5.6727e-12 1.7537e-11 -1.5368e-11 -1.3867e-12 1.0863e-11 1.6566e-12 1.0456e-12 3.8699e-12 -7.9878e-12 8.9017e-12 -1.2461e-12 -5.0666e-12 -5.4233e-15 -3.0431e-15 -7.3711e-15 -2.1927e-15 -7.3372e-15 3.8658e-14 -1.9255e-15 -1.9624e-12 -1.4156e-12 5.5401e-13 3.0402e-12 1.0000e+00 -1.2594e-12 9.7788e-13
-9.1500e-12 -1.0191e-11 -1.5103e-11 5.5167e-11 -5.2186e-11 -8.2689e-12 3.5399e-11 7.2937e-12 5.9575e-12 8.2999e-12 -2.2409e-11 2.4990e-11 -3.0068e-12 -1.4493e-11 -1.9245e-14 -1.4111e-14 -2.0531e-14 -3.4209e-15 -2.6204e-14 1.3716e-13 -1.4924e-14 -8.5724e-12 -6.0668e-12 1.5645e-12 1.2323e-11 -3.2618e-12 1.0000e+00 4.7766e-12
-8.0318e-12 -6.2291e-12 -7.3260e-12 3.3918e-11 -3.0391e-11 -4.4816e-12 2.0567e-11 3.8645e-12 3.1610e-12 7.0981e-12 -1.3488e-11 1.5378e-11 -1.6548e-12 -8.3258e-12 -1.0798e-14 -7.3257e-15 -1.6175e-14 -9.8324e-15 -1.6168e-14 8.1937e-14 -1.1156e-14 -4.8879e-12 -3.4423e-12 8.7740e-13 7.2360e-12 -1.9028e-12 -2.8276e-12 1.0000e+00
-------------- next part --------------
#include <armadillo>
#include <iostream>
using namespace std;
int main(int argc, char** argv) {
cout << "Program start." << endl;
// orthogonal matrix B
arma::mat B = {
{4.4143233708559126e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02, 5.2073847122169718e+02, 3.7152475267862627e+02,
7.2780748992934855e+02, 3.3080575372348960e+02, 2.8314623993465835e+02,
8.8281007483849464e+02, 5.2067395105565447e+02, 2.4141143908559127e+02,
3.7152475267862627e+02, 7.2780748992934855e+02, 3.3080575372348960e+02,
2.8314623993465835e+02, 8.8281007483849464e+02, 5.2067395105565447e+02,
1.6208440694948538e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02},
{3.7152475267862627e+02, 6.2101084268392844e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02, 3.7152475267862627e+02, 7.0031453182003440e+02,
9.9030414001983547e+02, 3.6861421044114815e+02, 3.4005724170194975e+02,
1.2395375465611073e+03, 7.6604313005847553e+02, 3.7152475267862627e+02,
4.2098749968392849e+02, 9.9030414001983547e+02, 3.6861421044114815e+02,
3.4005724170194975e+02, 1.2395375465611073e+03, 7.6604313005847553e+02,
3.7152475267862627e+02, 3.4166046754782263e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02},
{7.2780748992934855e+02, 9.9030414001983547e+02, 2.5045845237611029e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03, 7.2780748992934855e+02, 9.9030414001983547e+02,
2.5838737408972088e+03, 6.8488997533484269e+02, 6.1454067400014367e+02,
3.1970393777315066e+03, 1.7912477710073911e+03, 7.2780748992934855e+02,
9.9030414001983547e+02, 2.3045467087611028e+03, 6.8488997533484269e+02,
6.1454067400014367e+02, 3.1970393777315066e+03, 1.7912477710073911e+03,
7.2780748992934855e+02, 9.9030414001983547e+02, 2.2252196766249967e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03},
{3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
4.2406647996299085e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02, 3.3080575372348960e+02, 3.6861421044114815e+02,
6.8488997533484269e+02, 5.0337276309909669e+02, 2.8444507405058255e+02,
8.2970199478326299e+02, 5.0392512378152617e+02, 3.3080575372348960e+02,
3.6861421044114815e+02, 6.8488997533484269e+02, 2.2404573096299083e+02,
2.8444507405058255e+02, 8.2970199478326299e+02, 5.0392512378152617e+02,
3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
1.4471869882688495e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02},
{2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 4.3532000801648621e+02, 9.3167339337363285e+02,
6.1263109070298322e+02, 2.8314623993465835e+02, 3.4005724170194969e+02,
6.1454067400014367e+02, 2.8444507405058255e+02, 5.1462622415259204e+02,
9.3167339337363285e+02, 6.1263109070298322e+02, 2.8314623993465835e+02,
3.4005724170194969e+02, 6.1454067400014367e+02, 2.8444507405058255e+02,
2.3529919201648619e+02, 9.3167339337363285e+02, 6.1263109070298322e+02,
2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 1.5597215988038030e+02, 9.3167339337363285e+02,
6.1263109070298322e+02},
{8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.8470512303564547e+03,
2.7217094006090792e+03, 8.8281007483849464e+02, 1.2395375465611073e+03,
3.1970393777315066e+03, 8.2970199478326299e+02, 9.3167339337363285e+02,
4.9263330664925606e+03, 2.7217094006090792e+03, 8.8281007483849464e+02,
1.2395375465611073e+03, 3.1970393777315066e+03, 8.2970199478326299e+02,
9.3167339337363285e+02, 4.6470060343564546e+03, 2.7217094006090792e+03,
8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.5676790022203486e+03,
2.7217094006090792e+03},
{5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.7360181038689775e+03, 5.2067395105565436e+02, 7.6604313005847553e+02,
1.7912477710073908e+03, 5.0392512378152611e+02, 6.1263109070298322e+02,
2.7217094006090788e+03, 1.8153116640050832e+03, 5.2067395105565436e+02,
7.6604313005847553e+02, 1.7912477710073908e+03, 5.0392512378152611e+02,
6.1263109070298322e+02, 2.7217094006090788e+03, 1.5359846318689774e+03,
5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.4566575997328716e+03},
{5.2073847122169718e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02, 6.6301240963303087e+02, 3.7152475267862627e+02,
7.2780748992934855e+02, 3.3080575372348960e+02, 2.8314623993465835e+02,
8.8281007483849464e+02, 5.2067395105565447e+02, 1.6208440694948538e+02,
3.7152475267862627e+02, 7.2780748992934855e+02, 3.3080575372348960e+02,
2.8314623993465835e+02, 8.8281007483849464e+02, 5.2067395105565447e+02,
1.9829594538151753e+01, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02},
{3.7152475267862627e+02, 7.0031453182003440e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02, 3.7152475267862627e+02, 8.4259049323136810e+02,
9.9030414001983547e+02, 3.6861421044114815e+02, 3.4005724170194975e+02,
1.2395375465611073e+03, 7.6604313005847553e+02, 3.7152475267862627e+02,
3.4166046754782258e+02, 9.9030414001983547e+02, 3.6861421044114815e+02,
3.4005724170194975e+02, 1.2395375465611073e+03, 7.6604313005847553e+02,
3.7152475267862627e+02, 1.9940565513648897e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02},
{7.2780748992934855e+02, 9.9030414001983547e+02, 2.5838737408972088e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03, 7.2780748992934855e+02, 9.9030414001983547e+02,
2.7261616853085425e+03, 6.8488997533484269e+02, 6.1454067400014367e+02,
3.1970393777315066e+03, 1.7912477710073911e+03, 7.2780748992934855e+02,
9.9030414001983547e+02, 2.2252196766249967e+03, 6.8488997533484269e+02,
6.1454067400014367e+02, 3.1970393777315066e+03, 1.7912477710073911e+03,
7.2780748992934855e+02, 9.9030414001983547e+02, 2.0829648642136631e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03},
{3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
5.0337276309909669e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02, 3.3080575372348960e+02, 3.6861421044114815e+02,
6.8488997533484269e+02, 6.4564657851043034e+02, 2.8444507405058255e+02,
8.2970199478326299e+02, 5.0392512378152617e+02, 3.3080575372348960e+02,
3.6861421044114815e+02, 6.8488997533484269e+02, 1.4471869882688495e+02,
2.8444507405058255e+02, 8.2970199478326299e+02, 5.0392512378152617e+02,
3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
2.4638864155513134e+00, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02},
{2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 5.1462622415259204e+02, 9.3167339337363285e+02,
6.1263109070298322e+02, 2.8314623993465835e+02, 3.4005724170194969e+02,
6.1454067400014367e+02, 2.8444507405058255e+02, 6.5690009456392568e+02,
9.3167339337363285e+02, 6.1263109070298322e+02, 2.8314623993465835e+02,
3.4005724170194969e+02, 6.1454067400014367e+02, 2.8444507405058255e+02,
1.5597215988038030e+02, 9.3167339337363285e+02, 6.1263109070298322e+02,
2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 1.3717347469046670e+01, 9.3167339337363285e+02,
6.1263109070298322e+02},
{8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.9263330664925606e+03,
2.7217094006090792e+03, 8.8281007483849464e+02, 1.2395375465611073e+03,
3.1970393777315066e+03, 8.2970199478326299e+02, 9.3167339337363285e+02,
5.0686271259038940e+03, 2.7217094006090792e+03, 8.8281007483849464e+02,
1.2395375465611073e+03, 3.1970393777315066e+03, 8.2970199478326299e+02,
9.3167339337363285e+02, 4.5676790022203486e+03, 2.7217094006090792e+03,
8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.4254241898090149e+03,
2.7217094006090792e+03},
{5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.8153116640050832e+03, 5.2067395105565436e+02, 7.6604313005847553e+02,
1.7912477710073908e+03, 5.0392512378152611e+02, 6.1263109070298322e+02,
2.7217094006090788e+03, 1.9575960134164168e+03, 5.2067395105565436e+02,
7.6604313005847553e+02, 1.7912477710073908e+03, 5.0392512378152611e+02,
6.1263109070298322e+02, 2.7217094006090788e+03, 1.4566575997328716e+03,
5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.3144027873215380e+03},
{2.4141143908559127e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02, 1.6208440694948538e+02, 3.7152475267862627e+02,
7.2780748992934855e+02, 3.3080575372348960e+02, 2.8314623993465835e+02,
8.8281007483849464e+02, 5.2067395105565447e+02, 9.0353825815650282e+02,
4.7152475267862627e+02, 8.2780748992934855e+02, 4.3080575372348960e+02,
3.8314623993465835e+02, 9.8281007483849464e+02, 6.2067395105565447e+02,
5.2073847122169718e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02},
{3.7152475267862627e+02, 4.2098749968392849e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02, 3.7152475267862627e+02, 3.4166046754782263e+02,
9.9030414001983547e+02, 3.6861421044114815e+02, 3.4005724170194975e+02,
1.2395375465611073e+03, 7.6604313005847553e+02, 4.7152475267862627e+02,
1.0831143187548400e+03, 1.0903041400198354e+03, 4.6861421044114815e+02,
4.4005724170194975e+02, 1.3395375465611073e+03, 8.6604313005847553e+02,
3.7152475267862627e+02, 7.0031453182003440e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02},
{7.2780748992934855e+02, 9.9030414001983547e+02, 2.3045467087611028e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03, 7.2780748992934855e+02, 9.9030414001983547e+02,
2.2252196766249967e+03, 6.8488997533484269e+02, 6.1454067400014367e+02,
3.1970393777315066e+03, 1.7912477710073911e+03, 8.2780748992934855e+02,
1.0903041400198354e+03, 2.9666735278320143e+03, 7.8488997533484269e+02,
7.1454067400014367e+02, 3.2970393777315066e+03, 1.8912477710073911e+03,
7.2780748992934855e+02, 9.9030414001983547e+02, 2.5838737408972088e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03},
{3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
2.2404573096299083e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02, 3.3080575372348960e+02, 3.6861421044114815e+02,
6.8488997533484269e+02, 1.4471869882688495e+02, 2.8444507405058255e+02,
8.2970199478326299e+02, 5.0392512378152617e+02, 4.3080575372348960e+02,
4.6861421044114815e+02, 7.8488997533484269e+02, 5.2405573096299076e+02,
3.8444507405058255e+02, 9.2970199478326299e+02, 6.0392512378152617e+02,
3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
5.0337276309909669e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02},
{2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 2.3529919201648619e+02, 9.3167339337363285e+02,
6.1263109070298322e+02, 2.8314623993465835e+02, 3.4005724170194969e+02,
6.1454067400014367e+02, 2.8444507405058255e+02, 1.5597215988038030e+02,
9.3167339337363285e+02, 6.1263109070298322e+02, 3.8314623993465835e+02,
4.4005724170194969e+02, 7.1454067400014367e+02, 3.8444507405058255e+02,
8.9742601108739768e+02, 1.0316733933736327e+03, 7.1263109070298322e+02,
2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 5.1462622415259204e+02, 9.3167339337363285e+02,
6.1263109070298322e+02},
{8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.6470060343564546e+03,
2.7217094006090792e+03, 8.8281007483849464e+02, 1.2395375465611073e+03,
3.1970393777315066e+03, 8.2970199478326299e+02, 9.3167339337363285e+02,
4.5676790022203486e+03, 2.7217094006090792e+03, 9.8281007483849464e+02,
1.3395375465611073e+03, 3.2970393777315066e+03, 9.2970199478326299e+02,
1.0316733933736327e+03, 5.3091328534273662e+03, 2.8217094006090792e+03,
8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.9263330664925606e+03,
2.7217094006090792e+03},
{5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.5359846318689774e+03, 5.2067395105565436e+02, 7.6604313005847553e+02,
1.7912477710073908e+03, 5.0392512378152611e+02, 6.1263109070298322e+02,
2.7217094006090788e+03, 1.4566575997328716e+03, 6.2067395105565436e+02,
8.6604313005847553e+02, 1.8912477710073908e+03, 6.0392512378152605e+02,
7.1263109070298322e+02, 2.8217094006090788e+03, 2.1981114509398890e+03,
5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.8153116640050832e+03},
{1.6208440694948538e+02, 3.7152475267862627e+02, 7.2780748992934855e+02,
3.3080575372348960e+02, 2.8314623993465835e+02, 8.8281007483849464e+02,
5.2067395105565447e+02, 1.9829594538151753e+01, 3.7152475267862627e+02,
7.2780748992934855e+02, 3.3080575372348960e+02, 2.8314623993465835e+02,
8.8281007483849464e+02, 5.2067395105565447e+02, 5.2073847122169718e+02,
3.7152475267862627e+02, 7.2780748992934855e+02, 3.3080575372348960e+02,
2.8314623993465835e+02, 8.8281007483849464e+02, 5.2067395105565447e+02,
1.4123766281385688e+03, 4.7152475267862627e+02, 8.2780748992934855e+02,
4.3080575372348960e+02, 3.8314623993465835e+02, 9.8281007483849464e+02,
6.2067395105565447e+02},
{3.7152475267862627e+02, 3.4166046754782258e+02, 9.9030414001983547e+02,
3.6861421044114815e+02, 3.4005724170194975e+02, 1.2395375465611073e+03,
7.6604313005847553e+02, 3.7152475267862627e+02, 1.9940565513648897e+02,
9.9030414001983547e+02, 3.6861421044114815e+02, 3.4005724170194975e+02,
1.2395375465611073e+03, 7.6604313005847553e+02, 3.7152475267862627e+02,
7.0031453182003440e+02, 9.9030414001983547e+02, 3.6861421044114815e+02,
3.4005724170194975e+02, 1.2395375465611073e+03, 7.6604313005847553e+02,
4.7152475267862627e+02, 1.5919526887369059e+03, 1.0903041400198354e+03,
4.6861421044114815e+02, 4.4005724170194975e+02, 1.3395375465611073e+03,
8.6604313005847553e+02},
{7.2780748992934855e+02, 9.9030414001983547e+02, 2.2252196766249967e+03,
6.8488997533484269e+02, 6.1454067400014367e+02, 3.1970393777315066e+03,
1.7912477710073911e+03, 7.2780748992934855e+02, 9.9030414001983547e+02,
2.0829648642136631e+03, 6.8488997533484269e+02, 6.1454067400014367e+02,
3.1970393777315066e+03, 1.7912477710073911e+03, 7.2780748992934855e+02,
9.9030414001983547e+02, 2.5838737408972088e+03, 6.8488997533484269e+02,
6.1454067400014367e+02, 3.1970393777315066e+03, 1.7912477710073911e+03,
8.2780748992934855e+02, 1.0903041400198354e+03, 3.4755118978140804e+03,
7.8488997533484269e+02, 7.1454067400014367e+02, 3.2970393777315066e+03,
1.8912477710073911e+03},
{3.3080575372348960e+02, 3.6861421044114815e+02, 6.8488997533484269e+02,
1.4471869882688495e+02, 2.8444507405058255e+02, 8.2970199478326299e+02,
5.0392512378152617e+02, 3.3080575372348960e+02, 3.6861421044114815e+02,
6.8488997533484269e+02, 2.4638864155513134e+00, 2.8444507405058255e+02,
8.2970199478326299e+02, 5.0392512378152617e+02, 3.3080575372348960e+02,
3.6861421044114815e+02, 6.8488997533484269e+02, 5.0337276309909669e+02,
2.8444507405058255e+02, 8.2970199478326299e+02, 5.0392512378152617e+02,
4.3080575372348960e+02, 4.6861421044114815e+02, 7.8488997533484269e+02,
7.4563757551043034e+02, 3.8444507405058255e+02, 9.2970199478326299e+02,
6.0392512378152617e+02},
{2.8314623993465835e+02, 3.4005724170194969e+02, 6.1454067400014367e+02,
2.8444507405058255e+02, 1.5597215988038030e+02, 9.3167339337363285e+02,
6.1263109070298322e+02, 2.8314623993465835e+02, 3.4005724170194969e+02,
6.1454067400014367e+02, 2.8444507405058255e+02, 1.3717347469046670e+01,
9.3167339337363285e+02, 6.1263109070298322e+02, 2.8314623993465835e+02,
3.4005724170194969e+02, 6.1454067400014367e+02, 2.8444507405058255e+02,
5.1462622415259204e+02, 9.3167339337363285e+02, 6.1263109070298322e+02,
3.8314623993465835e+02, 4.4005724170194969e+02, 7.1454067400014367e+02,
3.8444507405058255e+02, 1.4062643810694635e+03, 1.0316733933736327e+03,
7.1263109070298322e+02},
{8.8281007483849464e+02, 1.2395375465611073e+03, 3.1970393777315066e+03,
8.2970199478326299e+02, 9.3167339337363285e+02, 4.5676790022203486e+03,
2.7217094006090792e+03, 8.8281007483849464e+02, 1.2395375465611073e+03,
3.1970393777315066e+03, 8.2970199478326299e+02, 9.3167339337363285e+02,
4.4254241898090149e+03, 2.7217094006090792e+03, 8.8281007483849464e+02,
1.2395375465611073e+03, 3.1970393777315066e+03, 8.2970199478326299e+02,
9.3167339337363285e+02, 4.9263330664925606e+03, 2.7217094006090792e+03,
9.8281007483849464e+02, 1.3395375465611073e+03, 3.2970393777315066e+03,
9.2970199478326299e+02, 1.0316733933736327e+03, 5.8179712234094322e+03,
2.8217094006090792e+03},
{5.2067395105565436e+02, 7.6604313005847553e+02, 1.7912477710073908e+03,
5.0392512378152611e+02, 6.1263109070298322e+02, 2.7217094006090788e+03,
1.4566575997328716e+03, 5.2067395105565436e+02, 7.6604313005847553e+02,
1.7912477710073908e+03, 5.0392512378152611e+02, 6.1263109070298322e+02,
2.7217094006090788e+03, 1.3144027873215380e+03, 5.2067395105565436e+02,
7.6604313005847553e+02, 1.7912477710073908e+03, 5.0392512378152611e+02,
6.1263109070298322e+02, 2.7217094006090788e+03, 1.8153116640050832e+03,
6.2067395105565436e+02, 8.6604313005847553e+02, 1.8912477710073908e+03,
6.0392512378152605e+02, 7.1263109070298322e+02, 2.8217094006090788e+03,
2.7069498209219551e+03}};
cout << "Test armadillo matrix inversion" << endl;
arma::arma_version ver;
cout << "ARMA version: " << ver.as_string() << endl;
cout << "B:\n" << B << endl;
// Inverse of the Matrix
arma::mat B_inv = arma::inv(B);
cout << "inv(B):\n" << B_inv << endl;
// Matrix times its inverse should give the unit matrix
arma::mat I = B*B_inv;
cout << "B*inv(B):\n" << I << endl;
return 0;
}
More information about the debian-science-maintainers
mailing list