Closed spedygiorgio closed 9 years ago
It's not complicated to extract the code writtent by M Maechler. But the true question is : will we use other functions of expm in the future? for continuous time markov chains, expm and logm might be needed?
I do intend to implement functionality for CTMCs (for eg, computing P(t) - states probability at time t) that will require matrix exponentials.
So it's safer not to reimplement in markovchain what M Maechler has done in expm.
Agreed
Inviato da iPhone
Il giorno 03/lug/2015, alle ore 10:33, Christophe Dutang notifications@github.com ha scritto:
So it's safer not to reimplement in markovchain what M Maechler has done in expm.
— Reply to this email directly or view it on GitHub.
Currently the markovchain package relies on expm package for powering of transition matrix (it loads the operator %^%). Check whether it were possible to implement the power internally