I have a problem Matthew...this is a Fortran code..which I don't think this routine is compatible with. Is there any other way?
Matthew Knepley wrote: > On Dec 3, 2007 11:49 AM, Tim Stitt <timothy.stitt at ichec.ie> wrote: > >> Hi all, >> >> Is there a quick way to map a global index for a parallel vector to a >> local mapping tuple (p,i) were 'p' represents the process containing the >> value and 'i' is the local index number on that process? >> > > PetscMapGetGlobalRange(&v->map,const &range); > for(p = 0; p < numProcs; ++p) if (range[p+1] > globalInd) break; > localInd = globalInd - range[p]; > > Matt > > >> As always, thanks in advance for any information provided. >> >> Tim. >> >> -- >> Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie> >> HPC Application Consultant - ICHEC (www.ichec.ie) >> >> Dublin Institute for Advanced Studies >> 5 Merrion Square - Dublin 2 - Ireland >> >> +353-1-6621333 (tel) / +353-1-6621477 (fax) >> >> >> > > > > -- Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie> HPC Application Consultant - ICHEC (www.ichec.ie) Dublin Institute for Advanced Studies 5 Merrion Square - Dublin 2 - Ireland +353-1-6621333 (tel) / +353-1-6621477 (fax)
