intersect()

static2 void intersect	(	char *	uplo,
		char *	diag,
		Int	j,
		Int	start,
		Int	end,
		Int	action,
		Int *	ptrsizebuff,
		double **	pptrbuff,
		double *	ptrblock,
		Int	m,
		Int	n,
		MDESC *	ma,
		Int	ia,
		Int	ja,
		Int	templateheight0,
		Int	templatewidth0,
		MDESC *	mb,
		Int	ib,
		Int	jb,
		Int	templateheight1,
		Int	templatewidth1
	)

Definition at line 173 of file pdtrmr2.c.

{
  /* int       un = 1; only when we use dcopy instead of memcpy */
  double *ptrstart;
  Int   offset, nbline;
  Int   intervalsize;
  assert(start < end);
  assert(j >= 0 && j < n);
  nbline =
        insidemat(uplo, diag, start, j, m, n, &offset);
  if (nbline <= 0)
    return;
  start += offset;
  if (start >= end)
    return;
  intervalsize = min(end - start, nbline);
  (*ptrsizebuff) += intervalsize;
  switch (action) {
  case SENDBUFF:        /* fill buff with local elements to be sent */
    ptrstart = ptrblock + localindice(start + ia, j + ja,
                                      templateheight0, templatewidth0, ma);
    memcpy((char *) (*pptrbuff), (char *) ptrstart,
           intervalsize * sizeof(double));
    /* dcopy_(&intervalsize, (char *) (ptrstart), &un, (char *) (*pptrbuff),
     * &un); */
    (*pptrbuff) += intervalsize;
    break;
  case RECVBUFF:        /* fill local memory with the values received */
    ptrstart = ptrblock + localindice(start + ib, j + jb,
                                      templateheight1, templatewidth1, mb);
    memcpy((char *) ptrstart, (char *) (*pptrbuff),
           intervalsize * sizeof(double));
    /* dcopy_(&intervalsize, (char *) (*pptrbuff), &un, (char *) (ptrstart),
     * &un); */
    (*pptrbuff) += intervalsize;
    break;
  case SIZEBUFF:        /* computation of sizebuff */
    break;
  default:
    printf("action is  %d outside the scope of the case [0..2] !! \n ", action);
    exit(0);
    break;
  };    /* switch (action) */
}/* intersect() */

Here is the call graph for this function:

Here is the caller graph for this function: