Displaying 2 results from an estimated 2 matches for "do_it_here".
2010 Aug 23
1
Speeding up matrix multiplies
...------------------------------------------------------------
static void matprod(double *x, int nrx, int ncx,
double *y, int nry, int ncy, double *z)
{
char *transa = "N", *transb = "N";
int i, j, k;
double one = 1.0, zero = 0.0;
LDOUBLE sum;
Rboolean do_it_here;
/* NOTE: ncx must be equal to nry. */
if (nrx==1 && ncy==1) {
/* Do dot product quickly. */
sum = 0.0;
for (i = 0; i<ncx; i++)
sum += x[i] * y[i];
z[0] = sum;
} else if (nrx > 0 && ncx > 0 && nry > 0 &...
2010 Aug 23
1
Internal state indicating if a data object has NAs/no NAs/not sure (Was: Re: Speeding up matrix multiplies)
...--------
>
> static void matprod(double *x, int nrx, int ncx,
> ? ? ? ? ? ? ? ? ? ?double *y, int nry, int ncy, double *z)
> {
> ? ?char *transa = "N", *transb = "N";
> ? ?int i, ?j, k;
> ? ?double one = 1.0, zero = 0.0;
> ? ?LDOUBLE sum;
> ? ?Rboolean do_it_here;
>
> ? ?/* NOTE: ncx must be equal to nry. */
>
> ? ?if (nrx==1 && ncy==1) {
>
> ? ? ? ?/* Do dot product quickly. */
>
> ? ? ? ?sum = 0.0;
> ? ? ? ?for (i = 0; i<ncx; i++)
> ? ? ? ? ? ?sum += x[i] * y[i];
> ? ? ? ?z[0] = sum;
>
> ? ?} else if (nrx...