[petsc-users] DMDAVecGetArrayF90

[Carol.Brickley at awe.co.uk]

Hi,

I am using DMDAVecGetArrayF90 with a pointer in petsc 3.4.3. Is there anyway of finding out what the start and finish indices of the pointer coming back from this routine are per processor?

I have declared:


DM                              ::   da
Vec                             ::  X
PetscScalar,pointer      :: xv(:,:)=>null()


DMDAVecGetArrayF90(da,X,xv,ierr)


Do j = 1, m
  Do i = 1, n
    xv(i-1.j-1) = rhs2d(i,j)
  enddo
enddo


DMDAVecRestoreArrayF90(da,X,xv,ierr)

When we have xv(0:n-1,0:m-1) from DMDAVecGetArrayF90 it works fine but other processors may give back eg; xv(0:n-1,1:m-1) or something else and it will fail as the second index
now goes below 0 to 1. Any ideas of what I can do next?


Thanks,

Carol


Dr Carol Brickley
BSc,PhD,ARCS,DIC,MBCS

Senior Software Engineer
Applied Computer Science
DS+T,
AWE
Aldermaston
Reading
Berkshire
RG7 4PR

Direct: 0118 9855035



___________________________________________________
____________________________

The information in this email and in any attachment(s) is 
commercial in confidence. If you are not the named addressee(s) 
or 
if you receive this email in error then any distribution, copying or 
use of this communication or the information in it is strictly 
prohibited. Please notify us immediately by email at 
admin.internet(at)awe.co.uk, and then delete this message from 
your computer. While attachments are virus checked, AWE plc 
does not accept any liability in respect of any virus which is not 
detected.

AWE Plc
Registered in England and Wales
Registration No 02763902
AWE, Aldermaston, Reading, RG7 4PR