Sure I'd be happy to contribute it. I'll need to do some work to convert it from "quick hack" into something people would actually want to use. I'm still trying to make it faster too - I think a lot of the dot product calls in the current version can be expressed in terms of calls to gemm(), which might give quite a nice speedup if I can get the details right.
On Thu, Jun 5, 2014 at 7:20 AM, Miles Lubin <[email protected]> wrote: >> Any idea why type inference is failing on full(chol(A, :U)) ? > > > It can't be determined from the types of the arguments what the result will > be, because you could give different symbols to access different parts of > the factorization. I think this was a design choice to make it easier to > access. > > If you'd like to contribute your chol() implementation to DualNumbers, we'd > be happy to have it. It'll need to be restructured a bit to work with > cholfact().
