Re: [Rd] Problem using F77_CALL(dgemm) in a package

2011-02-22 Thread Jason Rudy
Thanks for the tip. That API could make my work considerably easier. Jason On Tue, Feb 22, 2011 at 7:37 AM, Douglas Bates wrote: > On Sun, Feb 20, 2011 at 4:39 PM, Jason Rudy wrote: >> I've never used C++ before, so for this project I think I will stick >> with just using

Re: [Rd] Problem using F77_CALL(dgemm) in a package

2011-02-22 Thread Jason Rudy
B)) user system elapsed 2.787 0.029 1.586 Even so, it seems the .Call interface has a lot of advantages. I think I will use it for this project and see how I like it. Jason On Tue, Feb 22, 2011 at 7:27 AM, Simon Urbanek wrote: > On Feb 22, 2011, at 1:55 AM, Jason Rudy wrote: >

Re: [Rd] Problem using F77_CALL(dgemm) in a package

2011-02-22 Thread Jason Rudy
0.658   0.004   0.665 > > (it's not just a measurement error - it's consistent for more replications > etc. - but it's really negligible - possibly just due to dispatch of %*%) > > Cheers, > Simon > > > On Feb 20, 2011, at 5:23 PM, Jason Rudy wrote: > >

Re: [Rd] Problem using F77_CALL(dgemm) in a package

2011-02-20 Thread Jason Rudy
I've never used C++ before, so for this project I think I will stick with just using the BLAS and LAPACK routines directly. Another issue is that I will need to do some sparse matrix computations, for which I am planning to use CSPARSE, at least to begin with. I am interested by RcppArmadillo, an

Re: [Rd] Problem using F77_CALL(dgemm) in a package

2011-02-20 Thread Jason Rudy
completing this project. -Jason On Sun, Feb 20, 2011 at 7:42 AM, Prof Brian Ripley wrote: > Look a close look at matprod in src/main/array in the R sources. > Hint: it is the other dimensions you have wrong. > > And as BLAS is Fortran, counts do start at 1. > > On Sat, 19 Feb 20

[Rd] Problem using F77_CALL(dgemm) in a package

2011-02-20 Thread Jason Rudy
m missing something simple or obvious, as I have never done this before and am proceeding with only google and the R docs as my guide. I am wondering if anybody can see what I'm doing wrong here, or perhaps something I could do to try to fix it. Any assistance would be greatly appreciated.