Skip to content

Instantly share code, notes, and snippets.

@evanbiederstedt
Created July 2, 2015 16:51
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save evanbiederstedt/857e6480f9e374aeeffa to your computer and use it in GitHub Desktop.
Save evanbiederstedt/857e6480f9e374aeeffa to your computer and use it in GitHub Desktop.
np.linalg.inv() numerical errors
testmattt = Sij[13]
testmattt_inv = np.linalg.inv(Sij[13])
print testmattt
OUTPUT
[[ 3.60611595e-10 3.54999009e-10 3.49422963e-10 ..., -3.54999009e-10
-3.60611595e-10 -3.54999009e-10]
[ 3.54999009e-10 3.60611595e-10 3.54999009e-10 ..., -3.49422963e-10
-3.54999009e-10 -3.60611595e-10]
[ 3.49422963e-10 3.54999009e-10 3.60611595e-10 ..., -3.54999009e-10
-3.49422963e-10 -3.54999009e-10]
...,
[ -3.54999009e-10 -3.49422963e-10 -3.54999009e-10 ..., 3.60611595e-10
3.54999009e-10 3.49422963e-10]
[ -3.60611595e-10 -3.54999009e-10 -3.49422963e-10 ..., 3.54999009e-10
3.60611595e-10 3.54999009e-10]
[ -3.54999009e-10 -3.60611595e-10 -3.54999009e-10 ..., 3.49422963e-10
3.54999009e-10 3.60611595e-10]]
print testmattt_inv
OUTPUT
[[ 5.90171114e+23 -2.97719733e+23 3.29856106e+23 ..., -2.46714190e+23
2.11331633e+22 -3.16779300e+23]
[ -1.02998439e+24 -3.45833499e+23 -5.59535211e+23 ..., 2.60455420e+23
6.80907569e+23 6.68450292e+21]
[ -4.99010695e+23 -2.36591710e+23 5.93246262e+21 ..., 6.66918182e+22
1.95742816e+23 -1.61205246e+23]
...,
[ 4.69633660e+23 -2.46248642e+22 2.33103810e+23 ..., -1.43035832e+23
-1.87582879e+22 -4.42910477e+22]
[ -2.41150378e+23 3.97836968e+22 1.62738017e+22 ..., 1.76385611e+23
-3.67537289e+23 -7.58324328e+22]
[ 1.40915869e+24 -1.07070275e+23 8.14514766e+23 ..., -2.60085893e+23
-6.93957565e+23 -3.74134341e+23]]
print np.dot(testmattt, testmattt_inv)
[[ 0.25 0.25 -0.71875 ..., -0.25 0. 0.75 ]
[ 3.75 -0.75 -1.78125 ..., -0.625 0.5 -0.625 ]
[-2.125 -0.5 -2.5546875 ..., -0.375 0.25 0.5 ]
...,
[-3.125 -2. -0.875 ..., 2.125 2.5 1.5 ]
[ 1.125 0.25 0.1875 ..., 0.625 1.25 0.25 ]
[-2.5 0.5 3.546875 ..., -0.875 -1.75 0.75 ]]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment