Update calculation script to match what we do internally

tools/calc_cat_matrices.py

@@ -66,8 +66,8 @@ def pre_calculate_cat(src_white, dest_white, m):
     src = alg.dot(m, src_white)
     dest = alg.dot(m, dest_white)
     m2 = alg.diag(alg.divide(dest, src))
-    to_d50 = alg.solve(m, alg.dot(m2, m))
-    to_d65 = alg.inv(to_d50)
+    to_d65 = alg.dot(alg.solve(m, m2), m)
+    to_d50 = alg.inv(to_d65)
     return to_d65, to_d50
 
 
@@ -78,7 +78,7 @@ def pre_calculate_cat(src_white, dest_white, m):
     # afterwards, we just retrieve the previously calculated matrices from the cache.
     # This allows us to check what the actual values we are using.
 
-    to_d65, to_d50 = pre_calculate_cat(white_d65, white_d50, bradford_m)
+    to_d65, to_d50 = pre_calculate_cat(white_d50, white_d65, bradford_m)
     print('===== Bradford To D65 =====')
     alg.pprint(to_d65)
     print('===== Bradford To D50 =====')