On Behalf Of Martin Maechler
Sent: Saturday, April 05, 2008 1:52 PM
To: Rory Winston
Cc: r-devel@r-project.org
Subject: Re: [Rd] Adding a Matrix Exponentiation Operator
>>>>> "RW" == Rory Winston <[EMAIL PROTECTED]>
>>>>> on Sat, 5 Apr 2008 14:4
Orange
-Original Message-
From: Martin Maechler <[EMAIL PROTECTED]>
Date: Tue, 8 Apr 2008 15:57:18
To:Vincent Goulet <[EMAIL PROTECTED]>
Cc:Rory Winston <[EMAIL PROTECTED]>, r-devel@r-project.org
Subject: Re: [Rd] Adding a Matrix Exponentiation Operator
>>>>&g
> "rw" == rory winston <[EMAIL PROTECTED]>
> on Tue, 8 Apr 2008 14:06:26 + writes:
rw> Hi martin
rw> Actually I did read your email. If you had actually read *my* email you
would have seen that I am not using the windows version, which is what you
explicitly referred to e
> "VG" == Vincent Goulet <[EMAIL PROTECTED]>
> on Tue, 08 Apr 2008 09:28:00 -0400 writes:
VG> Le dim. 6 avr. à 07:01, Rory Winston a écrit :
>> Hi Martin
>>
>> Thanks for the detailed reply. I had a look at the matrix power
>> implementation in the actuar package a
Le dim. 6 avr. à 07:01, Rory Winston a écrit :
> Hi Martin
>
> Thanks for the detailed reply. I had a look at the matrix power
> implementation in the actuar package and the modified version in the
> expm
> package. I have a couple of questions/comments:
>
> 1. Firstly, I seem to have trouble lo
Hi Martin
I have changed the implementation slightly so that it now handles complex
matrices as well. Take a look and see what you think. I realise that a lot
of the if/else mode checking could be merged.
Cheers
Rory
SEXP do_matexp(SEXP call, SEXP op, SEXP args, SEXP rho)
{
int nrows, ncols;
Hi Martin
Thanks for the detailed reply. I had a look at the matrix power
implementation in the actuar package and the modified version in the expm
package. I have a couple of questions/comments:
1. Firstly, I seem to have trouble loading expm.
> install.packages("expm",repos="http://R-Forge.R-p
Thanks Antonio, thats a good suggestion.
On Sat, Apr 5, 2008 at 5:51 PM, Antonio, Fabio Di Narzo <
[EMAIL PROTECTED]> wrote:
> 2008/4/5, Rory Winston <[EMAIL PROTECTED]>:
>
> >
> > /* Convenience function */
> > static void copyMatrixData(SEXP a, SEXP b, int nrows, int ncols, int
> mode) {
> >
> "RW" == Rory Winston <[EMAIL PROTECTED]>
> on Sat, 5 Apr 2008 14:44:44 +0100 writes:
RW> Hi all I recently started to write a matrix
RW> exponentiation operator for R (by adding a new operator
RW> definition to names.c, and adding the following code to
RW> arrays.c).
2008/4/5, Rory Winston <[EMAIL PROTECTED]>:
>
> /* Convenience function */
> static void copyMatrixData(SEXP a, SEXP b, int nrows, int ncols, int mode) {
> for (int i=0; i < ncols; ++i)
> for (int j=0; j < nrows; ++j)
> REAL(b)[i * nrows + j] = REAL(a)[i * nrows + j];
> }
I wou
Hi all
I recently started to write a matrix exponentiation operator for R (by
adding a new operator definition to names.c, and adding the following code
to arrays.c). It is not finished yet, but I would like to solicit some
comments, as there are a few areas of R's internals that I am still feelin
11 matches
Mail list logo
