(Ted Harding) <[EMAIL PROTECTED]> writes: > Sorry! Missed a trick here: > > At <- t(A) > Bt <- t(B) > E <- B%*%A > Et <- t(E) > A%*%solve(Et%*%E + C)%*%Et%*%D > > (saves 2 multiplications at the relatively cheap cost of 1 transpose)
Well, you might consider getting rid of the first two transposes since you're not actually using them for anything.... -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html