Hi, I need to adapt typhon.oem.retrieval_gain_matrix to allow broadcasting. That would require two changes:
- replace every occurrence of .T by .swapaxes(-1, -2) - replace scipy.linalg.inv by numpy.linalg.inv, as the former does not broadcast but the latter does. The former is fully backward compatible. I don't know if that is the case for the latter. Has there been a conscious choice between the two? As far as I can see this code is not covered by unit tests so changing it may carry a risk. regards, Gerrit. _______________________________________________ typhon.mi mailing list typhon.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/typhon.mi
